数学建模实验四
数学建模实验报告

湖南城市学院数学与计算科学学院《数学建模》实验报告专业:学号:姓名:指导教师:成绩:年月日目录实验一 初等模型........................................................................ 错误!未定义书签。
实验二 优化模型........................................................................ 错误!未定义书签。
实验三 微分方程模型................................................................ 错误!未定义书签。
实验四 稳定性模型.................................................................... 错误!未定义书签。
实验五 差分方程模型................................................................ 错误!未定义书签。
实验六 离散模型........................................................................ 错误!未定义书签。
实验七 数据处理........................................................................ 错误!未定义书签。
实验八 回归分析模型................................................................ 错误!未定义书签。
实验一 初等模型实验目的:掌握数学建模的基本步骤,会用初等数学知识分析和解决实际问题。
实验内容:A 、B 两题选作一题,撰写实验报告,包括问题分析、模型假设、模型构建、模型求解和结果分析与解释五个步骤。
数学建模实验报告4

数学建模实验报告班级:姓名:学号:元件可靠性问题一、实验问题:给出3种不同情况的元件连接方式,分别求解他们的正常运行概率。
其中每个元件的正常运行概率均为p。
元件数为N,方式2与方式3用到了与A元件相同的N个B元件。
连接方式如图:方式1:方式2:方式3:二、问题分析:N个元件的连接方式,相当于电阻的串并联,所以可以用电阻串并联的关系去分析各无件之间的关系:对于方式一来说,相当于电阻的串联。
所以,他的正常运行的概率为p^n.对于方式二来说,相当于电阻先串联再并联。
所以,他的正常运行的概率为:1-(1-P^n)(1-P^n)=2P^n-P^2n.对于方式三来说,相当于电阻先并联再串联。
所以,他的正常运行的概率为:(1-(1-P^n)^2)^n=(2p-p^2)^n现在再比较三个系统正常工作概率大小P1- P2= p^n–(2p^n-p^2n )= p^2n–p^n 由于0<p<1,所以易知P^2n-P^n<0。
所以有P1< P2P2- P3=(2p^n- p^2n)- (2p-p^2)^n= p^n[(2- p^n)-(2-p)^n]因为p^n>0,所以只要比较[(2- p^n)-(2-p)^n]大小即可。
对此式求导有-n[p^(n-1)-(2-p)^n-1]可见此式恒大于零,所以函数单调递增。
当p=1时,[(2- p^n)-(2-p)^n]=0.所以P2- P3 <0,再由上求导可知所以P2<P3所以P3最大。
即其的可靠性最高。
理发店问题一、实验题目:某单人理发店有4反椅子接待顾客排队理发,当4把椅子都坐满人时,后来的顾客就不进店而离去。
顾客平均到达速率为4人/H,理发时间平均10min/人。
设到达过程为泊松流,服务时间服从负指数颁布。
求:(1)顾客一到达就能理发的概率;(2)系统中顾客数的期望值和排队等待顾客数的期望值;(3)顾客在理发店内逗留的全部时间的期望值;(4)在可能到达的顾客中因客满离开的概率。
数学建模基础实验报告(3篇)

第1篇一、实验目的本次实验旨在让学生掌握数学建模的基本步骤,学会运用数学知识分析和解决实际问题。
通过本次实验,培养学生主动探索、努力进取的学风,增强学生的应用意识和创新能力,为今后从事科研工作打下初步的基础。
二、实验内容本次实验选取了一道实际问题进行建模与分析,具体如下:题目:某公司想用全行业的销售额作为自变量来预测公司的销售量。
表中给出了1977—1981年公司的销售额和行业销售额的分季度数据(单位:百万元)。
1. 数据准备:将数据整理成表格形式,并输入到计算机中。
2. 数据分析:观察数据分布情况,初步判断是否适合使用线性回归模型进行拟合。
3. 模型建立:利用统计软件(如MATLAB、SPSS等)进行线性回归分析,建立公司销售额对全行业的回归模型。
4. 模型检验:对模型进行检验,包括残差分析、DW检验等,以判断模型的拟合效果。
5. 结果分析:分析模型的拟合效果,并对公司销售量的预测进行评估。
三、实验步骤1. 数据准备将数据整理成表格形式,包括年份、季度、公司销售额和行业销售额。
将数据输入到计算机中,为后续分析做准备。
2. 数据分析观察数据分布情况,绘制散点图,初步判断是否适合使用线性回归模型进行拟合。
3. 模型建立利用统计软件进行线性回归分析,建立公司销售额对全行业的回归模型。
具体步骤如下:(1)选择合适的统计软件,如MATLAB。
(2)输入数据,进行数据预处理。
(3)编写线性回归分析程序,计算回归系数。
(4)输出回归系数、截距等参数。
4. 模型检验对模型进行检验,包括残差分析、DW检验等。
(1)残差分析:计算残差,绘制残差图,观察残差的分布情况。
(2)DW检验:计算DW值,判断随机误差项是否存在自相关性。
5. 结果分析分析模型的拟合效果,并对公司销售量的预测进行评估。
四、实验结果与分析1. 数据分析通过绘制散点图,观察数据分布情况,初步判断数据适合使用线性回归模型进行拟合。
2. 模型建立利用MATLAB进行线性回归分析,得到回归模型如下:公司销售额 = 0.9656 行业销售额 + 0.01143. 模型检验(1)残差分析:绘制残差图,观察残差的分布情况,发现残差基本呈随机分布,说明模型拟合效果较好。
数学建模实验报告

