function?sample（）
t0=10;???????????	%定義時間長度
ts=0.001;????fs=1/ts;??
t=[-t0/2:ts:t0/2];?	%定義時間序列
df=0.5;???????		%定義頻率分辨率
x=sin（200*t）;?m=x./（200*t+eps）;
w=t0/（2*ts）+1;???	%確定t=0的點
m（w）=1;??????????%修正t=0點的信號值
m=m.*m;?
[Mmndfy]=fft_seq（mtsdf）;?	?%傅立葉變換
M=M/fs;
f=[0:dfy:dfy*length（mn）-dfy]-fs/2;?%定義頻率序列
figure（1）?
subplot（211）;?plot（tm）;
xlabel（‘時間‘）;ylabel（‘幅值‘）;title（‘原始信號（fh=200/2piHz）的波形‘）;
axis（[-0.150.1501.5]）;
subplot（212）;
plot（fabs（fftshift（M）））;
xlabel（‘頻率‘）;ylabel（‘幅值‘）;
axis（[-50050000.03]）;title（‘原始信號的頻譜‘）;
t0=10;?????		?%信號持續的時間
ts1=0.005;???		?%滿足抽樣條件的抽樣間隔
fs1=1/ts1;
t1=[-t0/2:ts1:t0/2];	?%定義滿足抽樣條件的時間序列
x1=sin（200*t1）;??m1=x1./（200*t1+eps）;??w1=t0/（2*ts1）+1;????
m1（w1）=1;?????	?%修正t=0時的信號值
m1=m1.*m1;???	?%定義信號
[M1mn1df1]=fft_seq（m1ts1df）;?%對滿抽樣條件的信號進行傅立葉變換
M1=M1/fs1;N1=[M1M1M1M1M1M1M1M1M1M1M1M1M1];
f1=[-7*df1*length（mn1）:df1:6*df1*length（mn1）-df1]-fs1/2;
figure（2）???
subplot（211）;?stem（t1m1）;
xlabel（‘時間‘）;ylabel（‘幅值‘）;
title（‘抽樣正常（fs=200Hz）時的信號波形‘）;axis（[-0.150.1501]）;
subplot（212）?
plot（f1abs（fftshift（N1）））;
xlabel（‘頻率‘）;ylabel（‘幅值‘）;axis（[-50050000.05]）;
title（‘抽樣正常時的信號頻譜‘）;axis（[-500500-0.010.03]）;
t0=10;?????	?	%信號持續的時間
ts2=0.01;???		%不滿足抽樣條件的抽樣間隔
fs2=1/ts2;
t2=[-t0/2:ts2:t0/2];??%定義不滿足抽樣條件的時間序列
x2=sin（200*t2）;?m2=x2./（200*t2+eps）;??w2=t0/（2*ts2）+1;???
m2（w2）=1;?????	%修正t=0時的信號值
m2=m2.*m2;??????%定義信號
[M2mn2df2]=fft_seq（m2ts2df）;%對不滿足抽樣條件的信號進行傅立葉變換
M2=M2/fs2;N2=[M2M2M2M2M2M2M2M2M2M2M2M2M2];
f2=[-7*df2*length（mn2）:df2:6*df2*length（mn2）-df2]-fs2/2;
figure（3）?
subplot（211）;?stem（t2m2）;
xlabel（‘時間‘）;ylabel（‘幅值‘）;title（‘抽樣失真（fs=100Hz）時的信號波形‘）;
axis（[-0.150.1501]）;subplot（212）
plot（f2abs（fftshift（N2）））;
xlabel（‘頻率‘）;ylabel（‘幅值‘）;axis（[-50050000.02]）;
title（‘抽樣失真時的信號頻譜‘）;axis（[-5005000.0050.02]）;
function?[Mmdf]=fft_seq（mtsdf）
fs=1/ts;
if?nargin==2
?????n1=0
else
?????n1=fs/df
end
n2=length（m）;n=2^（max（nextpow2（n1）nextpow2（n2）））;
M=fft（mn）;m=[mzeros（1n-n2）];df=fs/n

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????2298??2019-01-20?10:03??信號自學報告\e.m

?????文件??????78112??2019-01-20?10:46??信號自學報告\《信號與系統》自學報告.docx

?????目錄??????????0??2019-05-03?18:53??信號自學報告

-----------?---------??----------?-----??----

????????????????80410????????????????????3