差分方程求解

1、差分方程求解差分方程求解實驗步驟：主界面下進入實驗五的“差分方程求解”子系統(tǒng)輸入希望看到的輸出樣點數(shù)1)輸入差分方程系數(shù)向量輸入順序為：，。其中N1為差分方程兩邊)(10Naaa?)(10Nbbb?系數(shù)最大數(shù)目，如果有一邊輸入系數(shù)個數(shù)小于N1，將按不足系數(shù)為零計算。2)輸入系統(tǒng)初始狀態(tài)向量輸入順序為)1()1()0(???Nyyyzizizi?鼠標(biāo)單擊確定按鈕，以數(shù)值和圖形兩種方式顯示單位沖激響應(yīng)和輸出響應(yīng)3）確定差分方程形式：y(k)

2、y(k1)0.35y(k2)=2x(k)x(k1)，4）實驗程序:%y(k)y(k1)0.35y(k2)=2x(k)x(k1)N=2%差分方程階數(shù)a=[110.35]%差分方程系數(shù)向量a(0)~a(N)b=[21]%差分方程系數(shù)向量b(0)~b(N)k=20%輸出樣點數(shù)目zi=[21]%初始狀態(tài)N個yzi=[0ones(1kN1)]h=yziyzs=yzifn=1:Nyzi(n)=zi(Nn1)endy=yzin=[N:k]x=imps

3、eq(0Nk)zic=filtic(bazi)h(N1:end)=filter(bax(N1:end))%或由命令dimpulse實現(xiàn)x=stepseq(0Nk)yzs(N1:end)=filter(bax(N1:end))yzi(N1:end)=filter([00]ax(N1:end)zic)y(N1:end)=filter(bax(N1:end)zic)figuresubplot(311)stem(nx)title(激勵x(n))