系统辨识相关分析法脉冲响应模型参数-system identify.m

系统辨识相关分析法脉冲响应模型参数-system identify.m
摘要：这是一个应用举例，A=【1 -1.5  0.7】，B=【1 0.5】，nk=1  ， e为噪声项，分为三种情况，在这三种情况下，采用伪随机信号作为输入，应用相关分析法辨识系统的脉冲响应函数，并辨识出系统的模型参数。
关键词：系统辨识，相关分析法，脉冲响应函数，最小二乘

%---------------------------理想系统的阶跃响应---------------------------------
A  =  [1  -1.5   0.7];
B  =  [0   1    0.5];
Model  =  idpoly（AB）;       %理想系统模型
Step（Model [0 100]）;        %理想系统模型的阶越响应
Grid;

%---------------------------相关法辨识脉冲函数---------------------------------
%获得辨识所需数据
%1、白噪声
N=1000;          		%设置实验长度
A = [1  -1.5  0.7];	%建立系统模型
B= [0 1 0.5];
Model= idpoly（AB）;       %理想系统模型
Model_white_error= idpoly（AB1）;      %模型中加入白噪声
U= iddata（[]idinput（N‘prbs‘））;	%设置输入信号
E= iddata（[]idinput（N‘rgs‘））;		%设置噪声信号为白噪声
Y= sim（Model_white_error[U E]）;  	%获得输出数据
%2、有色噪声
C = [1 0.5];
D = [1 0.5];
Model_1_error= idpoly（ABC）;      %模型中的噪声为有色噪声1
Model_2_error= idpoly（AB1D）;    %模型中的噪声为有色噪声2
Y1 = sim（Model_1_error[U E]）;  		%获得输出数据
Y2 = sim（Model_2_error[U E]）;