【正文】 % 子圖4：作出序列圖f[k+2]subplot(2,3,4)。filled39。 axis([7,7, 0, 20])。f=f[k+2]39。xlabel(39。) 。f39。grid on % 子圖5：作出序列圖f[k2]*u[k2]subplot(2,3,5)。filled39。 axis([7,7, 0, 20])。f=f[k2]*u[k2]39。xlabel(39。) 。f39。grid on7（2）執(zhí)行結(jié)果如上圖，繪出了指定的序列波形若，計(jì)算。f1=[1 1 1 1 0 0] 。 %求卷積和k=0:length(f)1。filled39。 axis([1,12, 0, 6])。f=x[k]*x[k]39。xlabel(39。) 。f39。grid on 8（2）執(zhí)行結(jié)果如上圖，作出了離散序列卷積的波形圖實(shí)驗(yàn)三 連續(xù)時(shí)間系統(tǒng)時(shí)域分析一、實(shí)驗(yàn)?zāi)康恼莆站矸e計(jì)算方法。運(yùn)用課堂上學(xué)到的理論知識(shí)，從RC、RL一階電路的響應(yīng)中正確區(qū)分零輸入響應(yīng)、零狀態(tài)響應(yīng)、自由響應(yīng)與受迫響應(yīng)。其中L=1，R=2。te=5。t=ts:dt:te。a=[1 2]。y=lsim(b,a,f,t)。 %畫(huà)圖xlabel(39。) %橫坐標(biāo)定義ylabel(39。) %縱坐標(biāo)定義title(39。)。syms t xa=[1 2]。e=exp(x)。 %手算系統(tǒng)的單位沖激響應(yīng)i=int(e.*h,x,0,t)。 %符號(hào)函數(shù)作圖xlabel(39。)。y(t)39。title(39。)。（1）源程序%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 求上述系統(tǒng)的沖激響應(yīng) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%ts=0。dt=。a=[1 2]。impulse(b,a,5)xlabel(39。)。h(t)39。title(39。)。te=5。t=ts:dt:te。b=[1]。Time(sec)39。ylabel(39。)。系統(tǒng)的階躍響應(yīng)y(t)39。grid on（2）執(zhí)行結(jié)果。二、實(shí)驗(yàn)內(nèi)容假設(shè)某系統(tǒng)的單位函數(shù)響應(yīng)，系統(tǒng)激勵(lì)信號(hào)，求系統(tǒng)的零狀態(tài)響應(yīng)，并畫(huà)圖表示。 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% close。 k=0:3。e=[1 1 1 1]。 t=0:6。filled39。 title(39。)。采用P61 計(jì)算系統(tǒng)響應(yīng)：畫(huà)出時(shí)的系統(tǒng)響應(yīng)(1)源程序12%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 采用P61 計(jì)算系統(tǒng)響應(yīng)：% 畫(huà)出 時(shí)的系統(tǒng)響應(yīng)%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% T4_2a=[ ]。y0=[1 0]。n=0:10。% recur function beginN=length(a)。y=[y0 zeros(1,length(n))]。a1=a(length(a):1:1)。for i=N+1:N+length(n) y(i)= a1*y(iN:i1)39。endy=y(N+1:N+length(n))。filled39。xlabel(39。)。y[n]39。(2)執(zhí)行結(jié)果(3)分析:，利用遞歸函數(shù)計(jì)算并作圖，完成了10階差分方程的求解，得到系統(tǒng)響應(yīng)。 掌握系統(tǒng)無(wú)失真?zhèn)鬏數(shù)幕緱l件。圖10（1）源程序：13%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 求如圖10所示周期矩形脈沖信號(hào)%的Fourier級(jí)數(shù)表達(dá)式，畫(huà)出頻譜圖，%并用前N次諧波合成的信號(hào)近似。T=2。 %函數(shù)表達(dá)式為 f(t)=[u(t0)u(t1)][u(t1)u(t2)]，兩個(gè)門t2=linspace(1,2,101)。 %采用9次諧波擬合原信號(hào)gt1=exp(i*pi*n39。 %被積函數(shù)gt2=exp(i*pi*n39。Fn=*trapz(gt1,2)*t1(2)*trapz(gt2,2)*t2(2) %求出傅里葉級(jí)數(shù)的系數(shù)（復(fù)指數(shù)形式） %梯形積分法，積分運(yùn)算近似成求和計(jì)算,trapz(被積函數(shù)，維數(shù)) %代表每一點(diǎn)t的函數(shù)值乘以delta t的疊加，這里gt1為10行100列的二維數(shù)組， %deltat用t1(2)表示 ，t1數(shù)組的第二個(gè)元素，即為(10)/100=sig=abs(Fn)。figure。stem(n,abs(Fn))。n39。ylabel(39。)。振幅譜39。grid on。stem(n,angle(Fn))。n39。ylabel(39。)。相位譜39。grid on。%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%clear。y=0Mark=input(39。) %輸入疊加的諧波次數(shù)for k=1:Mark y=y+2/(k*pi)*(1cos(k*pi))*sin(k*pi*t)。plot(t,y)。t39。ylabel(39。)。14（2）執(zhí)行結(jié)果：（3）分析： ，操作窗口求出了周期矩形脈沖信號(hào)的Fourier級(jí)數(shù)（復(fù)指數(shù)形式）的系數(shù)，（其中n=9，9次諧波）（復(fù)指數(shù)形式）系數(shù)的振幅譜和相位譜圖像，（其中n=9，9次諧波）。 Fw=fourier(t*exp(3*t)*heaviside(t)) %求符號(hào)函數(shù)的Fourier變換FFw=abs(Fw)。 %符號(hào)函數(shù)作圖xlabel(39。)。|F(w)|39。grid on %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% (2)求下列信號(hào)的傅里葉變換 ,并作圖 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% syms t。 %符號(hào)函數(shù)可以表示成2*u(t)1FFw=abs(Fw)。 %符號(hào)函數(shù)作圖函數(shù)xlabel(39。)。|F(w)|39。grid （2）執(zhí)行結(jié)果：（3）分析：。設(shè)載波信號(hào)的頻率為100Hz。w1=500*pi:*pi:500*pi。 %調(diào)制信號(hào)(原信號(hào)) Sa(t) x1=subs(x,t,t1)。 %畫(huà)原信號(hào)的時(shí)域圖像plot(t1,x1)。t/s39。ylabel(39。)。原信號(hào)波形 x=sinc(t)39。 Fx=fourier(x)。ezplot(Fx,[2*pi,2*pi])。w (rad/s)39。ylabel(39。)。原信號(hào)頻譜 39。 f=cos(2*pi*100*t)。 % 調(diào)制y1=subs(y,t,t1)。 % 畫(huà)調(diào)制后信號(hào)的時(shí)域圖像plot(t1,y1)。t/s39。ylabel(39。)。 調(diào)制后的信號(hào) y(t)=x*cos(wt) 39。 Fy=fourier(y)。subplot(4,2,4)。plot(w1,FFya)。w (rad/s)39。ylabel(39。)。 調(diào)制后的信號(hào)頻譜 39。 f1=y*cos(2*pi*100*t)。 %化為頻域 ff1=subs(f1,t,t1)。 plot(t1,ff1)。t/s39。ylabel(39。)。解調(diào)信號(hào) f1=y*cos(2*pi*100*t)39。 Ff1a=abs(Ff1)。FFf1=subs(Ff1a,w,w1)。 xlabel(39。)。FFf139。title(39。)。 %低通濾波器，F(xiàn)f2 = Ff1 * Fd。 %低通后取出的信號(hào)，時(shí)域subplot(4,2,7)。 %低通后信號(hào)的時(shí)域表示plot(t1,ff2)。t/s39。ylabel(39。)。低通取出的信號(hào)波形 39。 subplot(4,2,8)。 %低通后信號(hào)的頻譜xlabel(39。)。Ff239。title(39。)。（子圖4）。然后通過(guò)一個(gè)幅值為 2 截止頻率為 ws（wm ws w0wm）的低通濾波器后取出原信號(hào)頻譜，再轉(zhuǎn)化成時(shí)域信號(hào)，即完成了信號(hào)重建。 掌握系統(tǒng)無(wú)失真?zhèn)鬏數(shù)幕緱l件。三、實(shí)驗(yàn)內(nèi)容如圖11所示系統(tǒng)：圖10 (a) 對(duì)不同的RC值，用freqs函數(shù)畫(huà)出系統(tǒng)的幅頻曲線。提示：為最大值的處對(duì)應(yīng)的頻率為通帶截止頻率，首先求取并找到和RC關(guān)系，然后根據(jù)題意選定即可確定RC值。for RC=1:1:5。 b=[0 , 1]。hold on。 %輸入信號(hào)Fx=simplify(fourier(fx))。 %低通濾波 Fy=Fx*H。 fx1=subs(fx,t,t1)。 subplot(2,1,1)。xlabel(39。)。fx139。title(39。)。plot(t1,fy1)。t/s39。ylabel(39。)。 輸出信號(hào) y(t) 39。19(2) 執(zhí)行結(jié)果:（3）分析：，表示出了RC=1，2,3,4,5，時(shí)的幅頻曲線和相頻曲線，可以觀察到RC越大，曲線開(kāi)始下降的頻率越低，即RC與截止頻率成反比。信號(hào)的高頻分量被大程度衰減，低頻分量保留。提示：利用fourier求取信號(hào)的傅立葉變換E(jω)，然后設(shè)計(jì)使之滿足不同條件，計(jì)算R( jω)= E(jω) H(jω)并畫(huà)頻譜圖。 close 。e=exp(2*abs(t))。 figure(1)。 Fe=fourier(e,v)。 ezplot(Fe,[3,3])title(39。)。 i=1。 C(i)=angle(R11)。end b=3::3。相位譜39。 axis([3,3,1,1])。 %系統(tǒng)函數(shù)(1)