頻率提取的MATLAB實驗

我寫的壹部分，能畫出原信號圖和fft變換後的頻譜圖：

fs=1000;%設定采樣頻率

N=1024;

n=0:N-1;

t=n/(fs*50);

sf=100;%設定正弦信號頻率

%生成信號

x=sin(2*pi*sf*t)+sin(6*sf*pi*t);

figure;

subplot(221);

plot(t,x);%作信號的時域波形

xlabel('t');

ylabel('y');

title('時域波形');

grid;

subplot(222);

T=1;

t=linspace(0,T,N);

x=sin(2*pi*sf*t)+cos(6*sf*pi*t);

dt=t(2)-t(1);

f=1/dt;

X=fft(x);

F=X(1:N/2+1);

f=f*(0:N/2)/N;

plot(f,abs(F)*2/N)

xlabel('頻率(Hz)');

ylabel('幅值');

title('幅頻譜圖N=128');

grid;

上一篇:開發人員工具顯示的代碼為什麽跟右鍵查看源代碼顯示的代碼不壹樣呢？

下一篇:java中定義HashSet時能否指定其容量？若能，如何定義？