信號與系統(tǒng)實(shí)驗(yàn)三new

實(shí)驗(yàn)內(nèi)容：驗(yàn)證性實(shí)驗(yàn)：編寫實(shí)現(xiàn)信號的傅里葉變換和傅里葉逆變換傅里葉變換：已知連續(xù)時(shí)間信號f=exp((-2)*abs(t))，通過程序完成f（t）的傅里葉變換。Matlab程序：clc;clearall;symst;f=fourier(exp((-2)*abs(t)));ezplot(f);信號的傅里葉變換如圖所示。試畫出f=2/3*exp(-3*t)*sym('Heaviside(t)');的波形及其幅頻特性曲線。Matlab程序：clearall;clc;symstvwff=2/3*exp(-3*t)*sym('Heaviside(t)');F=fourier(f);subplot(2,1,1);ezplot(f);subplot(2,1,2);ezplot(abs(F));信號的幅頻特性曲線如圖所示。傅里葉逆變換：已知f(jw)=1/(1+w^2)，求信號的傅里葉逆變換：Matlab程序：symstw;ifourier(1/(1+w^2),t)結(jié)果：ans=((pi*heaviside(t))/exp(t)+pi*heaviside(-t)*exp(t))/(2*pi)傅里葉變換的數(shù)值計(jì)算：已知門函數(shù)，試采用數(shù)值計(jì)算的方法確定信號的傅里葉變換F（jw）。Matlab程序：clc;clearall;R=0.02;t=-2:R:2;f=stepfun(t,-1)-stepfun(t,1);W1=2*pi*5;N=500;k=0:N;W=k*W1/N;F=f*exp(-j*t'*W)*R;F=real(F);W=[-fliplr(W),W(2:501)];F=[fliplr(F),F(2:501)];subplot(2,1,1);plot(t,f);xlabel('t');ylabel('f(t)');axis([-2,2,-0.5,2]);title('f(t)=U(t+1)-U(t-1)');subplot(2,1,2);plot(W,F);xlabel('w');ylabel('F(W)');title('f（t）的傅里葉變換');信號的傅里葉變換如圖所示。(4)連續(xù)函數(shù)的傅立葉變換Matlab程序：clf;symsw;dt=2*pi/8;w=linspace(-2*pi,2*pi,2000)/dt;k=-2:2;f=ones(-1,5);F=f*exp(-j*k'*w)fl=abs(F);plot(f,F);grid;連續(xù)函數(shù)的傅立葉變換如圖所示。程序設(shè)計(jì)實(shí)驗(yàn)：（1）是確定下列信號的傅里葉變換的數(shù)學(xué)表達(dá)式：（a）symst;f=fourier(sym('Heaviside(t+1)')-sym('Heaviside(t-1)'))結(jié)果：f=-(1/exp(w*i))*transform::fourier(Heaviside(t),t,-w)+exp(w*i)*transform::fourier(Heaviside(t),t,-w)（b）symst;f=fourier(exp(-3*t)*sym('Heaviside(t)'))結(jié)果：f=transform::fourier(Heaviside(t)/exp(3*t),t,-w)（c）symst;f=fourier(exp(-t)*sym('Heaviside(t)'))結(jié)果：f=transform::fourier(Heaviside(t)/exp(t),t,-w)2）試畫出信號f（t）=exp(-3*t)*U（t）及信號f（t-4），f（t）*exp(-j*4*t)的頻譜圖：clc;clearall;r=0.01;t=-10:r:10;N=200;Y=2*pi;k=-N:N;w=k*Y/N;f1=exp(-3*t).*stepfun(t,0);F1=r*f1*exp(-j*t'*w);subplot(2,3,1);plot(w,abs(F1));grid;xlabel('w');ylabel('F(jw)');subplot(2,3,4);plot(w,angle(F1));grid;xlabel('w');ylabel('?à??');f2=exp(-3*(t-4)).*stepfun(t,4);F2=r*f2*exp(-j*t'*w);subplot(2,3,2);plot(w,abs(F2));grid;xlabel('w');ylabel('F(jw)');subplot(2,3,5);plot(w,angle(F2));grid;xlabel('w');ylabel('?à??');f3=f1.*exp(-j*4*t);F3=r*f3*exp(-j*t'*w);subplot(2,3,3);plot(w,ab