重庆大学数学实验第六次
重庆市实验中学2023-2024学年八年级数学第一学期期末达标检测试题【含解析】
重庆市实验中学2023-2024学年八年级数学第一学期期末达标检测试题检测试题注意事项1.考试结束后,请将本试卷和答题卡一并交回.2.答题前,请务必将自己的姓名、准考证号用0.5毫米黑色墨水的签字笔填写在试卷及答题卡的规定位置.3.请认真核对监考员在答题卡上所粘贴的条形码上的姓名、准考证号与本人是否相符.4.作答选择题,必须用2B 铅笔将答题卡上对应选项的方框涂满、涂黑;如需改动,请用橡皮擦干净后,再选涂其他答案.作答非选择题,必须用05毫米黑色墨水的签字笔在答题卡上的指定位置作答,在其他位置作答一律无效.5.如需作图,须用2B 铅笔绘、写清楚,线条、符号等须加黑、加粗.一、选择题(每题4分,共48分)1.如图为张晓亮的答卷,每个小题判断正确得20分,他的得分应是()A .100分B .80分C .60分D .40分2.甲种防腐药水含药30%,乙种防腐药水含药75%,现用这两种防腐药水配制含药50%的防腐药水18千克,两种药水各需要多少千克?设甲种药水需要x 千克,乙种药水需要y 千克,则所列方程组正确的是()A .1830%75%1850%x y x y +=⎧⎨+=⨯⎩B .1830%75%18x y x y +=⎧⎨+=⎩C .1875%30%1850%x y x y +=⎧⎨+=⨯⎩D .1875%30%18x y x y +=⎧⎨+=⎩3.下列线段长能构成三角形的是()A .3、4、8B .2、3、6C .5、6、11D .5、6、104.下列四组数据,能组成三角形的是()A .2,2,6B .3,4,5C .359,,D .5,8,135.如图,Rt ABC ∆中,90C ∠=︒,AC BC =,60ADC ∠=︒,则BAD ∠的度数等于()A .10︒B .15︒C .30°D .45︒6.如图,在ABC ∆中,AB 的垂直平分线交AB 于点D ,交BC 于点E .ABC ∆的周长为19,ACE ∆的周长为13,则AB 的长为()A .3B .6C .12D .167.已知一个等腰三角形的两边长a 、b 满足方程组23{3a b a b -=+=则此等腰三角形的周长为()A .5B .4C .3D .5或48.下列说法错误的是()A .0.350是精确到0.001的近似数B .3.80万是精确到百位的近似数C .近似数26.9与26.90表示的意义相同D .近似数2.20是由数a 四会五入得到的,那么数a 的取值范围是2.195 2.205a < 9.4的平方根是()A .4B .4±C .2±D .210.用科学记数法表示0.0000000052为()A .105210-⨯B .95.210-⨯C .105.210-⨯D .115.210-⨯11.下列语句不属于命题的是()A .直角都等于90°B .两点之间线段最短C .作线段ABD .若a=b ,则a 2=b 212.如图,点A ,D ,C ,F 在一条直线上,AB=DE ,∠A=∠EDF ,下列条件不能判定△ABC ≌△DEF 的是()A .AD=CFB .∠BCA=∠FC .∠B=∠ED .BC=EF二、填空题(每题4分,共24分)13.如图,在扇形BCD 中,∠BCD=150°,以点B 为圆心,BC 长为半径画弧交BD 于点A ,连接AC ,若BC=8,则图中阴影部分的面积为________14.已知点P (1﹣a ,a+2)关于y 轴的对称点在第二象限,则a 的取值范围是______.15.若一组数据2,3,4,5,x 的方差与另一组数据5,6,7,8,9的方差相等,则x=_______________.16.若8m a =,2n a =,,m n 为正整数,则2m n a +=___________.17.如图,在平行四边形ABCD 中,10,8,AB m AD m AC BC ==⊥,则平行四边形ABCD 的面积为____________.18.已知2,3m n a a ==,则3m n a +=____.三、解答题(共78分)19.(8分)2019年11月30日上午符离大道正式开通,同时宿州至徐州的K902路城际公交开通试运营,小明先乘K902路城际公交车到五柳站下车,再步行到五柳景区游玩,从出发地到五柳景区全程31千米,共用了1个小时,已知步行的速度每小时4千米,K902路城际公交的速度是步行速度的10倍,求小明乘公交车所行驶的路程和步行的路程.20.(8分)(1)已知a +b =7,ab =10,求a 2+b 2,(a -b)2的值;(2)先化简(22221a a a +--2221a a a a --+)÷1a a +,并回答:原代数式的值可以等于-1吗?为什么?21.(8分)数轴上点表示,点关于原点的对称点为,设点所表示的数为,(1)求的值;(2)求的值.22.(10分)如图,已知四边形ABCD ,AB =DC ,AC 、BD 交于点O ,要使AOB DOC △≌△,还需添加一个条件.请从条件:(1)OB =OC ;(2)AC =DB 中选择一个合适的条件,并证明你的结论.解:我选择添加的条件是____,证明如下:23.(10分)利用乘法公式计算16827816878⨯-⨯:24.(10分)先化简,再求值:2112111x x x x +⎛⎫-÷⎪-+-⎝⎭,其中x 满足240x -=.25.(12分)分解因式:16n 4﹣126.计算:[(x 2+y 2)﹣(x ﹣y )2+2y (x ﹣y )]÷4y .参考答案一、选择题(每题4分,共48分)1、B【解析】解:c ac b++≠ab,1判断正确;227是有理数,2判断正确;﹣0.6,3判断错误;∵2<3,∴1﹣1<2,4判断正确;数轴上有无理数,5判断正确;张晓亮的答卷,判断正确的有4个,得80分.故选B.【点睛】本题主要考查了实数的大小比较,实数的分类等知识点,属于基础知识,同学们要熟练掌握.2、A【解析】根据等量关系:甲种防腐药水+乙种防腐药水=18千克,甲种防腐药+乙种防腐药=18×50%千克,可得出关于x、y的二元一次方程组,解之即可得出结论.【详解】由题意得:18 30%75%1850%x yx y+=⎧⎨+=⨯⎩.故选A.【点睛】本题考查由实际问题抽象出二元一次方程组,根据数量关系找出关于x、y的二元一次方程是解题关键.3、D【分析】根据三角形任意两边之和都大于第三边逐个判断即可.【详解】解:A、3+4<8,不符合三角形三边关系定理,故本选项错误;B、2+3<6,不符合三角形三边关系定理,故本选项错误;C、5+6=11,不符合三角形三边关系定理,故本选项错误;D、5+6>10,6+10>5,5+10>6,符合三角形三边关系定理,故本选项正确;故选D.【点睛】本题考查了三角形的三边关系定理的应用,主要考查学生对三角形的三边关系定理的理解能力,注意:三角形的两边之和大于第三边,三角形的两边之差小于第三边.4、B【分析】根据三角形三条边的关系计算即可,三角形任意两边之和大于第三边,任意两边之差小于第三边.【详解】A.∵2+2<6,∴2,2,6不能组成三角形;B.∵3+4>5,∴3,4,5能组成三角形;C.∵3+5<9,∴3,5,9不能组成三角形;D.∵5+8=13,∴5,8,13不能组成三角形;故选B.【点睛】本题考查了三角形三条边的关系,熟练掌握三角形三条边的关系是解答本题的关键.5、B【分析】先根据等腰三角形的性质可求出B Ð的度数,再根据三角形的外角性质即可得.【详解】90,C AC BC∠=︒=1(18090)452B BAC ∠=∠=︒-︒=∴︒60,ADC ADC B BAD∠=︒∠=∠+∠604515BAD ADC B ∴∠=∠-∠=︒-︒=︒故选:B .【点睛】本题考查了等腰三角形的性质、三角形的外角性质,熟记各性质是解题关键.6、B【分析】根据线段垂直平分线的性质和等腰三角形的性质即可得到结论.【详解】∵AB 的垂直平分线交AB 于点D ,∴AE=BE ,∵△ACE 的周长=AC+AE+CE=AC+BC=13,△ABC 的周长=AC+BC+AB=19,∴AB=△ABC 的周长-△ACE 的周长=19-13=6,故答案为:B .【点睛】本题考查了线段垂直平分线的性质:垂直平分线垂直且平分其所在线段;垂直平分线上任意一点,到线段两端点的距离相等.7、A【解析】试题分析:解方程组23{3a b a b -=+=得:21a b =⎧⎨=⎩所以,等腰三角形的两边长为2,1.若腰长为1,底边长为2,由1+1=2知,这样的三角形不存在.若腰长为2,底边长为1,则三角形的周长为2.所以这个等腰三角形的周长为2.故选A.考点:1.等腰三角形的性质;2.解二元一次方程组.8、C【分析】根据近似数的精确度对各项进行判断选择即可.【详解】A.0.350是精确到0.001的近似数,正确;B.3.80万是精确到百位的近似数,正确;C.近似数26.9精确到十分位,26.90精确到百分位,表示的意义不相同,所以错误;D.近似数2.20是由数a 四会五入得到的,那么数a 的取值范围是2.195 2.205a < ,正确;综上,选C.【点睛】本题考查了近似数,精确到第几位是精确度常用的表示形式,熟知此知识点是解题的关键.9、C【分析】根据平方根的性质,正数有两个平方根且互为相反数,开方求解即可.【详解】∵一个正数有两个平方根且互为相反数∴4的平方根是2±故选:C .【点睛】本题主要考查平方根的性质,熟知一个正数有两个平方根并互为相反数是解题的关键,区分平方根与算术平方根是易错点.10、B【分析】绝对值小于1的正数也可以利用科学记数法表示,一般形式为10n a -⨯,与较大数的科学记数法不同的是其所使用的是负指数幂,指数由原数左边起第一个不为零的数字前面的0的个数所决定.【详解】0.0000000052=95.210-⨯.故选:B .【点睛】本题考查用科学记数法表示较小的数,一般形式为10n a -⨯,其中110a ≤<,n 为由原数左边起第一个不为零的数字前面的0的个数所决定.11、C【分析】根据命题的定义对四个选项进行逐一分析即可.【详解】解:A 、正确,对直角的性质作出了判断,故不符合题意;B 、正确,两点之间,线段最短,作出了判断,故不符合题意;C 、错误,是叙述一件事,没作出任何判断,故符合题意;D 、正确,对a 2和b 2的关系作了判断,故不符合题意;故选C.【点睛】本题考查的是命题的定义,即判断一件事情的语句叫命题.12、D【解析】根据全等三角形的判定方法分别进行分析即可.【详解】AD=CF ,可用SAS 证明△ABC ≌△DEF ,故A 选项不符合题意,∠BCA=∠F ,可用AAS 证明△ABC ≌△DEF ,故B 选项不符合题意,∠B=∠E ,可用ASA 证明△ABC ≌△DEF ,故C 选项不符合题意,BC=EF ,不能证明△ABC ≌△DEF ,故D 选项符合题意,故选D.【点睛】本题考查三角形全等的判定方法,判定两个三角形全等的一般方法有:SSS 、SAS 、ASA 、AAS 、HL .但是AAA 、SSA 不能判定两个三角形全等,判定两个三角形全等时,必须有边的参与,若有两边一角对应相等时,角必须是两边的夹角.二、填空题(每题4分,共24分)13、16π+【分析】连接AB ,判断出ABC 是等边三角形,然后根据扇形及三角形的面积公式,即可求得阴影部分的面积为:()S S S S ABCBCD ABC =--阴影扇形扇形.【详解】解:连接AB ,∵BC AC AB 8===,∴ABC 是等边三角形,∴SABC 182=⨯⨯=,ABC 60∠=,∴()ABCBCD ABC S S S S=--阴影扇形扇形22150π860π8360360⎛⨯⨯=-- ⎝16π=+故答案为:16π+.【点睛】本题考察扇形中不规则图形面积的求解,掌握扇形的面积公式是解题的关键.14、21a -<<.【解析】试题分析:点P (1,2)a a -+关于y 轴的对称点在第二象限,在P 在第一象限,则10{,20a a ->+>2 1.a ∴-<<考点:关于x 轴、y 轴对称的点的坐标.15、1或1【解析】∵一组数据2,3,4,5,x 的方差与另一组数据5,1,7,8,9的方差相等,∴这组数据可能是2,3,4,5,1或1,2,3,4,5,∴x=1或1,故答案是:1或1.16、1【分析】根据同底数幂的乘法及幂的乘方的逆运算即可解答.【详解】解:222()m n m n m n a a a a a +=⋅=⋅∵8m a =,2n a =∴22()8232m n a a ⋅=⨯=,故答案为:1.【点睛】本题考查了同底数幂的乘法及幂的乘方的逆运算,解题的关键是熟练掌握同底数幂的乘法及幂的乘方的逆运算.17、48m 1【分析】由平行四边形的性质可得BC=AD=8m ,然后利用勾股定理求出AC ,根据底乘高即可得出面积.【详解】∵四边形ABCD 为平行四边形∴BC=AD=8m ∵AC ⊥BC∴△ABC 为直角三角形∴平行四边形ABCD 的面积=BC AC=86=48⋅⨯m 1故答案为:48m 1.【点睛】本题考查了平行四边形的性质与勾股定理,题目较简单,根据平行四边形的性质找到直角三角形的边长是解题的关键.18、1【分析】根据幂的乘方以及同底数幂乘法的逆用进行计算即可.【详解】解:∵2,3m n a a ==,∴()33332354m n m n m n a a a a a +=⋅=⋅=⨯=,故答案为:1.【点睛】本题主要考查了幂的乘方以及同底数幂的乘法,熟练掌握幂的运算性质是解答本题的关键.三、解答题(共78分)19、30千米;1千米【分析】设小明行驶的路程为x 千米,步行的路程y 千米,根据题意可得等量关系:①步行的路程+行驶的路程=31千米;②公交车行驶x 千米时间+步行y 千米的时间=1小时,根据题意列出方程组即可.【详解】解:设小明乘车路程为x 千米,步行的路程y 千米,∵公交的速度是步行速度的10倍,步行的速度每小时4千米,∴公交的速度是每小时40千米,由题意得:311404x y x y +=⎧⎪⎨+=⎪⎩,解得:301x y =⎧⎨=⎩,∴小明乘公交车所行驶的路程为30千米,步行的路程为1千米.【点睛】本题主要考查了二元一次方程组的应用,关键是正确理解题意,找出题目中的等量关系.20、(1)a 2+b 2=29,(a -b)2=9;(2)原代数式的值不能等于-1,理由见解析.【解析】试题分析:(1)根据完全平方公式,即可解答;(2)原式括号中两项约分后,利用乘法分配律化简,约分后利用同分母分式的减法法则计算得到最简结果,令原式的值为-1,求出x 的值,代入原式检验即可得到结果.试题解析:(1)a 2+b 2=(a +b)2-2ab =72-2×10=49-20=29,(a -b)2=(a +b)2-4ab =72-4×10=49-40=9.(2)原式=()()()()()22111[]•111a a a a a a a a a +-+-+--=()21111a a a a ++---=11a a +-,原式的值为-1,即11a a +-=-1,去分母得:a+1=-a+1,解得:a=0,代入原式检验,分母为0,不合题意,则原式的值不可能为-1.21、(1);(2)1.【解析】由对称性求出点B 表示的数,即为x 的值将x 的值代入原式计算即可得到结果.【详解】解:(1)∵数轴上点A 表示,点A 关于原点的对称点为B ,∴数轴上表示点B 表示-,即x=-(2)由(1)得,x=-将x=-代入原式,则=(-2)2+=8-2=1.【点睛】此题考查了实数与数轴,熟练掌握运算法则是解本题的关键.22、条件是(2)AC =DB ,证明见解析【分析】根据三角形全等的条件进行选择判断,先证明ABC DCB ∆≅∆,可以得到=BAC CDB ∠∠,从而可以证明出AOB DOC ∆≅∆.【详解】解:选择的条件是(2)AC DB =,证明如下:在ABC DCB ∆∆和中,∵AB DC AC DB BC CB =⎧⎪=⎨⎪=⎩,∴ABC DCB∆≅∆∴=BAC CDB∠∠在AOB DOB ∆∆和中,∵AOB DOC BAC CDB AB DC ∠=∠⎧⎪∠=∠⎨⎪=⎩,∴AOB DOC∆≅∆【点睛】本题考查了全等三角形的判定,在全等三角形的5种判定方法中,选用合适的方法进行判定是解题的关键.23、33600【分析】根据乘法分配律的逆运算进行计算,即可得到答案.【详解】解:16827816878⨯-⨯=)(16827878⨯-=168200⨯=33600;【点睛】本题考查了有理数的混合运算,解题的关键是熟练掌握运算法则进行解题.24、22x +,12.【分析】根据分式混合运算的法则把原式进行化简,再求出x 的值代入进行计算即可.【详解】原式11(1)(1)()112x x x x x +-=-⨯-++1122x x x x +-=-++22x =+因为:240x -=2x =当2x =时,原式12=.【点睛】本题考查分式的化简求值,熟练掌握计算法则是解题关键.25、(4n 2+1)(2n +1)(2n -1)【分析】根据公式法,利用平方差公式,即可分解因式.【详解】解:原式=(4n 2+1)(4n 2-1)=(4n 2+1)(2n+1)(2n-1).【点睛】本题考查分解因式,较容易,熟练掌握公式法分解因式,即可顺利解题.26、x ﹣12y 【分析】首先利用完全平方公式计算小括号,然后再去括号,合并同类项,最后再计算除法即可.【详解】解:原式=(x 2+y 2﹣x 2+2xy ﹣y 2+2xy ﹣2y 2)÷4y ,=(4xy ﹣2y 2)÷4y ,=x ﹣12y .【点睛】此题主要考查了整式的混合运算,关键是掌握计算顺序:有乘方、乘除的混合运算中,要按照先乘方后乘除的顺序运算,其运算顺序和有理数的混合运算顺序相似.。
数学实验_重庆大学中国大学mooc课后章节答案期末考试题库2023年
数学实验_重庆大学中国大学mooc课后章节答案期末考试题库2023年1.无向图中边的端点地位是平等的、边是无序点对。
而有向图中边的端点的地位不平等,边是有序点对,不可以交换。
参考答案:正确2.人口数量与下列因素都有关,人口基数、出生率、死亡率、年龄结构、性别比例、医疗水平、工农业生产水平、环境、生育政策等等。
参考答案:正确3.一元5次代数方程在复数范围内有多少个根?参考答案:54.任何贪心算法都能求出最优解。
参考答案:错误5.二维插值函数z=interp2(x0,y0,z0,x,y,’method’)中,method的缺省值是()参考答案:linear6.在当前文件夹和搜索路径中都有文件ex1.m,在命令行窗口输入ex1时,则执行的文件是当前文件夹中的ex1.m参考答案:正确7.下列关于Dijkstra算法的哪些说法正确参考答案:Dijkstra算法是求加权图G中从某固定起点到其余各点最短路径的有效算法;_Dijkstra算法的时间复杂度为O(n2),其中n为顶点数;_Dijkstra算法可用于求解无向图、有向图和混合图的最短路径问题;8.如果x=1: 2 : 10,则x(1)和x(5)分别是( )参考答案:1,99.人口是按指数规律无限增长的。
参考答案:错误10.在包汤圆问题的整个建模过程,包括了如下几个步骤(1)找出问题涉及的主要因素(变量),重新梳理问题使之更明确(2)作出简化、合理的假设(3)用数学的语言来描述问题(4)用几何的知识解决问题(5)模型应用参考答案:正确11.下面程序所解的微分方程组,对应的方程和初始条件为:(1)函数M文件weif.m:function xdot=weif(t, x)xdot=[3*x(1)+x(3);2*x(1)+6;-3*x(2)^2+2*x(3)];(2)脚本M文件main.m:x0=[1,2,3] ;[t,x]=ode23(‘weif’,[0,1],x0),plot(t,x’),figure(2),plot3(x( :,1),x( :,2),x( :,3)参考答案:___12.某公司投资2000万元建成一条生产线。
工程数学实验
工程数学实验报告成绩:2015—2016—2学期学部:班级:姓名:学号:电话:Ⅰ展示图形之美篇要求:涉及到的文字用中文宋体五号字,Mathematica程序中的字体用Times New Roamn。
【数学实验一】题目:利用Mathematica制作如下图形(1),,其中k的取值为自己学号的后三位。
(2),其中k的取值为自己学号的后三位。
Mathematica程序:(1)ParametricPlot[{423Sin[t],423Sin[2t]},{t,0,2Pi}](2)x=Sin[u]Cos[423v]y=Sin[u]Cos[v]z=Cos[u]ParametricPlot3D[{Sin[u]Cos[423v],Sin[u]Cos[v],Cos[u]},{u,0,Pi},{v,0,2Pi}]运行结果:【数学实验二】题目:请用Mathematica制作五个形态各异三维立体图形,图形函数自选,也可以由几个函数构成更美观、更复杂的图形;并用简短的语言说明选择该图形的理由和意义.Mathematica程序:x[u_,v_]:=Sin[u]Cos[v];y[u_,v_]:=Sin[u]Sin[v];z[u_,v_]:=Cos[u];ParametricPlot3D[{x[u,v],y[u,v],z[u,v]},{u,—Pi/12,Pi/12Pi},{v,0,4Pi},Boxed—〉False,BoxRatios{1,1,1}]运行结果:图片像一个窝窝头,粮食是人类的生存之本Mathematica程序:ParametricPlot3D[{r,Exp[—r^2Cos[4r]^2]*Cos[t],Exp[—r^2Cos[4r]^2]Sin[t]},{r,-1.2,1.2},{t,0,2Pi}]运行结果:图片像一块奶糖Mathematica程序:ContourPlot3D[(2x^2+y^2+z^2-1)^3-(x^2z^3)/10—y^2*z^30,{x,—1。
数学实验报告-非线性规划与多目标规划实验
1)建立函数M文件:
functionf=fun(x)
f=-20*exp(-0.2*(0.5*(x(1)^2+x(2)^2)^0.5))-exp(0.5*(cos(2*pi*x(1))+cos(2*pi*x(2))))+22.713;
2)
x0=[0,0];
options=optimset('display', 'iter', 'tolfun',1e-10);
1982
0.117
0.465
0.215
0.187
0.213
0.311
-0.019
0.084
1983
0.092
-0.015
0.224
0.235
0.217
0.08
0.237
-0.128
1984
0.103
0.159
0.061
0.03
-0.097
0.15
0.074
-0.175
1985
0.08
0.366
0.316
[3]熟悉MATLAB软件求解非线性规划模型的基本命令;
[4]通过范例学习,了解建立非线性规划模型的全过程,与线性规划比较其难点何在。
本实验包括基础实验、应用实验和创新实验,基础实验和应用实验要求独立完成,创新实验要求合作完成。通过该实验的学习,使学生掌握最优化技术,认识面对什么样的实际问题,提出假设和建立优化模型,并且使学生学会使用MATLAB软件和Lingo软件求解非线性规划模型,注意初始解的选择不同会导致软件求出的解的变化(是局部最优解还是整体最优解)。解决现实生活中的最优化问题是本科生学习阶段中一门重要的课程,因此,本实验对学生的学习尤为重要。
最新数学实验报告
重庆大学学生实验报告实验课程名称数学实验开课实验室DS1407学院自动化年级2013 专业班自动化02班学生姓名侯刚学号20134615开课时间2014 至2015 学年第二学期数学与统计学院制开课学院、实验室:数统学院DS1407实验时间:2014年4月3日课程名称数学实验实验项目名称种群数量的状态转移——微分方程实验项目类型验证演示综合设计其他指导教师龚劬成绩√实验目的[1] 归纳和学习求解常微分方程(组)的基本原理和方法;[2] 掌握解析、数值解法,并学会用图形观察解的形态和进行解的定性分析;[3] 熟悉MATLAB软件关于微分方程求解的各种命令;[4] 通过范例学习建立微分方程方面的数学模型以及求解全过程;基础实验一、实验内容1.微分方程及方程组的解析求解法;2.微分方程及方程组的数值求解法——欧拉、欧拉改进算法;3.直接使用MATLAB命令对微分方程(组)进行求解(包括解析解、数值解);4.利用图形对解的特征作定性分析;5.建立微分方程方面的数学模型,并了解建立数学模型的全过程。
二、实验过程1.求微分方程的解析解, 并画出它们的图形,y’= y + 2x, y(0) = 1, 0<x<1;(1)求解:输入:dsolve('Dy=y+2*x','y(0)=1','x')输出:ans=-2*x-2+3*exp(x)(3)作图:输入:>> x=0:0.1:1;>> y2=-2*x-2+3*exp(x);>> plot(x,y2)输出:图表 1 方程特解图形分析:注意dsolve的用法。
2.用向前欧拉公式和改进的欧拉公式求方程y’= y - 2x/y, y(0) = 1 (0≤x≤1,h = 0.1) 的数值解,要求编写程序,并比较两种方法的计算结果,说明了什么问题?(1)求解析解输入: dsolve('Dy=y-2*x/y','y(0)=1','x')输出: ans =(2*x+1)^(1/2)(2)用向前欧拉公式和改进的欧拉公式求方程的数值解并与解析解作图比较程序:x1(1)=0;y1(1)=1;y2(1)=1;h=0.1;for k=1:10x1(k+1)=x1(k)+h;y1(k+1)=y1(k)+h*(y1(k)-2*x1(k)/y1(k));k1=y2(k)-2*x1(k)/y2(k);k2=y2(k)+h*k1-2*x1(k+1)/(y2(k)+h*k1);y2(k+1)=y2(k)+h*(k1+k2)/2;endx1,y1,y2x=0:0.1:1;y=(2*x+1).^(1/2);plot(x,y,x,y1,'o',x,y2,'+')结果:x1 =0 0.1000 0.2000 0.3000 0.4000 0.5000 0.6000 0.7000 0.8000 0.9000 1.0000y1 =1.0000 1.1000 1.1918 1.2774 1.3582 1.4351 1.5090 1.58031.6498 1.7178 1.7848y2 =1.0000 1.0959 1.1841 1.2662 1.3434 1.4164 1.4860 1.55251.6165 1.6782 1.7379图表 2 向前欧拉公式和改进的欧拉公式所求方程数值解与解析解的比较由图可得,改进后的欧拉公式求得的数值解更贴合解析解。
数学实验综合实验报告
数学实验综合实验报告《数学实验综合实验报告》摘要:本实验旨在通过数学实验的方式,探索和验证数学理论,并通过实验数据的分析和处理,得出结论和结论。
本实验涉及到数学的多个领域,包括代数、几何、概率统计等。
通过实验,我们得出了一些有趣的结论和发现,验证了数学理论的正确性,并对数学知识有了更深入的理解。
一、实验目的1. 验证代数公式的正确性2. 探索几何图形的性质3. 分析概率统计的实验数据4. 探讨数学理论的应用二、实验方法1. 代数公式验证实验:通过代数运算和数值计算,验证代数公式的正确性。
2. 几何图形性质探索实验:通过几何构造和图形分析,探索几何图形的性质。
3. 概率统计数据分析实验:通过实验数据的收集和处理,分析概率统计的规律和特性。
4. 数学理论应用实验:通过实际问题的分析和解决,探讨数学理论在实际中的应用。
三、实验结果与分析1. 代数公式验证实验结果表明,代数公式在特定条件下成立,验证了代数理论的正确性。
2. 几何图形性质探索实验发现,某些几何图形具有特定的性质和规律,进一步加深了对几何学的理解。
3. 概率统计数据分析实验得出了一些概率统计的规律和结论,对概率统计理论有了更深入的认识。
4. 数学理论应用实验通过具体问题的分析和解决,验证了数学理论在实际中的应用性。
四、结论通过本次数学实验,我们验证了代数、几何、概率统计等数学理论的正确性,得出了一些有意义的结论和发现。
实验结果进一步加深了对数学知识的理解和应用,对数学理论的研究和发展具有一定的参考价值。
五、展望本次实验虽然取得了一些有意义的结果,但也存在一些不足之处,如实验方法的局限性、实验数据的局限性等。
未来可以进一步完善实验设计和方法,开展更深入的数学实验研究,为数学理论的发展和应用提供更多的支持和帮助。
重庆大学数学实验一 matlab的基本应用 参考答案
《数学实验》第一次上机实验1. 设有分块矩阵⎥⎦⎤⎢⎣⎡=⨯⨯⨯⨯22322333S O R E A ,其中E,R,O,S 分别为单位阵、随机阵、零阵和对角阵,试通过数值计算验证⎥⎦⎤⎢⎣⎡+=22S 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.632.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.632.740 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.002.某零售店有9种商品的单件进价(元)、售价(元)及一周的销量如表1.1,问哪种商品的利润最大,哪种商品的利润最小;按收入由小到大,列出所有商品及其收入;求这一周该10种商品的总收入和总利润。
表1.11)程序: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];s=sum((b-a).*c)i=b.*cmax((b-a).*c)min((b-a).*c)[m,n]=sort(b.*c)2)运行结果:s =4.6052e+004i =1.0e+004 *0.6305 1.8075 0.4518 0.9425 0.3911 3.8398 3.1990 1.95621.0757ans =1.3087e+004ans =1.2719e+003m =1.0e+004 *0.3911 0.4518 0.6305 0.9425 1.0757 1.8075 1.9562 3.1990 3.8398n =5 3 1 4 9 2 8 7 63. 近景图将x的取值范围局限于较小的区间内可以画出函数的近景图,用于显示函数的局部特性。
重庆大学--数学实验-温罗生
0.6959 0.6999 0.8200 0.9686 0.7788 0.8754 0.9577 0.6951 0.8444 0.7805
ans =
0.9686
0.9436
0.6718
0.9577
0.8754
ans =
0.0942 0.0336 0.4076 0.1056 0.0908 0.2810 0.2407 0.0680 0.2240 0.0067
ans =
0.9686
0.9436
0.6718
0.9577
0.8754
ans =
2.2825 1.7604 3.1313 2.5417 1.7132 2.5810 3.2978 2.4001 2.3355 2.4515
ans =
0.0942 0.0336 0.4076 0.1056 0.0908 0.2810 0.2407 0.0680 0.2240 0.00670.4229 0.0688 0.5309 0.3251 0.1537 0.4401 0.5181 0.2891 0.2548 0.38680.4709 0.3196 0.6544 0.5313 0.2665 0.4574 0.6377 0.6718 0.3445 0.60220.5985 0.6385 0.7184 0.6110 0.4235 0.5271 0.9436 0.6761 0.6678 0.67530.6959 0.6999 0.8200 0.9686 0.7788 0.8754 0.9577 0.6951 0.8444 0.7805
图形
5.代码
[x,y]=meshgrid(-5:0.1:5);
z=sin(pi*sqrt(x.^2+4*y.^2));
重庆大学--数学模型--数学实验作业二
显然X^2=3,是x^2-3=0的一个零点。
5.编辑函数-M文件h.m(见附件),在MATLAB命令行窗口键入Untitled5.m(见附件),得到解。由于方程系数随机产生,故每次运行结果都不同。
6.编辑函数-M文件liu.m(见附件),在MATLAB命令行窗口键入Untitled6.m(见附件),得到解如下(仅列出部分):
a(1)*x(1)^2+a(2)*y(1)^2+a(3)*x(1)*y(1)+a(4)*x(1)+a(5)*y(1)=1
a(1)*x(2)^2+a(2)*y(2)^2+a(3)*x(2)*y(2)+a(4)*x(2)+a(5)*y(2)=1
a(1)*x(3)^2+a(2)*y(3)^2+a(3)*x(3)*y(3)+a(4)*x(3)+a(5)*y(3)=1
x=a;
17、近年来,我国积极推广“无车日”活动,以节约能源和保护环境。科学家也正在研制太阳能汽车和燃料电池汽车,减少对空气的污染。end
7、食盐、白糖、碱面、味精的颗粒都是有规则几何外形的固体,人们把这样的固体物质叫做晶体。自然界中的大部分固体物质都是晶体或由晶体组成。
5.h.m:
3、除了我们日常生活产生的家庭垃圾外,工厂、学校、医院、建筑工地等每天也在产生大量的垃圾。functionf=h(x)
应用实验或综合实验一问题重述小行星的运动轨道问题一天文学家要确定一颗小行星绕太阳运行的轨道他在轨道平面内建立以太阳为原点的直角坐标系其单位为天文测量单位
开课学院、实验室:数统学院实验时间:
课程
名称
数学实验
实验项目
名称
重庆大学数学实验课程实验1下
分析:
先用x=linspace(x1,x2)函数确定图像所画的范围,本次实验范围采用[-5,5]即x=linspace(-5,5);再分别用四个函数y1,y2,y3,y4分别表示所需要做图的函数,用plot(x,y,s)函数将四个函数表示在一起,为更好区分不同函数的图像,s分别区实线(‘-’)、点线(‘.’)、虚线(‘:’)、o线(‘o’)做图,最后用legend()函数将全部函数进行标注、区分,做出的图像如图二
由图像可知y1=e^x的图像在x取值越大的时候变化越快;y3=1+x+(1/2)x^2恒大于0;y2=1+x,'y4=1+x+(1/2)x^2+(1/6)x^3都是有正有负
使用subplot()函数将屏幕分成四块,并在每一个块中分别使用plot()函数绘制函数图像,使用title()函数为每一个块赋一个标题即函数解析式,做出的图像如图三
6.作出下列曲面的3维图形,
7.作出函数y=x4-4x3+3x+5(x[0,6])的图形,找出其在[0,6]之间的最小值点。
二、实验过程(一般应包括实验原理或问题分析,变量说明、程序、调试情况记录、图表等,实验结果及分析)
4:
程序:
functionf = sq(x)
a=x;
x=1;
y=(x+a/x)/2;
重庆大学
学生实验报告
实验课程名称数学实验
开课实验室
学院计算机年级2015专业班计科6
学生姓名王博韬学号********
2 实验指导书(2010.4)
页眉内容《数学实验》实验指导书龚劬重庆大学数学实验教学示范中心目录预备实验——桥梁分析.............................................................. 错误!未定义书签。
实验1 MATLAB软件入门.......................................................... 错误!未定义书签。
实验2 方程模型及其求解算法............................................... 错误!未定义书签。
实验3 收敛与混沌——迭代................................................... 错误!未定义书签。
实验4 微分方程模型、求解及稳定性分析........................... 错误!未定义书签。
实验5 插值方法....................................................................... 错误!未定义书签。
实验6 数据拟合及参数辨识方法........................................... 错误!未定义书签。
实验7 回归分析模型、求解及检验....................................... 错误!未定义书签。
实验8 连续系统与离散系统的计算机模拟........................... 错误!未定义书签。
实验9 线性规划模型、求解及灵敏度分析........................... 错误!未定义书签。
实验10 非线性规划与多目标规划模型及其求解................. 错误!未定义书签。
实验11 如何表示二元关系—图的模型及矩阵表示............. 错误!未定义书签。
重庆大学通识选修课
课程类别序号课程代码中文名称学分开课单位教学周次176100120中国文化之论语2高研院1月9日276100220中国文化之《资治通鉴》导读2高研院1月9日376100320莎士比亚戏剧选读2外国语1月9日476100420音乐与心理调适2艺术1月9日576100520电影美学2电影1月9日676100620中国古代思想与文化2高研院1月9日776100720中国文化史专题2高研院1月9日876100920中国文化之《老子》读解2高研院10月17日976101120中国传统文化经典解读2公共管理1月9日1076101220中国文化之《诗经》导读2高研院10月17日1176101320西方哲学史专题2公共管理1月9日1276101420柏拉图《理想国》导读2高研院1月9日1376101520汉语历史与文化2高研院1月9日1476101620中国美术史及名作鉴赏2艺术1-9;10-171576101720外国美术史及名作鉴赏2艺术1月9日1676101820中国电影艺术史2电影1月9日1776101920媒介文化2新闻1月9日1876102020计算机辅助古典诗词写作2外国语1月9日1976102120全球通识与跨文化交流2外国语1月9日2076102220爱与正义:法律与文学导读2高研院1月9日2176102320柏拉图对话讲疏:美诺篇2高研院1月9日2276102620《世说新语》导读2高研院1月9日2376102720乌托邦名著选读2高研院1月9日2476102820荷马史诗2高研院1月9日2576102920阅读城市与建筑2建筑城规10月17日2676103020“小说中国”:文学与批评导读2高研院10月17日2776103120文学视野下的中国近代史2高研院1-9;10-172876103220中国文学经典导读2高研院10月17日2976103320明清经典小说鉴赏2高研院10月17日3074510110音乐鉴赏1艺术10月17日3174530110播音与主持艺术概论1电影10月17日3274540110多媒体作品创作与鉴赏 1通信1-9;10-173374540210现代音响技术与音乐鉴赏1材料10月17日3474550410中国民间美术鉴赏1艺术1月9日3574560110唐代小说鉴赏1高研院1月9日3674570110摄影与艺术鉴赏1材料1-9;10-173774610710中国文化概论1公共管理1-9;10-173874620110俄罗斯语言与文化1外国语1-9;10-173974620210法国语言与文化1外国语1-9;10-174074630110德国语言与文化1外国语1-9;10-174174630210日语语言与文化1外国语1月9日176200220逻辑学2公共管理1月9日276200320法律与社会2法学院1月9日376200420创业与合同2法学院10月17日476200520经济学导论2经管10月17日576200720现代质量管理2机械1月9日676201020公益慈善与社会发展2公共管理10月17日776201220经济学原理2公共管理1月9日876201320信息素养2经管10月17日人文学科人文学科社会科学。
重庆大学理论力学教案考点
重庆大学理论力学教案考点一、引言理论力学是力学的基础课程之一,是研究物体运动规律的数学模型和方法的学科。
本教案旨在总结重庆大学理论力学课程的教学要点和考点,帮助学生更好地掌握和应用理论力学的知识。
二、课程概述理论力学课程主要包括牛顿力学和拉格朗日力学两个部分。
牛顿力学主要研究质点和刚体的运动规律,而拉格朗日力学则是一种更为普遍的力学方法,可以描述任意系统的运动。
三、教学要点1. 牛顿力学部分1.1 牛顿第一定律:惯性系、非惯性系、惯性力的概念及其应用。
1.2 牛顿第二定律:质点的运动方程、质点系的运动方程、力的合成与分解。
1.3 牛顿第三定律:作用力与反作用力的性质和应用。
1.4 力的合成与分解:力的合成定理、力的分解定理及其应用。
2. 拉格朗日力学部分2.1 广义坐标与广义速度:自由度的概念、广义坐标与广义速度的定义及其关系。
2.2 拉格朗日方程:拉格朗日方程的推导、广义力与约束力的概念。
2.3 拉格朗日方程的应用:保守系统与非保守系统的判定、稳定平衡与平衡条件、运动积分与守恒定律。
四、考点总结1. 牛顿力学考点1.1 质点的运动方程:根据给定的力和初始条件,求解质点的运动方程。
1.2 力的合成与分解:根据给定的力和角度,求解合力和分解力的大小和方向。
1.3 非惯性系中的运动:根据给定的非惯性系情况,求解惯性力和运动方程。
1.4 质点系的运动方程:根据给定的质点系情况,求解质点系的运动方程。
2. 拉格朗日力学考点2.1 广义坐标与广义速度:根据给定的系统,确定广义坐标和广义速度。
2.2 拉格朗日方程的推导:根据给定的系统和约束条件,推导出拉格朗日方程。
2.3 广义力与约束力:根据给定的系统和约束条件,确定广义力和约束力。
2.4 运动积分与守恒定律:根据给定的系统和势能函数,判断是否存在运动积分和守恒定律。
五、教学方法1. 理论讲解:通过讲解理论知识,帮助学生理解基本概念和原理。
2. 示例分析:通过实际例子,帮助学生掌握解题方法和思路。
利用矩阵翻转法求最佳H圈
第24卷第1期后 勤 工 程 学 院 学 报VoL.24No .1 2008年1月JOURNAL OF LOGI STI C AL E NGI N EER I N G UN I V ERSI TY Jan .2008 文章编号:1672-7843(2008)01-0102-05利用矩阵翻转法求最佳H 圈杨秀文1,陈振杰2,李爱玲3,田艳芳1(1.后勤工程学院基础部,重庆400016;2.后勤工程学院后勤信息工程系,重庆400042;3.湖南省金融培训学校,长沙410200)摘 要 利用矩阵翻转实现二边逐次修正法求最佳哈密尔顿圈(H 圈)。
首先构造完备加权图,并用距离矩阵表示之,使所选初始圈的顶点为矩阵主对角线的上方元素对应的顶点;然后对距离矩阵加边框并进行若干次“翻转”,直到矩阵不满足二边逐次修正法的修正原则,最后得到的矩阵主对角线的上方元素确定了最佳H 圈的权重及路线。
关键词 二边逐次修正法;最佳哈密尔顿圈;距离矩阵;矩阵翻转中图分类号:O157.5文献标识码:ASeeking the Best Ha m ilt on Cycle Thr ough Matrix Turning Y ANG Xiuwen 1,CHEN Zhengjie 2,L IA iling 3,TI A N Yanfang1(1.Dep t .of Foundati on Studies,LE U ,Chongqing 400016,China;2.Dep t .of Logistical I nf or mati on Engineering,LE U ,Chongqing 400042,China;3.Hunan Financial Training Acade my,Changsha 410200,China )ABSTRACT The paper utilizes the matrix turning t o realize that using the p rinci p le of one by one revisi on of t w o sides t o obtain the best m ilt on Circle .W e construct comp lete weighted graph at first,and rep resent it by matrix,such that the points of the initial Circle are the points corres ponding t o the upper ele ments of the matrix main diagonal .Then carrying on "turning"on matrix several ti m es until the matrix does not satisfy with the p rinci p le of one by one revisi on of t w o sides .A t last,the weight and r oute of the best Ha m ilt on Circle are confir med by the upper ele ments of the final matrix main diagonal .Keywords the p rinci p le of one by one revisi on of t w o sides;the best Ha m ilt on Circle;distance matrix;matrix turning求最佳哈密尔顿圈是一个NP 问题。
重庆大学数学与统计学院
学建模 竞赛组委会。
5 【 O 名人名言 】 教学 家就像 情人 ,给一 个数学 家最小的原理 , 就会从 中引出你必须承认的结果 .并且从 这个又引出另外一个。 他
— —
弗坦 内里 ( 法国 )
专栏 I ln u s C m o
大学 科 普
重 庆大 学 数 学 与 统 计 学 皖 师姿 队 伍 一 赀
大学 科普
- ̄Io n ' Cu s t - l m
重 庆 大 学 数 学 与 统 计 学 院
重庆 大学数学 与统计学 院 的前 身是始建 于 1 2 9 9年 的 重 庆大学 理学 院和 13 建立 的重庆大 学商学 院 ,理学 9 7年 院是重庆 大学 最早设立 的 个学 院之 一 ,首任 院长为数学 家何鲁先 生 ,著名数学家柯 召院士 、胡坤陛 、段调元 、潘 璞 、周雪 欧等数学名家早年都 曾任教于此 。商学 院首任 院 长是 马寅初 先生 ,保险和统计是商学 院最 早设立 的 个 系
动站培养 了一批计算数学专业人 才 ,2 0 年 三校合并后 , 00
实现强强 结合 、发展尤 为迅 速 ,成长起一批博士生导 师和 知名学者 ,在最优化数值解 、数 值逼 近与图论 、偏微分方
究工作,国内著名复分析专家顾永兴教授是上一任所长、
现任所 长穆春来 教授 ,研 究人员有李 江涛 f 博士 副教授 、 复函) 、舒 永 录 ( 博士 副教授 、动力 系统 )、黄小军 ( 博 士副教授 、复分 析 、圆填充 、几何 函数论 )、朱 长荣 ( 博 士副教授 、动力系统 、生物数学 )、蒋 卫生 ( 博士副教授 、 控制论 )、张万雄 ( 博士副教授 、动力 系统 )、张 良才 ( 博
一
、
数学实验(第二版)课后习题答案
贵州师范学院2012级数本一班李刚数学实验课后练习答案习题2.11. syms x y;>> x=-5:0.01:5;>> y=x.^1/2;>> plot(x,y)2. f plot('exp(-x.^2)',[-5,5])3. ezplot('x.^3+y.^3-3*x*y',[-5,5])4 . ezplot('y.^2-x.^3/(1-x)',[-5,5])5.t=0:0.1:2*pi;x=t-sin(t);y=2*(1-cos(t));plot(x,y)6. t=0:0.1:2*pi; x=cos(t).^3; >> y=sin(t).^3;>> plot(t,y)>>7: t=0:0.1:2*pi; x=cos(t); y=2*sin(t); z=3*t; plot3(x,y,z)8: x =0:0.1:2*pi; r=x; polar(x,r)9: x =0:0.1:2*pi; r=exp(x); polar(x,r)10: x=0:0.1:2*pi; r=sqrt(cos(2*x)); polar(x,r)11: x=0:0.1:2*pi; r=sqrt(sin(2*x)); polar(x,r)12: x =0:0.1:2*pi; r=1+cos(x); polar(x,r)练习2.2 1:(1)(2):syms n; limit('sqrt(n+2)-2*(sqrt(n+1))+sqrt(n)',n,inf)Ans= 0 (3):: (4):(5):(6):2:3:fplot('x.^2*sin(x.^2-x-2)',[-2,2])练习2.3 1:(2):2:练习2.4 1:(1)(2):(3)(4):2:(1):syms x;int(x^(-x),x,0,1)ans =int(x^(-x),x = 0 .. 1)vpa(ans,10)ans =1.291285997(2):syms x;int(exp(2*x)*cos(x)^3,x,0,2*pi)ans =-22/65+22/65*exp(4*pi)(3):syms x; int(exp(x^2/2)/sqrt(2*pi),x,0,1)ans =-1125899906842624/5644425081792261*i*erf(1/2*i*2^(1/2))*pi^(1/2)*2^(1/2) >> vpa(ans,10)ans =.4767191345(4):syms x;int(x*log(x^4)*asin(1/x^2),x,1,3)ans =int(x*log(x^4)*asin(1/x^2),x = 1 .. 3)>> vpa(ans,10)ans =2.459772128(5):syms x ;int(exp(x^2/2)/sqrt(2*pi),x,-inf,inf)ans =Inf(6):syms x ;int(sin(x)/x,x,0,inf)ans =1/2*pi(7):syms x ;int(tan(x)/sqrt(x),x,0,1)Warning: Explicit integral could not be found. > In sym.int at 58ans =int(tan(x)/x^(1/2),x = 0 .. 1)>> vpa(ans,10)ans =.7968288892(8):syms x ;int(exp(-x^2/2)/(1+x^4),x,-inf,inf)ans =1/4*pi^(3/2)*2^(1/2)*(AngerJ(1/2,1/2)-2/pi^(1/2)*sin(1/2)+2/pi^(1/2)*cos(1/2)-WeberE(1/2,1/2 ))>> vpa(ans,10)ans =1.696392536(9):syms x ;int(sin(x)/sqrt(1-x^2),x,0,1)ans =1/2*pi*StruveH(0,1)>> vpa(ans,10)ans =.8932437410练习2.5(1):syms n;symsum(1/n^2^n,n,1,inf)ans =sum(1/((n^2)^n),n = 1 .. Inf)(2):s yms n ;symsum(sin(1/n),n,1,inf)ans =sum(sin(1/n),n = 1 .. Inf)(3):syms n ;symsum(log(n)/n^3,n,1,inf) ans =-zeta(1,3)(4):syms n ;symsum(1/(log(n))^n,n,3,inf) ans =sum(1/(log(n)^n),n = 3 .. Inf)(5):syms n;symsum(1/(n*log(n)),n,2,inf) ans =sum(1/n/log(n),n = 2 .. Inf)(6):yms n;symsum((-1)^n*n/(n^2+1),n,1,inf)ans =-1/4*Psi(1-1/2*i)+1/4*Psi(1/2-1/2*i)-1/4*Psi(1+1/2*i)+1/4*Psi(1/2+1/2*i)第三章练习3.11:(1):a=-30:1:30;b=-30:1:30;[x,y]=meshgrid(a,b);z=10*sin(sqrt(x.^2+y.^2))./(sqrt(1+x.^2+y.^2)); meshc(x,y,z)(2):a=-30:1:30;b=-30:1:30;[x,y]=meshgrid(a,b);z=4*x.^2/9+y.^2;meshc(x,y,z)(3):(4):a=-30:1:30;b=-30:1:30;[x,y]=meshgrid(a,b); z=x.^2/3-y.^2/3; meshc(x,y,z)(5):a=-30:1:30;>> b=-30:1:30;>> [x,y]=meshgrid(a,b); >> z=x*y;>> meshc(x,y,z)(6):(7):a=-30:1:30;>> b=-30:1:30;>> [x,y]=meshgrid(a,b); >> z=sqrt(x.^2+y.^2); >> meshc(x,y,z)(8):(9):a=-30:1:30;>> b=-30:1:30;>> [x,y]=meshgrid(a,b);>> z=atan(x./y);>> meshc(x,y,z)练习3.21;a=-1:0.1:1;>> b=0:0.1:2;>> [x,y]=meshgrid(a,b);>> z=x.*exp(-x.^2-y.^2);>> [px,py]=gradient(z,0.1,0.1);>> contour(a,b,z)>> hold on>> quiver(a,b,px,py)2:a=-2:0.1:1;>> b=-7:0.1:1;>> [x,y]=meshgrid(a,b);>> z=y.^3/9+3*x.^2.*y+9*x.^2+y.^2+x.*y+9; >> plot3(x,y,z)>> grid on3:[x,y]=meshgrid(-2*pi:0.2:2*pi); z=x.^2+2*y.^2;plot3(x,y,z)hold onezplot('x^2+y^2-1',[-2*pi,2*pi]) ; grid on4:t=0:0.03:2*pi;>> s=[0:0.03:2*pi]';>> x=(0*s+1)*cos(t);y=(0*s+1)*sin(t);z=s*(0*t+1); >> mesh(x,y,z)>> hold on>> [x,y]=meshgrid(-1:0.1:1);>> z=1-x+y;>> mesh(x,y,z)5:syms x y z dx dyz=75-x^2-y^2+x*y;zx=diff(z,x),zy=diff(z,y)zx =-2*x+yzy =-2*y+x练习3.31:ezplot('x^2+y^2-2*x',[-2,2]);>> grid onsyms x y ;s=int(int(x+y+1,y,-sqrt(1-(x-1)^2),sqrt(1-(x-1)^2)),x,0,2)s =2*pi2:syms r t ;>> s=int(int(sqrt(1+r^2*sin(t)),r,0,1),t,0,2*pi)s =int(1/2*((1+sin(t))^(1/2)*sin(t)^(1/2)+log(sin(t)^(1/2)+(1+sin(t))^(1/2)))/sin(t)^(1/2),t = 0 .. 2*pi) 3:syms x y z ;>> s=int(int(int(1/(1+x+y+z)^3,z,0,1-x-y),y,0,1-x),x,0,1)s =-5/16+1/2*log(2)4:s=vpa(int(int(x*exp(-x^2-y^2),y,0,2),x,-1,10))s =0.16224980455070416645061789474030练习3.41:(1):y=dsolve('Dy=x+y','y(0)=1','x')得:y =-1-x+2*exp(x)(2):y=dsolve('Dy=2*x+y^2','y(0)=0')y =tan(t*x^(1/2)*2^(1/2))*x^(1/2)*2^(1/2)练习4.11:(1):p=[5 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -6 8 0 0 0 -5 0 0]; >> x=roots(p)x =0.97680.9388 + 0.2682i0.9388 - 0.2682i0.8554 + 0.5363i0.8554 - 0.5363i0.6615 + 0.8064i0.6615 - 0.8064i0.3516 + 0.9878i0.3516 - 0.9878i-0.0345 + 1.0150i-0.0345 - 1.0150i-0.4609 + 0.9458i-0.4609 - 0.9458i-0.1150 + 0.8340i-0.1150 - 0.8340i-0.7821 + 0.7376i-0.7821 - 0.7376i-0.9859 + 0.4106i-0.9859 - 0.4106i-1.0416-0.7927(2): p=[8 36 54 23];x=roots(p)x =-1.8969 + 0.6874i-1.8969 - 0.6874i-0.70632:p1=[1 0 -3 -2 -1];p2=[1 -2 5];[q2,r2]=deconv(p1,p2)q2 =1 2 -4r2 =0 0 0 -20 19 3:syms x;f=x^4+3*x^3-x^2-4*x-3;g=3*x^3+10*x^2+2*x-3;p1=factor(f),p2=factor(g)p1 =(x+3)*(x^3-x-1)p2 =(x+3)*(3*x^2+x-1)4:syms x ;f=x^12-1;p=factor(f)p =(-1+x)*(1+x^2+x)*(1+x)*(1-x+x^2)*(1+x^2)*(x^4-x^2+1)5: (1):p=[1 0 1];q=[1 0 0 0 1];[a,b,r]=residue(p,q)a =-0.0000 - 0.3536i-0.0000 + 0.3536i0.0000 - 0.3536i0.0000 + 0.3536ib =0.7071 + 0.7071i0.7071 - 0.7071i-0.7071 + 0.7071i-0.7071 - 0.7071ir =[](2):p=[1];q=[1 0 0 0 1];[a,b,r]=residue(p,q)a =-0.1768 - 0.1768i -0.1768 + 0.1768i0.1768 - 0.1768i0.1768 + 0.1768ib =0.7071 + 0.7071i0.7071 - 0.7071i -0.7071 + 0.7071i -0.7071 - 0.7071ir =[](3):p=[1 0 1];q=[1 1 -1 -1];[a,b,r]=residue(p,q)a =0.5000-1.00000.5000b =-1.0000-1.00001.0000r =[] (4): p=[1 1 0 0 0 -8];[a,b,r]=residue(p,q)a =-4-38b =-11r =1 1 1练习 4.21:(1):D=[2 1 3 1;3 -1 2 1;1 2 3 2;5 0 6 2];det(D)ans =6(2):syms a b c dD=[a 1 0 0 ;-1 b 1 0;0 -1 c 1;0 0 -1 d];det(D)ans =a*b*c*d+a*b+a*d+c*d+12:(1):D=[1 1 1 1; a b c d;a^2 b^2 c^2 d^2;a^3 b^3 c^3 d^3];det(D)ans =b*c^2*d^3-b*d^2*c^3-b^2*c*d^3+b^2*d*c^3+b^3*c*d^2-b^3*d*c^2-a*c^2*d^3+a*d^2*c^3+a *b^2*d^3-a*b^2*c^3-a*b^3*d^2+a*b^3*c^2+a^2*c*d^3-a^2*d*c^3-a^2*b*d^3+a^2*b*c^3+a^ 2*b^3*d-a^2*b^3*c-a^3*c*d^2+a^3*d*c^2+a^3*b*d^2-a^3*b*c^2-a^3*b^2*d+a^3*b^2*c(2): s yms a b x y zD=[a*x+b*y a*y+b*z a*z+b*x; a*y+b*z a*z+b*x a*x+b*y;a*z+b*x a*x+b*y a*y+b*z];det(D)ans =3*a^3*x*z*y+3*b^3*y*x*z-a^3*x^3-a^3*y^3-b^3*z^3-a^3*z^3-b^3*x^3-b^3*y^33: (1): D=[1 1 1 1;1 2 -1 4;2 -3 -1 -5;3 1 2 11];D1=[5 1 1 1;-2 2 -1 4;-2 -3 -1 -5;0 1 2 11];D2=[1 5 1 1;1 -2 -1 4;2 -2 -1 -5;3 0 2 11];D3=[1 1 5 1;1 2 -2 4;2 -3 -2 -5;3 1 0 11];D4=[1 1 1 5;1 2 -1 -2;2 -3 -1 -2;3 1 2 0];x1=det(D1)/det(D);x2=det(D2)/det(D);x3=det(D3)/det(D);x4=det(D4)/det(D);x1,x2,x3,x4x1 =1x2 =2x3 =3x4 =-1(2):D=[5 6 0 0 0;1 5 6 0 0;0 1 5 6 0;0 0 1 5 6;0 0 0 1 5]; D1=[1 6 0 0 0;0 5 6 0 0;0 1 5 6 0;0 0 1 5 6;1 0 0 1 5]; D2=[5 1 0 0 0;1 0 6 0 0;0 0 5 6 0;0 0 1 5 6;0 1 0 1 5]; D3=[5 6 1 0 0;1 5 0 0 0;0 1 0 6 0;0 0 0 5 6;0 0 1 1 5]; D4=[5 6 0 1 0;1 5 6 0 0;0 1 5 0 0;0 0 1 0 6;0 0 0 1 5]; D5=[5 6 0 0 1;1 5 6 0 0;0 1 5 6 0;0 0 1 5 0;0 0 0 1 1]; x1=det(D1)/det(D);x2=det(D2)/det(D);x3=det(D3)/det(D);x4=det(D4)/det(D);x5=det(D5)/det(D);x1,x2,x3,x4,x5x1 =2.2662x2 =-1.7218x3 =1.0571x4 =-0.5940x5 =0.3188练习 4.3 1:A=[1 2 0;3 4 -1; 1 1 -1];B=[1 2 3;-1 0 1;-2 4 -3];A',2+A,2*A-B,A*B,A^2,A^(-1)ans =1 3 12 4 10 -1 -1ans =3 4 25 6 13 3 1ans =1 2 -37 8 -34 -2 1ans =-1 2 51 2 162 -2 7ans =7 10 -214 21 -33 5 0ans =-3.0000 2.0000 -2.00002.0000 -1.0000 1.0000-1.0000 1.0000 -2.0000 2:(1):B=[2 4 3];B'ans =243(2):A=[1 2 3];B=[2 4 3];A.*B,B.*Aans =2 8 9ans =2 8 93:(1):A=[0 1 0;1 0 0;0 0 1];B=[1 0 0;0 0 1;0 1 0];C=[1 -4 3;2 0 -1;1 -2 0];A^(-1),B^(-1),X=A^(-1)*C*B^(-1) ans =0 1 01 0 00 0 1ans =1 0 00 0 10 1 0X =2 -1 01 3 -41 0 -2(2):>> A=[1 2 3;2 2 3;3 5 1];B=[1 0 0;2 0 0;3 0 0];A^(-1),x=A^(-1)*Bans =-1.0000 1.0000 0.00000.5385 -0.6154 0.23080.3077 0.0769 -0.1538x =1 0 00 0 00 0 0练习 4.41:(1):A=[4 2 -1;3 -1 2;11 3 0];b=[2;10;8];B=[A,b];rank(A),rank(B)ans =2ans =3(2):A=[2 1 -1 1;3 -2 1 -3;1 4 -3 5];b=[1;4;-2];B=[A,b];rank(A),rank(B)ans =2ans =2(3):A=[ 1 1 1 1; 1 2 -1 4;2 -3 -1 -5;3 1 2 11];b=[5;-2;-2;0];B=[A,b];rank(A),rank(B)ans =4ans =4(4):A=[ 1 1 2 -1; 2 1 1 -1;2 2 1 2];b=[0;0;0];B=[A,b];rank(A),rank(B)ans =3ans =32:syms a;A=[-2 1 1;1 -2 1;1 1 -2];b=[-2;a;a^2];B=[A,b];rank(A),rank(B)ans =2ans =3练习4.51:(1):A=[0 1;-1 0];[a,b]=eig(A)a =0.7071 0.70710 + 0.7071i 0 - 0.7071ib =0 + 1.0000i 000 - 1.0000i(2):A=[0 0 1;0 1 0;1 0 0];[a,b]=eig(A)a =0.7071 0.7071 00 0 -1.0000-0.7071 0.7071 0b =-1 0 00 1 00 0 1(3):A=[4 1 -1;3 2 -6;1 -5 3];[a,b]=eig(A)a =0.0185 -0.9009 -0.3066-0.7693 -0.1240 -0.7248-0.6386 -0.4158 0.6170b =-3.0527 0 00 3.6760 00 0 8.3766(4):A=[1 1 1 1;1 1 -1 -1;1 -1 1 -1;1 1 -1 1];[a,b]=eig(A)a =0.5615 0.3366 0.2673 -0.7683-0.5615 -0.3366 0.0000 -0.0000-0.5615 -0.3366 -0.5345 -0.6236-0.2326 0.8125 0.8018 -0.1447b =-1.4142 0 0 00 1.4142 0 00 0 2.0000 00 0 0 2.0000(5):A=[5 7 6 5;7 10 8 7;6 8 10 9;5 7 9 10];[a,b]=eig(A)a =0.8304 0.0933 0.3963 0.3803-0.5016 -0.3017 0.6149 0.5286-0.2086 0.7603 -0.2716 0.55200.1237 -0.5676 -0.6254 0.5209b =0.0102 0 0 00 0.8431 0 00 0 3.8581 00 0 0 30.2887(6):A=[5 6 0 0 0;1 5 6 0 0 ;0 1 5 6 0 ;0 0 1 5 6; 0 0 0 1 5 ]; [a,b]=eig(A)a =0.7843 -0.7843 -0.9860 -0.9237 -0.92370.5546 0.5546 0.0000 0.3771 -0.37710.2614 -0.2614 0.1643 -0.0000 0.00000.0924 0.0924 0.0000 -0.0628 0.06280.0218 -0.0218 -0.0274 0.0257 0.02579.2426 0 0 0 00 0.7574 0 0 00 0 5.0000 0 00 0 0 2.5505 00 0 0 0 7.4495 2:(1):A=[0 1;-1 0];[a,b]=eig(A)a =0.7071 0.70710 + 0.7071i 0 - 0.7071ib =0 + 1.0000i 00 0 - 1.0000i>> P=orth(a),B=P'*A*P,P*P'P =-0.7071 -0.70710 - 0.7071i 0 + 0.7071iB =0 + 1.0000i 0 - 0.0000i0 - 0.0000i 0 - 1.0000ians =1.0000 0 + 0.0000i0 - 0.0000i 1.0000>> inv(a)*A*a0 + 1.0000i 000 - 1.0000i3:(1):A=[2 0 0;0 3 2;0 2 3]; [a,b]=eig(A)a =0 1.0000 0-0.7071 0 0.70710.7071 0 0.7071b =1.0000 0 00 2.0000 00 0 5.0000>> P=orth(a),B=P'*A*P,P*P'P =-1.0000 0 -0.00000.0000 0.7071 0.7071-0.0000 -0.7071 0.7071B =2.0000 0.0000 0.00000.0000 1.0000 00.0000 0 5.0000ans =1.0000 -0.0000 0.0000-0.0000 1.0000 -0.00000.0000 -0.0000 1.0000(2):A=[1 1 0 -1;1 1 -1 0;0 -1 1 1;-1 0 1 1];[a,b]=eig(A)a =-0.5000 0.7071 0.0000 0.50000.5000 -0.0000 0.7071 0.50000.5000 0.7071 0.0000 -0.5000-0.5000 0 0.7071 -0.5000 b =-1.0000 0 0 00 1.0000 0 00 0 1.0000 00 0 0 3.0000 >> P=orth(a),B=P'*A*P,P*P'P =-0.5000 -0.4998 -0.4783 -0.52100.5000 -0.4822 0.5212 -0.49580.5000 0.4998 -0.4964 -0.5037-0.5000 0.5175 0.5031 -0.4786 B =-1.0000 0.0000 0.0000 0.00000.0000 2.9988 -0.0362 0.03440.0000 -0.0362 1.0007 -0.00060.0000 0.0344 -0.0006 1.0006 ans =1.0000 0.0000 0.0000 -0.00000.0000 1.0000 -0.0000 00.0000 -0.0000 1.0000 0.0000-0.0000 0 0.0000 1.0000练习5.3 1: [m,v]=unifstat(1,11)m =6v =8.33332:[m,v]=normstat(0,16)m =v =256>> s=sqrt(v)s =163:x=randn(200,6);s=std(x)s =0.9094 0.9757 0.9702 0.9393 0.9272 1.09824: x=normrnd(0,16,300,1);hist(x,10)练习 5.61:x=[352 373 411 441 462 490 529 577 641 692 743];y=[166 153 177 201 216 208 227 238 268 268 274];plot(x,y,'*')4:(1):x=[10 10 10 15 15 15 20 20 20 25 25 25 30 30 30];y=[25.2 27.3 28.7 29.8 31.1 27.8 31.2 32.6 29.7 31.7 30.1 32.3 29.4 30.8 32.8]; plot(x,y,'*')。
重庆大学数学实验实验报告(参考)
( x) x ( x) 1 ( x)
2 x 5x
4 x 2 x 5 x 3 1 5 2 15 x 2 5 ,算法设计为: 4 1 3 2 1 2 x 5 x 1 5 2 15 x 5 3
1 5
1
x1=0; x2=((2*x1-5*x1^3-1)^(1/5)-1/5*(2*x1-5*x1^3-1)^(-4/5)*(2*x1-15*x1^3))/(1-1/5*(2*x1-5*x1^3-1)^(4/5)*(2-15*x1^2)); for k=1:100 x1=x2; x2=((2*x1-5*x1^3-1)^(1/5)-1/5*(2*x1-5*x1^3-1)^(-4/5)*(2*x1-15*x1^3))/(1-1/5*(2*x1-5*x1^31)^(-4/5)*(2-15*x1^2)); end x1 输出结果为: x1 = -0.1483 + 0.7585i
2 2 1 2 3 2 2 x x 2 3 3 x ( x) x ( x) x x x x ,算法设计为: x h( x ) 2 1 ( x) 2 5 2 3 3x 2 x x
1
2
x1=0; x2=((0.4*x1-0.2*x1^5-0.2)^(1/3)-1/15*(0.4*x1-0.2*x1^5-0.2)^(-2/3)*(2*x1-5*x1^5))/(1-(1/15*(0.4*x 1-0.2*x1^5-0.2)^(-2/3)*(2-5*x1^4))); while abs(x1-x2)>10^(-5) x1=x2; x2=((0.4*x1-0.2*x1^5-0.2)^(1/3)-1/15*(0.4*x1-0.2*x1^5-0.2)^(-2/3)*(2*x1-5*x1^5))/(1-(1/15*(0 .4*x1-0.2*x1^5-0.2)^(-2/3)*(2-5*x1^4))); end x1 输出结果为: x1 = 0.4004 + 0.2860i ③迭代函数为 ( x)
重庆大学数学与统计学院岗位聘任结果
重庆大学数学与统计学院岗位聘任结果(经院设岗聘任委员会11.9.29讨论通过)院长、副院长全部有高岗安排http://202.202.1.211/gangw.asp数学与统计学院第五轮岗位聘任设置关键岗位:教学科研岗位9级2个,8级8个,7级10个;教学岗位8级2个,7级7个。
所有待聘关键岗位欢迎海内外学者前来应聘(Email:yh@)由此进入数统学院岗位聘任在线申请、考核、公示专板,该板为认证版面,仅限本院教师使用,新注册的教师请向学院申请认证●第五轮岗位聘任关键岗位(21人)团队岗位:数学与统计创新团队兼职岗位:8人杨虎(院长)、何传江(副院长)、穆春来(副院长)、李声杰、黄小军、朱长荣、张志民、王浩(成教毕业的院办主任)(全院1个院长、两个副院长,实力怎么样,请到网上查)教学科研九级:2数学一级学科博士点负责人:李声杰统计学一级学科博士点负责人:杨虎(院长,这么大个远,到底有多少时间在搞科研,不过凭借院长资源笼络了一批优秀学生,其他人羡慕)教学科研八级:8基础数学博士点负责人:待聘计算数学博士点负责人:何传江(副院长)应用数学博士点负责人:穆春来(副院长)运筹学与控制论博士点负责人:待聘数理统计博士点建设岗:待聘生物卫生统计博士点建设岗:待聘金融统计与风险管理和精算博士点建设岗:待聘应用统计博士点建设岗:待聘教学科研七级:10基础数学硕士点负责人:黄小军(低聘六级)计算数学硕士点负责人:曾理应用数学硕士点负责人:舒永录(低聘六级)运筹学与控制论硕士点负责人:伍俊良数理统计硕士点建设岗:待聘生物卫生统计硕士点建设岗:待聘金融统计与风险管理和精算硕士点建设岗:张志民(低聘五级)(院长的学生,刚毕业的毛头小子)应用统计硕士点建设岗:黎雅莲(低聘五级)(院长的学生)应用统计专业硕士点负责人:易正俊(低聘六级)软件工程硕士点负责人:何光辉(低聘五级)教学八级:2数学实验国家级精品课程负责人:龚劬工程数学市级教学团队负责人(兼概率论与数理统计课程负责人):刘琼荪教学七级:7高等数学课程第一负责人:叶仲泉线性代数课程第一负责人:段正敏(软件工程硕士,居然在数学拿七级岗)复变函数第一负责人:李江涛研究生课程第一负责人:钟波研究生课程第二负责人:王开荣数学教学实验中心负责人:谢德政信息与计算科学市级优势与特色专业负责人:胡小兵●第五轮岗位聘任一般岗位(59人)教学六级:3高等数学课程第二负责人:王新质数学与应用数学专业负责人:张谋统计学专业负责人:阴文革(低聘五级)教学与专技5级:18刘思泽、田玉芳、蒋卫生、肖志祥、于光磊、张平、温罗生、潘明勇、黄薇、黄光辉朱长荣、蒲学科、魏曙光、张良才、王晓宏、杨木洪、方延洪、曹蓓(院长夫人、实验管理员,居然和副教授、博士同一级,很多博士副教授还在4级中)教学与专技4级:21张敏、王克金、王海鹰、张万雄、王汉明、谭宏、肖剑、颜军、刘霞、荣腾中蔡薇、徐建文、李寒宇、周云华、陈纯荣、刘德强、张应应、刘朝林、刘琼芳、李曼曼王坤教学与专技3级:7邓林、胥斌、罗广萍、董海云、李小娅、潘致锋、彭智军(都是老革命了,还3级)管理六级:1王浩管理五级:2胡小兰赖大伟管理四级:7李黎黎李慧许菲戴娟罗强段曦盛刘忠英。
重庆大学--数学模型--数学实验作业四
开课学院、实验室:数统学院实验时间:2015年10月28日当固定参数b=2, c=4时,试讨论随参数a 由小到大变化(如a ∈(0,0.65))而方程解的变化情况,并且画出空间曲线图形,观察空间曲线是否形成混沌状?4.Apollo 卫星的运动轨迹的绘制二、实验过程1.编辑程序代码Untitle1:s=dsolve('Dy=y+2*x','y(0)=1','x')ezplot(s,[0,1])运行结果如下:s =3*exp(x) - 2*x – 2图形为:1133121331212222121()()2,2,1/82.45,1,(),()(0) 1.2,(0)0,(0)0,(0) 1.04935751x x x y x r r y y y x y r r r x y r x y x x y y μμμμμμμμμμμ+-=+--=-+--==-=++=-+====-⎪⎩⎪⎨⎧-+=+=--=)('''c x z b z ayx y z y x2.编写程序代码Untitle2:clcy=dsolve('Dy=y-2*x/y','y(0)=1','x')ezplot(y,[0,1])hold onx=[];x(1)=0;y1(1)=1;y2(1)=1;for n=1:10x(n+1)=x(n)+0.1;y1(n+1)=1.1*y1(n)-0.2*x(n)/y1(n);k1=y2(n)-2*x(n)/y2(n);k2=y2(n)+0.1*k1-2*x(n)/(y2(n)+0.1*k1);y2(n+1)=y2(n)+0.05*(k1+k2);endplot(x,y1,'k:',x,y2,'k-.')运行得到y =(2*x + 1)^(1/2),这是解析解。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
2、提交的纸质稿要求双面打印,中途提交批改不需要封面,但最后一次需将该课程所有实验项目内页与封面一起装订成册提交。
通过该实验的学习,使学生掌握最优化技术,认识面对什么样的实际问题,提出假设和建立优化模型,并且使学生学会使用MATLAB软件进行线性规划模型求解的基本命令,并进行灵敏度分析。解决现实生活中的最优化问题是本科生学习阶段中一门重要的课程,因此,本实验对学生的学习尤为重要。
实验内容
1.最优化问题的提出,提出不同的假设可以建立不同的最优化模型;
X41 6.521739 0.000000
X12 3.913043 0.000000
X42 0.000000 0.3036000E-01
X13 0.000000 0.000000
X43 0.000000 0.000000
X14 4.500000 0.000000
X44 0.000000 0.2640000E-01
总的收回:x11(1+15%) +x42(1+6%)
第三年初:可投资金额:x11(1+15%) +x42(1+6%)可投资项目:A、B、D项目
A的投资金额:x13(将在第四年末收回)B的投资金额:x2(将在第五年末收回且x2<4)
D的投资金额:x43
则必有x13+x2+x43= x11(1+15%) +x42(1+6%)
Max y= x3(1+40%)+ x2(1+25%)+ x14(1+15%)+[ x13(1+15%) +x44(1+6%)]*(1+6%)
整理得Max y=1.4x3+1.25x2+1.15x14+1.219x13+1.1236x44
约束条件整理为:
x11+x41=10
x12+x3+x42=x41(1+6%)
7 0.000000 0.000000
8 0.000000 0.000000
9 0.000000 0.000000
10 4.500000 0.000000
11 0.000000 0.000000
即最有投资策略为:第一年分别投资A项目3.4783万元D项目6.5217万元
第二年分别投资A项目3.9130万元C项目3万元
Row Slack or Surplus Dual Price
1 14.37500 1.000000
2 0.000000 0.7750000E-01
3 0.000000 0.3100000E-01
4 0.000000 1.150000
5 6.521739 0.000000
6 3.913043 0.000000
第三年末:收回第二年A项目的本利息:x12(1+15%)第三年D项目的本利息:x43(1+6%)总的收回x12(1+15%) +x43(1+6%)
第四年初:可投资金额:x12(1+15%) +x43(1+6%)可投资项目:A、D项目
A的投资金额:x14(第五年末收回)D的投资金额:x44
则必有x14+x44= x12(1+15%) +x43(1+6%)
第三年只投资B项目4万元第四年只投资A项目4.5万元
第五年获得的最大本利息为:14.37500
总结与体会
设计记录表格,包括碰到的问题汇总及解决情况
注行距:选最小值16磅,每一图应有简短确切的题名,连同图号置于图下。每一表应有简短确切的题名,连同表号置于表上。图表的题名及其中的文字采用小5号宋体。公式应该有编号,编号靠右端。
x41>0;
x12>0;
x42>0;
x13>0;
x43>0;
x14>0;
x44>0;
end
求解结果为:
Global optimal solution found.
Objective value: 14.37500
Total solver iterations: 2
Variable Value Reduced Cost
3、综合实验要求3人合作完成,请在实验报告上注明合作者的姓名。
则必有x11+x41=10
第一年末:收回D项目的本利息:x41*(1+6%)
第二年初:可投资金额:x41*(1+6%)可投资项目:A、C、D项目
A的投资金额:x12(将在第三年末收回)C的投资金额:x3(将在第五年末收回且x3<3)
D的投资金额:x42
则必有x12+x3+x42=x41(1+6%)
第二年末:收回第一年A项目的本利息:x11(1+15%)第二年D的本利息:x42(1+6%)
开课学院、实验室:实验时间:年月日
课程
名称
实验项目
名称
实验项目类型
验证
演示
综合
设计
其他
指导
教师
成绩
实验目的
[1]学习最优化技术和基本原理,了解最优化问题的分类;
[2]掌握线性规划的建模技巧和求解方法;
[3]学习灵敏度分析问题的思维方法;
[4]熟悉MATLAB软件求解线性规划模型的基本命令;
[5]通过范例学习,熟悉建立线性规划模型的基本要素和求解方法。
model:
max= 0.0465*x41-1.4*x12-0.075*x42-0.031*x13-1.25*x43+1.15*x14+1.1236*x44+14.3750;
1.06*x41-x12-x42-3<0;
1.06*x42-1.15*x41-x43-x13+7.5<0;
x14+x44-1.15*x12-1.06*x43=0;
[x,f]=linprog(c,A,b,[],[],zeros(2,1))
输出:
Optimization terminated.
x =
0.0000
15.0000
3.0000
f =
-78.0000
应用实验(或综合实验)
3.投资策略
某部门现有资金10万元,五年内有以下投资项目可供选择:
项目A:从第一年到第四年每年初投资,次年末收回本金且获利15%;
第四年A项目的本利息:x14(1+15%)第五年D项目的本利息;[ x13(1+15%) +x44(1+6%)]*(1+6%)
则五年的本息总额为:
y= x3(1+40%)+ x2(1+25%)+ x14(1+15%)+[ x13(1+15%) +x44(1+6%)]*(1+6%)
由以上分析整个模型的目标函数为:
项目B:第三年初投资,第五年末收回本金且获利25%,最大投资额为4万元;
项目C:第二年初投资,第五年末收回本金且获利40%,最大投资额为3万元;
项目D:每年初投资,年末收回本金且获利6%;
问如何确定投资策略使第五年末本息总额达最大?
第一年初:可投资金额:10万可投资项目:A、D项目
A的投资金额:x11(将在第二年末收回)D的投资金额:x41
第四年末:收回第三年A项目的本利息:x13(1+15%)第四年D项目的本利息:x44(1+6%)总的收回x13(1+15%) +x44(1+6%)
第五年初:可投资金额:x13(1+15%) +x44(1+6%)可投资项目:D项目
D的投资金额:x45=x13(1+15%) +x44(1+6%)
第五年末:收回C项目的本利息:x3(1+40%) B项目本利息:x2(1+25%)
x14+x44-1.15x12-1.06x43=0
将 带入 得x2=1.06x42-1.15x41-x43-x13+11.5<4
将 式带入目标函数得
Maxy= 0.0465x41-1.4x12-0.075x42-0.031x13-1.25x43+1.15x14+1.1236x44+14.3750;
程序:
2.建立线性规划模型的基本要素和步骤;
3.使用MATLAB命令对线性规划模型进行计算与灵敏度分析;
4.利用优化数值解与图形解对最优化特征作定性与定量分析;
基础实验
1.求解下述线性规划问题
min
s.t.
程序:
c=[-5 -4 -6];
A=[1 -1 1;3 2 4;3 2 0];
b=[20;42;30];
x13+x2+x43= x11(1+15%) +x42(1+6%)
x14+x44= x12(1+15%) +x43(1+6%)
x3<3
x2<4 (所有变量均大于0)
首先将上述约束条件转化为x11=10-x41
x3=1.06x41-x12-x42<3
x2=1.15x11+1.06x42-x13-x43<4