MATLAB系統(tǒng)辨識(shí)仿真程序
5000多字 15頁(yè)


作業(yè)：用Matlab實(shí)現(xiàn)一個(gè)仿真例子的系統(tǒng)辨識(shí)，分別使用最小二乘整批算法、遞推算法和廣義最小二乘法實(shí)現(xiàn)，具有演示性。

要求：
1.對(duì)圖1、圖2分別用整批、遞推、廣義最小二乘算法實(shí)現(xiàn)，并對(duì)比這些算法對(duì)不同噪聲模型的辨識(shí)精度（λ取相同的值），并對(duì)結(jié)果給予合理的解釋；（兩種模型分別代表辨識(shí)模型類中殘差是白噪聲和有色噪聲的情況）
2.改變 值，比較辨識(shí)結(jié)果的精度，說明信噪比對(duì)各種算法辨識(shí)精度的影響，并解釋仿真結(jié)果；
3.顯示輸入輸出數(shù)據(jù)；
4.顯示遞推/迭代的收斂過程；
5.驗(yàn)證數(shù)據(jù)飽和現(xiàn)象（可能有兩種情況：1. 誤差不發(fā)散。找出臨界點(diǎn)，在圖中及其報(bào)告中加以說明。2. 誤差發(fā)散，要求如1）；
6.顯示預(yù)報(bào)誤差（數(shù)據(jù)分成兩段，前一部分用來辨識(shí)系統(tǒng)，后一部分用來檢驗(yàn)?zāi)Ｐ捅孀R(shí)結(jié)果）；
7.確定最佳階數(shù)，可以參考Matlab中的命令arxstruc、ivstruc、selstruc實(shí)現(xiàn)；
8.寫一份報(bào)告，對(duì)應(yīng)的結(jié)果和圖在報(bào)告中要有解釋說明；
9.7月9日前同時(shí)提交文字報(bào)告和相應(yīng)的程序，并在報(bào)告中簡(jiǎn)要說明程序的功能。
注：程序可以直接調(diào)用Matlab中的相應(yīng)命令實(shí)現(xiàn)，但提倡自己編程實(shí)現(xiàn)；要求具有良好的演示效果，在figure上信息要全。

統(tǒng)辨識(shí)MATLAB仿真程序，共分三個(gè)部分，數(shù)據(jù)產(chǎn)生，數(shù)據(jù)處理和數(shù)據(jù)顯示。數(shù)據(jù)產(chǎn)生，根據(jù)所選模型和參數(shù)，生成系統(tǒng)輸入、輸出數(shù)據(jù)；數(shù)據(jù)處理，根據(jù)觀測(cè)數(shù)據(jù)，進(jìn)行整批、遞推和廣義最小二乘法運(yùn)算，還可以預(yù)報(bào)誤差，定階處理。數(shù)據(jù)顯示，以圖表和文字方式顯示數(shù)據(jù)處理結(jié)果。
1，該仿真程序已在MATLAB 6.5.X上測(cè)試運(yùn)行，對(duì)其他版本兼容性未知；
2，每次在“輸入信號(hào)”框內(nèi)生成輸入序列后，以后系統(tǒng)的每步操作都是基于該輸入序列，除非重新點(diǎn)擊“畫圖”生成新的輸入序列；輸入序列的生成是隨機(jī)的；
3，輸入序列生成后必須在“輸出信號(hào)”框內(nèi)點(diǎn)擊“畫圖”按鈕生成輸出序列才可進(jìn)行系統(tǒng)下一步的辨識(shí)。輸入序列、信噪比、辨識(shí)模型改變后必須重新生成輸出序列。