搜索
Bouc-Wen模型因数字处理方便简单而得到较为广泛的应用，力可以表示为：利用遗传算法工具箱（GA）对Bouc-Wen模型进行参数识别。实验数据来源于对磁流变阻尼器（MRdamper）进行性能测试，试验获得的数据包括力F，位移x，采用频率已知，速度和加速度可以由位移求导得出。参数识别出现程序如下：（文件名：Copy_0_of_BoucWen）functionj=myfung(x)y0=[0];yy=y0;tspan=[]';s=[]';v=[]';Ft=[]';rr=max(size(s));%计算数据个数i=1;while(irr)&(~isnan(y0(1,1)))&(~(max(abs(y0))1e5))%%判断是否出现奇异点，具体忘了。。[ty]=ode45(@uubird,[tspan(i),tspan(i+1)],y0,[],v(i),x);%参考论坛的y0=y(end,:);yy=[yy;y0];i=i+1;kk=max(size(y));ifkk150%微分方程计算，停止是有条件的（具体没去研究），这边设置150次，不管有没有收敛，都停止，不然整个程序运行的实际太久，你也可以改成其他的，慢慢研究break;endendif(i==rr)&(~isnan(yy(1,1)))==1%判断是否出现奇异点（就是NAN），如果没有出现，就是正常的F=x(:,4)*yy(:,1)+x(:,5)*(s-ones(size(s))*x(:,6))+x(:,7)*v;%x(:,4)代表alpha5代表k0,6代表s07代表c0位移s就是公式中的xj=sum((F-Ft).*(F-Ft));i=i+1;elsei(rr-1)%出现奇异点（ＮＡＮ）j=1e10;%因为出现奇异值，所以随便给一个目标函数值（这个要足够大），目的是排除这组优化值（也就是这个种群不要）i=rr;endfunctiondy=uubird(t,y,v,x)dy=-x(:,1)*abs(v)*y*abs(y)-x(:,2)*(v)*y^2+x(:,3)*v;%1r2beita3A取n=2运行结果拟合效果对比程序（文件名BoucWenjianyan247）：clccleardatasorce=load('247-0.txt');II=1;a=2000;b=3000;%为了拟合。。datasorce=datasorce(a:b,:);%去一个周期进行对比s=datasorce(:,1)*0.7857-0.027860;%位移换算（你不用管）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;fori=1:(b-a+1)A(i)=i*0.001;endT=A';p=polyfit(T,s(:),10);%拟合，用于求解y10=polyval(p,T);symsxy=p(1)*x^10+p(2)*x^9+p(3)*x^8+p(4)*x^7+p(5)*x^6+p(6)*x^5+p(7)*x^4+p(8)*x^3+p(9)*x^2+p(10)*x+p(11);f=diff(y);ac=diff(y,2);fori=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.76512.2913.1084.9942.1639253.002];%%参数识别出出来的结果while(irr)[ty]=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;%4alpha5k06x07c0s=s(100:2:800,:);v=v(100:2:800,:);F=F(100:2:800,:);Ft=Ft(100:2:800,:);figure(1),plot(v,Ft,'k',v,F,'r--'),holdonfigure(2),plot(s,Ft,'k',s,F,'r--'),holdon实线是试验测得的结果。[attach]188480[/attach][attach]188480[/attach][attach]188480[/attach]如果觉得效果不够好，可以增加迭代次数。