squarefft2.m

利用FFT分析方波信号的频谱

Fs=618;
N=310;
n=0:N-1;
t=n/Fs;

x=square(2*pi*200*t)+square(2*pi*202*t)+square(2*pi*206*t);
y=square(2*pi*200*t)+square(2*pi*205*t)+square(2*pi*210*t)+randn(1,length(t));
%信号最高频率f0=206hz，频率分辨率F=2hz
%tp=0.5；
%T取样周期T《=1/2*f0
%N 》=2*f0/F=206

fx=fft(y);
%N个点有N-1个间隔；


ff=(0:N/2-1)*Fs/N;
f=(0:length(fx)-1)'*Fs/length(fx);%进行对应的频率转换
subplot(411);
plot(t,y);grid; title('受噪声干扰方波信号');
subplot(412);
plot(f(1:N/2),abs(fx(1:N/2))); grid;title('受噪声干扰，取样点310,采样率为618hz');axis([0,250,0,300]);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Fs=1236;
N=310;
n=0:N-1;
t=n/Fs;
x=square(2*pi*200*t)+square(2*pi*205*t)+square(2*pi*210*t);
y=square(2*pi*200*t)+square(2*pi*205*t)+square(2*pi*210*t)+randn(1,length(t));
%信号最高频率f0=206hz，频率分辨率F=2hz
%tp=0.5；
%T取样周期T《=1/2*f0
%N 》=2*f0/F=206

fx=fft(y);
%N个点有N-1个间隔；
ff=(0:N/2-1)*Fs/N;
subplot(413);
f=(0:length(fx)-1)'*Fs/length(fx);%进行对应的频率转换

plot(f(1:N/2),abs(fx(1:N/2))); grid;title('受噪声干扰，取样点310,采样率为1236hz');axis([0,250,0,300]);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Fs=618;
N=600
n=0:N-1;
t=n/Fs;

x=square(2*pi*200*t)+square(2*pi*205*t)+square(2*pi*210*t);
y=square(2*pi*200*t)+square(2*pi*205*t)+square(2*pi*210*t)+randn(1,length(t));


%信号最高频率f0=206hz，频率分辨率F=2hz
%tp=0.5；
%T取样周期T《=1/2*f0
%N 》=2*f0/F=206

fx=fft(y);
%N个点有N-1个间隔；

f=(0:length(fx)-1)'*Fs/length(fx);%进行对应的频率转换
ff=(0:N/2-1)*Fs/N;
subplot(414);
plot(f(1:N/2),abs(fx(1:N/2))); grid;title('受噪声干扰，增加取样点500');axis([0,250,0,300]);