电力变压器铁心柱截面的优化设计之程序实现(matlab)

《电力变压器铁心柱截面的优化设计之程序实现(matlab)》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、电力变压器铁心柱截面的优化设计之程序实现（MATLAB）附录7.1附录1funf='f=-1*(x(1)*sqrt(325^2-x(1)^2)+x(2)*(sqrt(325^2-x(2)^2)-sqrt(325^2-x(1)^2))+x(3)*(sqrt(325^2-x(3)^2)-sqrt(325^2-x(2)^2))+x(4)*(sqrt(325^2-x(4)^2)-sqrt(325^2-x(3)^2))+x(5)*(sqrt(325^2-x(5)^2)-sqrt(325^2-x(4)^2))+x(6)

2、*(sqrt(325^2-x(6)^2)-sqrt(325^2-x(5)^2))+x(7)*(sqrt(325^2-x(7)^2)-sqrt(325^2-x(6)^2))+x(8)*(sqrt(325^2-x(8)^2)-sqrt(325^2-x(7)^2))+x(9)*(sqrt(325^2-x(9)^2)-sqrt(325^2-x(8)^2))+x(10)*(sqrt(325^2-x(10)^2)-sqrt(325^2-x(9)^2))+x(11)*(sqrt(325^2-x(11)^2)-sqrt(3

3、25^2-x(10)^2))+x(12)*(sqrt(325^2-x(12)^2)-sqrt(325^2-x(11)^2))+x(13)*(sqrt(325^2-x(13)^2)-sqrt(325^2-x(12)^2))+x(14)*(sqrt(325^2-x(14)^2)-sqrt(325^2-x(13)^2)));';%最大面积的目标函数fung='g=[x(2)-x(1)+5;x(3)-x(2)+5;x(4)-x(3)+5;x(5)-x(4)+5;x(6)-x(5)+5;x(7)-x(6)+5;x(8

4、)-x(7)+5;x(9)-x(8)+5;x(10)-x(9)+5;x(11)-x(10)+5;x(12)-x(11)+5;x(13)-x(12)+5;x(14)-x(13)+5];';%宽度逐级递减的约束条件fun=[funffung];x0=[180185190195200205210215225225230235240245];%初始值options=[];vlb=[8580757065605550454035302520];%下界vub=[395390385380375370365360355350

5、345340335330];%上界[x,options]=constr(fun,x0,options,vlb,vub);y=zeros(1,14);x=x.*2;x=(round(x./10)).*10;%以10为倍数的宽度fori=1:14yy=sum(y',1);y(i)=sqrt(325^2-x(i)^2)-yy;endy=round(y);xys=x.*y;sum(s',1)(ans)/((325^2)*pi)7.2附录2max=0;hh1=395;yy1=(325^2-hh1^2)^(1/2);s

6、s1=hh1*yy1;forb=1:6hh2=385+(b-3)*5;yy2=(325^2-(hh2)^2)^(1/2)-yy1;ss2=hh2*yy2;forc=1:6hh3=375+(c-3)*5;yy3=(325^2-(hh3)^2)^(1/2)-yy1-yy2;ss3=hh3*yy3;ford=1:6hh4=365+(d-3)*5;yy4=(325^2-(hh4)^2)^(1/2)-yy1-yy2-yy3;ss4=hh4*yy4;fore=1:6hh5=355+(e-3)*5;yy5=(325^2-

7、(hh5)^2)^(1/2)-yy1-yy2-yy3-yy4;ss5=hh5*yy5;forf=1:6hh6=345+(f-3)*5;yy6=(325^2-(hh6)^2)^(1/2)-yy1-yy2-yy3-yy4-yy5;ss6=hh6*yy6;forg=1:6hh7=335+(g-3)*5;yy7=(325^2-(hh7)^2)^(1/2)-yy1-yy2-yy3-yy4-yy5-yy6;ss7=hh7*yy7;forh=1:6hh8=325+(h-3)*5;yy8=(325^2-(hh8)^2)^(

8、1/2)-yy1-yy2-yy3-yy4-yy5-yy6-yy7;ss8=hh8*yy8;fori=1:6hh9=315+(i-3)*5;yy9=(325^2-(hh9)^2)^(1/2)-yy1-yy2-yy3-yy4-yy5-yy6-yy7-yy8;ss9=hh9*yy9;forj=1:6hh10=305+(j-3)*5;yy10=(325^2-(hh10)^2)^(1/2)-yy1-yy2-yy3-yy4-