大学数学实验八线性规划
实验报告
【实验目的】
1、掌握用MATLAB优化工具箱解线性规划的方法。
2、练习建立实际问题的线性规划模型。
【实验内容】
6某银行经理计划用一笔资金进行证券投资，可供购进的证券以及其信用等级、到期年限、收益如下表所示。按照规定，市政证券的收益可以免税，其他证券的收益需按50％的税率纳税，此外还有以下限制：
①政府及代办机构的证券总共至少要购进400万元；
②所购证券的平均信用等级不超过1.4（信用等级数字越小，信用程度越高）；
③所购证券的平均到期年限不超过5年。
试问：
①若该经理有1000万元资金，应如何投资？
②如果能够以2.75%的利率借到不超过100万元资金，该经理应如何操作？并考虑利率在什么范围内变化时，投资方案不改变？
③在1000万元资金情况下，若证券A的税前收益增加为4.5%,投资应否改变？若证券C 的税前收益减少为4.8%，投资应否改变？
证券名称证券种类信用等级到期年限/年到期税前收益/%
A市政 2 9 4.3
B代办机构 2 15 5.4
C政府 1 4 5.0
D政府 1 3 4.4
E市政 5 2 4.5
一、线性规划的模型建立
这很明显是一个约束优化的问题。其模型应该包括决策变量、目标函数和约束条件三个部分。
（一）决策变量
具体到这个问题，决策变量应该指的是这个银行经理的投资意向。设他对这5个证券分别投资了、、、、万元。
（二）目标函数
毫无疑问，这个银行经理进行投资活动，一定是希望能够利用手头的资金获得最大的收益，所以目标函数是最终的总收益。最终的总收益即每种证券的收益之和，再扣除非市政证券的收益税收。具体到这个问题，非市政证券的所得税率是50%，相当于非市政证券的收益减小了一半。故有：
（三）约束条件
这个问题有来自几方面的约束。
（1）政府和代办机构证券的最低购入量限制：
（2）所购证券的信用等级限制：
即
（3）所购证券的到期年限限制：
即
（4）购买资金总额限制（按照常理，要让手头的资金发挥最大的价值就必须全部投资）：
（5）决策变量皆为非负：
、、、、
综上所述，可得如下基本模型。
、、、、
二、利用MATLAB分析投资方案
利用MATLAB的linprog程序求解这个问题的最优解。有几个问题需要说明：（1）程序默认求解最小值，故将目标函数取负。（2）约束条件中有大于号的应两边取负变为小于号。具体程序如下。
运行上面的程序后得到的最优解和最优值如下。
可见，这个银行经理手中的1000万元应该这样投资：购买A 证券218.1818万元，购买C 证券736.3636万元，购买E 证券45.4545万元。最终的收益为29.8364万元。
输出的lag.ineqlin 如下：
[0.0000, 0.0062, 0.0024]
这表明第一项约束式并没有起作用，而第二个和第三个约束条件起了作用，即在这个最优解下，第2、3个约束条件是等式约束。
下面尝试使用不同的算法进行计算，通过比较它们的迭代次数来反映其优劣。前面使用的方法是内点法，程序输出的迭代次数为6。根据教材上对内点法的介绍，它在解决特别大规模的问题时有明显优势，本题属于小规模的线性规划问题，内点法是否有优势不能确定。用有效集方法和单纯形算法计算的程序如下（加在上面的程序之后即可）。
运行程序后发现，输出的最优解和最优值完全一样，而有效集方法的迭代次数为4，单纯形算法的迭代次数仅为1。可见，在像本题一样的小规模线性规划问题求解过程中，内点法并没有明显优势。
三、分析可贷款情况下的投资方案
在上述程序运行完毕的基础上，让MATLAB 输出lag.eqlin ，输出量为0.0298。这个数字的意义是，总投资金额每增加1万元，最终收益会增加0.0298万元。也就是说投资金额每增加100万元，最终收益会增加2.98万元。扣除贷款利率的2.75万元，最终收益能够增加0.23万元。这样来看，最终收益约为30.07万元。
将上述程序中的1000改为1100再运行，得到的最优解是
得到的最优值为
= 32.8200
简单计算一下，发现对证券A 、C 、E 的投资分别增加了21.8182、73.6364、4.5455。也就是说，新增加的
100万元投资也是按照原来对三个证券的投资比例进行投资的，可见投资方案
并没有改变。而最终的收益为
32.8200 – 2.75 = 30.0700 （万元）
与之前的分析吻合。
根据上述分析，只要贷款利率小于2.98%，将贷款来的钱再投资进去，最终收益就会增加。一旦利率等于或大于这个值，投资收益的增加不能抵过贷款利息的偿还，最终收益反而会减小，这样的话就不应贷款。
四、分析收益率扰动对投资方案的影响
将程序做如下修改，以观察证券A 、C 的收益率的扰动对投资方案会不会产生影响。
当A 的收益率增加到4.5%时，输出的最优解与之前的结果完全一样，说明投资方案不用改变。而当C 的收益率减小到4.8%时，输出的最优解发生了很大的变化，结果为
虽然变化量都是0.2%，但是产生的结果却完全不一样，这说明了最优投资计划对证券C 收益率的敏感程度高于对证券A 收益率的敏感程度。
上面所作的，实际上是线性规划的敏感性分析，就是要求最优解不变条件下目标函数系数的允许的变化范围。遗憾的是，MATLAB 没有给出这种敏感性分析的结果。事实上，经过进一步的枚举试验，发现证券A 收益率在3.1%~4.6%范围内变动，最优解均不变，而证券C 收益率的变化范围是4.9%~8.4%。实际问题中的“波动”是不可能预知是增加还是减小的，证券C 的收益率变化范围虽然比较大，但是一旦降低超过0.1%最优解就发生了变化，相比之下，证券A 的收益率上下变动0.3%范围内最优解都不会发生变化，可见最优投资计划对证券C 收益率的敏感程度高于对证券A
收益率的敏感程度。
8某牧场主知道，对于一匹平均年龄的马来说，最低的营养需求为：40磅蛋白质，20磅碳水化合物，45磅粗饲料。这些营养成分是从不同的饲料中得到的，饲料及其价格在下表中列出。建立数学模型，确定如何以最低的成本满足最低的营养需求。
蛋白质/磅碳水化合物/磅粗饲料/磅价格/美元干草/捆0.5 2.0 5.0 1.80
燕麦片/袋 1.0 4.0 2.0 3.50
饲料块/块 2.0 0.5 1.0 0.40
高蛋白浓缩料/袋 6.0 1.0 2.5 1.00
每匹马的需求/天40.0 20.0 45.0 -
一、线性规划的模型建立
这个约束优化的问题的模型也应该包括决策变量、目标函数和约束条件三项内容。（一）决策变量
这个问题的决策变量是干草、燕麦片、饲料快和高蛋白浓缩料这4种饲料的用量，设为、、、，单位为上表中第一列中的单位。
（二）目标函数
这个问题的核心在于如何在满足马的基本营养需求的条件下，把饲料的成本压得最低，故目标函数应为饲料的总价格，设为，且有
（三）约束条件
这个问题比较简单，只有两个约束条件。
（1）马的基本营养需求限制：
（2）决策变量皆为非负：
、、、
综上所述，可得到如下模型：
、、、
二、利用MATLAB分析饲料购买方案
这里还是利用MATLAB里的linprog程序求解。
运行这个程序，输出的最优解为
根据题意，决策变量应该取整数，可是这样的话，前三个约束条件中的等号就不能严格满足。实际情况是，将等号改为大于等于号，也可以确保马的最低营养需要被满足。程序用上面的即可，只是在调用linprog命令时，需要改变一下参数。具体如下。
再运行这个程序，得到的最优解为
最优值为
= 17
输出的lag.ineqlin如下：
[0.0000, 0.4000, 0.2000]
这表明第一项约束式并没有起作用，而第二个和第三个约束条件起了作用，即在这个最优解下，第2、3个约束条件是等式约束。
可见，这个牧场主应该给他的马配备的饲料是：5捆干草和20块饲料块，这样饲料的
在上一题中，经过分析发现内点法并没有明显优势。本题亦是一道小规模的线性规划问题，内点法是不是依然没有优势呢？下面比较一下用不同算法计算时的迭代次数。上面用的算法是内点法，输出的迭代次数是10。
在之前程序的基础上运行上面的程序。首先，不同算法得到的最优解和最优值是完全一样的。其次比较它们的迭代次数。有效集方法的迭代次数为5，单纯形算法的迭代次数为2。这再次说明，在小规模的线性规划问题中，内点法没有明显优势。
【实验总结】
这次的实验内容比以往任何一次都来得简单、易于操作。首先，线性规划的模型建立遵循着约束优化模型建立的一般步骤，根据题目的叙述很容易建立相应模型。其次，MATLAB软件提供了求解线性规划问题的专门程序，还提供了不同的算法，可以进行准确计算。另外，这已经是第六次数学实验了，MATLAB软件的应用已经不成问题。
虽然这次实验很简单，但是它的重要性却很大，因为生活中也好，其他学科领域也好，许多问题最后要归结到线性规划问题上。通过本次实验，我了解了线性规划问题求解的算法，掌握了：（1）建立实际问题的线性规划模型；（2）利用MATLAB提供的程序求解线性规划问题；（3）对MATLAB提供的结果进行简单分析。
谢谢老师和助教的指导。
北京交通大学海滨学院考试试题 课程名称：数学实验2010-2011第一学期出题教师：数学组适用专业: 09机械, 物流, 土木, 自动化 班级：学号：姓名： 选做题目序号： 1.一对刚出生的幼兔经过一个月可以长成成兔, 成兔再经过一个月后可以 繁殖出一对幼兔. 如果不计算兔子的死亡数, 请用Matlab程序给出在未来24个月中每个月的兔子对数。 解: 由题意每月的成兔与幼兔的数量如下表所示： 1 2 3 4 5 6 ··· 成兔0 1 1 2 3 5··· 幼兔 1 0 1 1 2 3··· 运用Matlab程序： x=zeros(1,24); x(1)=1;x(2)=1； for i=2:24 x(i+1)=x(i)+x(i-1); end x 结果为x = 1 1 2 3 5 8 13 21 3 4 5 5 89 144 233 377 610 987 1597 2584 4181 6765 1094 6 7711 2865 7 46368 2.定积分的过程可以分为分割、求和、取极限三部分, 以1 x e dx 为例, 利用
已学过的Matlab 命令, 通过作图演示计算积分的过程, 并与使用命令int() 直接积分的结果进行比较. 解：根据求积分的过程，我们先对区间[0,1]进行n 等分， 然后针对函数x e 取和，取和的形式为10 1 i n x i e e dx n ξ=≈ ∑ ? ，其中1[ ,]i i i n n ξ-?。这里取i ξ为区间的右端点，则当10n =时，1 x e dx ?可用10 101 1.805610 i i e ==∑ 来近似计算， 当10n =0时，100 100 1 01 =1.7269100 i x i e e dx =≈ ∑?，当10n =000时，10000 10000 1 1 =1.718410000 i x i e e dx =≈ ∑ ?. 示意图如下图，Matlab 命令如下： x=linspace (0,1,21); y=exp(x); y1=y(1:20); s1=sum(y1)/20 y2=y(2:21); s2=sum(y2)/20 plot(x,y); hold on for i=1:20 fill([x(i),x(i+1),x(i+1),x(i),x(i)],[0,0,y(i),y(i),0],'b') end syms k;symsum(exp(k/10)/10,k,1,10)；%n=10 symsum(exp(k/100)/100,k,1,100)；%n=100 symsum(exp(k/10000)/10000,k,1,10000)；%n=10000
一、实验内容 P206第六题 function f=wuyan2(c) y=[3.9 5.3 7.2 9.6 12.9 17.1 23.2 31.41 38.6 50.2 62.9 76.0 92.0 106.5 123.2 131.7 150.7 179.3 204.0 226.5 251.4 281.4] t=[0 10 20 30 40 50 60 70 80 90 100 110 120 130 140 150 160 170 180 190 200 210] f=y-c(1)/(1+c(1)/3.9-1)*exp^(-c(2)*t) c0=[1 1] c=lsqnonlin('wuyan2',c0) P206第七题 function f=wuyan1(c) q=[0.4518 0.4862 0.5295 0.5934 0.7171 0.8964 1.0202 1.1963 1.4928 1.6909 1.8548 2.1618 2.6638 3.4634 4.6759 5.8478 6.7885 7.4463 7.8345 8.2068 8.9468 9.7315 10.5172 11.7390 13.6876 ]; k=[0.0911 0.0961 0.1230 0.1430 0.1860 0.2543 0.3121 0.3792 0.4754 0.4410 0.4517 0.5595 0.8080 1.3072 1.7042 2.0019 2.2914 2.4941 2.8406 2.9855 3.2918 3.7214 4.3500 5.5567 7.0477]; l=[4.2361 4.3725 4.5295 4.6436 4.8179 4.9873 5.1282 5.2783 5.4334 5.5329 6.4749 6.5491 6.6152 6.6808 6.7455 6.8065 6.8950 6.9820 7.0637 7.1394 7.2085 7.3025 7.3470 7.4432 7.5200]; f=q-c(1)*k.^c(2).*l.^c(3) c0=[1 1 1] c=lsqnonlin('wuyan1',c0) c = 0.4091 0.6401 1.1446 a=0.4091 α=0.6401 β=1.1446 P239第五题 c=[-20 -30]; A=[1 2;5 4]; b=[20 70]; v1=[0 0]; [x,f,ef,out,lag]=linprog(c,A,b,[],[],v1) z=-f x = 10.0000 5.0000
实验2 方程模型及其求解算法 一、实验目的及意义 [1] 复习求解方程及方程组的基本原理和方法； [2] 掌握迭代算法； [3] 熟悉MATLAB软件编程环境；掌握MATLAB编程语句(特别是循环、条件、控制等语句)； [4] 通过范例展现求解实际问题的初步建模过程； 通过该实验的学习，复习和归纳方程求解或方程组求解的各种数值解法（简单迭代法、二分法、牛顿法、割线法等），初步了解数学建模过程。这对于学生深入理解数学概念，掌握数学的思维方法，熟悉处理大量的工程计算问题的方法具有十分重要的意义。 二、实验内容 1．方程求解和方程组的各种数值解法练习 2．直接使用MATLAB命令对方程和方程组进行求解练习 3．针对实际问题，试建立数学模型，并求解。 三、实验步骤 1．开启软件平台——MATLAB，开启MATLAB编辑窗口； 2．根据各种数值解法步骤编写M文件 3．保存文件并运行； 4．观察运行结果(数值或图形)； 5．根据观察到的结果写出实验报告，并浅谈学习心得体会。 四、实验要求与任务 基础实验 1．用图形放大法求解方程x sin(x) = 1. 并观察该方程有多少个根。 画出图形程序： x=-10:0.01:10; y=x.*sin(x)-1; y1=zeros(size(x)); plot(x,y,x,y1) MATLAB运行结果：
-10-8-6-4-20246810 -8-6 -4 -2 2 4 6 8 扩大区间画图程序： x=-50:0.01:50; y=x.*sin(x)-1; y1=zeros(size(x)); plot(x,y,x,y1) MATLAB 运行结果： -50-40-30-20-1001020304050 由上图可知，该方程有偶数个无数的根。
《管理运筹学》实验报告 实验日期： 2016年 04月 21日—— 2016 年 05 月 18 日 班级2014级04班姓名杨艺玲学号56 实验 管理运筹学问题的计算机求解 名称 实验目的： 通过实验学生应该熟练掌握“管理运筹学”软件的使用，并能利用“管理运筹学”对具体问题进行问题处理，且能对软件处理结果进行解释和说明。 实验所用软件及版本： 管理运筹学 实验过程：（含基本步骤及异常情况记录等) 一、实验步骤（以P31页习题1 为例） 1.打开软件“管理运筹学” 2.在主菜单中选择线性规划模型，屏幕中会出现线性规划页面
3.在点击“新建”按钮以后，按软件的要求输入目标函数个数和约束条件个数，输入目标函数级约束条件的歌变量的系数和b值，并选择好“≤”、“≥”或“＝”，如图二所示，最后点击解决 4.注意事项： （1）输入的系数可以是整数、小数，但不能是分数，要把分数化为小数再输入。（2）输入前要合并同类项。 当约束条件输入完毕后，请点击“解决”按钮，屏幕上讲显现线性规划问题的结果，如图所示
5.输出结果如下
5.课后习题： 一、P31习题1 某家具公司生产甲、乙两种型号的组合柜，每种组合柜需要两种工艺（制白坯和油漆）.甲型号组合柜需要制白坯6工时，油漆8工时：乙型号组合柜需要制白坯12工时，油漆4工时.已知制白坯工艺的生产能力为120工时/天，油漆工艺的生产能力为64工时/天，甲型号组合柜单位利润200元，乙型号组合柜单位利润为240元. 约束条件： 问题： （1）甲、乙两种柜的日产量是多少这时最大利润是多少 答：由实验过程中的输出结果得甲组合柜的日产量是4个，乙的事8个。 （2）图中的对偶价格的含义是什么 答： 对偶价格的含义是约束条件2中，每增加一个工时的油漆工作，利润会增加元。 （3）对图中的常数项范围的上、下限的含义给予具体说明，并阐述如何使用这些信息。 答：当约束条件1的常数项在48~192范围内变化，且其他约束条件不变时，约束条件1的对偶价格不变，仍为；当约束条件2的常数项在40~180范围内变化，而其他约束条件的常数项不变时，约束条件2的对偶价格不然，仍为。 （4）若甲组合柜的利润变为300，最优解不变为什么 . 0,0,6448,120126; 240200 z max ≥≥≤+≤++=y x y x y x y x
东华大学M A T L A B数学实验第二版答案(胡良 剑) -CAL-FENGHAI-(2020YEAR-YICAI)_JINGBIAN
数学实验答案 Chapter 1 Page20,ex1 (5) 等于[exp(1),exp(2);exp(3),exp(4)] (7) 3=1*3, 8=2*4 (8) a为各列最小值，b为最小值所在的行号 (10) 1>=4,false, 2>=3,false, 3>=2, ture, 4>=1,ture (11) 答案表明：编址第2元素满足不等式(30>=20)和编址第4元素满足不等式(40>=10) (12) 答案表明：编址第2行第1列元素满足不等式(30>=20)和编址第2行第2列元素满足不等式(40>=10) Page20, ex2 (1)a, b, c的值尽管都是1，但数据类型分别为数值，字符，逻辑，注意a与c 相等，但他们不等于b (2)double(fun)输出的分别是字符a,b,s,(,x,)的ASCII码 Page20,ex3 >> r=2;p=0.5;n=12; >> T=log(r)/n/log(1+0.01*p) Page20,ex4 >> x=-2:0.05:2;f=x.^4-2.^x; >> [fmin,min_index]=min(f) 最小值最小值点编址 >> x(min_index) ans = 0.6500 最小值点 >> [f1,x1_index]=min(abs(f)) 求近似根--绝对值最小的点 f1 = 0.0328 x1_index = 24 >> x(x1_index) ans = -0.8500 >> x(x1_index)=[];f=x.^4-2.^x; 删去绝对值最小的点以求函数绝对值次小的点>> [f2,x2_index]=min(abs(f)) 求另一近似根--函数绝对值次小的点 f2 = 0.0630 x2_index = 65 >> x(x2_index) ans =
《管理运筹学》实验报告 实验日期：2016年04月21日——2016年05月18日 实验目的： 通过实验学生应该熟练掌握“管理运筹学 3.0”软件的使用，并能利用“管理运筹学 3.0” 对具体问题进行问题处理，且能对软件处理结果进行解释和说明。实验所用软件及版本：管理运筹学3.0 实验过程：（含基本步骤及异常情况记录等―） 一、实验步骤（以P31页习题1为例） 1?打开软件“管理运筹学3.0” 2?在主菜单中选择线性规划模型，屏幕中会出现线性规划页面 3?在点击“新建”按钮以后，按软件的要求输入目标函数个数和约束条件个数，输入目标函数级约束条件的歌变量的系数和b值，并选择好“w”、“》”或“二”, 如图二所示，最后点击解决 班级2014级04班姓名杨艺玲学号2014190456实验 名称 管理运筹学问题的计算机求解 n 幵 目标的数 娈童个数约束条件个数 芙 遇出 保存解决关于
X 4?注意事项： (1)输入的系数可以是整数、小数，但不能是分数，要把分数化为小数再输入。 (2)输入前要合并同类项。 当约束条件输入完毕后，请点击“解决”按钮，屏幕上讲显现线性规划问题的结果, 如 图所示 D tiff 0% 关于遇出 变童个数约朿条件个数F目标的数3V 标淮北结杲: 上一曲
5.输出结果如下 me車最优解如下***#尊1林*祜除目标函数最优值知2?20 变1 最优解相差値 XI 4.00 0.00 X2 8.00 0100 釣束松弛颅11余变量对偶价格 01. 00 16. 5€ 0.00 13.33 目标函数系数范園: 娈1下限当前值上限 XI 120. 30 200.00430. 00 X2 100. 0D 240.00400.00 常数【页范園； 的束T眼当前值上限 143.00120 00152.00 240.00 64.00 160.00 5.课后习题： 一、P31习题1 某家具公司生产甲、乙两种型号的组合柜，每种组合柜需要两种工艺（制白坯和油漆）.甲型号组合柜需要制白坯6工时，油漆8工时：乙型号组合柜需要制白坯12工时，油漆4工时.已知制白坯工艺的生产能力为120工时/天，油漆工艺的生产能力为64工时/天，甲型号组合柜单位利润200元，乙型号组合柜单位利润为240 元. max z = 200x 240y; 约束条件：6x，12心2°, 8x +4y 兰64, x 一0, y -0. 问题： （1）甲、乙两种柜的日产量是多少？这时最大利润是多少？ 答：由实验过程中的输出结果得甲组合柜的日产量是4个，乙的事8个
实验7：线性规划与非线性规划 班级：2015级电科班，学号：222015333210187，姓名：吴京宣，第1组 ====================================================================== 一、实验目的： 1. 了解线性规划的基本内容。 2. 直观了解非线性规划的基本内容。 3. 掌握用数学软件求解优化问题。 二、实验内容 1. 两个引例. 2. 用数学软件包MATLAB求解线性规划与非线性规划问题. 3. 用数学软件包LINDO、LINGO求解线性规划问题. 4. 建模案例：投资的收益与风险. 5. 非线性规划的基本理论 6. 钢管订购及运输优化模型． 三、实验步骤 对以下问题,编写M文件: 1.某厂生产甲乙两种口味的饮料,每百箱甲饮料需用原料6千克,工人10名,可获利10万元;每百箱乙饮料需用原料5千克,工人20名,可获利9万元.今工厂共有原料60千克,工人150名,又由于其他条件所限甲饮料产量不超过800箱.问如何安排生产计划,即两种饮料各生产多少使获利最大.进一步讨论: 1)若投资0.8万元可增加原料1千克,问应否作这项投资. 2)若每100箱甲饮料获利可增加1万元,问应否改变生产计划. 2.某厂向用户提供发动机，合同规定，第一、二、三季度末分别交货40台、60 台、80台．每季度的生产费用为（单位:元）, 其中x 是该季度生产的台数．若交货后有剩余，可用于下季度交货，但需支付存储费，每台每季度c元．已知工厂每季度最大生产能力为100台，第一季度开始时无存货，设a=50、b=0.2、c=4，问:工厂应如何安排生产计划，才能既满足合同又使总费用最低．讨论a、b、c变化对计划的影响，并作出合理的解释．
数学建模与实验课程实验报告 实验名称三、线性规划与整数规划实验地点日期2014-10-28 姓名班级学号成绩 【实验目的及意义】 [1] 学习最优化技术和基本原理，了解最优化问题的分类； [2] 掌握规划的建模技巧和求解方法； [3] 学习灵敏度分析问题的思维方法； [4] 熟悉MATLAB软件求解规划模型的基本命令； [5] 通过范例学习，熟悉建立规划模型的基本要素和求解方法。 通过该实验的学习，使学生掌握最优化技术，认识面对什么样的实际问题，提出假设和 建立优化模型，并且使学生学会使用MATLAB、Lingo软件进行规划模型求解的基本命令， 并进行灵敏度分析。解决现实生活中的最优化问题是本科生学习阶段中一门重要的课程，因 此，本实验对学生的学习尤为重要。 【实验要求与任务】 根据实验内容和步骤，完成以下实验，要求写出实验报告（符号说明—模型的建立—模型 的求解（程序）—结论） A组 高校资金投资问题 高校现有一笔资金100万元，现有4个投资项目可供投资。 项目A：从第一年到底四年年初需要投资，并于次年年末回收本利115%。 项目B：从第三年年初需要投资，并于第5年末才回收本利135%，但是规定最大投资总 额不超过40万元。 项目C：从第二年年初需要投资，并于第5年末才回收本利M%，但是规定最大投资总 额不超过30万元。(其中M为你学号的后三位+10) 项目D：五年内每年年初可以买公债，并于当年年末归还，并可获得6%的利息。 试为该校确定投资方案，使得第5年末他拥有的资金本利总额最大。 该校在第3年有个校庆，学校准备拿出8万元来筹办，又应该如何安排投资方案，使得 第5年末他拥有的资金本利总额最大。 B组题 1)最短路问题, 图1中弧上的数字为相邻2点之间的路程，求从1到7的最短路。 图1 图 2 r为你的学号后2位+10 其中 1 2)最大车流量, 图1中弧上的数字为相邻2点之间每小时的最大车流量。求每小时1到7最大
二、判断题（“对”在题号前（）中打√×）（10分） （√）1、误差是指测量值与真值之差，即误差=测量值－真值，如此定义的误差反映的是测量值偏离真值的大小和方向，既有大小又有正负符号。 （×）2、残差（偏差）是指测量值与其算术平均值之差，它与误差定义一样。（√）3、精密度是指重复测量所得结果相互接近程度，反映的是随机误差大小的程度。 （√）4、测量不确定度是评价测量质量的一个重要指标，是指测量误差可能出现的范围。 （×）7、分光计设计了两个角游标是为了消除视差。 （×）9、调节气垫导轨水平时发现在滑块运动方向上不水平，应该先调节单脚螺钉再调节双脚螺钉。 （×）10、用一级千分尺测量某一长度（Δ仪=0.004mm），单次测量结果为N=8.000mm，用不确定度评定测量结果为N=（8.000±0.004）mm。 三、简答题（共15分） 1．示波器实验中，（1）CH1（x）输入信号频率为50Hz，CH2（y）输入信号频率为100Hz；（2）CH1（x）输入信号频率为150Hz，CH2（y）输入信号频率为50Hz；画出这两种情况下，示波器上显示的李萨如图形。（8分）
差法处理数据的优点是什么？（7分） 答：自变量应满足等间距变化的要求，且满足分组要求。（4分） 优点：充分利用数据；消除部分定值系统误差 四、计算题（20分，每题10分） 1、用1/50游标卡尺，测得某金属板的长和宽数据如下表所示，求金属板的面 解：（1）金属块长度平均值：)(02.10mm L = 长度不确定度： )(01.03/02.0mm u L == 金属块长度为：mm L 01.002.10±= %10.0=B （2分） （2）金属块宽度平均值：)(05.4mm d = 宽度不确定度： )(01.03/02.0mm u d == 金属块宽度是：mm d 01.005.4±= %20.0=B （2分） （3）面积最佳估计值：258.40mm d L S =?= 不确定度：2222222 221.0mm L d d s L s d L d L S =+=??? ????+??? ????=σσσσσ 相对百分误差：B ＝%100?S s σ=0.25％ （4分） （4）结果表达：21.06.40mm S ±= B ＝0.25％ （2分） 注：注意有效数字位数，有误者酌情扣 5、测量中的千分尺的零点误差属于已定系统误差；米尺刻度不均匀的误差属于未
第一次练习 教学要求：熟练掌握Matlab 软件的基本命令和操作，会作二维、三维几何图形，能够用Matlab 软件解决微积分、线性代数与解析几何中的计算问题。 补充命令 vpa(x,n) 显示x 的n 位有效数字，教材102页 fplot(‘f(x)’,[a,b]) 函数作图命令，画出f(x)在区间[a,b]上的图形 在下面的题目中m 为你的学号的后3位（1-9班）或4位（10班以上） 计算30sin lim x mx mx x →-与3 sin lim x mx mx x →∞- 程序： syms x limit((1001*x-sin(1001*x))/x^3,x,0) 结果： 程序： syms x limit((1001*x-sin(1001*x))/x^3,x,inf) 结果： 0 cos 1000 x mx y e =，求''y 程序： syms x diff(exp(x)*cos(1001*x/1000),2) 结果： -2001/1000000*exp(x)*cos(1001/1000*x)-1001/500*exp(x)*sin(1001/1000*x)
计算 2 2 11 00 x y e dxdy +?? 程序： dblquad(@(x,y) exp(x.^2+y.^2),0,1,0,1) 结果： 计算4 2 2 4x dx m x +? 程序： syms x int(x^4/(1000^2+4*x^2)) 结果： (10)cos , x y e mx y =求 程序： syms x diff(exp(x)*cos(1000*x),10) 结果： 给出 0x =的泰勒展式(最高次幂为4). 程序： syms x taylor(sqrt(1001/1000+x),5) 结果： Fibonacci 数列{}n x 的定义是121,1x x ==， 12,(3,4,)n n n x x x n --=+=L 用循环语句编程给出该数列的前20项（要求将结果用向量的形式给出）。 程序： x=[1,1]; for n=3:20 x(n)=x(n-1)+x(n-2); end x 结果： Columns 1 through 10 1 1 2 3 5 8 13 21 3 4 5 5 Columns 11 through 20 89 144 233 377 610 987 1597 2584 4181 6765
实验8 线性规划 实验目的： 1）掌握用matlab优化工具箱解线性规划的方法 2）练习建立实际问题的线性规划模型 实验内容： 6，信息如下表所示。市政证券的收益可以免税，其他证券的收益需按50%的税率纳税。此 （1）政府及代办机构的证券总共至少要购进400万元； （2）所购证券的平均信用等级不超过1.4 （3）所购证券的平均到期年限不超过5年。 ①若该经理有1000万元资金，应该如何投资？ ②如果能够以2.75%的李律借到不超过100万元资金，该经理应如何操作？ ③在1000万元资金的情况下，若证券A的税前收益增加为4.5%，投资应否改变？若证券C的税前收益减少为4.8%，投资应否改变？ 解： ① 问题的建模： 令总收益为z，每项投资各为x1,x2,x3,x4,x5 则有： z=0.043x1+0.054*0.5x2+0.050*0.5x3+0.044*0.5x4+0.045x5 约束为： x2+x3+x4≥400 0.6x1+0.6x2-0.4x3-0.4x4+3.6x5≤0 (平均信用等级不超过1.4) 4x1+10x2-x3-2x4-3x5≤0 （平均到期年限不超过5年） x1+x2+x3+x4+x5=1000 x1,x2,x3,x4,x5≥0 模型的求解： matlab代码如下： c=[0.043 0.054*0.5 0.050*0.5 0.044*0.5 0.045]; A1=[0 -1 -1 -1 0;0.6 0.6 -0.4 -0.4 3.6;4 10 -1 -2 -3]; b1=[-400 0 0] A2=[1 1 1 1 1]; b2=[1000]; v1=[0 0 0 0 0];
《管理运筹学》实验报告实验日期： 2016年 04月 21日—— 2016 年 05 月 18 日
3.在点击“新建”按钮以后，按软件的要求输入目标函数个数和约束条件个数，输入目标函数级约束条件的歌变量的系数和b值，并选择好“≤”、“≥”或“＝”，如图二所示，最后点击解决
4.注意事项： （1）输入的系数可以是整数、小数，但不能是分数，要把分数化为小数再输入。（2）输入前要合并同类项。 当约束条件输入完毕后，请点击“解决”按钮，屏幕上讲显现线性规划问题的结果，如图所示
5.输出结果如下
5.课后习题： 一、P31习题1 某家具公司生产甲、乙两种型号的组合柜，每种组合柜需要两种工艺（制白坯和油漆）.甲型号组合柜需要制白坯6工时，油漆8工时：乙型号组合柜需要制白坯12工时，油漆4工时.已知制白坯工艺的生产能力为120工时/天，油漆工艺的生产能力为64工时/天，甲型号组合柜单位利润200元，乙型号组合柜单位利润为240元. 约束条件： 问题： （1）甲、乙两种柜的日产量是多少？这时最大利润是多少？ 答：由实验过程中的输出结果得甲组合柜的日产量是4个，乙的事8个。 . 0,0,6448,120126;240200 z max ≥≥≤+≤++=y x y x y x y x
（2）图中的对偶价格13.333的含义是什么？ 答： 对偶价格13.333的含义是约束条件2中，每增加一个工时的油漆工作，利润会增加13.33元。 （3）对图中的常数项围的上、下限的含义给予具体说明，并阐述如何使用这些信息。 答：当约束条件1的常数项在48~192围变化，且其他约束条件不变时，约束条件1的对偶价格不变，仍为15.56；当约束条件2的常数项在40~180围变化，而其他约束条件的常数项不变时，约束条件2的对偶价格不然，仍为13.333。 （4）若甲组合柜的利润变为300，最优解不变？为什么？ 答：目标函数的最优值会变，因为甲组合柜的利润增加，所以总利润和对偶价格增加；甲、乙的工艺耗时不变，所以甲、乙的生产安排不变。 二、学号题 约束条件： 无约束条件 （学号）学号43214321432143214321 0 0,30 9991285376)(53432max x x x x x x x x x x x x x x x x x x x x z ≤≥≤-+-+≥-+-+=-++-+++=??????????????-≥?-?-?-?-?-7606165060~5154050~414 )30(40~313)20(30~21210 20~11 10~1）（学号）（学号）（学号学号学号）（学号不变学号规则
《数学实验》第一次上机实验 1． 设有分块矩阵?? ? ???= ????22322333S O R E A ，其中E,R,O,S 分别为单位阵、随机阵、零阵和对角阵，试通过数值计算验证?? ????+= 22 S 0RS R E A 。 程序及结果： E=eye(3); %创建单位矩阵E% R=rand(3,2); %创建随机矩阵R% O=zeros(2,3); %创建0矩阵% S=diag(1:2); %创建对角矩阵% A=[E,R;O,S]; %创建A 矩阵% B=[E,(R+R*S);zeros(2,3),S^2] %计算等号右边的值% A^2 %计算等号左边的值% 运行结果： B = 1.00 0 0 1.63 2.74 0 1.00 0 1.81 1.90 0 0 1.00 0.25 0.29 0 0 0 1.00 0 0 0 0 0 4.00 ans = 1.00 0 0 1.63 2.74 0 1.00 0 1.81 1.90 0 0 1.00 0.25 0.29 0 0 0 1.00 0 0 0 0 0 4.00 2．某零售店有9种商品的单件进价（元）、售价（元）及一周的销量如表1.1，问哪种商品的利润最大，哪种商品的利润最小；按收入由小到大，列出所有商品及其收入；求这一周该10种商品的总收入和总利润。 表1.1 1)程序： a=[7.15 8.25 3.20 10.30 6.68 12.03 16.85 17.51 9.30]; b=[11.10 15.00 6.00 16.25 9.90 18.25 20.80 24.15 15.50]; c=[568 1205 753 580 395 2104 1538 810 694];
线性规划综合性实验参考选题 1.某工厂生产A、B两种产品，均需经过两道工序，每生产一吨产品A需要经第一道工序加工2小时，第二道工序加工3小时；每生产一吨产品B需要经第一道工序加工3小时，第二道工序加工4小时。可供利用的第一道工序为12小时，第二道工序为24小时。生产产品B的同时产出副产品C，每生产一吨产品B，可同时得到2吨产品C而毋需外加任何费用；副产品C一部分可以盈利，剩下的只能报废。出售产品A每吨能盈利400元、产品B每吨能盈利1000元，每销售一吨副产品C能盈利300元，而剩余要报废的则每吨损失200元。经市场预测，在计划期内产品C最大销量为5吨。 根据以上资料该工厂应如何制定生产方案，使工厂总的利润最大。 2.某厂接受了一批加工定货，客户要求加工100套钢架，每套由长2.9米、2.1米和1.5米的圆钢各一根组成。现在仅有一批长7.4米的棒料毛坯，问应如何下料，使所用的棒料根数最少？ 3.某公司在5年内考虑下列投资，已知：项目A可从第一年至第四年的年初投资，并于次年末收回本利共115%；项目B在第三年的年初投资，到第五年的年末收回本利135%,但规定投资额不能超过4万元；项目C在第二年的年初投资，到第五年的年末收回本利145%,但规定投资额不能超过3万元；项目D每年年初购买债券，年底归还，利息是0.06。公司现有资金10万元，问如何投资，才能使第五年年末拥有的资金最多？ 4.某企业在今后三年内有四种投资机会。第一种是在三年内每年年初投资，年底可回收本利和120％；第二种是在第一年年初投资，第二年年底可回收本利和150％，但该项投资不得超过2万元；第三种是在第二年年初投资，第三年年底回收本利和160％，但该项投资不得超过1.5万元；第四种是在第三年年初投资，该年年底可回收本利和140％，该项投资不得超过1万元。现在该企业准备拿出3万元资金，问如何制订投资计划，使到第三年年末本利和最大？ 5. D&D Corporation是一家专门从事艺术品买卖业务的公司。最近，D&D以低价收购了AT&T，Bell，Cisco，Dell，Epson公司的一些艺术品。这些艺术品可分为五类，不妨称其为A类，B类，C类，D类和E类。在D&D的广告宣传下，很多顾客来D&D购买这些艺术品，每个顾客都给D&D留下了要求购买的艺术品的数量，并提供了愿意出的价格。有关数据资料如下：设A类，B类，C类，D类和E类艺术品数量分别为3 件、3件、3件、1件和1件；设有5个顾客分别为Alan、Betty、Carl、David和Elton，他们需要艺术品的最多数量分别为5件、5件、2件、1件和1件。顾客Alan对五类艺术品愿意出的价格分别为10,10,10,30,50；顾客Betty对五类艺术品愿意出的价格分别为20,5,18,40,20；顾客Carl对五类艺术品愿意出的价格分别为15,20,20,20,20；顾客David对五类艺术品愿意出的价格分别为40,40,40,60,60；顾客Elton 对五类艺术品愿意出的价格分别为25,25,25,55,55. 现在任命你为D&D的销售部经理，要求你制定一个艺术品销售方案(即向上述五位顾客如何销售艺术品)，将所有艺术品全部售出，并使D&D的收入最大。 6.某公司有钢材、铝材、铜材1200吨，800吨和650吨，拟调往物资紧张的地区甲、乙、丙。已知甲、乙、丙对上述物资的总需求为：900吨，800吨和1000吨，各种物资在各地销售每吨的获利如下表所示。
大学物理实验试题（一） 一、单项选择题（每小题3分，共10小题） （1）.在光栅测量波长的实验中，所用的实验方法是[ ] （A）模拟法（B）干涉法（C）稳态法（D）补偿法 （2）.用箱式惠斯登电桥测电阻时，若被测电阻值约为4700欧姆，则倍率选[ ] （A）0.01 （B） 0.1 （C）10 （D） 1 （3）.用某尺子对一物体的长度进行15次重复测量，计算得A类不确定度为0.01mm，B类不确定度是0.6mm，如果用该尺子测量类似长度，应选择的合理测量次数为 （A）1次（B）6次（C）15次（D） 30次 （4）.用惠斯登电桥测电阻时，如果出现下列情况，试选择出仍能正常测 量的情况[ ] （A）有一个桥臂电阻恒为零（B）有一个桥臂电阻恒为无穷大 （C）检流计支路不通（断线）（D）电源与检流计位置互换 （5）.研究二极管伏安特性曲线时，正确的接线方法是[ ] （A）测量正向伏安特性曲线时用外接法；测量反向伏安特性曲线时用内接法（B）测量正向伏安特性曲线时用内接法；测量反向伏安特性曲线时用外接法（C）测量正向伏安特性曲线时用内接法；测量反向伏安特性曲线时用内接法（D）测量正向伏安特性曲线时用外接法；测量反向伏安特性曲线时用外接法（6）.在测量钢丝的杨氏模量实验中，预加1Kg砝码的目的是[ ] （A）消除摩擦力（B）使系统稳定 （C）拉直钢丝（D）增大钢丝伸长量 （7）.调节气垫导轨水平时发现在滑块运动方向上不水平，应该[ ] （A）只调节单脚螺钉（B）先调节单脚螺钉再调节双脚螺钉（C）只调节双脚螺钉（D）先调节双脚螺钉再调节单脚螺钉（8）.示波管的主要组成部分包括[ ] （A）磁聚集系统、偏转系统、显示屏（B）电子枪、偏转系统、显示屏（C）电聚集系统、偏转系统、显示屏（D）控制极、偏转系统、显示屏（9）.分光计设计了两个角游标是为了消除[ ]
实验二___线性规划灵敏度分析
实验二线性规划模型及灵敏度分析 （一）实验目的：掌握使用Excel软件进行灵敏度分析的操作方法。 （二）实验内容和要求：用Excel软件完成案例。 （三）实例操作： （1）建立电子表格模型； （2）使用Excel规划求解功能求解问题并生成“敏感性报告”； （3）结果分析：哪些问题可以直接利用“敏感性报告”中的信息求解，哪些问题需要重新规划求解，并对结果提出你的看法； （4）在Word文档中书写实验报告，包括线性规划模型、电子表格模型、敏感性报告和结果分析等。 案例1 市场调查问题 某市场调查公司受某厂的委托，调查消费者对某种新产品的了解和反应情况。该厂对市场调查公司提出了以下要求： （1）共对500个家庭进行调查；
（2）在被调查家庭中，至少有200个是没有孩子的家庭，同时至少有200个是有孩子的家庭； （3）至少对300个被调查家庭采用问卷式书面调查，对其余家庭可采用口头调查； （4）在有孩子的被调查家庭中，至少对50%的家庭采用问卷式书面调查； （5）在没有孩子的被调查家庭中，至少对60%的家庭采用问卷式书面调查。 对不同家庭采用不同调查方式的费用如下表所示： 市场调查费用表 家庭类型调查费用（元） 问卷式书面调查口头调查 有孩子的家庭50 30 没有孩子的家庭40 25 问：市场调查公司应如何进行调查，使得在
满足厂方要求的条件下，使得总调查费用最少？ 案例2 经理会议建议的分析 某公司生产三种产品A1，A2，A3，它们在B1，B2两种设备上加工，并耗用C1，C2两种原材料，已知生产单位产品耗用的工时和原材料以及设备和原材料的每天最多可使用量如下表所示： 生产三种产品的有关数据 资源产品A1 产品A2 产品A3 每天最多可使用量 设备B1(min) 1 2 1 430 设备B2(min) 3 0 2 460 原料C1(kg) 1 4 0 420 原料C2(kg) 1 1 1 300 每件利润(元) 30 20 50
实验三 线性规划 程序： linprog c=[-0.4 -0.28 -0.32 -0.72 -0.64 -0.6]; A=[0.01 0.01 0.01 0.03 0.03 0.03;0.02 0 0 0.05 0 0;0 0.02 0 0 0.05 0;0 0 0.03 0 0 0.08]; b=[850;700;100;900]; Aeq=[]; beq=[]; vlb=[0;0;0;0;0;0]; vub=[]; [x,fval]=linprog(c,A,b,Aeq,beq,vlb,vub) Exam5: function f=fun3(x); f=-x(1)-2*x(2)+(1/2)*x(1)^2+(1/2)*x(2)^2 实验目的 2、掌握用数学软件包求解线性规划问题。 1、了解线性规划的基本内容。 例1 max 6543216.064.072.032.028.04.0x x x x x x z +++++= 85003.003.003.001.001.001.0..654321≤+++++x x x x x x t s 70005.002.041≤+x x 10005.002.052≤+x x 90008.003.063≤+x x 6,2,10 =≥j x j
x0=[1;1]; A=[2 3 ;1 4]; b=[6;5]; Aeq=[];beq=[]; VLB=[0;0]; VUB=[]; [x,fval]=fmincon('fun3',x0,A,b,Aeq,beq,VLB,VUB) 书 求下列非线性规划 2221232212322 1232 12223123min 8020 ..2023,,0x x x x x x x x x s t x x x x x x x +++?-+≥?++≤??--+=??+=? ?≥? 在Matlab 2013软件中输入如下程序： （i ）编写M 文件fun1.m 定义目标函数 function f=fun1(x); f=sum(x.^2)+8; （ii ）编写M 文件fun2.m 定义非线性约束条件 function [g,h]=fun2(x); g=[-x(1)^2+x(2)-x(3)^2 x(1)+x(2)^2+x(3)^3-20]; %非线性不等式约束 h=[-x(1)-x(2)^2+2 x(2)+2*x(3)^2-3]; %非线性等式约束 （iii ）编写主程序文件example2.m 如下： options=optimset('largescale','off'); [x,y]=fmincon('fun1',rand(3,1),[],[],[],[],zeros(3,1),[], ... 'fun2', options) 就可以求得当1230.5522 1.2033,,0.9478x x x ===时，最小值y =10.6511。 4. 选修课的策略 决策目标为选修的课程总数最少，即 921min x x x +++ 约束条件： (1) 满足课程要求：（至少２门数学课程，３门运筹学课程和２门计算机课程）
最优化算法实验指导书 1.线性规划求解 1.1 生产销售计划 问题 一奶制品加工厂用牛奶生产A 1、A 2两种普通奶制品，以及B 1、B 2两种高级奶制品，分别是由A 1、A 2深加工开发得到的，已知每1桶牛奶可以在甲类设备上用12h 加工成3kg A 1，或者在乙类设备上用8h 加工成4kg A 2；深加工时，用2h 并花1.5元加工费，可将1kg A 1加工成0.8kg B 1，也可将1kg A 2加工成0.75kg B 2，根据市场需求，生产的4种奶制品全部能售出，且每公斤A 1、A 2、 B 1、B 2获利分别为12元、8元、22元、16元。 现在加工厂每天能得到50桶牛奶的供应，每天正式工人总的劳动时间最多为480h ，并且乙类设备和深加工设备的加工能力没有限制，但甲类设备的数量相对较少，每天至多能加工100kg A 1，试为该厂制定一个生产销售计划，使每天的净利润最大，并讨论以下问题： （1）若投资15元可以增加供应1桶牛奶，应否作这项投资； （2）若可以聘用临时工人以增加劳动时间，支付给临时工人的工资最多是每小时几 元？ （3）如果B 1、B 2的获利经常有10%的波动，波动后是否需要制定新的生产销售计划？ 模型 这是一个有约束的优化问题，其模型应包含决策变量、目标函数和约束条件。 决策变量用以表述生产销售计划，它并不是唯一的，设A 1、A 2、 B 1、B 2每天的销售量分别为1234,,,x x x x （kg ），34,x x 也是B 1、B 2的产量，设工厂用5x （kg ）A 1加工B 1，6x （kg ）A 2加工B 2（增设决策变量5x 、6x 可以使模型表达更清晰）。 目标函数是工厂每天的净利润z ，即A 1、A 2、 B 1、B 2的获利之和扣除深加工费，容易写出1234561282216 1.5 1.5z x x x x x x =+++--(元)。 约束条件 原料供应：A 1每天的产量为15x x +（kg ），用牛奶13()/3x x +（桶），A 2的每天产量为26x x +（kg ），用牛奶26()/4x x +（桶），二者之和不得超过每天的供应量50（桶）。 劳动时间：每天生产A 1、A 2的时间分别为154()x x +和262()x x +，加工B 1、B 2的时间分别为52x 和62x ，二者之和不得超过总的劳动时间480h 。 设备能力：A 1每天的产量15x x +，不得超过甲类设备的加工能力100（kg ）。 加工约束：1（kg ）A 1加工成0.8（kg ）B 1，故350.8x x =；类似的460.75x x =。 非负约束：123456,,,,,x x x x x x 均为非负。 由此得如下基本模型： 123456max 1282216 1.5 1.5z x x x x x x =+++--
实验5 线性规划 分1 黄浩 43 一、实验目的 1.掌握用MATLAB工具箱求解线性规划的方法 2.练习建立实际问题的线性规划模型 二、实验内容 1.《数学实验》第二版（问题6） 问题叙述： 某银行经理计划用一笔资金进行有价证券的投资，可供购进的证券以及其信用等级、到期年限、收益如下表所示。按照规定，市政证券的收益可以免税，其他证券的收益需按50%的税率纳税。此外还有如下限制： (1).政府及代办机构的证券总共至少要购进400万元； (2).所购证券的平均信用等级不超过1.4(信用等级数字越小，信用程度越高)； (3).所购证券的平均到期年限不超过5年 I.若该经理有1000万元资金，该如何投资？ II.如果能够以2.75%的利率借到不超过100万元资金，该经理应如何操作？ III.在1000万元资金情况下，若证券A的税前收益增加为4.5%，投资应否改变？若证券C的税前收益减少为4.8%，投资应否改变？ 模型转换及实验过程： I. 设经理对于上述五种证券A、B、C、D、E的投资额分别为：、、、、(万
元)，全部到期后的总收益为z万元。 由题目中的已知条件，可以列出约束条件为： 而决策变量的上下界约束为： 目标函数 将上述条件转变为matlab的要求形式： 使用matlab解上述的线性规划问题（程序见四.1），并整理成表格： 得出结论： 当经理对A、B、C、D、E五种证券分别投资218.18、0、736.36、0、45.45万元时，在全部收回时可得到29.836万元的税后收益，而且这种投资方式所得收益是最大的。 讨论： 尝试输出该约束条件下的拉格朗日乘子： 该乘子表示，第一个约束条件对目标函数的取值不起作用，而剩余三个约束条件取严格等号的时候，目标函数达到最优解。下面验证之： 由解得的x值，代入四个约束条件中，得：