在MatlabSimulink构造一缓慢时变线性系统如下图所示.docVIP

在Matlab/Simulink构造一缓慢时变线性系统。试根据系统的输入生产数据分别用带遗忘因子最小二乘法和广义最小二乘法辨识系统的参数。 （1）带遗忘因子最小二乘法 t=[1:1:300]; w=sin(t); %干扰信号 x=idinput(300,prbs); %输入信号 a1=[1 0.5];b1=[0 0.4]; th0=arx2th(a1,b1,1,1); %线性系统1 y1=idsim([x w],th0); %系统1输出 a2=[1 0.7];b2=[0 0.5]; th0=arx2th(a2,b2,1,1); %线性系统1 y2=idsim([x w],th0); %系统1输出 for k=1:300 y(k,1)=((300-k)/300)*y1(k,1)+(k/300)*y2(k,1); end %时变系统输出 z=[y x]; na=1;nb=1;nk=1; [thm,yhat]=rarx(z,[na nb nk],ff,0.97); % 带遗忘因子最小二乘法 plot(thm(:,1),r-); grid; figure(2); plot(thm(:,2),:);grid; 结果如下： （2）广义最小二乘法 t=[1:1:300]; w=sin(t); %干扰信号 x=idinput(300,prbs); %输入信号 a1=[1 0.5];b1=[0 0.4]; th0=arx2th(a1,b1,1,1); %线性系统1 y1=idsim([x w],th0); %线性1输出 a2=[1 0.7];b2=[0 0.5]; th0=arx2th(a2,b2,1,1); %线性系统2 y2=idsim([x w],th0); %线性2输出 for k=1:300 y(k,1)=((300-k)/300)*y1(k,1)+(k/300)*y2(k,1); end %时变系统输出 z=[y x]; na=1;nb=1;nc=1;nk=1; [thm,yhat]=rarmax(z,[na nb nc nk],ff,0.97); %广义最小二乘法 plot(thm(:,1),-); grid; figure(2); plot(thm(:,2),r:);grid; 结果如下： 3－8.在Matlab/Simulink构造一缓慢时变线性系统如下图所示。试根据系统的输入生产数据分别用带遗忘因子最小二乘法和广义最小二乘法辨识系统的参数。 （1）带遗忘因子最小二乘法 t=[1:1:300]; w=sin(t); %干扰信号 x=idinput(300,prbs); %输入信号 a1=[1 0.5];b1=[0 0.4]; th0=arx2th(a1,b1,1,1); %线性系统1 y1=idsim([x w],th0); %系统1输出 a2=[1 0.7];b2=[0 0.5]; th0=arx2th(a2,b2,1,1); %线性系统1 y2=idsim([x w],th0); %系统1输出 for k=1:300 y(k,1)=((300-k)/300)*y1(k,1)+(k/300)*y2(k,1); end %时变系统输出 z=[y x]; na=1;nb=1;nk=1; [thm,yhat]=rarx(z,[na nb nk],ff,0.97); % 带遗忘因子最小二乘法 plot(thm(:,1),r-); grid; figure(2); plot(thm(:,2),:);grid; 结果如下： （2）广义最小二乘法 t=[1:1:300]; w=sin(t); %干扰信号 x=idinput(300,prbs); %输入信号 a1=[1 0.5];b1=[0 0.4]; th0=arx2th(a1,b1,1,1); %线性系统1 y1=idsim([x w],th0); %线性1输出 a2=[1 0.7];b2=[0 0.5]; th0=arx2th(a2,b2,1,1); %线性系统2 y2=idsim([x w],th0); %线性2输出 for k=1:300 y(k,1)=((300-k)/3