十二章matlab答案.docVIP

实验十二：连续时间系统的频域分析 例题： 1、使用MATLAB函数实现下列信号的傅里叶变换，并画出变换后的曲线 求出e-|2t|的傅里叶变换，并画出变换后的曲线 clear all; syms t f; f=fourier(exp((-2)*abs(t))); ezplot(f); 使用MATLAB函数实现下列信号的傅里叶逆变换 已知F（jw）=1/1+w2,求信号的逆傅里叶变换 clear all; syms t w; ifourier(1/(1+(w^2)),t); ans 1/2*exp(-t)*heaviside(t)+1/2*exp(t)*heaviside(-t) 3、使用MATLAB函数实现傅里叶的时移特性 画出f（t）=1/2e-2tu(t)和f(t-1)的频谱图，观察信号时移对频谱的影响 clear all; r=0.02; 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)); subplot(3,1,3); plot(w,P1*180/pi); grid; xlabel(w); ylabel(相位度); 使用MATLAB函数实现下列信号的频移变换 已知f（t）为门函数，求f1（t）= f（t）e-j5t以及f2（t）ej5t的频谱图 clear all; R=0.02; t=-2:R:2; f=stepfun(t,-1)-stepfun(t,1); f1=f.*exp(-j*5*t); f2=f.*exp(j*5*t); N=500; W1=5*pi; k=-N:N; W=k*W1/N; F1=f1*exp(-j*t*W)*R; F2=f1*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（1）求下列函数的傅里叶变换 ① clear all; syms t f; fourier(sin(2*t)*sin(2*t)) ans ans = 1/2*pi*(-dirac(w-4)+2*dirac(w)-dirac(w+4)) ② clear all; syms t f; fourier(cos(pi*t)+cos(2*pi*t)); ans ans = pi*(dirac(w-pi)+dirac(w+2*pi)+dirac(w-2*pi)+dirac(w+pi)) ③ clear all; syms t f; fourier(cos(2*t)*sin(3*t)) ans ans = 1/2*i*pi*(dirac(w+5)-dirac(w-5)+dirac(w+1)-dirac(w-1)) 5、(2)求下列信号的频谱图 ①j将第3题的f1等式改为f1=exp(-3*t+2).*stepfun(t,-1);即可 ② 将第3题的f1等式改为f1=exp(-abs(t)).*cos(t); ③ 将第3题的f1等式改为f1=exp(-t).*sin(2*t).*stepfun(t,0); ④将第3题的f1等式改为f1=sin(t)./t; 同时在t=-5:r:5;后面加以句t(251)=exp(-20); （3）求下列信号的傅里叶逆变换 ① clear all; syms t w; ifourier(4*sin(w)/w.*cos(2*w)); ans ans = heaviside(x+3)-heaviside(x-3)-heaviside(x+1)+heaviside(x-1) ② clear all; syms t w; ifourier(3/(-w*w+j*w-2)); ans ans = -exp(-2*x)*heaviside(x)-exp(x)*heaviside(-x) ③ clear all; syms t w; ifourier(sin(w/4)/(w/4).*sin(w/4)/(w/4)); ans （