数学实验六
06实验六 空间图形的画法
图6-4
绘制二元函数图形也可用简捷绘制的ezsurf指令，它的使 用格式为: ezsurf(f(x,y),[a,b,u,v]) 即可绘制函数在区域[a,b]×[u,v]上的图形。当省略区域 时，默认区间是[-2 ,2 ]×[- 2 , 2 ]。例如输入: ezsurf('x*exp(-x^2-y^2)') 则输出如图6-5所示。
图6-12
【例7】可以证明:函数z=xy的图形是双曲抛物面。在区 域-2≤x≤2，-2≤y≤2上作出它的图形。
输入: x=-2:0.1:2; y=-2:0.1:2; [xx,yy]=meshgrid(x,y); zz=xx.*yy; surf(xx,yy,zz) 输出如图6-13所示。
图6-13
例如，画出曲面 z x y 的图形。输入: x=-2:0.1:2; y=-2:0.1:2; [x,y]=meshgrid(x,y); z=x.^2+y.^2; surf(x,y,z) z x 2 y 2，见图6-2。 得到曲面
2 2
图6-2
执行下面的程序: x=-2:0.015:2; y=-2:0.015:2; [x,y]=meshgrid(x,y); z=x.^2+y.^2; i=find(x.^2+y.^2>4); z(i)=NaN; surf(x,y,z) 同样得到曲面(见图6-3)。 由于自变量的取值范围不同，图形也不同。不过，后者比 较好地反映了旋转曲面的特点，因此是常用的方法。
图6-3
又如，参数方程: x 2sin cos , y 2sin sin , z 2cos 是以原 点为中心、2为半径的球面，其中 0 , 0 2 因此只要输入: t=0:0.1:pi; r=0:0.1:2*pi; [r,t]=meshgrid(r,t); x=2*sin(t).*cos(r); y=2*sin(t).*sin(r); z=2*cos(t); surf(x,y,z) 2 2 2 2 便作出了方程为 x y z 2 的球面(见图6-4)。
实验六 古典密码与破译
定义在 Z26上的方阵 A 模 26 可逆的充要条件是：
det(A) 不能被 2 和 13 整除
模 26 可逆
Z26 中具有模 26 倒数的整数及其模 26 倒数表2
a a-1 1 1 3 5 7 9 11 15 17 19 21 23 25 9 21 15 3 19 7 23 11 5 17 25
N O P Q R S T U V W X Y Z 14 15 16 17 18 19 20 21 22 23 24 25 0
Hill2 解密举例
① 将密文字母分组，通过查表列出每组字母对应的向量 ② 求出加密矩阵 A 的 模 26 逆矩阵
1 8 B 0 9
③ 用 B 左乘每组密文字母组成的向量，然后再反查字母
Hill2 加密举例
反查字母表值得每个向量对应的字母组为：
16 , 12
1 15 , , 12 20
20 20
PL
AL
OT
TT
HDSDSXX
Hill2 加密
K L M
PLALOTTT
A B C D E F G H I J
③ 令 = B mod(m) ，由 的分量反查字母表值表，
得到相应的明文字母
Hill2 解密举例
甲方收到乙方（己方）的一个密文信息，内容为： WKVACPEAOCIXGWIZUROQWAB ALOHDKCEAFCLWWCVLEMIMCC 按照甲方与乙方的约定，他们之间采用 Hill2密码，密钥 1 2 为 A ，字母表值见下表，问这段密文的原文 0 3 是什么？ A 1 B 2 C 3 D 4 E 5 F 6 G 7 H 8 I 9 J K L M 10 11 12 13
数学实验完整
数学实验一1， 立方曲线y=x 3>> x=-10:0.1:10; >> y=x.^3; >> plot(x,y)-10-8-6-4-2246810-1000-800-600-400-200020040060080010002．立方抛物线y=x 1/3>> x=-10:0.1:10; >> y=x.^(1/3); >> plot(x,y);-10-8-6-4-2024681000.511.522.53,高斯曲线y=2x e以参数方程表示的曲线；>> x=-10:0.1:10; >> y=exp(-x.^2); >> plot(x,y);>>-10-8-6-4-224681000.10.20.30.40.50.60.70.80.91以参数方程表示的曲线： 4，奈尔抛物线 x=t 3y=t 2(y=x 2/3)>> x=-10:0.1:10; >> y=x.^(2/3); >> plot(x,y);-10-8-6-4-20246810-3-2-10123455,半立方抛物线 x=t 2y=t 3(y 2=x 3)>> x=-10:0.1:10; >> y=x.^(3/2); >> plot(x,y);-10-8-6-4-20246810-5051015202530356,笛卡尔曲线 x=3at 2) , y=3at 2/(1+t 2) (x 3+y 3-3axy=0)>> t=-10:0.1:10; >> a=8;>> x=3*a*t./(1+t.^2);y=3*a*t.^2./(1+t.^2); >> plot(x,y);-15-10-505101505101520257,蔓叶线 x=at 2/(1+t 2),y=at 3/(1+t 2) (y 2=x 3/(a-x)) t=-10:0.1:10; >> a=2;>> x=a*t.^2./(1+t.^2);y=a*t.^3./(1+t.^2); >> plot(x,y);00.20.40.60.81 1.2 1.4 1.6 1.82-20-15-10-5051015208,摆线 x=a(t-sint), y=b(1-cost) >> a=10; >> b=15; >> t=-10:0.1:10;>> x=a*(t-sin(t));y=a*(1-cos(t)); >> plot(x,y);-150-100-5050100150024681012141618209,内摆线（星形线） x=acos 3t, y=asin 3t (x 2/3+y2/3=a2/3)>> t=-10:0.1:10; >> a=5;>> x=a*(cos(t)).^3; y=a*(sin(t)).^3; >> plot(x,y);-5-4-3-2-101234510,圆的渐伸线（渐开线）x=a(cost+tsint) , y=a(sint- tcost )t=-10:0.1:10;>> a=5;>> x=a.*(cos(t)+t.*sin(t));y=a.*(sin(t)-t.*cos(t));>> plot(x,y);-40-30-20-1010203040-50-40-30-20-100102030405011,空间螺线 x=acost , y=bsint , z=ct 。
小学数学 六年级上 北师大数学实验班(六上)测试(一)
北师大数学实验班（六上）测试（一）学号：_______ 班级：_________ 姓名：_________ 成绩：_________一、直接写出得数。
（每题1分，共8分）3.14×3= 15×3.14= 8×3.14= 52×3.14=20×3.14= π×22= 72×π= 162=二、填空。
（每题2分，共26分）1、在同一个圆里，所有的直径都（），所有的半径都（），半径的长度是直径的（）。
2、一个圆的半径是5厘米，它的直径是（），周长是（），面积是（）。
3、用圆规画一个周长为31.4厘米的圆，两脚尖张开的距离是（）厘米。
4、一个圆的半径扩大2倍，它的直径扩大（）倍，周长扩大（）倍，面积扩大（）倍。
5、一座石英挂钟，它的分钟长8厘米，1小时分针尖端走过（）厘米。
6、圆是轴对称图形，任何一条（）都是圆的对称轴。
圆有（）条对称轴。
三、判断。
（每题2分，共8分）1、同一个圆的直径都相等。
（）2、半圆的周长等于整个圆周长的一半。
（）3、π是一个无限不循环小数。
（）4、半径是2厘米的圆，它的周长与面积相等。
（）四、选择。
（每题2分，共6分）1、用12米长的绳子围成三种图形，其中面积最小的是（）。
A、圆B、长方形C、正方形2、把一根长10米的铜丝，在一个圆盘上绕了3圈，还多0.58米，这个圆盘的半径是（ ）。
A 、0.5米B 、1米C 、1.5米3、大圆周长是小圆周长的4倍，小圆直径是大圆直径的（ ）。
A 、4倍B 、41 C 、81 D 、8倍 五、求下面各圆的周长和面积。
（每题6分，共12分）d=4.2分米 C=28.26米六、解决问题。
（每题10分，共40分）1、学校操场的跑道如下图所示，跑道两边是两个半圆，沿跑道跑一周的长度是多少米？2、育苗小学有一个圆形景点，它的直径是10米，景点中央有一个直径是5米的圆形喷水池，环绕水池的部分种植花草，种花草部分的面积是多少平方米？3、在一个周长是80厘米的正方形木板上，锯下一个最大的圆，这个圆的面积是多少平方厘米？4、计算下图阴影部分的面积。
数学实验六插值问题
实验6 插值问题\一、实验目的：1．理解插值方法的基本原理；2．掌握用MATLAB进行各种插值的方法；3．能够利用插值方法解决实际问题。
二、实验内容：1．选择n（5≤n≤10）个节点，利用不同的插值方法计算下列函数在m（50≤m≤100）个插值点的函数值。
通过数值和图形输出，将插值结果与精确值比较。
适当增大n，再作比较，并作出分析结论。
(1) y=sin x , (0≤x ≤2pi);>> x=linspace(0,2*pi,5);>> y=sin(x);>> xi=linspace(0,2*pi,80);>> yi=sin(xi);>> y1=interp1(x,y,xi);>> y2=interp1(x,y,xi,'nearest');>> y3=interp1(x,y,xi,'spline');>> y4=interp1(x,y,xi,'cubic');>> plot(xi,yi,'o',xi,y1,'b',xi,y2,'g',xi,y3,'r',xi,y4,'y')(2) y=(1-x2)1/2 , (-1≤x ≤1);>> x=linspace(-1,1,5);y=(1-x).^(1/2);xi=linspace(-1,1 ,80);yi= (1-xi).^(1/2);y1=interp1(x,y,xi);y2=interp1(x,y,xi,'nearest');y3=interp1(x,y,xi,'spline');y4=interp1(x,y,xi,'cubic');>> plot(xi,yi,'o',xi,y1,'b',xi,y2,'g',xi,y3,'r',xi,y4,'y')(3) y=cos10x , (-2≤x ≤2);x=linspace(-2,2,5);y=cos(10*x);xi=linspace(-2,2,80);yi= cos(10*xi);y1=interp1(x,y,xi);y2=interp1(x,y,xi,'nearest');y3=interp1(x,y,xi,'spline');y4=interp1(x,y,xi,'cubic');plot(xi,yi,'o',xi,y1,'b',xi,y2,'g',xi,y3,'r',xi,y4,'y')(4)y=exp(-x2) , (-2≤x ≤2)x=linspace(-2,2,5);y=exp(-(x.^2));xi=linspace(-2,2,80);yi= exp(-(xi.^2));y1=interp1(x,y,xi);y2=interp1(x,y,xi,'nearest');y3=interp1(x,y,xi,'spline');y4=interp1(x,y,xi,'cubic');plot(xi,yi,'o',xi,y1,'b',xi,y2,'g',xi,y3,'r',xi,y4,'y')2．为了计算瑞士的国土面积，以由西向东方向为x轴，由南向北方向为y 轴，最西边界点到最东边界点在x轴上对应的区间为[0, 151]，将此区间适当划分为若干段，在每个分点的y方向上测出南边界点y1和北边界点y2的坐标数据（精确值为41288km2）。
清华大学数学实验报告6
实验六非线性方程求解实验目的1. 掌握用matlab软件求解非线性方程和方程组的基本用法, 并对结果做初步分析.2. 练习用非线性方程和方程组建立实际问题的模型并进行求解.实验内容题目3（1）小张夫妇以按揭方式贷款买了1 套价值20 万元的房子，首付了5 万元，每月还款1000 元，15 年还清。
问贷款利率是多少？（2）某人欲贷款50 万元购房，他咨询了两家银行，第一家银行开出的条件是每月还4500 元，15 年还清；第二家银行开出的条件是每年还450000 元，20 年还清。
从利率方面看，哪家银行较优惠（简单地假设年利率=月利率×12）？建立模型：设房价为b，首付款为b0，银行按照月利率(复利)来计算，月利率为r，月付款(月末支付)为a，共需要支付的月数为n。
根据经济学中资金的时间价值概念，可以得到：房价在n个月之后的实际价值为：b(1+r)n按揭购房期间交的所有款项在第n个月末的实际价值为：b0(1+r)n+a(1+r)n−1+(1+r)n−2+⋯+1=b0(1+r)n+a×(1+r)n−1由于在第n个月末还清了贷款，因此上述两个时间价值相等，则得到下面的关系式，即为解答此问题的方程：b(1+r)n=b0(1+r)n+a×(1+r)n−1即：(b−b0)(1+r)n−a×(1+r)n−1=0(1)代入已知条件：b=200000，b0=50000，a=1000，n=180，利用MATLAB解此非线性方程，经过简单的估测之后，给定初始值为r0=0.001，得到结果为：r=0.0020812，即贷款月利率为0.20812%。
（2）I.第一家银行相应的已知条件为：b=500000，b0=0，a=4500，n=180，利用MATLAB计算，经过简单的估测之后，给定初始值为r0=0.005，得到结果为：r=0.0058508，即这家银行的贷款月利率为0.58508%。
人教版数学六年级下册实验报告
第（）组
蹬一圈，自行车走（）
实验报告（一）
第（）组
蹬一圈，自行车走（）
实验报告（一）
第（）组
蹬一圈，自行车走（）
实验报告（一）
第（）组
蹬一圈，自行车走（）
实验报告（一）
第（）组
蹬一圈，自行车走（车走（）
实验报告（一）
第（）组
蹬一圈，自行车走（）
实验报告
第（）组
实验课题：蹬一圈，自行车能走多远？
1、车轮的直径是30CM。
2、前齿轮油（）个齿，后齿轮有（）个齿。
3、前齿轮的齿数×前齿轮转的圈数=后齿轮的齿数×（）
根据以上规律得出：
4、前齿轮转一圈，后齿轮转（）圈，后车轮转（）圈。
5、结论：
蹬一圈自行车走的路程（）×
6、蹬一圈自行车走（）CM。
实验报告（二）
第（）组
实验课题：蹬一圈，自行车能走多远？
1、车轮的直径是30CM。
2、前齿轮油（）个齿，后齿轮有（）个齿。
3、前齿轮的齿数×前齿轮转的圈数=后齿轮的齿数×（）
根据以上规律得出：
4、前齿轮转一圈，后齿轮转（）圈，后车轮转（）圈。
5、结论：
蹬一圈自行车走的路程（）×
6、蹬一圈自行车走（）CM。
实验小学六年级数学(1)
安丘市实验小学六年级数学第一单元测试题班级：姓名：等级：书写：一、填空。
1、如果下降5米，记作－5米，那么上升4米记作（）米；如果＋2千克表示增加2千克，那么－3千克表示（）。
2、二月份，妈妈在银行存入5000元，存折上应记作（）元。
三月一日妈妈又取出1000元，存折上应记作（）元。
3、＋8.7读作（），－25读作（）。
4、海平面的海拔高度记作0m，海拔高度为＋450米，表示（），海拔高度为－102米，表示（）。
5、如果把平均成绩记为0分，＋9分表示比平均成绩（），－18分表示（），比平均成绩少2分，记作（）。
6、数轴上所有的负数都在0的（）边，所有正数都在0的（）边。
7、在数轴上，从表示0的点出发，向右移动3个单位长度到A 点，A点表示的数是（）；从表示0的点出发向左移动6个单位长度到B点，B点表示的数是（）。
8、比较大小。
-10○0.1 -0.3○0.02 －25○-52-99○3－7○－5 1.5○520○－2.5 －3○3二、判断对错。
1、零上12℃（＋12℃）和零下12℃（－12℃）是两种相反意义的量。
（）2、0是正数。
（）3、数轴上左边的数比右边的数小。
（）4、死海低于海平面400米，记作＋400米。
（）5、在8.2、－4、0、6、－27中，负数有3个。
（）6、数轴是一条直线。
（）7、在我们学的数中，既没有最大的数也没有最小的数。
（）8、因为8大于6，所以-8大于-6。
（）9、在0和-1之间没有负数。
（）10、带有“-”的数一定是负数。
（）三、选择正确答案的序号填在括号里。
1、低于正常水位0.16米记为－0.16，高于正常水位0.02米记作（）。
A、＋0.02B、－0.02C、＋0.18D、－0.142、以明明家为起点，向东走为正，向西走为负。
如果明明从家走了＋30米，又走了－30米，这时明明离家的距离是（）米。
A、30B、－30C、60D、03、数轴上，－12在－18的（）边。
A、左B、右C、北D、无法确定4、规定10吨记为0吨，11吨记为＋1吨，则下列说法错误的是（）。
实验小学六年级上册数学期末测试卷及答案
六年级数学期末测试卷同学们，一学期结束了，你一定有许多收获，快来挑战吧，只要认真审题，仔细答卷，就能取得好成绩哦！一、填空 （第9题2分，其余每空1分，共25分） 1、%40 3915.172 ⨯=⨯=⨯=÷）（）（）（）（ 2、一张半圆形纸片的半径为1分米，它的周长为（ ）分米，要剪一个这样的半圆形，至少需要一张面积为（ ）平方分米的长方形纸片。
3、六年级(3)班有男生24人，女生30人，女生人数是男生人数的（ ）%，女生人数比男生人数多（ ）%，男生人数比女生人数少（ ）%。
4、出勤率、出粉率、发芽率、合格率中，不可能达到100%的是（ ）。
5、盐占盐水的113，那么盐与水的质量比是（ ）。
6、工程队修一段路，4天修了全长的52，这时还有30米没修，这条路全长（ ）米。
7、合唱队有男生25人，女生20人。
男生人数是女生的（ ）倍，女生人数与男生人数的最简单的整数比是（ ），女生人数占总人数的）（）（ 。
8、比80米多25%是（ ），（ ）的51等于24的32。
9、丽丽面向北站立，向右转400后所面对的方向是（ ）偏（ ）（ ）0；丁丁面向西站立，向左转400后所面对的方向是（ ）。
10、如果要画一个周长是12.56厘米的圆,，圆规两脚之间的距离应该是（ ）厘米，这个圆的面积是（ ）平方厘米。
11、用109千克小麦可以磨出43千克面粉，每千克小麦可以磨面粉（ ）千克，要磨1千克面粉需要小麦（ ）千克。
12、小敏把一个数除以32错算为乘以32，得出结果6。
那么正确的计算结果应是（ ）。
13、）（ =+++++641321161814121。
二、判断 （对的在括号里面画“√”，错的在括号里画“X"）（共5分） 1、两根同样长的铁管，第一根用去了53米，第二根用去53，第二根用去的多。
（ ） 2、假分数的倒数一定比这个假分数小。
（ ）3、小明和小军进行投篮练习，小明的命中率是65%，小军的命中率是70%，小军投中的次数多一些。
东明县刘楼镇实验学校六年级数学实践性作业
3.在设计、测量、整理等实践活动中感受数学与生活的密切联系，进一步提高
学习兴趣，发展自我反思能力。
设计方案
想一想，绘制家到学校的路线图前，要先做哪些方面的准备?我家到学校的 路线图中，需要收集哪些数据?如何收集这些数据?
绘制路线：
语言描述路线：
东明县刘楼镇实验学校六年级数学实践作业
活动内容 绘制家到学校的路线图 班级
姓名
1.通过“绘制家到学校的路线图”的实践操作活动，进一步理解并综合运用
图形位置、测量、数据收集等知识，积累“从头到尾思考问题的经验
活动目的
2.经历设计方案、动手实践、及时反思的活动过程，发展统筹规划和按方案操 作等综合实践能力。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
实验六1415100076陈田实验项目名称：优化实验时间：2016－５-27、201６-6—３实验地点：理学楼52５实验目得:1、掌握Matlaｂ优化工具箱得基本用法,利用优化工具包解线性规划与非线性规划得问题；对不同算法作初步分析、比较；2、练习实际问题得非线性最小二乘拟合.3、初步了解实际问题得线性规划与非线性规划得模型建立。
实验内容:1、教材2０5页习题2；funcｔion f=exａm０702fun（x)f=exｐ（x(1)）＊(4*x（１)、^２+2＊x（2)、^2＋４＊x（1)*ｘ（2)＋2＊ｘ（2）+1）；enｄM文件:formａtｓhoｒtｅx０＝［—1,１］;’—－-ｃase1：bfgs-——＇ｏpt１=optimseｔ('larges’，’off'，’mａｘfｕnevals',1００0）;［ｘ1，v1,exit1，out1]＝fminunc(＇ｅxａm0７02fun’,x０,ｏｐt1) pause＇—－—case2:dfｐ-－—＇fopt=ｏpｔimseｔ（oｐt1，'hｅssｕpdate'，’dｆp');[x２，ｖ2,exiｔ2,oｕt2］=ｆmiｎunc（'ｅｘam070２fｕn',x０，fopt）pausｅ’——-ｃａse３：steｅp－－-＇fopt=optimset（opt1,'heｓsupｄａte’，’stｅepdｅsc’);［x３,v3,eｘｉｔ3，out3]=fminｕnc(＇exam0702fｕn’，x０，fopt）pause＇＋+＋rｅｓｕlt ｏfｓolutions ++＋’ｓoｌutiｏns=［ｘ1；x2；x3];funｖalues=[v1；v2；v３］;iterationｓ=[oｕt1、funｃCouｎｔ；ｏut2、funcCｏunt;out3、fun ｃCoｕnｔ］；[soｌutions，ｆunvalｕes，iteratｉoｎｓ]运行结果：＞〉ｃ6２ans=－—-case1：bｆgs—--Ｌocａl mｉnｉｍuｍfｏund、Optiｍizationｐlｅteｄｂecause ｔhe siｚe of the gradiｅnt ｉs less thanthe defauｌt valｕｅof tｈe optｉmalｉｔy tolｅrａｎｃe、〈stopp ｉnｇｃｒitｅrｉa deｔaiｌs>x1＝5、0000e-0１－1、0０00e+00v１= 3、6609e—15exit1 ＝ 1out1 = ｉteratiｏｎs：8funcCount：6６stepｓize: 6、3３６1e-07lssteｐleｎgｔｈ: １firstoｒderｏpt: １、2284e－０7aｌgorｉthm：'qｕasi－newton＇message：'Ｌｏｃal mｉnimum ｆouｎd、…＇ans =－——case2:dfp—-－Ｌｏｃａｌmｉnimum ｆoｕｎd、Optiｍization ｐlｅｔｅd bｅｃauｓe the sizｅoｆｔhe gradｉeｎt is lｅｓs thaｎｔhe defａuｌｔｖalue of tｈe opｔiｍalitｙｔolerancｅ、＜sto ｐping ｃrｉterｉa details＞ｘ２= 5、0000e-01 －1、00０0e+０0v2 = 1、2４47e—14ｅxiｔ２= 1ouｔ2 = itｅｒaｔiｏnｓ: 8funｃＣouｎｔ: ５1steｐsｉzｅ：３、7463e-06lｓsteｐｌength: 1fｉrｓｔｏrderopt：5、8963e—07algoｒithm：'quasｉ—ｎewtｏn＇ｍeｓsａｇe: ＇Lｏcal miｎimuｍfｏunｄ、…’anｓ=-－－ｃａsｅ3：sｔeep—－－Local ｍinimum possiｂle、fminuｎc stoppedｂecaｕｓe ｔｈe siｚe ｏf tｈe curｒent steｐiｓless thanthe dｅfaｕｌｔvaｌue of theｓtep ｓizｅｔoｌerancｅ、〈ｓt ｏppｉng cｒiteria detａiｌs>x３= 5、0000e－01－1、0000e+00v3 = 9、0５30e—12ｅxit3 = 2out３= ｉｔeratｉons：49fuｎcＣouｎt: ３06stepsｉzｅ：1、1892ｅ-06lssteｐlengtｈ: 1、0０03e-01fｉｒstｏrｄeropt: 8、4759e－0６ａlｇｏriｔhｍ: 'ｑｕasｉ—neｗton'messａge: ’Local minｉmｕｍpｏsｓible、…'anｓ=+++rｅsult oｆsｏlutｉons ＋＋+ans=5、0000e-01-1、00０0ｅ+0０3、6609e-1５6、6000e+015、00００e—０1 -1、０000e+00 １、2447e—1４5、1０00e＋015、000０e-01-1、00０0e+００9、０530e—12 3、0600e+022、教材２05页习题４，要求：实现非线性最小二乘拟合即可,对“用GN与LM两种方法求解、对y作一扰动后观察迭代收敛变慢得情况”不做；fｕnction y=F(x，t,c）y=ｘ(1）+x（２)＊exｐ（—ｘ（4）*t)+ｘ(3)＊exｐ(－x（5)*t)－ｃ；ｅnｄM文件：ｘ0=［、5 1、5 －1 、0１、０2];ｔ＝zｅｒos（1，33);fｏｒi=1：33t（ｉ）=10＊(i－1）；ｅnｄc１=［、８4４、908 、９３2 、936 、９25 、908 、881、８50、81８、784 、751 、718、6８5 、658 、６２8 、６03 、580 、５58 、538 、522、506 、490 、478 、4６7 、457 、４４8 、438 、43１、424 、４20、41４、4１1 、406］;c2＝raｎｄ(1,3３)/1０—０、０５；＇--—caｓe1：LM-—－’opt１=ｏptimset(＇LarｇeScａle','oｆf','Maｘfuｎevａls＇，１000）;[ｘ１1,noｒm１1,res1１]=lsqｎoｎliｎ（’F’，ｘ0，［］，［],oｐt1，ｔ，c１)＇-－-casｅ2:GM-－-'opｔ2=optimset(opt1，'Ａlｇｏrithm＇,＇Ｌｅveｎberg—Marｑuaｒdt’)；［x12，norm1２，reｓ１２］=lsｑｎｏnlin(’F’,x０，［]，［]，opt2,t，ｃ１）’—－－cａse3：ＬM2--—’opt3＝oｐtiｍset(oｐt1，’Disp＇，'iter’);［ｘ21,norｍ２1，rｅs2１］=lsqnonlin('F'，x0,［］,［],opt３，t,c2）＇---casｅ4：ＧＭ2—－—'oｐｔ4=opｔimｓｅt（opt2,'Diｓｐ’，’itｅr'）；［x22,ｎorm22，ｒes2２]=lsqnｏｎlｉn（'F’,x0，[]，［],ｏｐt4，t，c2)运行结果：〉> c6４ans ＝——-case1：LＭ－－-x11=3、7541e—01 1、９358e＋０0 -1、４647e＋00 1、286８e—０2 ２、2123e—02norm１１=5、４649e-05aｎs =—-—caｓｅ２：ＧM－-－ｘ12 =3、7５41e－011、9３５8e＋０0 —1、464７e+0０1、2868e－0２2、212３e-02norm12 =5、４6４9e-05ans =－-—caｓe３:LM2－－—x2１=—５、8167e＋01 3、62１７ｅ-01 5、8127e+01 1、004９e－01 －3、651３ｅ—06norm2１=１、3721e-01anｓ=-——cａse４:ＧM2－-－x２2 ＝-3、90５6e+0１6、09４2e+01-2、1883ｅ+011、2573e-0４３、６１０２e－04ｎorm２2 =3、073３e-0２3、教材23８页习题１;M文件:c＝[1 -2 －1］;A1=［-1 12；1 1 -1;1 2 0；0 1 1];A2＝［0 1 —1];b1=[2４ 4 4］;ｂ2=[1］;v１＝［00 0];[x，f，eｘitflag,outpuｔ，lag]=ｌiｎｐrog(c，A1，b１,A２，b2，v1)lag、iｎｅqlｉn运行结果:〉>c81Oｐtimizatioｎteｒminatｅｄ、x ＝7、9９96e-011、6000ｅ+０05、９9９9e－0１ｆ=－3、００00ｅ+００ｅｘｉtfｌag ＝1outｐut＝iｔerａtiｏｎｓ: 6algorｉｔhm：'inｔerｉｏr-pｏinｔ—lｅgaｃy'cgiterａtions：0messaｇe: ’Ｏｐｔiｍizatiｏnｔｅrｍiｎated、'cｏｎstrｖiolatｉon：2、64２３e—14fｉｒstorderｏｐt：2、47７9e—13laｇ=ineqlin: [4x１dｏｕbｌe]eｑlin：1、００00e+00uｐper：[３x1double］lｏwer：［3x１doubｌｅ］>>lａg、ｉnｅqlinans＝1、0000e+00１、１569e－１4１、7４55ｅ－131、863６e-14由运行结果可知所求最大值为—3、00０0e＋00，起作用约束为第一个约束4、教材2３９页习题５。
M文件:ｃ1=[-20 -３0］；A=［1 ２；５ 4］;b＝［2０ 70]；v1=[００］；[x１,ｆ１,ｅxitｆlａｇ1,ｏutpｕt1,laｇ1］=linprｏg（c1,A，b，［],[］，v1）y1=ｌag１、ineqｌin运行结果:〉> c85Oｐtimｉzａtioｎtermｉnated、x1 =1、０000ｅ+01５、0００0e+0０f1=-３、5000ｅ+02ｅxitflａｇ1 =1outｐut1 =ｉｔeraｔions: 5ａlｇｏritｈm：’iｎterior-poiｎｔ—lｅgaｃｙ'ｃgiteｒatｉons：0mｅssaｇe：’Optｉmiｚaｔion ｔerｍｉnａted、＇ｃonstrvｉｏlatｉｏn：0firstｏｒderopt：６、4７31e—１1ｌag１=ineqlin: ［２x1 double］eqliｎ: ［０x1 double］upper：［２x１double］lｏwer: ［2x1 doublｅ]y１=1、16６7e+０１１、６66７e+０0由运行结果可知以甲乙分别2０、30元出售时安排甲乙分别生产１0、5件可使收入最大，当A、B增加1公斤时收入分别增加1１、667、1、６66７,考虑成本后由11、６６7-６=5、6６7,1、６667－2＝-1、6667＜0可知每进1公斤买Ａ原料可增加利润5、６６7元，买B原料则会减少利润－1、６６7元,故应全部买进A原料。