《数学建模实验》实验报告学院名称数学与信息学院专业名称提交日期课程教师实验一:数学规划模型AMPL求解实验内容1. 用AMPL求解下列问题并作灵敏度分析:一奶制品加工厂用牛奶生产A1和A2两种奶制品,1桶牛奶可以在甲类设备上用12小时加工成3公斤A1或者在乙类设备上用8小时加工成4公斤A2,且都能全部售出,且每公斤A1获利24元,每公斤A2获利16元。
先加工厂每天能得到50桶牛奶的供应,每天工人总的劳动时间为480小时,并且甲类设备每天至多加工100公斤A1,乙类设备的加工能力没有限制,试为该厂制定一个计划,使每天的获利最大。
(1)建立模型文件:milk.modset Products ordered;param Time{i in Products }>0;param Quan{i in Products}>0;param Profit{i in Products}>0;var x{i in Products}>=0;maximize profit: sum{i in Products} Profit [i]* Quan [i]*x[i];subject to raw: sum{i in Products}x[i] <=50;subject to time:sum{i in Products}Time[i]*x[i]<=480;subject to capacity: Quan[first(Products)]*x[first(Products)]<=100;(2)建立数据文件milk.datset Products:=A1 A2;param Time:=A1 12 A2 8;param Quan:=A1 3 A2 4;param Profit:=A1 24 A2 16;(3) 建立批处理文件milk.runmodel milk.mod;data milk.dat;option solver cplex;solve;display x;(4)运行运行结果:CPLEX 11.0.0: optimal solution; objective 33602 dual simplex iterations (1 in phase I)x [*] :=A1 20A2 30;(5)灵敏度分析:model milk.mod;data milk.dat;option solver cplex;option cplex_options 'sensitivity';solve;display x;display x.rc, x.down, x.up;display raw, time, capacity;display raw.down, raw.up,raw.current, raw.slack;得到结果:【灵敏度分析】: x.rc x.down x.up:=A1 -3.55271e-15 64 96A2 0 48 72;raw = 48time = 2capacity = 0raw.down = 43.3333raw.up = 60raw.current = 50raw.slack = 0某公司有6个建筑工地,位置坐标为(a i, b i)(单位:公里),水泥日用量d i (单位:吨)1) 现有j j j吨,制定每天的供应计划,即从A, B两料场分别向各工地运送多少吨水泥,使总的吨公里数最小。
数学建模——混合整数规划

