典型函数的频谱.doc

典型函数的频谱 典型函数的频谱〔矩形窗函数, 汉宁窗函数，直线，阶跃函数，δ函数，方波，三角波等〕，如图13~18所示 21.510.50矩形窗函数的时域波形图050100150202350300矩形窗函数频域波形图200幅值1000050频率100150 图13 10.80.60.40.20230__0150202350300δ函数的时域波形图21.5δ函数的频域波形图幅值10.50050频率100150 图 14 方波的时域波形图10.5001500.020.040.060.080.10.120.140.160.180.2方波的频域波形图100幅值500050频率100150 图 15 汉宁窗函数的时域波形图10.80.60.40.20230__0150202350300150汉宁窗函数频域波形图100幅值500050频率100150 图 16 21.510.50阶跃函数的时域波形图050100150202350300300阶跃函数的频域波形图200幅值1000050频率100150 图 17 10.50-0.5-1三角波的时域波形图00.020.040.060.080.10.120.140.160.180.28060三角波的频域波形图幅值40200050频率100150 图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); ?T运算 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); ?T运算 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); ?T运算 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); ?T运算 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);?T运算 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); ?T运算 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); ?T运算 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('幅值'); 第 页 共 页。