連續(xù)和離散系統(tǒng)分析

1、實(shí)驗(yàn)一連續(xù)和離散系統(tǒng)分析一、實(shí)驗(yàn)?zāi)康膶W(xué)習(xí)連續(xù)系統(tǒng)和離散系統(tǒng)響應(yīng)的matlab 求解方法；二、實(shí)驗(yàn)主要儀器設(shè)備和材料計(jì)算機(jī)三、實(shí)驗(yàn)方法、步驟及結(jié)果測(cè)試實(shí)驗(yàn)方法：編程，上機(jī)調(diào)試，分析實(shí)驗(yàn)結(jié)果；步驟：編程實(shí)現(xiàn)上述各實(shí)驗(yàn)容四、實(shí)驗(yàn)結(jié)果1、某系統(tǒng)的傳遞函數(shù)為：1Y(s) / X ( x)(s1)( s2)試求系統(tǒng)的沖激響應(yīng)和階躍響應(yīng)。2、編制程序求解下列兩個(gè)系統(tǒng)的單位沖激響應(yīng)和階躍響應(yīng)，并繪出其圖形。要求分別用 filter、conv、impz 三種函數(shù)完成。給出理論計(jì)算結(jié)果和程序計(jì)算結(jié)果并討論。（I ） y n0.75yn1 0.125yn 2xn x n1理論計(jì)算結(jié)果：?jiǎn)挝粵_激響應(yīng) :n012345

2、h(n)1-1.751.19-0.670.355-0.18單位階躍響應(yīng) :n012345y(n)1-0.750.44-0.2340.12-0.06程序計(jì)算結(jié)果：A：?jiǎn)挝粵_激響應(yīng)（ 1）用 Filter函數(shù)（ 2）用Conv函數(shù)（ 3）用 impz 函數(shù)B：?jiǎn)挝浑A躍響應(yīng)(1)用 Filter 函數(shù)（ 2）用 Conv 函數(shù)（ 3）用 Impz 函數(shù)（II ） y n0.25 xn1x n2x n3xn4理論計(jì)算結(jié)果：?jiǎn)挝粵_激響應(yīng):n012345h(n)00.250.250.250.250單位階躍響應(yīng) :N012345y(n)00.250.50.7511程序計(jì)算結(jié)果：A ：?jiǎn)挝粵_激響應(yīng)（1）用 f

3、ilter 函數(shù)（ 2）用 Conv 函數(shù)（3）用 Impz 函數(shù)B：?jiǎn)挝浑A躍響應(yīng)（1）用 filter 函數(shù)（ 2）用 Conv 函數(shù)（ 3）用 Impz 函數(shù)結(jié)論：y=filter(p,d,x) 用來(lái)實(shí)現(xiàn)差分方程， d 表示差分方程輸出 y 的系數(shù)， p 表示輸入 x 的系數(shù)，而 x 表示輸入序列。輸出結(jié)果長(zhǎng)度數(shù)等于 x 的長(zhǎng)度；而 y=conv(x,h)是用來(lái)實(shí)現(xiàn)卷積的，對(duì)x 序列和h 序列進(jìn)行卷積，輸出的結(jié)果個(gè)數(shù)等于x 的長(zhǎng)度與h 的長(zhǎng)度之和減去1；y=impz(p,d,N)是用來(lái)實(shí)現(xiàn)沖擊響應(yīng)的，d 和p 的定義同filter ， N表示沖擊響應(yīng)輸出的序列個(gè)數(shù)。3 已知描述某連續(xù)系統(tǒng)的

4、微分方程為y'' (t )5 y ' (t )6 y(t )2x ' (t)8x(t)求該系統(tǒng)的單位沖激響應(yīng)。程序：b=2 8;a=1 5 6;sys=tf(b,a);t=0:0.1:10;y=impulse(sys,t);plot(t,y);xlabel( 時(shí)間 (t) );ylabel( y(t)單位沖);tltle(激響應(yīng) );4 計(jì)算上述系統(tǒng)在輸入為x(t )e t u(t) 時(shí)的零狀態(tài)響應(yīng)。程序：ts=0;te=10;dt=0.01;sys=tf(2 8,1 5 6);t=ts:dt:te;f=exp(-t);y=lsim(sys,f,t);plot(

5、t,y)xlabel('time(sec)')ylabel('y(t)')12z 15已知系統(tǒng)函數(shù) H ( z)10.4z 1 0.12 z 2 ，求1）離散系統(tǒng)的單位沖激響應(yīng)h(n);2) 輸入為 x(n) u(n) ，求系統(tǒng)的零狀態(tài)響應(yīng)。程序：num=1 2;den=1 0.4 0.12;% system modely1=impz(num,den,20);subplot(121); stem(y1); title('impulse response');n=0:20;x2=ones(1,21);y1filter=filter(num,den,

6、x2);subplot(122);stem(n,y1filter);title('filter_step');xlabel('x');ylabel('y');6 已知某離散因果系統(tǒng)的系統(tǒng)函數(shù) H ( z)1z 12 ，試分析該系統(tǒng)的幅頻特z10.5z1性。程序：num=1 1; den=1 -1 0.5;sys=tf(num,den);w=-10:0.001:10;H=freqs(num,den,w);subplot(2,1,1);plot(w,abs(H);title('幅頻特性 ');subplot(2,1,2);plot(w

7、,angle(H);title(' 相頻特性 ' );五、思考題1）連續(xù)系統(tǒng)響應(yīng)的計(jì)算機(jī)求解可以分為哪些方法？各是什么原理？答：連續(xù)時(shí)間系統(tǒng)零狀態(tài)響應(yīng)的求解ylsim ( sys , f , t )t 表示計(jì)算系統(tǒng)響應(yīng)的抽樣點(diǎn)向量， f 是系統(tǒng)輸入信號(hào)向量 sys 是 LTI 系統(tǒng)模型，借助 tf 函數(shù)獲得。連續(xù)系統(tǒng)沖激響應(yīng)和階躍響應(yīng)的求解連續(xù)時(shí)間系統(tǒng)沖激響應(yīng)可用impulse 函數(shù)直接求出，其調(diào)用形式為：y=impulse(sys,t)連續(xù)時(shí)間系統(tǒng)階躍響應(yīng)可用step 函數(shù)直接求出，其調(diào)用形式為：y=step(sys,t)t 表示計(jì)算系統(tǒng)響應(yīng)的抽樣點(diǎn)向量，sys 是 LTI 系統(tǒng)模型 .2） matlab 中用于離散系統(tǒng)求解的命令有哪些？各基于什么求解方法？答：y=filter(p,d,x) 用來(lái)實(shí)現(xiàn)差分方程， d 表示差分方程輸出 y 的系數(shù)， p 表示輸入 x 的系數(shù)，而 x 表示輸入序列。輸出結(jié)果長(zhǎng)度數(shù)等于 x 的長(zhǎng)度；y=co