电力变压器铁心柱截面的优化设计之程序实现(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(3 25^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(3 25^2-x(13)^2)-sqrt(325^2-x(12)^2))+x(14)*(sqrt(325^2-x(14)^2)-sqrt(32 5^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;