plot(t0,f0,'r',t0,y,'b--'); figure(3) for i=1:WaveNum F(i)=sqrt(A_List(i).^2+B_List(i).^2)./2; end F=[A0 F]; xx=0:length(F)-1; xx=[-fliplr(xx) xx(2:end)]; F=[fliplr(F) F(2:end)]; stem(xx,F); title('三角波-幅度谱'); figure(4) beta=-atan(B_List./A_List); beta=[0 beta]; xx=-length(beta)+1:length(beta)-1; beta=[fliplr(beta),beta(2:end)]; stem(xx,beta); title('三角波-相位谱'); end function y=testFunction2(t) global T; y=abs(2-mod(t,T))-1; end function y=testFunction1 global h; global a; global T; a=1; T=4; t=-3*a:0.01:T-a; e1=1/2+1/2.*sign(t); e2=1/2+1/2.*sign(t-a); y=h.*(e1-e2); end 三、实验总结： （1）通过MATLAB实现周期信号傅立叶级数分解与综合方法，掌握基本方法，对傅里叶级数有更深的认识 （2）通过利用MATLAB 观察周期信号频谱特点，对周期信号频谱特点进行了理解
实验题目 一、实验目的： 实验四 连续时间信号的傅立叶变换 （1）掌握连续信号傅立叶变换与逆变换的计算方法 （2）掌握利用MATLAB 实现连续时间信号傅立叶变换的方法 二、实验内容： 1、利用 fourier()命令求解如下信号的傅立叶变换，给出f(t)的波形图以及F(j?)的表达式和幅度频谱图： (1) 钟形脉冲：f(t)?e (2) 符号函数： 实验程序： (1)钟形脉冲： t)2?(2,???t??? 实验结果： syms t w ft=exp(-(t/2)^2); Fw=fourier(ft); subplot(2,1,1); ezplot(ft); subplot(2,1,2) ezplot(Fw); (2) 符号函数： 实验结果： syms t w ft=sym('heaviside(t)-heaviside(-t)'); Fw=fourier(ft) ezplot(ft,[-,2]); 实验画图由于在t=0处傅里叶变换有冲激函数，无