MATLAB数据处理与应用.docVIP

1.求解微分方程组，画出解函数图。 0t30 1. 源程序： （1） 编写Matlab函数M文件qimo_1如下function?du=qimo_1(t,u)?u(1)代表x,u(2)代表ydu=[-1*u(1)^3?-?u(2);u(1)?-?u(2)^3];（2）在命令窗口调用M文件qimo_1(t,u)[t,u]?=?ode23(qimo_1,[0,20],[1,0.5]);subplot(2,1,1)plot(t,u(:,1));grid;title(x),xlabel(t);subplot(2,1,2);plot(t,u(:,2)),grid; 2.运行结果及函数图像： 2.用生成一组稀疏数据，并用一维插值对该组数据进行插值，并与理论曲线进行比较，画出相应曲线。 （1）先做函数取0,0.5,1……9.5,10,每隔0.5取一个值，做散点图。（2）对散点分别用线性插值，三次多项式和样条插值并做出图像。 A．线性插值 1.源程序 2.运行结果及函数图像： B三次多项式插值 1.源程序 2.运行结果及函数图像： c.样条插值 1.源程序 2.运行结果及函数图像： 3. 厂向用户提供发动机，合同规定，第一、二、三季度末分别交货40台、60台、80台．每季度的生产费用为 （单位:元）, 其中x是该季度生产的台数．若交货后有剩余，可用于下季度交货，但需支付存储费，每台每季度c元．已知工厂每季度最大生产能力为100台，第一季度开始时无存货，设a=50、b=0.2、c=4，问:工厂应如何安排生产计划，才能既满足合同又使总费用最低．讨论a、b、c变化对计划的影响，并作出合理的解释。 ⒈问题的分析和假设：??⑴分析：问题的关键在于由于工厂的生产能力足以满足每个季度用户的需求，但是为了使总费用最少，那么利用每个季度生产费用的不同，可用利用上个生产费用低的季度多生产来为下个季度进行准备，前提是本月节省下的费用减去总的发动机存储费用还有剩余，这样生产才有价值，才可能满足合同的同时又能使总费用最低。?⑵基本假设：①工厂的生产能力不受外界环境因素影响。②为使总费用最低，又能满足合同要求，各个季度之间的生产数量之间是有联系的。③第一季度开始时无存货。④工厂每季度的生关费用与本季度生产的发动机台数有关。⑤生产要按定单的数量来进行，生产的数量应和订单的数量相同，以避免生产出无用的机器。? ⒉符号规定：X1―――第一季度生产发动机的数量?X2―――第二季度生产发动机的数量?X3―――第三季度生产发动机的数量⒊建模：1三个季度发动机的总的生产量为180台。?????????2每个季度的生产量和库存机器的数量之和要大于等于本季度的交货数量。?????3每个月的生产数量要符合工厂的生产能力。?????????4将实际问题转化为非线性规划问题，建立非线性规划模型。目标函数??min?f(x)=50(x1+x2+x3)+0.2(x12+x22+x32)+4(x1-40)+4(x1+x2-100)?整理，得???min?f(x)=50(x1+x2+x3)+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程序代码：⑴编写Matlab函数M文件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)?⑵在命令窗口调用M文件fun?(x)主程序?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.0000fval?=???????11280x?=?50.0000??60.0000??70.0000?fval?=?11280???⑵问题分析讨论：?