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

电力变压器铁心柱截面的优化设计之程序实现（MATLAB） 附录 7.1 附录 1 funf=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=[funf fung]; x0=[180 185 190 195 200 205 210 215 225 225 230 235 240 245]; %初始值 options=[]; vlb=[85 80 75 70 65 60 55 50 45 40 35 30 25 20]; %下界 vub=[395 390 385 380 375 370 365 360 355 350 345 340 335 330]; %上界 [x,options]=constr(fun,x0,options,vlb,vub); y=zeros(1,14); x=x.*2; x=(round(x./10)).*10; %以10 为倍数的宽度 for i=1:14 yy=sum(y,1); y(i)=sqrt(325^2-x(i)^2)-yy; end y=round(y); x y s=x.*y; sum(s,1) (ans)/((325^2)*pi) 7.2 附录 2 max=0; hh1=395; yy1=(325^2-hh1^2)^(1/2); ss1=hh1*yy1; for b=1:6 hh2=385+(b-3)*5; yy2=(325^2-(hh2)^2)^(1/2)-yy1; ss2=hh2*yy2; for c=1:6 hh3=375+(c-3)*5; yy3=(325^2-(hh3)^2)^(1/2)-yy1-yy2; ss3=hh3*yy3; for d=1:6 hh4=365+(d-3)*5; yy4=(325^2-(hh4)^2)^(1/2)-yy1-yy2-yy3; ss4=hh4*yy4; for e=1:6 hh5=355+(e-3)*5; yy5=(325^2-(hh5)^2)^(1/2)-yy1-yy2-yy3-yy4; ss5=hh5*yy5; for f=1:6 hh6=345+(f-3)*5; yy6=(325^2-(hh6)^2)^(1/2)-yy1-yy2-yy3-yy4-yy5; ss6=hh6*yy6; for g=1:6 hh7=335+(g-3)*5; yy7=(325^2-(hh7)^2)^(1/2)-yy1-yy2-yy3-yy4-yy5-yy6; ss7=hh7*yy7; for h=1:6 hh8=325+(h-3)*5; yy8=(325^2-(hh8)^2)^(1/2)-yy1-yy2-yy3-yy4-yy5-yy6-yy7; ss8=hh8*yy8; for i=1:6 hh9=315+(i-3)*5; yy9=(325^2-(hh9)