资源描述:
《连续LTI系统的频域分析》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库。
1、-.验证性试验1.傅里叶变换>>symst;>>f=fourier(exp(-2*abs(t)));>>ezplot(f);2.F(t)=2/3e-3tU(t)symstvwf;f二sym('(2/3咤xp(-3*t))*(heaviside(t))‘);>>F=fourier(f);>>subplot(2,l,l);>>ezplot(f);>>subplot(2丄2);>>ezplot(abs(F));3.傅里叶反变换>>symstw;>>ifourier(l/(l+wA2)/t)>>simplify(ans)ans=heaviside(t)/(2*exp(t))-(
2、exp(t)*(heaviside(t)-1))/24.傅里叶数值计算f(t)=U(t+l)-U(t-l)t=-2:0.02:2;>>f=stepfun(t,-1)-stepfun(t,l);>>wl=2*pi*5;>>n=500;>>k=0:n;>>w=k*wl/n;HFigure1V凹>F=f*exp(-j*t'*w)*0.02;>>F=real(F);>>w=[-fliplr(w)/w(2:501)];>>F=[-fliplr(F),F(2:5Ol)];>>subplot(2,l,l);>>plot(t,f);>>xlabel('t');>>ylabelCfft)
3、1);>>axis(卜22・0.5,2]);>>titleCf(t)=U(t+l)-U(t-l)');>>subplot(2,l,2);>>plot(w,F);>>xlabel('w');ylabel('F(w)');>>title('f(t)的傅里叶变换');1.连续函数的傅里叶变换>>dt=2*pi/8;>>w=linspace(-2*piz2*pi,2000)/dt;>>k=-2:2;>>f二ones(lz5);>>F=f*exp(-j*k'*w);>>fl=abs(F);plot(w,fl);>>grid;2.连续周期信号的傅里叶级数clearall;N二&nl
4、=-N:-l;p=nl.A2;cl=-4*j*sin(nl*pi/2)/(piA2)*(l./p);c0=0;n2=l:N;pl=n2.A2;c2=-4*j*sin(n2*pi/2)/(piA2)*(l./pl);cn=[clcOc2];n=-N:N;subplot(2/l/l);stem(n,abs(cn));ylabel('AmofCN');subplot(2/l,2);stem(n,angle(cn));ylabel('phaseofCN');xlabel('omegaAomegaO');1.F(t)=l/Ze-2tU(t)r=0.02;>>t=-5:r:5;N
5、=200;Y=2*pi;k=-N:N;w=k*Y/N;>>fl=l/2*exp(-2*t).*stepfun(t,0);>>F=r*fl*exp(-j*t**w);>>Fl=abs(F);Pl=angle(F);subplot(3,l,l);>>plot(t,fl);grid;>>xlabel('t');ylabel('f(t)');title('f(t)‘);>>subplot(3/l,2);plot(w,Fl);>>xlabel(W);ylabel('f(jw)');grid>>subplot(3,l,3);plot(w/Pl*180/pi);>>xlabel('w
6、');ylabel(相位(度)*);grid2.F(t・l)的频谱r=0.02;t=-5:r:5;N=200;Y=2*pi;k=-N:N;w=k*Y/N;fl=l/2*exp(-2*(t-l)).*stepfun(t,l);F=r*fl*exp(-j*t'*w);Fl=abs(F);Pl=angle(F);subplot(3/l/l);plot(t,fl);grid;xlabel('t');ylabel('f(t)');titlG('f(t)');subplot(3/l,2);plot(w/Fl);xlabelCw');ylabelCf(jw)的模');gridsubp
7、lot(3/l/3);plot(w,Pl*180/pi);xlabelCw');ylabel(相位(度)');grid3.傅里叶变换的频移特性R=0.02;t=-2:R:2;f=stepfun(t,-l)-stepfun(t,1);fl=f.*exp(-j*10*t);f2=f.*exp(j*10*t);Wl=10*pi;N=500;k=-N:N;W=k*Wl/N;Kt)Fl=fl*exp(-j*t**W)*R;F2=f2*exp(-j*t**W)*R;Fl=real(Fl);F2=real(F2);subplot(2/l,l);plot