clear?all;
bt=0;???????????????????????????%開始時間
dt=0.001;???????????????????????%時間間隔
N=512;??????????????????????????%傅立葉變換點數
et=bt+N*dt-dt;??????????????????%結束時間
t=bt:dt:et;?????????????????????%時間域
TT=et-bt;????????????????????????%總的時間
y?=?cos（2*pi*80*t）;?????????????%待分析的信號
figure（1）;
subplot（321）
plot（ty）;
title（‘信號‘）;
xlabel（‘時間（s）‘）;????
axis（[0?0.5?-1.2?1.2]）;?

Y?=fftshift（fft（yN））;
df=1/TT;?????????????????????????%頻率間隔
Tf=N*df;?????????????????????????%分析的頻寬?????????????????????????
f=-Tf/2+df:df:Tf/2;?????????????%頻率域
subplot（322）
plot（fabs（Y））;
title（‘頻譜‘）;
xlabel（‘頻率（Hz）‘）;
axis（[-100?100?0?max（abs（Y））]）;

Pyy=abs（Y/Tf）.^2/TT;???%計算功率譜密度
subplot（323）
plot（fPyy）;
title（‘根據定義求出的功率譜‘）;
xlabel（‘頻率（Hz）‘）;
axis（[-100?100?0?max（abs（Pyy））]）;

M=N;
Ryy=zeros（1M）;
for?m=1:M
????for?n=1:N-m+1
????????Ryy（m）=Ryy（m）+y（n）*y（n+m-1）;
????end
????Ryy（m）=Ryy（m）/（N-m+1）;
end
subplot（324）
plot（tRyy（1:M））;
title（‘信號自相關‘）;
xlabel（‘時間（s）‘）;??????
axis（[0?0.5?-1?1]）;

Pyy2=fftshift（fft（RyyM））/Tf;
subplot（325）
plot（fabs（Pyy2））;
title（‘自編自相關函數求出的功率譜‘）;
xlabel（‘頻率（Hz）‘）;
axis（[-100?100?0?max（abs（Pyy2））]）;

Ryy2=xcorr（yy‘unbiased‘）;
L=length（Ryy2）;
Pyy3=fftshift（fft（Ryy2（L-N+1:L）N））/Tf;
subplot（326）
plot（fabs（Pyy3））;
title（‘系統自相關函數求出的功率譜‘）;
xlabel（‘頻率（Hz）‘）;
axis（[-100?100?0?max（abs（Pyy3））]）;

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

?????文件???????1498??2009-11-03?09:28??實驗一?信號的時頻分析實驗\t1.m

?????文件???????1494??2009-11-03?09:28??實驗一?信號的時頻分析實驗\t2.m

?????文件???????1495??2009-11-03?09:28??實驗一?信號的時頻分析實驗\t3.m

?????文件???????1495??2009-11-03?09:28??實驗一?信號的時頻分析實驗\t4.m

?????文件???????1504??2009-11-03?09:08??實驗一?信號的時頻分析實驗\t5.m

?????文件???????1507??2009-11-03?09:21??實驗一?信號的時頻分析實驗\t6.m

?????文件???????1508??2009-11-03?09:21??實驗一?信號的時頻分析實驗\t7.m

?????文件???????1509??2009-11-03?09:21??實驗一?信號的時頻分析實驗\t8.m

?????文件?????732160??2009-11-03?09:35??實驗一?信號的時頻分析實驗\實驗一?信號的時頻分析實驗.doc

?????目錄??????????0??2010-03-02?21:51??實驗一?信號的時頻分析實驗

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

???????????????744170????????????????????10