隨機(jī)過程實(shí)驗(yàn)四 線性系統(tǒng)參數(shù)估計(jì)及隨機(jī)過程預(yù)測

1、HUNANHUNANUNIVERSITYUNIVERSITY課程實(shí)驗(yàn)報(bào)告題目：線性系統(tǒng)參數(shù)估計(jì)及隨機(jī)過程預(yù)測u=nmrnd(0a1N)%nmrnd為正態(tài)分布隨機(jī)數(shù)產(chǎn)生函數(shù)figure(1)%畫出第一個圖stem([1:N]u.r)%用紅色畫出1到N范圍內(nèi)的隨機(jī)數(shù)分布圖gridon%畫出網(wǎng)格title(u(n)波形)%(2)畫出x(n)的波形圖x=zeros(11998)x(1)=u(1)x(2)=0.9x(1)u(2)fn=3:Nx(n

2、)=0.9x(n1)0.2x(n2)u(n)endfigure(2)stem([1:N]x.r)gridontitle(x(n)波形)%(3)以x(n)(n=34…1500)為已知數(shù)據(jù)，估計(jì)系統(tǒng)參數(shù)m1=0m2=0m3=0m4=0n1=0n2=0%賦初值fn=3:1500m1=m1x(n1)^2m2=m2x(n1)x(n2)m3=m3x(n1)x(n2)m4=m4x(n2)^2%第一個矩陣的四個表達(dá)式n1=n1x(n)x(n1)n2=n

3、2x(n)x(n2)%第二個矩陣的兩個表達(dá)式endA=[m1m2m3m4]A1=inv(A)%求A的逆矩陣B=[n1n2]C=A1B%將兩個矩陣相乘a=C(1)b=C(2)%將所求得值賦值給ab%(4)在x(n)的波形圖上用不同的顏色畫出y(n)的波形圖y=zeros(12000)n=[1501:2000]y(n)=ax(n1)bx(n2)figure(3)stem(x.r)gridonholdonstem(y.b)%在x(n)的波形圖