離散系統(tǒng)的時(shí)域分析實(shí)驗(yàn)報(bào)告參考模板_第1頁(yè)
離散系統(tǒng)的時(shí)域分析實(shí)驗(yàn)報(bào)告參考模板_第2頁(yè)
離散系統(tǒng)的時(shí)域分析實(shí)驗(yàn)報(bào)告參考模板_第3頁(yè)
離散系統(tǒng)的時(shí)域分析實(shí)驗(yàn)報(bào)告參考模板_第4頁(yè)
離散系統(tǒng)的時(shí)域分析實(shí)驗(yàn)報(bào)告參考模板_第5頁(yè)
已閱讀5頁(yè),還剩2頁(yè)未讀, 繼續(xù)免費(fèi)閱讀

下載本文檔

版權(quán)說(shuō)明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請(qǐng)進(jìn)行舉報(bào)或認(rèn)領(lǐng)

文檔簡(jiǎn)介

1、實(shí)驗(yàn)2 離散系統(tǒng)的時(shí)域分析一、實(shí)驗(yàn)?zāi)康?、熟悉并掌握離散系統(tǒng)的差分方程表示法;2、加深對(duì)沖激響應(yīng)和卷積分析方法的理解。 二、實(shí)驗(yàn)原理在時(shí)域中,離散時(shí)間系統(tǒng)對(duì)輸入信號(hào)或者延遲信號(hào)進(jìn)行運(yùn)算處理,生成具有所需特性的輸出信號(hào),具體框圖如下:其輸入、輸出關(guān)系可用以下差分方程描述:輸入信號(hào)分解為沖激信號(hào),記系統(tǒng)單位沖激響應(yīng),則系統(tǒng)響應(yīng)為如下的卷積計(jì)算式:當(dāng)時(shí),hn是有限長(zhǎng)度的(),稱系統(tǒng)為FIR系統(tǒng);反之,稱系統(tǒng)為IIR系統(tǒng)。三、實(shí)驗(yàn)內(nèi)容1 / 71、用MATLAB求系統(tǒng)響應(yīng)1) 卷積的實(shí)現(xiàn)線性移不變系統(tǒng)可由它的單位脈沖響應(yīng)來(lái)表征。若已知了單位脈沖響應(yīng)和系統(tǒng)激勵(lì)就可通過(guò)卷積運(yùn)算來(lái)求取系統(tǒng)響應(yīng),即程序:x=

2、input(Type in the input sequence=); %輸入xh=input(Type in the impulse response sequence=); %輸入hy=conv(x,h); % 對(duì)x,h進(jìn)行卷積N=length(y)-1; %求出N的值n=0:1:N; %n從0開(kāi)始,間隔為1的取值取到N為止disp(output sequence=); disp(y); %輸出ystem(n,y); %畫出n為橫軸,y為縱軸的離散圖xlabel(Time index n); ylable(Amplitude); % 規(guī)定x軸y軸的標(biāo)簽輸入為: x=-2 0 1 -1 3

3、h=1 2 0 -1 圖形:2) 單位脈沖響應(yīng)的求取線性時(shí)不變因果系統(tǒng)可用MATLAB的函數(shù)filter來(lái)仿真 y=filter(b,a,x); 其中,x和y是長(zhǎng)度相等的兩個(gè)矢量。矢量x表示激勵(lì),矢量a,b表示系統(tǒng)函數(shù)形式濾波器的分子和分母系數(shù),得到的響應(yīng)為矢量y。例如計(jì)算以下系統(tǒng)的單位脈沖響應(yīng)y(n)+0.7y(n-1)-0.45y(y-2)-0.6y(y-3)=0.8x(n)-0.44x(n-1)+0.36x(n-2)+0.02x(n-3) 程序: N=input(Desired impuse response length=);b=input(Type in the vector b=)

