fbh.m

使用FFT分析方波信号的频谱

fs=100;%设定采样频率
N=16;
n=0:N-1;
t=n/fs;
f0=10;%设定方波信号频率
%生成方波信号
x=square(2*pi*f0*t);
figure(1);
subplot(221);
stem(t,x);%作方波信号的时域波形
xlabel('t');
ylabel('y');
axis([0 0.2 -1 1]);
title('10hz方波信号的时域波形');
grid;
%进行FFT变换并做频谱图
y=fft(x,N);%进行fft变换
mag=abs(y);%求幅值
f=(0:length(y)-1)'*fs/length(y);%进行对应的频率转换
figure(1);
subplot(222);
stem(f,mag);%做频谱图
axis([0,100,0,120]);
xlabel('频率(Hz)');
ylabel('幅值');
title('方波信号幅频谱图');
set(gca,'XTickMode','manual','XTick',[0,10,100]);
grid;
%求功率谱
sq=abs(y);
power=sq.^2;
figure(1);
subplot(223);
stem(f,power);
xlabel('频率(Hz)');
ylabel('功率谱');
title('方波信号功率谱');
set(gca,'XTickMode','manual','XTick',[0,10,100]);
grid;
%求相位谱
w=angle(x);
figure(1);
subplot(224);
stem(f,w);
xlabel('频率(Hz)');
ylabel('相位谱');
title('方波信号的相位谱');
axis([0,100,-4,4]);
grid;