实验四 混合整数规划一、问题重述某开放式基金现有总额为15亿元的资金可用于投资,目前共有8个项目可供投资者选择,每个项目可重复投资。
根据专家经验,对每个项目投资总额不能太高,应有上限。
这些项目所需要的投资额已知,一般情况下投资一年后各项目所得利润也可估算出来,如表1所示。
请帮该公司解决以下问题:(1) 就表1提供的数据,应该投资哪些项目,使得第一年所得利润最高?(2) 在具体投资这些项目时,实际还会出现项目之间互相影响的情况。
公司咨询有关专家后,得到以下可靠信息:同时投资项目A 1,A 3,它们的年利润分别是1005万元,1018.5万元;同时投资项目A 4,A 5,它们的年利润分别是1045万元,1276万元;同时投资项目A 2,A 6,A 7,A 8,它们的年利润分别是1353万元,840万元,1610万元,1350万元,该基金应如何投资? 其中M 为你的学号后3位乘以10。
(3) 如果考虑投资风险,则应如何投资,使收益尽可能大,而风险尽可能小。
投资项目总体风险可用投资项目中最大的一个风险来衡量。
专家预测出各项目的风险率,如表2所示。
二、符号说明i A ::投资额;i b :i A 个项目所获得的年利润;i C :第i A 个项目投资所获得的利润; 'i C :第i A 个项目同时投资所获得的利润;i m :投资i A 的上限; i y :表示0—1变量;i p :投资第i A 个项目的投资风险;三、模型的建立 对于问题一目标函数:81max i i i c x ==∑s.t. 150000i i i i i ib x b x m ⎧≤⎪⎨⎪≤⎩∑对于问题二 设定0—1变量131130...,1...,A A y A A ⎧⎨⎩项目不同时投资项目同时投资 452450...,1...,A A y A A ⎧⎨⎩项目不同时投资项目同时投资 2678326780...,,1...,,A A A A y A A A A ⎧⎨⎩,项目不同时投资,项目同时投资 目标函数:''''11133111332445524455''''322667788322667788max ()(1)()()(1)()()(1)()y x c x c y x c x c y x c x c y x c x c y x c x c x c x c y x c x c x c x c =++-++++-++++++-+++s.t. 11313124545232678267831500001000i i i i i ib x k y x xx x y ky x x x x y k y x x x x x x x x y kb x m ⎧≤⎪⎪=⎪⎪≤⎪⎪≥⎪⎪≤⎨⎪⎪≥⎪⎪≤⎪⎪≥⎪⎪≤⎩∑对于问题三:目标函数:max min max()i iii i i c x b x p =∑s.t. 150000i i i i i ib x b x m ⎧≤⎪⎨⎪≤⎩∑对于问题三模型的简化固定投资风险,优化收益,设a 为固定的最大风险。
数学建模的实验报告

数学建模实验报告姓名:学院:专业班级:学号:数学建模实验报告(一)——用最小二乘法进行数据拟合一.实验目的:1.学会用最小二乘法进行数据拟合。
2.熟悉掌握matlab软件的文件操作和命令环境。
3.掌握数据可视化的基本操作步骤。
4.通过matlab绘制二维图形以及三维图形。
二.实验任务:来自课本64页习题:用最小二乘法求一形如y=a+b x2的多项式,使之与下列数据拟合:三.实验过程:1.实验方法:用最小二乘法解决实际问题包含两个基本环节:先根据所给出数据点的变化趋势与问题的实际背景确定函数类;然后按照最小二乘法原则求最小二乘解来确定系数。
即要求出二次多项式: y=a+b x2的系数。
2.程序:x=[19 25 31 38 44]y=[19.0 32.3 49.0 73.3 97.8]ab=y/[ones(size(x));x.^2];a=ab(1),b=ab(2)xx=19:44;plot(xx,a+b*xx.^2,x,y,'.')3.上机调试得到结果如下:x = 19 25 31 38 44y=19.0000 32.3000 49.0000 73.3000 97.8000a = 0.9726b = 0.0500图形:四.心得体会通过本次的数学模型的建立与处理,我们学习并掌握了用最小二乘法进行数据拟合,及多项式数据拟合的方法,进一步学会了使用matlab软件,加深了我们的数学知识,提高了我们解决实际问题的能力,为以后深入学习数学建模打下了坚实的基础。
数学建模实验报告(二)——用Newton法求方程的解一.实验目的1.掌握Newton法求方程的解的原理和方法。
2.利用Matlab进行编程求近似解。
二.实验任务来自课本109页习题4-2:用Newton法求f(x)=x-cosx=0的近似解三.实验过程1.实验原理:把f(x)在x0点附近展开成泰勒级数f(x) = f(x0)+(x-x0)f'(x0)+(x-x0)^2*f''(x0)/2! +… 取其线性部分,作为非线性方程f(x) = 0的近似方程,即泰勒展开的前两项,则有f(x0)+f'(x0)(x-x0)=0 设f'(x0)≠0则其解为x1=x0-f(x0)/f'(x0) 这样,得到牛顿法的一个迭代序列:x(n+1)=x(n)-f(x(n))/f'(x(n))。
数学建模实验报告

