3)+0.2(x12+x22+x32)+4(2x1+x2-140)Р约束函数 s.t x1+x2≥100;Р X1+x2+x3=180;Р 40≤x1≤100;Р 0≤x2≤100;Р 0≤x3≤100;Р求解的Matlab程序代码:Р M-文件 fun.m:Р function f=fun (x);Рf=50*(x(1)+x(2)+x(3))+0.2*(x(1)^2+x(2)^2+x(3)^2)+4*(2*x(1) +x(2)-140)Р主程序fxxgh.m:Рx0=[60;60;60];РA=[-1 -1 0];b=[-100];РAeq=[1 1 1];beq=[180];Рvlb=[40;0;0];vub=[100;100;100];Р[x,fval]=fmincon('fun',x0,A,b,Aeq,beq,vlb,vub)Р计算结果与问题分析讨论:Р 计算结果:x =Р50.0000Р 60.0000Р 70.0000Рfval =Р11280 Р问题分析讨论:Р 由运算结果得:该厂第一季度、第二季度、第三季度的生产量分别是50台、60台和70台时,才能既满足合同又使总费用最低,费用最低为11280元。Р给出下列为微分方程数值解的求解程序Р1. (课本,P132页例4)Р2. (课本,P132页例5)Р3. (课本,P132页例4,系数“1000”被改为了“7”)Р4. Р解:(1) 微分方程M函数文件Рfunction dy=daodan01(x,y)Рdy=zeros(2,1);Рdy(1)=y(2);Рdy(2)=0.2*(1/(1-x))*sqrt(1+y(2)^2)Р(2)建立主程序求解Рclear all;Рt0=0;tf=10;Р[t,y]=ode45('daodan01',[t0,tf],[0,0]);Рplot(t,y,'r*');Рt , y
全文预览
