文档介绍：四、实验内容 1. 程序: clear all ; syms tf; f=fourier(exp((-2)*abs(t))); ezplot(f); 曲线: 2. 程序: clear all ; syms tw; ifourier(1/(1+(w^2)),t) 结果: ans = 1/2*exp(t)*heaviside(-t)+1/2*exp(-t)*heaviside(t) 3. 程序: clear all ; r=; t=-5:r:5; N=200; w=2*pi; k=-N:N; w=k*w/N; f1=1/2*exp(-2*t).*stepfun(t,0); F=r*f1*exp(-j*t'*w); F1=abs(F); P1=angle(F); subplot(3,1,1); plot(t,f1); grid on ; xlabel( 't' ); ylabel( 'f(t)' ); title( 'f(t)' ); subplot(3,1,2); plot(w,F1); xlabel( 'w' ); grid on ; ylabel( 'F(jw)' ); -6 -4 -2 0 2 4 6 0 0. 1 0. 2 0. 3 0. 4 0. 5 0. 6 0. 7 0. 8 0. 9 1 w 4 / (4+w 2) subplot(3,1,3); plot(w,P1*180/pi); grid; xlabel( 'w' ); ylabel( ' 相位(度) ' ); 频谱图: 附加: f(t-1) 程序: clear all ; r=; t=-5:r:5; N=200; w=2*pi; k=-N:N; w=k*w/N; f1=1/2*exp(-2*(t-1)).*stepfun(t-1,0); F=r*f1*exp(-j* (t -1) '*w); F1=abs(F); P1=angle(F); subplot(3,1,1); plot(t,f1); grid on ; xlabel( 't' ); ylabel( 'f(t)' ); title( 'f(t-1)' ); subplot(3,1,2); plot(w,F1); xlabel( 'w' ); grid on ; ylabel( 'F(jw)' ); subplot(3,1,3); plot(w,P1*180/pi); grid; xlabel( 'w' ); ylabel( ' 相位(度) ' ); 频谱图: -5 -4 -3 -2 -1 0 1 2 3 4 5 0 0. 5 t f (t ) f(t ) -8 -6 -4 -2 0 2 4 6 8 0 0. 2 0. 4 w F (j w) -8 -6 -4 -2 0 2 4 6 8 -100 0 100 w 相位( 度) 4. 程序: clear all ; R=; t=-2:R:2; f=stepfun(t,-1)-stepfun(t,1); f1=f.*exp(-j*5*t); f2=f.*exp(j*5*t); W1=pi*5; N=500; k=-N:N; W=k*W1/N; F1=f1*exp(-j*t'*W)*R; F2=f2*exp(-j*t'*W)*R; F1=real(F1); F2=real(F2); subplot(2,1,1); plot(W,F1); xlabel( 'w' ); ylabel( 'F1(jw)' ); title( ' 频谱 F1(jw)' ); subplot(2,1,2); plot(W,F2); xlabel( 'w' ); ylabel( 'F2(jw)' ); title( ' 频谱 F2(jw)' ); 频谱图: -5 -4 -3 -2 -1 0 1 2 3 4 5 0 0. 5 t f (t ) f(t -1) -8 -6 -4 -2 0 2 4 6 8 0 0. 2 0. 4 w F (j w) -8 -6 -4 -2 0 2 4 6 8 -100 0 100 w 相位( 度) 5. 使用 MATLAB 完成下列问题( 1 )求下列函数的傅立叶变换○ 1 sin 2 (2t) 程序: clear all ; syms tf; f=fourier((sin(2*t))^2); f 结果: f= 1/2*pi*(2*dirac(w)-dirac(w-4)-dirac(w+4)) ○ 2 cos( ∏ t)+cos(2 ∏ t) 程序: clear all ; syms tf; f=fourier(cos(pi*t)+cos(2*pi*t)); f 结果: f= pi*(dirac(w-pi)+dirac(w+pi)+d