fft1234.m

快速傅立叶算法,希望大家能够喜欢.我很喜欢这个网站.

f1=1;%f1=8;坐标调整在横轴0-32，每隔4个点去一个，效果最好
f2=12.375;
ts=1;
fs=1000;
t=0:1/fs:0.2;%t=0:1/fs:0.15;%不在一个完整周期内的采样，共取了150个点
%均值为0，方差为1的白色噪声
ww1=normrnd(0,1,1,size(t));%调节SNR数值的时候只需调方差即可
%带通滤波器，使白噪声通过后变成有色噪声
k=0.0201;
temp1=k*[1 0 -2 0 -1.23];
temp2=[1 -1.637 2.237 -1.307 0.641];
www1=filter(temp1,temp2,ww1);
%信噪比
snr1=1./(2.*var(www1)); 
disp('信噪比SNR(db)为:');
snr=-10*log10(snr1)%无分号
y=cos(2*pi*4.67*t)+cos(2*pi*9.34*t)+cos(2*pi*14.01*t)+cos(2*pi*23.35*t)+cos(2*pi*47.665*t)+www1;%randn(size(t));%调节信号幅值的时候，噪声会淹没信号
n=length(y);
r=fft(y);%r=fft(y)/n;
r=fftshift(r);
f=linspace(-fs/2,fs/2,n);
subplot(2,1,1);
plot(t,y);
subplot(2,1,2);
plot(f,abs(r));
xlabel('频率/Hz');
ylabel('dB');
set(gca,'XTick',-fs/2:50:fs/2);
axis tight;
grid on;