内江师范学院中学数学建模实验报告册编制数学建模组审定牟廉明专业:班级:级班学号:姓名:数学与信息科学学院2016年3月说明1.学生在做实验之前必须要准备实验,主要包括预习与本次实验相关的理论知识,熟练与本次实验相关的软件操作,收集整理相关的实验参考资料,要求学生在做实验时能带上充足的参考资料;若准备不充分,则学生不得参加本次实验,不得书写实验报告;2.要求学生要认真做实验,主要是指不得迟到、早退和旷课,在做实验过程中要严格遵守实验室规章制度,认真完成实验内容,极积主动地向实验教师提问等;若学生无故旷课,则本次实验成绩不合格;3.学生要认真工整地书写实验报告,实验报告的内容要紧扣实验的要求和目的,不得抄袭他人的实验报告;4.实验成绩评定分为优秀、合格、不合格,实验只是对学生的动手能力进行考核,跟据所做的的情况酌情给分。
根据实验准备、实验态度、实验报告的书写、实验报告的内容进行综合评定。
实验名称:数学规划模型(实验一)指导教师:实验时数: 4 实验设备:安装了VC++、mathematica、matlab的计算机实验日期:年月日实验地点:实验目的:掌握优化问题的建模思想和方法,熟悉优化问题的软件实现。
实验准备:1.在开始本实验之前,请回顾教科书的相关内容;2.需要一台准备安装Windows XP Professional操作系统和装有数学软件的计算机。
实验内容及要求原料钢管每根17米,客户需求4米50根,6米20根,8米15根,如何下料最节省?若客户增加需求:5米10根,由于采用不同切割模式太多,会增加生产和管理成本,规定切割模式不能超过3种,如何下料最节省?实验过程:摘要:生活中我们常常遇到对原材料进行加工、切割、裁剪的问题,将原材料加工成所需大小的过程,称为原料下料问题。
按工艺要求,确定下料方案,使用料最省,或利润最大是典型的优化问题。
以此次钢管下料问题我们采用数学中的线性规划模型.对模型进行了合理的理论证明和推导,然后借助于解决线性规划的专业软件Lingo 对题目所提供的数据进行计算从而得出最优解。
数学建模计算实验

