【正文】 xlabel(39。k39。) 。ylabel(39。f39。) 。grid on % 子圖3：作出序列圖f[k]subplot(2,3,3)。stem(k,f,39。filled39。)。 axis([7,7, 0, 20])。 title(39。f=f[k]39。)。xlabel(39。k39。) 。ylabel(39。f39。) 。grid on % 子圖4：作出序列圖f[k+2]subplot(2,3,4)。stem(k+2,f,39。filled39。)。 axis([7,7, 0, 20])。 title(39。f=f[k+2]39。)。xlabel(39。k39。) 。ylabel(39。f39。) 。grid on % 子圖5：作出序列圖f[k2]*u[k2]subplot(2,3,5)。stem(k2,f.*stepfun(k2,0),39。filled39。)。 axis([7,7, 0, 20])。 title(39。f=f[k2]*u[k2]39。)。xlabel(39。k39。) 。ylabel(39。f39。) 。grid on7（2）執(zhí)行結(jié)果如上圖，繪出了指定的序列波形若，計算。（1）源程序%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 求離散序列的卷積%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% clear。f1=[1 1 1 1 0 0] 。 %序列x[k]f=conv(f1,f1)。 %求卷積和k=0:length(f)1。 %起點的和是起點，終點的和是終點[0,5]+[0,5]=[0,10]stem(k,f,39。filled39。)。 axis([1,12, 0, 6])。 title(39。f=x[k]*x[k]39。)。xlabel(39。k39。) 。ylabel(39。f39。) 。grid on 8（2）執(zhí)行結(jié)果如上圖，作出了離散序列卷積的波形圖實驗三 連續(xù)時間系統(tǒng)時域分析一、實驗?zāi)康恼莆站矸e計算方法。掌握函數(shù)lsim，impulse，step的用法，lsim為求取零狀態(tài)響應(yīng)，impulse為求取單位脈沖響應(yīng)，step為求取單位階躍響應(yīng)。運用課堂上學(xué)到的理論知識，從RC、RL一階電路的響應(yīng)中正確區(qū)分零輸入響應(yīng)、零狀態(tài)響應(yīng)、自由響應(yīng)與受迫響應(yīng)。二、實驗內(nèi)容1. 分別用函數(shù)lsim和卷積積分兩種方法求如圖8所示系統(tǒng)的零狀態(tài)響應(yīng)。其中L=1，R=2。圖8（1）源程序9%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 用函數(shù)lsim求系統(tǒng)的零狀態(tài)響應(yīng) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% ts=0。te=5。dt=。t=ts:dt:te。f=exp(t).*stepfun(t,0)。a=[1 2]。b=[1]。y=lsim(b,a,f,t)。%求取零狀態(tài)響應(yīng)plot(t,y)。 %畫圖xlabel(39。Time(sec)39。) %橫坐標(biāo)定義ylabel(39。y(t)39。) %縱坐標(biāo)定義title(39。Yzs(t) 39。)。grid on%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 用卷積積分求系統(tǒng)的零狀態(tài)響應(yīng) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%clear。syms t xa=[1 2]。b=[1]。e=exp(x)。 %激勵h=exp(2.*(tx))。 %手算系統(tǒng)的單位沖激響應(yīng)i=int(e.*h,x,0,t)。 %用卷積的定義式計算系統(tǒng)的零狀態(tài)響應(yīng)ezplot(i,[0 5])。 %符號函數(shù)作圖xlabel(39。Time(sec)39。)。ylabel(39。y(t)39。)。title(39。Yzs(t) 39。)。grid on（2）執(zhí)行結(jié)果10 是函數(shù)lsim 求出的系統(tǒng)零狀態(tài)響應(yīng)，2. 求上述系統(tǒng)的沖激響應(yīng)與階躍響應(yīng)。（1）源程序%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 求上述系統(tǒng)的沖激響應(yīng) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%ts=0。te=5。dt=。t=ts:dt:te。a=[1 2]。b=[1]。impulse(b,a,5)xlabel(39。Time(sec)39。)。ylabel(39。h(t)39。)。title(39。系統(tǒng)的沖激響應(yīng)h(t)39。)。grid on%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 求上述系統(tǒng)的階躍響應(yīng) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%ts=0。te=5。dt=。t=ts:dt:te。a=[1 2]。b=[1]。step(b,a,5)xlabel(39。Time(sec)39。)。ylabel(39。y(t)39。)。title(39。系統(tǒng)的階躍響應(yīng)y(t)39。)。grid on（2）執(zhí)行結(jié)果。三、思考題1. 為什么連續(xù)時間系統(tǒng)的零狀態(tài)響應(yīng)為激勵與沖擊響應(yīng)的卷積？2. 利用卷積積分法計算系統(tǒng)響應(yīng)應(yīng)從幾個方面進(jìn)行?實驗四 離散時間系統(tǒng)時域分析一、實驗?zāi)康恼莆沼肕ATLAB計算離散時間系統(tǒng)響應(yīng)的方法。二、實驗內(nèi)容假設(shè)某系統(tǒng)的單位函數(shù)響應(yīng)，系統(tǒng)激勵信號，求系統(tǒng)的零狀態(tài)響應(yīng)，并畫圖表示。（1）源程序 11%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 假設(shè)某系統(tǒng)的單位函數(shù)響應(yīng) ，% 系統(tǒng)激勵信號 ，求系統(tǒng)的零狀態(tài)響應(yīng)，% 并畫圖表示。 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% close。 clear。 k=0:3。h=().^k.*stepfun(k,0)。e=[1 1 1 1]。y=conv(h,e)。 t=0:6。 stem(t,y,39。filled39。)。 title(39。零狀態(tài)響應(yīng)39。)。（2）執(zhí)行結(jié)果 （3）分析：，利用卷積和求出了離散時間系統(tǒng)的零狀態(tài)響應(yīng)。采用P61 計算系統(tǒng)響應(yīng)：畫出時的系統(tǒng)響應(yīng)(1)源程序12%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 采用P61 計算系統(tǒng)響應(yīng)：% 畫出 時的系統(tǒng)響應(yīng)%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% T4_2a=[ ]。 b=[ 0 1 0]。y0=[1 0]。 x0=[0 0]。n=0:10。x=ones(1,length(n))。% recur function beginN=length(a)。M=length(b)1。y=[y0 zeros(1,length(n))]。x=[x0 x]。a1=a(length(a):1:1)。b1=b(length(b):1:1)。for i=N+1:N+length(n) y(i)= a1*y(iN:i1)39。 + b1*x(iN:iN+M)39。endy=y(N+1:N+length(n))。 % end of recur stem(n, y, 39。filled39。)。xlabel(39。n39。)。ylabel(39。y[n]39。)。(2)執(zhí)行結(jié)果(3)分析:，利用遞歸函數(shù)計算并作圖，完成了10階差分方程的求解，得到系統(tǒng)響應(yīng)。12實驗五 連續(xù)時間信號頻域分析一、實驗?zāi)康恼莆者B續(xù)時間系統(tǒng)變換域分析的基本方法。 掌握系統(tǒng)無失真?zhèn)鬏數(shù)幕緱l件。二、實驗內(nèi)容求如圖10所示周期矩形脈沖信號的Fourier級數(shù)表達(dá)式，畫出頻譜圖，并用前N次諧波合成的信號近似。圖10（1）源程序：13%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 求如圖10所示周期矩形脈沖信號%的Fourier級數(shù)表達(dá)式，畫出頻譜圖，%并用前N次諧波合成的信號近似。%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%clear。T=2。 %令周期T為2t1=linspace(0,1,101)。 %函數(shù)表達(dá)式為 f(t)=[u(t0)u(t1)][u(t1)u(t2)]，兩個門t2=linspace(1,2,101)。 %劃分為兩個門信號疊加n=0:9。 %采用9次諧波擬合原信號gt1=exp(i*pi*n39。*t1)。 %被積函數(shù)gt2=exp(i*pi*n39。*t2)。Fn=*trapz(gt1,2)*t1(2)*trapz(gt2,2)*t2(2) %求出傅里葉級數(shù)的系數(shù)（復(fù)指數(shù)形式） %梯形積分法，積分運算近似成求和計算,trapz(被積函數(shù)，維數(shù)) %代表每一點t的函數(shù)值乘以delta t的疊加，這里gt1為10行100列的二維數(shù)組， %deltat用t1(2)表示 ，t1數(shù)組的第二個元素，即為(10)/100=sig=abs(Fn)。 %由于數(shù)值計算會造成真實0值以極小值體現(xiàn)，因此去掉無意義的值Fn=Fn.*sig。figure。subplot(2,1,1)。stem(n,abs(Fn))。xlabel(39。n39。)。ylabel(39。|Fn|39。)。title(39。振幅譜39。)。grid on。 subplot(2,1,2)。stem(n,angle(Fn))。xlabel(39。n39。)。ylabel(39。Theta39。)。title(39。相位譜39。)。grid on。 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%用前N次諧波合成的信號近似。%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%clear。t=0::10。y=0M