典型函数的频谱

典型函数的频谱（矩形窗函数, 汉宁窗函数，直线，阶跃函数，δ函数，方波，三角波等），如图13~18所示。

21.510.50矩形窗函数的时域波形图050100150200250300矩形窗函数频域波形图200幅值1000050频率100150 图13

10.80.60.40.20050100150200250300δ函数的时域波形图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.20050100150200250300150汉宁窗函数频域波形图100幅值500050频率100150 图 16

21.510.50阶跃函数的时域波形图050100150200250300300阶跃函数的频域波形图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('幅值');

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('幅值');

本文来源：https://www.bwwdw.com/article/oboa.html