典型函数的频谱

典型函数的频谱（矩形窗函数, 汉宁窗函数，直线，阶跃函数，函数，方波，三角波等），如图1318所示。图13图 14图 15图 16图 17图18此部分MATLAB代码如下：t=0:0.001:0.2;N=256;FS=300;w=boxcar(N); %产生信号figure;plot(w);title(矩形窗函数的时域波形图);axis(0,260,0,2);grid on;y=fft(w,N); %FFT运算mag=abs(y);%取幅值f=(0:length(y)-1)*FS/length(y);figure;plot(f(1:N/2),mag(1:N/2); %输出FS/2点幅频谱图title(矩形窗函数频域波形图);grid;xlabel(频率);ylabel(幅值);t=0:0.001:0.2;N=256;FS=300;w=hanning(N); %产生信号figure;plot(w);title(汉宁窗函数的时域波形图);grid on;y=fft(w,N); %FFT运算mag=abs(y); %取幅值f=(0:length(y)-1)*FS/length(y);figure;plot(f(1:N/2),mag(1:N/2); %输出FS/2点幅频谱图title(汉宁窗函数频域波形图);grid on;xlabel(频率);ylabel(幅值);t=0:0.001:0.2;N=256;FS=300;w=1; %产生信号y=fft(w,N); %FFT运算mag=abs(y); %取幅值f=(0:length(y)-1)*FS/length(y);figure;plot(f(1:N/2),mag(1:N/2); %输出FS/2点幅频谱图title(直线频域波形图);grid on;xlabel(Frequency(Hz);ylabel(Magnitude);%阶跃函数的频域波图clc;clf;t=0:0.001:0.2;N=256;FS=300;w=ones(1,N); %产生信号figure;plot(w);title(阶跃函数的时域波形图);grid on;y=fft(w,N); %FFT运算mag=abs(y);%取幅值f=(0:length(y)-1)*FS/length(y);figure;plot(f(1:N/2),mag(1:N/2); %输出FS/2点幅频谱图title(阶跃函数的频域波形图);grid on;xlabel(频率);ylabel(幅值);t=0:0.001:0.2;N=256;FS=300;w=zeros(1,N);w(1)=1; %产生信号figure;plot(w);grid on;title(函数的时域波形图);y=fft(w,N);%FFT运算mag=abs(y);%取幅值f=(0:length(y)-1)*FS/length(y);figure;plot(f(1:N/2),mag(1:N/2); %输出FS/2点幅频谱图title(函数的频域波形图);grid on;xlabel(频率);ylabel(幅值);t=0:0.001:0.2;N=256;FS=300;w=square(2*pi*50*t); %产生信号figure;plot(t,w); title(方波的时域波形图);axis(0,0.2,-0.2,1.2);grid on; y=fft(w,N); %FFT运算mag=abs(y); %取幅值f=(0:length(y)-1)*FS/length(y);figure;plot(f(1:N/2),mag(1:N/2); %输出FS/2点幅频谱图title(方波的频域波形图);grid on;xlabel(频率);ylabel(幅值);t=0:0.001:0.2;N=256;FS=300;w=sawtooth(2*pi*50*t,0.5); figure;plot(t,w); grid on;title(三角波的时域波形图);%产生信号y=fft(w,N); %FFT运算mag=abs(y); %取幅值f=(0:length(y)-1)*FS/length(y);f