典型函数的频谱

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

图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);

figure;

plot(f(1:N/2),mag(1:N/2)); %输出FS/2点幅频谱图

title('三角波的频域波形图');

grid on;

xlabel('频率');ylabel('幅值');

本文来源：https://www.2haoxitong.net/k/doc/3ef0b0a9700abb68a882fb4d.html