【正文】 %點擊鼠標，標注a2曲線gtext(39。)gtext(39。:r39。eb139。ed239。)subplot(414)i=1:30。\leftarrow39。\leftarrow39。g39。)gtext(39。) %點擊鼠標，標注左箭頭gtext(39。:39。b2=c(4,1:k)。 %產(chǎn)生不考慮噪聲時的模型輸出 ymd(k)=x139。 %設定誤差允許值c=[c0,zeros(7,29)]。)gtext(39。) %第三個圖形，繪制辨識的誤差title(39。)gtext(39。,i,a2,39。*p0。 %估計誤差矩陣賦初值for k=3:30 x1=[y(k1) y(k2) u(k1) u(k2)]39。c3=XL39。四、實驗原理 最小二乘一次完成算法 公式辨識參數(shù)上式中 程序流程圖圖 1最小二乘一次完成程序流程圖 遞推最小二乘算法 遞推公式公式為其中， 算法流程圖圖 2 遞推最小二乘法實現(xiàn)程序框圖 增廣最小二乘遞推算法 遞推公式公式為：其中， 算法流程圖圖 3 增廣最小二乘法算法流程圖五、實驗結(jié)果 最小二乘法一次完成實驗結(jié)果XL = 0 0 0 0 0 0 0 0 0 0 0 yL = 0 0 辨識參數(shù)矩陣c = a1 =；a2 =；b1 =；b2 =下圖為輸入、輸出矩陣的根徑圖圖 4最小二乘法一次實現(xiàn)輸入輸出根徑圖 遞推最小二乘法算法辨識結(jié)果系統(tǒng)輸出矩陣：y = Columns 1 through 13 0 0 Columns 14 through 26 Columns 27 through 30 辨識參數(shù)矩陣（辨識過程執(zhí)行26次即滿足了誤差要求）：c = Columns 1 through 13 0 0 0 0 Columns 14 through 26 辨識誤差矩陣：e = Columns 1 through 13 0 0 0 0 0 0 0 0 0 0 0 Columns 14 through 26 輸入輸出根徑圖圖 5 遞推最小二乘法輸入輸出根徑圖辨識參數(shù)過程圖 6 遞推最小二乘法辨識參數(shù)（辨識過程）辨識誤差：圖 7 遞推最小二乘法辨識過程中的誤差曲線 增光最小二乘法實驗結(jié)果隨機噪聲e0 = Columns 1 through 13 Columns 14 through 26 Columns 27 through 30 考慮噪聲的系統(tǒng)輸出矩陣y = Columns 1 through 13 0 0 Columns 14 through 26 Columns 27 through 30 不考慮噪聲的系統(tǒng)輸出矩陣ys = Columns 1 through 13 0 0 Columns 14 through 26 Columns 27 through 30 不考慮噪聲時的模型輸出ym = Columns 1 through 13 0 0 Columns 14 through 26 Columns 27 through 30 考慮噪聲時的模型輸出ymd = Columns 1 through 13 0 0 Columns 14 through 26 Columns 27 through 30 辨識參數(shù)矩陣：c = Columns 1 through 13 0 0 0 0 0 0 0 Columns 14 through 26