clear?all;
bt=0;???????????????????????????%開始時(shí)間
dt=0.001;???????????????????????%時(shí)間間隔
N=512;??????????????????????????%傅立葉變換點(diǎn)數(shù)
et=bt+N*dt-dt;??????????????????%結(jié)束時(shí)間
t=bt:dt:et;?????????????????????%時(shí)間域
TT=et-bt;????????????????????????%總的時(shí)間
y?=?cos（2*pi*80*t）;?????????????%待分析的信號(hào)
figure（1）;
subplot（321）
plot（ty）;
title（‘信號(hào)‘）;
xlabel（‘時(shí)間（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;???%計(jì)算功率譜密度
subplot（323）
plot（fPyy）;
title（‘根據(jù)定義求出的功率譜‘）;
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（‘信號(hào)自相關(guān)‘）;
xlabel（‘時(shí)間（s）‘）;??????
axis（[0?0.5?-1?1]）;

Pyy2=fftshift（fft（RyyM））/Tf;
subplot（325）
plot（fabs（Pyy2））;
title（‘自編自相關(guān)函數(shù)求出的功率譜‘）;
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（‘系統(tǒng)自相關(guān)函數(shù)求出的功率譜‘）;
xlabel（‘頻率（Hz）‘）;
axis（[-100?100?0?max（abs（Pyy3））]）;

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件???????1498??2009-11-03?09:28??實(shí)驗(yàn)一?信號(hào)的時(shí)頻分析實(shí)驗(yàn)\t1.m

?????文件???????1494??2009-11-03?09:28??實(shí)驗(yàn)一?信號(hào)的時(shí)頻分析實(shí)驗(yàn)\t2.m

?????文件???????1495??2009-11-03?09:28??實(shí)驗(yàn)一?信號(hào)的時(shí)頻分析實(shí)驗(yàn)\t3.m

?????文件???????1495??2009-11-03?09:28??實(shí)驗(yàn)一?信號(hào)的時(shí)頻分析實(shí)驗(yàn)\t4.m

?????文件???????1504??2009-11-03?09:08??實(shí)驗(yàn)一?信號(hào)的時(shí)頻分析實(shí)驗(yàn)\t5.m

?????文件???????1507??2009-11-03?09:21??實(shí)驗(yàn)一?信號(hào)的時(shí)頻分析實(shí)驗(yàn)\t6.m

?????文件???????1508??2009-11-03?09:21??實(shí)驗(yàn)一?信號(hào)的時(shí)頻分析實(shí)驗(yàn)\t7.m

?????文件???????1509??2009-11-03?09:21??實(shí)驗(yàn)一?信號(hào)的時(shí)頻分析實(shí)驗(yàn)\t8.m

?????文件?????732160??2009-11-03?09:35??實(shí)驗(yàn)一?信號(hào)的時(shí)頻分析實(shí)驗(yàn)\實(shí)驗(yàn)一?信號(hào)的時(shí)頻分析實(shí)驗(yàn).doc

?????目錄??????????0??2010-03-02?21:51??實(shí)驗(yàn)一?信號(hào)的時(shí)頻分析實(shí)驗(yàn)

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

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