该楼层疑似违规已被系统折叠 隐藏此楼查看此楼

function optimization4()

A=[];b=[];

Aeq=[];beq=[];LB=[0.1;0.03;0.03;0.1;0.03;0.03];UB=[0.4;0.06;0.06;0.4;0.06;0.06];

nvars=6;

options=gaoptimset('TimeLimit',inf,'PlotFcns',{@gaplotbestf},'PopulationSize',10,'Generations',15,'PopInitRange',[0;0.6]);

[x,fval,exitflag]=ga(@fitnessfcn4,nvars,A,b,Aeq,beq,LB,UB,@nonlcon3,options)

function y=fitnessfcn4(x)

y=x(1)*x(2)+0.38*x(3)+x(4)*x(5)+0.2*x(6);

function [c ceq]=nonlcon3(x)

fid=fopen('nmsl.txt','w+');

fprintf(fid,'%6.2f%6.2f%6.2f%6.2f%6.2f%6.2f\n',x');

fclose(fid);

!SET KMP_STACKSIZE=2048k & C:\"Program Files"\"ANSYS Inc"\v150\ANSYS\bin\winx64\ansys150 -b -p struct -i C:\Users\Administrator\bett.mac -o C:\Users\Administrator\temp.txt

fid=fopen('C:\Users\Administrator\result.txt','r');

tline=fgetl(fid);

y=str2num(tline);

fclose(fid);

c(1)=y(1)-1.41e8;

c(2)=abs(y(2))-2e8;

c(3)=abs(y(3))-2.7e8;

c(4)=1.2e6-y(4);

ceq=[];