遺傳算法工具箱識別(GA)Bouc-Wen模型參數(shù)辨識-識別(共4頁)

1、精選優(yōu)質(zhì)文檔-傾情為你奉上Bouc-Wen模型因數(shù)字處理方便簡單而得到較為廣泛的應用，力可以表示為：                                                         

2、60;                                                        利用遺傳算法工具箱（GA）對Bouc-Wen模型進行參數(shù)識別。實驗數(shù)據(jù)來源于對磁流變阻尼器（MR damper）

3、進行性能測試，試驗獲得的數(shù)據(jù)包括力F，位移x，采用頻率已知，速度和加速度可以由位移求導得出。參數(shù)識別出現(xiàn)程序如下：（文件名：Copy_0_of_BoucWen）function j=myfung(x)y0=0;yy=y0;tspan='s='v='Ft='rr=max(size(s);%計算數(shù)據(jù)個數(shù)                i=1;               

4、;  while (i<rr)&(isnan(y0(1,1)&(max(abs(y0)>1e5)%判斷是否出現(xiàn)奇異點，具體忘了。                  t y=ode45(uubird,tspan(i),tspan(i+1),y0,v(i),x);%參考論壇的                   y

5、0=y(end,:);                   yy=yy;y0;                   i=i+1;                   kk=max(size(y);     &#

6、160;             if kk>150     %微分方程計算，停止是有條件的（具體沒去研究），這邊設(shè)置150次，不管有沒有收斂，都停止，不然整個程序運行的實際太久，你也可以改成其他的，慢慢研究                       break;     

7、0;             end                 end                 if  (i=rr)&(isnan(yy(1,1)=1%判斷是否出現(xiàn)奇異點（就是NAN），如果沒有出現(xiàn)，就是正常的   &

8、#160;                   F=x(:,4)*yy(:,1)+x(:,5)*(s-ones(size(s)*x(:,6)+x(:,7)*v;%x(:,4)代表alpha 5代表k0,6代表s0   7代表c0    位移s就是公式中的x                    &#

9、160;  j=sum(F-Ft).*(F-Ft);                       i=i+1;                                    

10、;else  i<(rr-1)%出現(xiàn)奇異點（）                       j=1e10;%因為出現(xiàn)奇異值，所以隨便給一個目標函數(shù)值（這個要足夠大），目的是排除這組優(yōu)化值（也就是這個種群不要）                       

11、i=rr;                             endfunction dy=uubird(t,y,v,x)dy=-x(:,1)*abs(v)*y*abs(y)-x(:,2)*(v)*y2+x(:,3)*v;% 1r 2beita 3A  取n=2運行結(jié)果擬合效果對比程序（文件名BoucWenjianyan247）：clccleardatasorce=loa

12、d('247-0.txt');II=1;a=2000;b=3000;                                          %為了擬合。datasorce=datasorce(a:b,:);%去一個周期進行對比s=datasorce(:,1)*0.7857-0.;%位移換算（

13、你不用管）Ft=datasorce(:,4)*1000/3.44;     %力換算sa=(max(s)-min(s)/2-max(s);s=s+sa;sFt=(max(Ft)-min(Ft)/2-max(Ft);Ft=Ft+sFt;for i=1:(b-a+1)    A(i)=i*0.001;endT=A'p=polyfit(T,s(:),10);%擬合，用于求解y10=polyval(p,T);syms xy= p(1)*x10 + p(2)*x9 +p(3)*x8 + p(4)*x7 + p(5)*x6 + p(6)*

14、x5 + p(7)*x4 + p(8)*x3 +p(9)*x2 + p(10)*x +p(11);f=diff(y);ac=diff(y,2);for i=1:(b-a+1)x=i*0.001;v(i)=eval(f);acc(i)=eval(ac);vm(i)=(v(i)*v(i)-s(i)*acc(i)0.5;I(i)=0;endv=v'T=T'y0=0;yy=y0;rr=max(size(s);tspan=0:0.001:(rr-1)*0.001;i=1;x=0.764        -0.7651 

15、      2.291        3.108        4.994        2.        3.002;%參數(shù)識別出出來的結(jié)果  while (i<rr)t y=ode45(uubird,tspan(i),tspan(i+1),y0,v(i),x);y0=y(end,:);yy=yy;y0;i=i+1;endF=x(:,4)*yy(:,1)+x(:,5)*(s-ones(size(s)*x(:,6)+x(:,7)*v;%4alpha 5k0 6x0 7c0s=s(100:2:800,:);v=