電力變壓器鐵心柱截面的優(yōu)化設計之程序實現(xiàn)(MATLAB)

1、電力變壓器鐵心柱截面的優(yōu)化設計之程序實現(xiàn)（MATLAB）附錄7.1 附錄 1funf=f=-1*(x(1)*sqrt(3252-x(1)2)+x(2)*(sqrt(3252-x(2)2)-sqrt(3252-x(1)2)+x(3)*(sqrt(3252-x(3)2)-sqrt(3252-x(2)2)+x(4)*(sqrt(3252-x(4)2)-sqrt(3252-x(3)2)+x(5)*(sqrt(3252-x(5)2)-sqrt(3252-x(4)2)+x(6)*(sqrt(3252-x(6)2)-sqrt(3252-x(5)2)+x(7)*(sqrt(3252-x(7)2)-sqrt(3

2、252-x(6)2)+x(8)*(sqrt(3252-x(8)2)-sqrt(3252-x(7)2)+x(9)*(sqrt(3252-x(9)2)-sqrt(3252-x(8)2)+x(10)*(sqrt(3252-x(10)2)-sqrt(3252-x(9)2)+x(11)*(sqrt(3252-x(11)2)-sqrt(3252-x(10)2)+x(12)*(sqrt(3252-x(12)2)-sqrt(3252-x(11)2)+x(13)*(sqrt(3252-x(13)2)-sqrt(3252-x(12)2)+x(14)*(sqrt(3252-x(14)2)-sqrt(3252-x(1

3、3)2); ; %最大面積的目標函數(shù)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

4、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 為倍數(shù)的寬度for i=1:14yy=sum(y,1);y(i)=sqrt(3252-x(i)2)-yy;endy=round(y);xys=x.*y;sum(s,1)(ans)/(3252)*pi)7.2 附錄 2

5、max=0;hh1=395;yy1=(3252-hh12)(1/2);ss1=hh1*yy1;for b=1:6hh2=385+(b-3)*5;yy2=(3252-(hh2)2)(1/2)-yy1;ss2=hh2*yy2;for c=1:6hh3=375+(c-3)*5;yy3=(3252-(hh3)2)(1/2)-yy1-yy2;ss3=hh3*yy3;for d=1:6hh4=365+(d-3)*5;yy4=(3252-(hh4)2)(1/2)-yy1-yy2-yy3;ss4=hh4*yy4;for e=1:6hh5=355+(e-3)*5;yy5=(3252-(hh5)2)(1/2)-y

6、y1-yy2-yy3-yy4;ss5=hh5*yy5;for f=1:6hh6=345+(f-3)*5;yy6=(3252-(hh6)2)(1/2)-yy1-yy2-yy3-yy4-yy5;ss6=hh6*yy6;for g=1:6hh7=335+(g-3)*5;yy7=(3252-(hh7)2)(1/2)-yy1-yy2-yy3-yy4-yy5-yy6;ss7=hh7*yy7;for h=1:6hh8=325+(h-3)*5;yy8=(3252-(hh8)2)(1/2)-yy1-yy2-yy3-yy4-yy5-yy6-yy7;ss8=hh8*yy8;for i=1:6hh9=315+(i-3

7、)*5;yy9=(3252-(hh9)2)(1/2)-yy1-yy2-yy3-yy4-yy5-yy6-yy7-yy8;ss9=hh9*yy9;for j=1:6hh10=305+(j-3)*5;yy10=(3252-(hh10)2)(1/2)-yy1-yy2-yy3-yy4-yy5-yy6-yy7-yy8-yy9;ss10=hh10*yy10;for k=1:6hh11=285+(k-3)*5;yy11=(3252-(hh11)2)(1/2)-yy1-yy2-yy3-yy4-yy5-yy6-yy7-yy8-yy9-yy10;ss11=hh11*yy11;for l=1:6hh12=265+(l

8、-3)*5;yy12=(3252-(hh12)2)(1/2)-yy1-yy2-yy3-yy4-yy5-yy6-yy7-yy8-yy9-yy10-yy11;ss12=hh12*yy12;for m=1:6hh13=245+(m-3)*5;yy13=(3252-(hh13)2)(1/2)-yy1-yy2-yy3-yy4-yy5-yy6-yy7-yy8-yy9-yy10-yy11-yy12;ss13=hh13*yy13;for n=1:6hh14=225+(n-3)*5;yy14=(3252-(hh14)2)(1/2)-yy1-yy2-yy3-yy4-yy5-yy6-yy7-yy8-yy9-yy10

9、-yy11-yy12-yy13;ss14=hh14*yy14;ss=ss1+ss2+ss3+ss4+ss5+ss6+ss7+ss8+ss9+ss10+ss11+ss12+ss13+ss14;if maxssmax=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;for i=length(s)-1:-1:1t=t s(i); %將23 級面積全部合并為一數(shù)組endnn=squre/7;