matlab与傅里叶论文

老师让我们写一篇MATLAB与傅里叶的论文，可是小弟我既不懂傅里叶级数与变换，又没听MATLAB的课，所以很是头疼！哪位好心的大善人帮帮忙，帮我写篇MATLAB与傅里叶的论文，小弟我感激不尽啊！即使是给点材料也可以！老师要求每句程序要加注释，最后还要在机子上运行，我很惆怅啊！不要太难，简单些就可以！谢谢各位啊！小弟感激涕零......

以正弦波为例

1．用Matlab产生正弦波,并显示时域波形图。

2. 进行FFT傅里叶变换，显示频谱图，其中采样率，频率、数据长度自选。

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

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

%***************正弦波****************%
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);
figure(1);
subplot(233);
plot(f,sq);
xlabel('频率(Hz)');
ylabel('均方根谱');
title('正弦信号y=2*pi*10t均方根谱');
grid;
%求功率谱
power=sq.^2;
figure(1);
subplot(234);
plot(f,power);
xlabel('频率(Hz)');
ylabel('功率谱');
title('正弦信号y=2*pi*10t功率谱');
grid;
%求对数谱
ln=log(sq);
figure(1);
subplot(235);
plot(f,ln);
xlabel('频率(Hz)');
ylabel('对数谱');
title('正弦信号y=2*pi*10t对数谱');
grid;
%用IFFT做傅里叶反变换恢复原始信号
xifft=ifft(y);
magx=real(xifft);
ti=[0:length(xifft)-1]/fs;
figure(1);
subplot(236);
plot(ti,magx);
xlabel('t');
ylabel('y');
title('通过IFFT转换的正弦信号波形');
grid;