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

电力变压器铁心柱截面的优化设计之程序实现（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)*(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(325^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)-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=[395390385380375370365360355350345340335330];%上界[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);ss1=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-(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)^(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-yy5-yy6-yy7-yy8-yy9;ss10=hh10*yy10;fork=1:6hh11=285+(k-3)*5;yy11=(325^2-(hh11)^2)^(1/2)-yy1-yy2-yy3-yy4-yy5-yy6-yy7-yy8-yy9-yy10;ss11=hh11*yy11;forl=1:6hh12=265+(l-3)*5;yy12=(325^2-(hh12)^2)^(1/2)-yy1-yy2-yy3-yy4-yy5-yy6-yy7-yy8-yy9-yy10-yy11;ss12=hh12*yy12;form=1:6hh13=245+(m-3)*5;yy13=(325^2-(hh13)^2)^(1/2)-yy1-yy2-yy3-yy4-yy5-yy6-yy7-yy8-yy9-yy10-yy11-yy12;ss13=hh13*yy13;forn=1:6hh14=225+(n-3)*5;yy14=(325^2-(hh14)^2)^(1/2)-yy1-yy2-yy3-yy4-yy5-yy6-yy7-yy8-yy9-yy10-yy11-yy12-yy13;ss14=hh14*yy14;ss=ss1+ss2+ss3+ss4+ss5+ss6+ss7+ss8+ss9+ss10+ss11+ss12+ss13+ss14;ifmaxssmax=ss;hh=[hh1,hh2,hh3,hh4,hh5,hh6,hh7,hh8,hh9,hh10,hh11,hh12,hh13,hh14];yy=[yy1,yy2,yy3,yy4,yy5,yy6,yy7,yy8,yy9,yy10,yy11,yy12,yy13,yy14];sss=[ss1,ss2,ss3,ss4,ss5,ss6,ss7,ss8,ss9,ss10,ss11,ss12,ss13,ss14];endend，end，end，end，end，end，end，end，end，end，endendend7.3附录3t=s;fori=length(s)-1:-1:1t=[ts(i)];%将23级面积全部合并为一数组endnn=squre/7;%平均分割时的面积ss=0;w=30000;%误差设置的初始值sss=[];whilelength(sss)~=7w=w-1000;fori=1:23ss=ss+t(i);ifnn-w=ss=nn+wsss=union(sss,ss);ss=0;endendendsss%油道分割的各级面积