4、;a=input(Type in the vector a=);x=1 zeros(1,N-1);y=filter(b,a,x);k=0:1:N-1;stem(k,y);xlabel(Time index n); ylable(Amplitude);輸入:N=41b=0.8 -0.44 0.36 0.02a=1 0.7 -0.45 -0.6圖形: 2、以下程序中分別使用conv和filter函數(shù)計(jì)算h和x的卷積y和y1,運(yùn)行程序,并分析y和y1是否有差別,為什么要使用xn補(bǔ)零后的x1來(lái)產(chǎn)生y1;具體分析當(dāng)hn有i個(gè)值,xn有j個(gè)值,使用filter完成卷積功能,需要如何補(bǔ)零?程序:clf;h

5、= 3 2 1 -2 1 0 -4 0 3; %impulse responsex = 1 -2 3 -4 3 2 1; %input sequencey = conv(h,x);n = 0:14;subplot(2,1,1);stem(n,y);xlabel('Time index n'); ylabel('Amplitude');title('Output Obtained by Convolution'); grid;x1 = x zeros(1,8);y1 = filter(h,1,x1);subplot(2,1,2);stem(n,y1

6、);xlabel('Time index n'); ylabel('Amplitude');title('Output Generated by Filtering'); grid; 圖形: 因?yàn)樵趛=filter(b,a,x)中,利用給定矢量a和b對(duì)x中的數(shù)據(jù)進(jìn)行濾波,結(jié)果放入y矢量中,y與x長(zhǎng)度要相等,所以要使用xn補(bǔ)零后的x1來(lái)產(chǎn)生y1。若hn有i個(gè)值,xn有j個(gè)值,則x1 = x zeros(1,i-1)3、編制程序求解下列兩個(gè)系統(tǒng)的單位沖激響應(yīng),分別用filter 和 impz實(shí)現(xiàn),并繪出其圖形。給出理論計(jì)算結(jié)果和程序計(jì)算結(jié)果并討論。第

7、一題:filter實(shí)現(xiàn):程序:N=input('Desired impuse response length=');b=input('Type in the vector b=');a=input('Type in the vector a=');x=1 zeros(1,N-1);y=filter(b,a,x);k=0:1:N-1;stem(k,y);xlabel('Time index n');ylabel('Amplitude');圖形: impz實(shí)現(xiàn):程序:b=input('Type in the v

8、ector b=');a=input('Type in the vector a='); N=25y=impz(b,a,N);k=0:1:N-1;stem(k,y);xlabel('Time index n'); ylabel('Amplitude');圖形: 第二題:filter實(shí)現(xiàn):程序:N=input('Desired impuse response length=');b=input('Type in the vector b=');a=input('Type in the vector a

9、=');x=1 zeros(1,N-1);y=filter(b,a,x);k=0:1:N-1;stem(k,y);xlabel('Time index n');ylabel('Amplitude');輸入:Type in the vector b=0.25 0.25 0.25 0.25Type in the vector a=1N =25圖形: impz實(shí)現(xiàn):程序: b=input('Type in the vector b=');a=input('Type in the vector a='); N=30y=impz(b,a,N);k=0:1:N-1;stem(k,y);xlabel('Time index n'); ylabel('A

溫馨提示

  • 1. 本站所有資源如無(wú)特殊說(shuō)明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請(qǐng)下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請(qǐng)聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁(yè)內(nèi)容里面會(huì)有圖紙預(yù)覽,若沒(méi)有圖紙預(yù)覽就沒(méi)有圖紙。
  • 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
  • 5. 人人文庫(kù)網(wǎng)僅提供信息存儲(chǔ)空間,僅對(duì)用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對(duì)用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對(duì)任何下載內(nèi)容負(fù)責(zé)。
  • 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請(qǐng)與我們聯(lián)系,我們立即糾正。
  • 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶因使用這些下載資源對(duì)自己和他人造成任何形式的傷害或損失。

評(píng)論

0/150

提交評(píng)論