最小二乘法辨識(shí)參數(shù)

1、姓名：專(zhuān)偉學(xué)號(hào)：201221014368專(zhuān)業(yè)：控制工程增廣型遞推最小二乘法仿真第一種模型的仿貞程序:%選擇的模型結(jié)構(gòu)(遞推最小二乘法)：Y=Fai*theta+E%輸入U(xiǎn)為一個(gè)偽隨機(jī)序列L=10;yl=l;y2=l;y3=l;y4=0;for i=l:L;xl=xor(y3,y4);x2=yl;x3=y2;x4=y3;y(i)=y4;if y(i)>0.5.u(i)=-0.1;else u(i)=0.1;endyl=xl ;y2=x2 ;y 3=x3 ;y4=x4;endv=randn(l JO);y(l )=0;y(2)=0;%設(shè)定輸出的初始值y(O),y( 1)theta0=0.01

2、;0;0.01;0;0.01;%給出待辨識(shí)參數(shù)的初始值PO=10A6*eye(5);%生成初始矩陣7x7的單位陣E=0.00005;%E為遞推結(jié)束的條件%下而進(jìn)行遞推運(yùn)算for k=3:10y(k)= 1.5*y(k-l)-0.7*y(k-2)+1.0*u(k-l)+0.5*u(k-2)+v(k);%產(chǎn)生輸出 yF£l=-y(k-l)-y(k-2) u(k-l) u(k-2) v(k)'%生成實(shí)測(cè)數(shù)據(jù)組FailK l=P0*Fa 訂 /(1 +Fair*PO*Fail);theta 1 =thetaO+K 1 *(y(k)-Fa 訂 Fhc(aO);el=thetal-the

3、taO;%e2=el/theta 1;Pl=PO-Kl*Fair*P():thetaO=thetal;%供下次遞推使用PO=P1:%供下次遞推使用if abs(el)<E,break;% 循壞結(jié)束條件endendtheta 1辨識(shí)的參數(shù)的最后結(jié)果為：theta 1 =-1.50000.70001.00000.50001.0000可以從結(jié)果看岀，辨識(shí)的效果很好，基本沒(méi)有誤差。第二種模型的仿真程序:%選擇的模型結(jié)構(gòu)(增廣最小二乘法)：Y=Fai*theta+E%輸入U(xiǎn)為一個(gè)偽隨機(jī)序列L=10;yl=l;y2=l;y3=l;y4=0;for i=l:L;xl=xor(y3,y4);x2=yl;

4、x3=y2;x4=y3;y(i)=y4;if y(i)>0.5.u(i)=-0.1;else u(i)=0.1;endyl=xl ;y2=x2;y3=x3 ;y4=x4;end v=randn(L10);%噪聲干擾隨機(jī)數(shù)序列y( 1 )=0;y(2)=0;%設(shè)定輸出的初始值y(O),y( 1)theta0=0.01 ;0;0.01 ;0;0.01 ;0;0;%給出待辨識(shí)參數(shù)的初始值PO=10A6*eye(7);%生成初始矩陣7x7的單位陣E=0.00005;%E為遞推結(jié)束的條件%下而進(jìn)行遞推運(yùn)算for k=3:10y(k)= 1.5*y(k-l )-0.7*y(k-2)+ 1.0*u(k

5、-l )+0.5*u(k-2)+v(k)-v(k-1 )+0.2*v(k-2);%產(chǎn)生輸出 yF£l=-y(k-l)-y(k-2) u(k-l)u(k-2) v(k) v(k-l) v(k-2)'；%生成實(shí)測(cè)數(shù)據(jù)組FailK l=P0*Fa il/(l +Fair*PO*Fail);theta 1 =thetaO+K 1 *(y(k)-Fa 訂 FhctaO);el=thetal-thetaO;%e2=e 1/theta 1;Pl=PO-Kl*Fair*PO;thetaO=theta 1;%供下次遞推使用PO=P1;%供下次遞推使用if abs(el )<E,break;% 循環(huán)結(jié)束條件endendtheta 1辨識(shí)的參