搜索
典型函数的频谱（矩形窗函数,汉宁窗函数，直线，阶跃函数，δ函数，方波，三角波等），如图13~18所示。05010015020025000.511.52矩形窗函数的时域波形图0501001500100200300矩形窗函数频域波形图频率幅值图1305010015020025030000.20.40.60.81δ函数的时域波形图05010015000.511.52δ函数的频域波形图频率幅值图1400.020.040.060.080.10.120.140.160.180.200.51方波的时域波形图050100150050100150方波的频域波形图频率幅值图1505010015020025030000.20.40.60.81汉宁窗函数的时域波形图050100150050100150汉宁窗函数频域波形图频率幅值图1605010015020025030000.511.52阶跃函数的时域波形图0501001500100200300阶跃函数的频域波形图频率幅值图1700.020.040.060.080.10.120.140.160.180.2-1-0.500.51三角波的时域波形图050100150020406080三角波的频域波形图频率幅值图18此部分MATLAB代码如下：t=0::;N=256;FS=300;w=boxcar(N);%产生信号figure;plot(w);title('矩形窗函数的时域波形图');axis([0,260,0,2]);gridon;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::;N=256;FS=300;w=hanning(N);%产生信号figure;plot(w);title('汉宁窗函数的时域波形图');gridon;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('汉宁窗函数频域波形图');gridon;xlabel('频率');ylabel('幅值');t=0::;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('直线频域波形图');gridon;xlabel('Frequency(Hz)');ylabel('Magnitude');%阶跃函数的频域波图clc;clf;t=0::;N=256;FS=300;w=ones(1,N);%产生信号figure;plot(w);title('阶跃函数的时域波形图');gridon;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('阶跃函数的频域波形图');gridon;xlabel('频率');ylabel('幅值');t=0::;N=256;FS=300;w=zeros(1,N);w(1)=1;%产生信号figure;plot(w);gridon;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('δ函数的频域波形图');gridon;xlabel('频率');ylabel('幅值');t=0::;N=256;FS=300;w=square(2*pi*50*t);%产生信号figure;plot(t,w);title('方波的时域波形图');axis([0,,,]);gridon;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('方波的频域波形图');gridon;xlabel('频率');ylabel('幅值');t=0::;N=256;FS=300;w=sawtooth(2*pi*50*t,;figure;plot(t,w);gridon;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('三角波的频域波形图');gridon;xlabel('频率');ylabel('幅值');