信号的Fourier分析的matlab编程

一、实验目的：

1、掌握周期信号的频谱-----Fourier级数的分析方法。

2、深入理解信号频谱的概念，掌握典型信号的频谱以及Fourier变换的主要性质。

（1）原函数：
t=-1:0.001:3;
T=2;
y=rectpuls(t-0.5,T*0.5)-rectpuls(t-0.75*T,0.5*T);
plot(t,y)

合成：方波信号的幅值与相位
N=9; 
n=-N:2:N; 
cn=4./n./pi; 
subplot(2,1,1); 
stem(n,abs(cn));
ylabel('Cn的幅值'); 
subplot(2,1,2); 
stem(n,angle(cn)); 
ylabel('Cn的相位');
xlabel('\omega/\omega0')

（2）吉布斯函数近似
写函数文件：
function s=f(m)
s=0;p=0;
t=0:0.01:10;
for n=1:2:m
    p=p+1/n*(sin(n*t));  %取Ω=1
end
s=(4/3.14)*p;
plot(t,s);
保存为f.m，然后在指令窗口执行：
>> clear
>> f(100)  %这里取N=100

2、求信号te-2tε(t)的幅度谱。
>> t=-10:0.01:10; f=abs(1./((2+j*t).^2)); plot(t,f)；

另：
N=8;
n1=-N:-1;
c1=1./((2+n1.*i).^2);
c0=0;
n2=1:N;
c2=1./((2+n2.*i).^2);
cn=[c1 c0 c2];
n=-N:N;
subplot(1,2,1);
stem(n,abs(cn));
ylabel('Cn的幅度');
xlabel('\Omerga');
subplot(1,2,2);
stem(n,angle(cn));
ylabel('Cn的相位');
xlabel('Omega');