学时:4学时 实验目的:掌握用Lindo求解线性规划问题的方法,能够阅读Lindo结果 报告。
实验内容:
解:
实例2:求解书本上P130的习题1。列出线性规划模型,然后用
Lindo求解,根据结果报告得出解决方案。
投资规划问题
某银行经理计划用一笔资金进行有价证券的投资,可供购进的证券
以及其信用等级、到期年限、收益如下表所示。按照规定,市政证券的
解:设投资证券A,B,C,D的金额分别为
(百万元),按照
规定限制1000万元的资金约束,则线性规划模型为:
0.043 +0.054*0.5 +0.050*0.5 +0.044*0.5 +0.045
实验三:用Lingo求解非线性规划问题
学时:2学时 实验目的:掌握用Lingo求解非线性规划问题的方法。 实验内容:
考虑如下的在线DVD租赁问题。顾客缴纳一定数量的月费成为会 员,订购DVD租赁服务。会员对哪些DVD有兴趣,只要在线提交订 单,网站就会通过快递的方式尽可能满足要求。会员提交的订单包括多 张DVD,这些DVD是基于其偏爱程度排序的。网站会根据手头现有的 DVD数量和会员的订单进行分发。每个会员每个月租赁次数不得超过2 次,每次获得3张DVD。会员看完3张DVD之后,只需要将DVDa放进网 站提供的信封里寄回(邮费由网站承担),就可以继续下次租赁。请考 虑以下问题:
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
Integer variables: 9
Total constraints: 10
Nonlinear constraints: 0
Total nonzeros: 38
Nonlinear nonzeros: 0
Variable Value Reduced Cost
问题三求解结果为目标函数值是16,此处6门课程分别是x1=x2=x3=x5=x7=x9=1,其余为零,学分最高为22,答案同问题二。
附录:
问题一程序
model:
min=x1+x2+x3+x4+x5+x6+x7+x8+x9;
x1+x2+x3+x4+x5>2;
x3+x5+x6+x8+x9>3;
x4+x6+x7+x9>2;
1 、以选课门数最少为目标建立选课策略数学规划模型。
2 、在选课门数最少的前提下,以学分最多为目标建立选课策略数学规划模型。
3 、既要选课门数少,又要学分多,建立2个目标的多目标规划模型,并用线性加权组合的方式将多个目标处理成一个新目标, 化为单目标规划进行求解。
问题分析与假设:
上述问题为数学规划中的多目标规划问题,可以建立0-1规划模型进行求解,设变量Xi表示课号,令Xi=1表示选课,Xi=0表示不选,然后根据题目建立目标函数和约束条件,最后用lingo软件求解。
对于第一个问题: 以选课门数最少为目标建立选课策略数学规划模型,可以直接列目标函数记为最小值求解;
对于第二个问题:在选课门数最少的前提下,以学分最多为目标建立选课策略数学规划模型,可以把问题一选课最少求解结果作为已知量再去建立学分最多目标函数记为最大值;
对于第三个个目标的多目标规划模型,并用线性加权组合的方式将多个目标处理成一个新目标, 化为单目标规划进行求解。
x4+x6+x7+x9>2;
2*x3-x1-x2<0;
x4-x7<0;
2*x5-x1-x2<0;
x6-x7<0;
x8-x5<0;
2*x9-x1-x2<0;
x1+x2+x3+x4+x5+x6+x7+x8+x9<6;
@bin(x1);@bin(x2);@bin(x3);@bin(x4);@bin(x5);@bin(x6);@bin(x7);@bin(x8);@bin(x9);
将以上目标函数和约束条件所建立的模型输入到lingo中求解(注意加上xi为0-1约束)
对于第二个问题,在选课最少的前提下,即选课门数为6,使得学分最高。目标函数即:
约束条件除问题一以上约束外再增加选课门数约束,即:
再次输入lingo求解
对于第三个问题,既要选课门数少,又要学分多,建立2个目标的多目标规划模型,并用线性加权组合的方式将多个目标处理成一个新目标, 化为单目标规划进行求解。所以建立目标函数为
end
程序结果:
Global optimal solution found.
Objective value:
Objective bound:
Infeasibilities:
Extended solver steps: 0
Total solver iterations: 0
Elapsed runtime seconds:
约束条件如问题一,此处略。输入lingo后求解。
结果分析与检验:
问题一求解结果为x1=x2=x3=x6=x9=1,其他变量为0.对照编号,若使选课最少,可选择微积分,线性代数,最优化方法,计算机模拟,计算机编程,数学实验,共6门课程,总学分21.
问题二得到结果为:在选课门数最少的前提下,学分最高为22,此处6门课程分别是x1=x2=x3=x5=x7=x9=1,其余为零。
X1
X2
X3
X4
X5
X6
X7
X8
X9
Row Slack or Surplus Dual Price
1
2
3
4
5
6
7
8
9
10
问题二程序:
model:
max=5*x1+4*x2+4*x3+3*x4+4*x5+3*x6+2*x7+2*x8+3*x9;
x1+x2+x3+x4+x5>2;
x3+x5+x6+x8+x9>3;
模型的建立和求解:
用Xi=1表示选修表中按编号顺序的9门课程(Xi=0表示不选;i=1,2…,9)。问题的目标为选修的课程总数最少,即
约束条件包括两个方面:
第一,每人最少要学习2门数学课,3门运筹课和2门计算机课,根据表中对每门课程所属类别的划分,这一约束可以表示为:
第二,某些课程要有先修的要求,例如“数据结构”的先修课程“计算机编程”,意味着如果x4=1必须x7=1,可以表示为x4 x7,同里其他先修课程亦是如此。其中x1 x3,x2 x3可以用一个约束表示为2x3-x2-x1 0,故所有课程的先修课程要求可以表示为以下约束:
数学建模实验四
《数学建模》实验报告
实验序号:实验四日期:2015年6月3日
班级
遥感1402班
姓名
郑永杰
学号
3
实验名称 整数规划与lingo程序设计 选课策略
实验目的:
1)学会lingo软件的基本使用方法。
2)加强学生对0-1变量的认识和理解;培养学生建模的能力。
3)了解多目标规划模型求解基本原理。
实验内容及要求:
Model Class: PILP
Total variables: 9
Nonlinear variables: 0
Integer variables: 9
Total constraints: 11
Nonlinear constraints: 0
Global optimal solution found.
Objective value:
Objective bound:
Infeasibilities:
Extended solver steps: 0
Total solver iterations: 0
Model Class: PILP
Total variables: 9
2*x3-x1-x2<0;
x4-x7<0;
2*x5-x1-x2<0;
x6-x7<0;
x8-x5<0;
2*x9-x1-x2<0;
@bin(x1);@bin(x2);@bin(x3);@bin(x4);@bin(x5);@bin(x6);@bin(x7);@bin(x8);@bin(x9);
end
程序结果