第三章数值计算.pptVIP

第三章 数值运算基础 定积分的近似计算 微分方程的求解 线性规划问题的求解 矩形法 定积分的定义： 矩形法 定积分的近似： 左点法、右点法和中点法 矩形法举例 矩形法举例 定积分几何意义 梯形法 梯形法 梯形法举例 抛物线法 Matlab 计算定积分函数介绍 梯形法：trapz trapz 举例 抛物线法 quad 举例 二重积分的计算 抛物线法计算二重积分： dblquad dblquad 举例 dblquad 举例 数值实验 数值实验 Matlab 微分方程 dsolve 求解析解 dsolve 的使用 dsolve 举例 dsolve 举例 dsolve 举例 Matlab函数数值求解 Matlab提供的ODE求解器 参数说明 数值求解举例 数值求解举例 odefun 为显式常微分方程，可以用命令 inline 定义，或在函数文件中定义，然后通过函数句柄调用。 fun=inline(-2*y+2*x^2+2*x,x,y); [x,y]=ode23(fun,[0,0.5],1); 注：也可以在 tspan 中指定对求解区间的分割，如： [x,y]=ode23(fun,[0:0.1:0.5],1); %此时 x=[0:0.1:0.5] [T,Y] = solver(odefun,tspan,y0) 求初值问题 的数值解，求解范围为 [0,0.5] 例 4： 如果需求解的问题是高阶常微分方程，则需将其化为一阶常微分方程组，此时需用函数文件来定义该常微分方程组。 令 ，则原方程可化为 求解 Ver der Pol 初值问题 例 5： 先编写函数文件 verderpol.m function xprime=verderpol(t,x) global mu; xprime=[x(2); mu*(1-x(1)^2)*x(2) - x(1)]; 再编写脚本文件 vdpl.m，在命令窗口直接运行该文件。 clear; global mu;mu=7; y0=[1;0]; [t,x]=ode45(verderpol,[0,40],y0); plot(t,x(:,1),r-); 问题一 ： 任务分配问题：某车间有甲、乙两台机床，可用于加工三种工件。假定这两台车床的可用台时数分别为800和900，三种工件的数量分别为400、600和500，且已知用三种不同车床加工单位数量不同工件所需的台时数和加工费用如下表。问怎样分配车床的加工任务，才能既满足加工工件的要求，又使加工费用最低？ 三、线性规划问题 解 设在甲车床上加工工件1、2、3的数量分别为x1、x2、x3，在乙车床上加工工件1、2、3的数量分别为x4、x5、x6。可建立以下线性规划模型： 问题二： 某厂每日8小时的产量不低于1800件。为了进行质量控制，计划聘请两种不同水平的检验员。一级检验员的标准为：速度25件/小时，正确率98%，计时工资4元/小时；二级检验员的标准为：速度15件/小时，正确率95%，计时工资3元/小时。检验员每错检一次，工厂要损失2元。为使总检验费用最省，该工厂应聘一级、二级检验员各几名？ 解 设需要一级和二级检验员的人数分别为x1、x2人, 则应付检验员的工资为： 因检验员错检而造成的损失为： 故目标函数为： 约束条件为： 线性规划模型： 用MATLAB优化工具箱解线性规划 min z=cX 1、模型： 命令：x=linprog（c，A，b） 2、模型：min z=cX 命令：x=linprog（c，A，b，Aeq,beq） 注意：若没有不等式： 存在，则令A=[ ]，b=[ ]. 3、模型：min z=cX VLB≤X≤VUB 命令：[1] x=linprog（c，A，b，Aeq,beq, VLB，VUB） [2] x=linprog（c，A，b，Aeq,beq, VLB，VUB, X0） 注意：[1] 若没有等式约束: , 则令Aeq=[ ], beq=[ ]. [2]其中X0表示初始点 4、命令：[x,fval]=linprog(…) 返回最优解ｘ及ｘ处的目标函数值fval. * * 一、定积分的近似计算 定积分计算的基本公式是牛顿－莱布尼兹公式。但当被积函数的原函数不知道时，如何计算？这时就需要利用近似计算。特别是在许