clear??all
clc
t=0:0.01:3;						%時間間隔為0.01說明采樣頻率為100?Hz
x=sin（2*pi*30*t）+sin（2*pi*500*t）;	%產生主要頻率為30Hz和500Hz的信號
f=x+3.4*randn（1length（t））;			%在信號中加入白噪聲
subplot（121）;
plot（f）;							%畫出原始信號的波形圖
ylabel（‘幅值‘）;
xlabel（‘時間‘）;
title（‘原始信號‘）;
y=fft（f1024）;	%對原始信號進行離散傅里葉變換，參加DFT的采樣點個數為1024
p=y.*conj（y）/1024;					%計算功率譜密度
ff=100*（0:511）/1024;				%計算變換后不同點所對應的頻率值
subplot（122）;
plot（ffp（1:512））;					%畫出信號的頻譜圖
ylabel（‘功率譜密度‘）;
xlabel（‘頻率‘）;
title（‘信號功率譜圖‘）;

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

?????文件????????588??2018-07-25?19:57??8\s8_1.m

?????文件???????3119??2018-07-25?19:57??8\s8_10.m

?????文件????????828??2018-07-25?19:57??8\s8_11.m

?????文件????????582??2018-07-25?19:57??8\s8_12.m

?????文件????????783??2018-07-25?19:57??8\s8_13.m

?????文件????????788??2018-07-25?19:57??8\s8_2.m

?????文件???????1797??2018-07-25?19:57??8\s8_3.m

?????文件????????404??2018-07-25?19:57??8\s8_4.m

?????文件???????1102??2018-07-25?19:57??8\s8_5.m

?????文件????????292??2018-07-25?19:57??8\s8_6.m

?????文件????????540??2018-07-25?19:57??8\s8_7.m

?????文件????????601??2018-07-25?19:57??8\s8_8.m

?????文件???????1112??2018-07-25?19:57??8\s8_9.m

?????目錄??????????0??2018-09-09?12:02??8

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

????????????????12536????????????????????14