频谱分析matlab编程,Matlab编程实现FFT实践及频谱分析

Matlab编程实现FFT实践及频谱分析 by li343

2007-12-17 08:38

1．用Matlab产生正弦波,矩形波,以及白噪声信号，并显示各自时域波形图2．进行FFT变换，显示各自频谱图，其中采样率，频率、数据长度自选

3．做出上述三种信号的均方根图谱,功率图谱,以及对数均方根图谱

4．用IFFT傅立叶反变换恢复信号，并显示恢复的正弦信号时域波形图

源程序

%******************************************************************** *****%

% FFT实践及频谱分

析 %

%******************************************************************** *****%

%******************************************************************** *****%

%***************1.正弦波****************%

fs=100;%设定采样频率

N=128;

n=0:N-1;

t=n/fs;

f0=10;%设定正弦信号频率

%生成正弦信号

x=sin(2*pi*f0*t);

figure(1);

subplot(231);

plot(t,x);%作正弦信号的时域波形

xlabel('t');

ylabel('y');

title('正弦信号y=2*pi*10t时域波形');

grid;

%进行FFT变换并做频谱图

y=fft(x,N);%进行fft变换

mag=abs(y);%求幅值

f=(0:length(y)-1)'*fs/length(y);%进行对应的频率转换

figure(1);

subplot(232);

plot(f,mag);%做频谱图

axis([0,100,0,80]);

xlabel('频率(Hz)');

ylabel('幅值');

title('正弦信号y=2*pi*10t幅频谱图N=128');

grid;

%求均方根谱

sq=abs(y);