数学建模作业43508

合集下载

08年全国大学生数学建模A题文档

R 旋转矩阵t 平移矢量K 相机标定矩阵AB（粗体）向量AB，粗体表示向量3尺寸、形状等信息，我们就可以确定靶标在三维空间中与相机的位置关系。

确定了靶标与相机的位置关系后，就可以很容易的将靶标上的点投射到靶标像平面上，当然也包括五个圆心。

在求取出五个圆心的空间坐标之后，将其投射到像平面上，就得到了第一问需要求的坐标。

在求解时涉及到一个问题，就是像平面上怎样确定切线。

因为像平面上的图形是不规则的，所以很难确定这些形状的切线。

因此我们考虑另外的方法，使用搜索的办法，利用模拟退火算法求解。

在如何检验模型的问题上，需要分两方面进行检验，一是精度，而是稳定性。

按照以上的方法求圆心在像平面上的坐标，并没有充分利用像平面上所有轮廓点的信息，因此可以利用这些点来检验模型的精度。

对于稳定性问题，可以采用计算机模拟的方法，随机修改图形的轮廓，并用以上的方法再次进行求解，通过比较修改前后的结果来分析模型的稳定性。

最后，考虑另外一台相机的定位相对位置问题。

根据前面模型，我们应能够对任意一台相机确定靶标相对它的位置，因此可以以这个靶标作为参照物，建立一个世界坐标系，将这两台相机的位置在这个坐标系里面表示出来，以此确定两台相机的相对位置。

四、模型假设1、假设靶标像的中心恰好在光轴上2、假设数码相机中图像平面与光轴垂直3、假设相机两个方向上焦距相等4、假设透镜的焦距很小，像距约等于焦距五、模型准备(一)靶图像矩阵表示首先将题目中的图片保存出来，得到的图像可以很方便的放到Matlab 里面进行处理。

但在处理之前还要进行进一步加工：(1) 将文件读入Matlab，使用imread()函数(2) 将矩阵变为0-1 矩阵对于用以上方式得到的矩阵，有两个值：0、15。

其中0 代表像素为白色的点，15代表像素为黑色的点。

为了方便下面处理，对需要把以上像素为15 的点值全部变为1。

以上两步的源代码见附录一。

4(二)图像轮廓的提取在提取图像轮廓时，首先要引入计算机图像处理技术中四邻域的概念。

数学建模案例分析-- 模糊数学方法建模5习题八

习题八
1、试举一个模糊综合评判的实际例子（由建立评判的因素集、评价集到综合评判结果）。

2、试对新式服装作模糊综合评判
设因素集U ={花色，式样，质料，舒适度，价格}
评价集V ={很欢迎，欢迎，不大欢迎，不欢迎} 权重分配)35.0,15.0,3.0,12.0,1.0(~
=A ，单因素评价矩阵
⎪⎪⎪⎪⎪⎪⎭
⎫
⎝
⎛=02
.03.05.03.06.01.001.05.04.00
1.05.03.02
.003.05.02.0~R 3、已知由于脾产生的肝病为模糊集合~
1A ，由于肝炎产生的肝病为模糊集合~
2A ，现有一病人的症状为模糊集合~
B ，试用模糊模式识别方法诊断该病人属于哪一种类型的肝病？
~
1A =1/(腹胀)+0.9/(乏力)+0.8/(肠呜)+0.1/(口干)+0.2/(苔薄黄)+0.85/(脉沉缓)
~
2A =0.7/(腹胀)+1/(乏力)+0.5/(肠呜)+0.9/(口干)+0.85/(苔薄黄)+1/(脉沉缓)
~
B =0.6/(腹胀)+1/(乏力)+0.2/(肠呜)+0.8/(口干)+0.5/(苔薄黄)+0.9/(脉沉缓)
4、松毛虫的每一个生态环境都具有一定的生态条件,它和气候、植被、土壤、地形、天敌等共同
构成自然地理景观。

从湖南省的考察资料中抽出8个地区的资料作为分类样本，选取6个主要因
试用模糊聚类法对这8个地区进行分类。

5、经多年摸索发现影响产生第一代玉米螟的主要因素是4月下旬至5月10日雨量及4月份雨日
12
等级划分如下表。

看看与实际情况是否吻合？。

数学建模答案(完整版)

1 建立一个命令M 文件:求数60.70.80,权数分别为1.1，1.3,1.2的加权平均数。

在指令窗口输入指令edit ，打开空白的M 文件编辑器；里面输入s=60*1.1+70*1.3+80*1.2;ave=s/3然后保存即可2 编写函数M 文件SQRT.M;函数()f x = x=567.889与0.0368处的近似值（保留有效数四位）在指令窗口输入指令edit ，打开空白的M 文件编辑器；里面输入syms x1 x2 s1 s2 zhi1 zhi2x1=567.889;x2=0.368;s1=sqrt(x1);s2=sqrt(x2);zhi1=vpa(s1,4)zhi2=vpa(s2,4)然后保存并命名为SQRT.M 即可3用matlab 计算()f x a b=-的值,其中a=2.3,b=4.89. >> syms a b>> a=2.3;b=4.89;>> sqrt(a^2+b^2)/abs(a-b)ans =2.08644用matlab 计算函数()f x =x=3π处的值. >> syms x>> x=pi/3;>> sqrt(sin(x)+cos(x))/abs(1-x^2)ans =12.09625用matlab 计算函数()arctan f x x =在x=1.23处的值.>> syms x>> x=1.23;>> atan(x)+sqrt(log(x+1))ans =1.78376 用matlab 计算函数222sin cos ()()31a b x x f x f x a b xπ++==--在x=-2.1处的值. >> syms x >> x=-2.1;>> 2-3^x*log(abs(x))ans =1.92617 用蓝色.点连线.叉号绘制函数y=2x 在[0,2]上步长为0.1的图像.>> syms x y>> x=0:0.2:2;y=2*sqrt(x);>> plot(x,y,'b.-')8 用紫色.叉号.实连线绘制函数ln 10y x =+在[20,15]--上步长为0.2的图像. >> syms x y>> x=-20:0.2:-15;y=log(abs(x+10));>> plot(x,y,'mx-')ln 10[20,15]y x x=+--9 用红色.加号连线虚线绘制函数sin()22x y π=-在[-10,10]上步长为0.2的图像. >> syms x y;>> x=-10:0.2:10;y=sin(x/2-pi/2);>> plot(x,y,'r+--')10用紫红色.圆圈.点连线绘制函数sin(2)3y x π=+在[0,4]π上步长为0.2的图像.sin(2)sin()[0,4]322x y x y πππ=+=- >> syms x y >> x=0:0.2:4*pi;y=sin(2*x+pi/3);>> plot(x,y,'mo-.')11 在同一坐标中,用分别青色.叉号.实连线与红色.星色.虚连线绘制y=cos3x 与y x =.>> syms x y1 y2>> x=0:pi/50:2*pi;y1=cos(3*sqrt(x));y2=3*cos(sqrt(x));>> plot(x,y1,'cx-',x,y2,'r*--')12 在同一坐标系中绘制函数234,,y x y x y x ===这三条曲线的图标,并要求用两种方法加各种标注.234,,y x y x y x ===>> syms x y1 y2 y3;>> x=-2:0.1:2;y1=x.^2;y2=x.^3;y3=x.^4;plot(x,y1,x,y2,x,y3);13 作曲线2sinx ty tz t⎧=⎪=⎨⎪=⎩的3维图像>> syms x y t z>> t=0:1/50:2*pi; >> x=t.^2;y=sin(t);z=t;>> stem3(x,y,z)14 作环面(1cos)cos(1cos)sinsinx u vy u vz u=+⎧⎪=+⎨⎪=⎩在(0,2)(0,2)ππ⨯上的3维图像>> syms x y u v z>> u=0:pi/50:2*pi;v=0:pi/50:2*pi;>>x=(1+cos(u)).*cos(v);y=(1+cos(u)).*sin(v);z=sin(u); >> plot3(x,y,z)15 求极限0lim x +→0lim x +→>> syms x y>> y=sin(2^0.5*x)/sqrt(1-cos(x));>> limit(y,x,0,'right')ans =216 求极限1201lim ()3x x +→ >> syms y x>> y=(1/3)^(1/(2*x));>> limit(y,x,0,'right')ans =17求极限lim x>> syms x y>> y=(x*cos(x))/sqrt(1+x^3);>> limit(y,x,+inf)ans =18 求极限21lim ()1x x x x →+∞+- >> syms x y>> y=((x+1)/(x-1))^(2*x);>> limit(y,x,+inf)ans =exp(4)19 求极限01cos 2lim sin x x x x→->> syms x y>> y=(1-cos(2*x))/(x*sin(x));>> limit(y,x,0)ans =220 求极限 0x → >> syms x y>> y=(sqrt(1+x)-sqrt(1-x))/x;>> limit(y,x,0)ans =121 求极限2221lim 2x x x x x →+∞++-+ >> syms x y>> y=(x^2+2*x+1)/(x^2-x+2);>> limit(y,x,+inf)ans =122 求函数y=5(21)arctan x x -+的导数>> syms x y>> y=(2*x-1)^5+atan(x);>> diff(y)ans =10*(2*x - 1)^4 + 1/(x^2 + 1)23 求函数y=2tan 1x x y x=+的导数 >> syms y x>> y=(x*tan(x))/(1+x^2);>> diff(y)ans =tan(x)/(x^2 + 1) + (x*(tan(x)^2 + 1))/(x^2 + 1) - (2*x^2*tan(x))/(x^2 + 1)^224 求函数3tan x y e x -=的导数>> syms y x>> y=exp^(-3*x)*tan(x)>> y=exp(-3*x)*tan(x)y =exp(-3*x)*tan(x)>> diff(y)ans =exp(-3*x)*(tan(x)^2 + 1) - 3*exp(-3*x)*tan(x)25 求函数y=22ln sin 2xx π+在x=1的导数>> syms x y>> y=(1-x)/(1+x);>> diff(y,x,2)ans =2/(x + 1)^2 - (2*(x - 1))/(x + 1)^3>> syms x y>> y=2*log(x)+sin(pi*x/2)^2;>> dxdy=diff(y)dxdy =2/x + pi*cos((pi*x)/2)*sin((pi*x)/2)zhi=subs(dxdy,1)zhi =226 求函数y=01cos 2limsin x x x x →-11x x-+的二阶导数 >> syms x y>> y=(1-x)/(1+x);>> diff(y,x,2)ans =2/(x + 1)^2 - (2*(x - 1))/(x + 1)^327 求函数; >> syms x y>> y=((x-1)^3*(3+2*x)^2/(1+x)^4)^0.2;>> diff(y)ans =(((8*x + 12)*(x - 1)^3)/(x + 1)^4 + (3*(2*x + 3)^2*(x - 1)^2)/(x + 1)^4 - (4*(2*x + 3)^2*(x -1)^3)/(x + 1)^5)/(5*(((2*x + 3)^2*(x - 1)^3)/(x + 1)^4)^(4/5))28在区间(,-∞+∞)内求函数43()341f x x x =-+的最值.>> f='-3*x^4+4*x^3-1';>> [x,y]=fminbnd(f,-inf,inf)x =NaNy =NaN>> f='3*x^4-4*x^3+1';>> [x,y]=fminbnd(f,-inf,inf)x =NaNy =NaN29在区间(-1,5)内求函数发()(f x x =-的最值.>> f='(x-1)*x^0.6';>> [x,y]=fminbnd(f,-1,5)x =0.3750y =-0.3470>>>> f='-(x-1)*x^0.6';>> [x,y]=fminbnd(f,-1,5)x =4.9999y =-10.505930 求不定积分(ln32sin )x x dx -⎰(ln32sin )x x dx -⎰ >> syms x y>> y=log(3*x)-2*sin(x);>> int(y)ans =2*cos(x) - x + x*log(3) + x*log(x)31求不定积分2sin x e xdx ⎰>> syms x y>> y=exp(x)*sin(x)^2;>> int(y)ans =-(exp(x)*(cos(2*x) + 2*sin(2*x) - 5))/1032. 求不定积分>> syms x y>> y=x*atan(x)/(1+x)^0.5;>> int(y)Warning: Explicit integral could not be found.ans =int((x*atan(x))/(x + 1)^(1/2), x)33.计算不定积分2(2cos )x x x e dx --⎰ >> syms x y>> y=1/exp(x^2)*(2*x-cos(x));>> int(y)Warning: Explicit integral could not be found.ans =int(exp(-x^2)*(2*x - cos(x)), x)34.计算定积分10(32)x e x dx -+⎰>> syms x y>> y=exp(-x)*(3*x+2);>> int(y,0,1)ans =5 - 8*exp(-1)10(32)x e x dx -+⎰35.计算定积分0lim x x→120(1)cos x arc xdx +⎰>> syms y x>> y=(x^2+1)*acos(x);>> int(y,0,1)ans =11/936.计算定积分10cos ln(1)x x dx +⎰>> syms x y>> y=(cos(x)*log(x+1));>> int(y,0,1)Warning: Explicit integral could not be found.ans =int(log(x + 1)*cos(x), x == 0..1)37计算广义积分2122x x dx +∞++-∞⎰;>> syms y x>> y=(1/(x^2+2*x+2));>> int(y,-inf,inf)ans =pi38.计算广义积分20x dx x e +∞-⎰； >> syms x y>> y=x^2*exp(-x);>> int(y,0,+inf)ans =2。

数学建模题目及答案-数学建模100题

资料范本本资料为word版本，可以直接编辑和打印，感谢您的下载数学建模题目及答案-数学建模100题地点：__________________时间：__________________说明：本资料适用于约定双方经过谈判，协商而共同承认，共同遵守的责任与义务，仅供参考，文档可直接下载或修改，不需要的部分可直接删除，使用时请详细阅读内容09级数模试题1. 把四只脚的连线呈长方形的椅子往不平的地面上一放，通常只有三只脚着地，放不稳，然后稍微挪动几次，就可以使四只脚同时着地，放稳了。

试作合理的假设并建立数学模型说明这个现象。

（15分）解：对于此题，如果不用任何假设很难证明，结果很可能是否定的。

因此对这个问题我们假设：（1）地面为连续曲面（2）长方形桌的四条腿长度相同（3）相对于地面的弯曲程度而言，方桌的腿是足够长的（4）方桌的腿只要有一点接触地面就算着地。

那么，总可以让桌子的三条腿是同时接触到地面。

现在，我们来证明：如果上述假设条件成立，那么答案是肯定的。

以长方桌的中心为坐标原点作直角坐标系如图所示，方桌的四条腿分别在A、B、C、D 处，A、B,C、D的初始位置在与x轴平行，再假设有一条在x轴上的线ab,则ab也与A、B，C、D平行。

当方桌绕中心0旋转时，对角线 ab与x轴的夹角记为。

容易看出，当四条腿尚未全部着地时，腿到地面的距离是不确定的。

为消除这一不确定性，令为A、B离地距离之和，为C、D离地距离之和，它们的值由唯一确定。

由假设（1），,均为的连续函数。

又由假设（3），三条腿总能同时着地，故=0必成立（）。

不妨设,g（若也为0，则初始时刻已四条腿着地，不必再旋转），于是问题归结为：已知,均为的连续函数，,且对任意有，求证存在某一，使。

证明：当θ=π时，AB与CD互换位置，故，。

作，显然，也是的连续函数，而，由连续函数的取零值定理，存在，，使得，即。

又由于，故必有，证毕。

2.学校共1000名学生，235人住在A宿舍，333人住在B宿舍，432人住在C宿舍。

数学建模作业43508

数学建模作业1、在甲乙双方的一场战争中，部分甲方部队被乙方部队包围长达4个月，乙方封锁了所有水陆交通通道，因此被包围的甲方只能依靠空中交通维持补给，运送4个月的供给依此分别需要2次、3次、3次、4次飞行，每次飞行编队由50架飞机组成，每架飞机都需要3名飞行员，每架飞机每月只能飞行一次，每名飞行员每月也只能飞行一次，每次执行完运输飞行任务后的返回途中有20%的飞机被乙方部队击落，导致机上的飞行员也牺牲或失踪。

在第一个月开始时，甲方拥有110架飞机和330名熟练的飞行员，每个月开始时，甲方可以招聘新飞行员和购买新飞机，新飞机必须经过一个月的检查磨合后才可以投入使用，新飞行员也必须在熟练飞行员的指导下经过一个月的训练才能成为熟练飞行员而投入飞行（作为教练的熟练飞行员本月不能参与飞行任务），每名熟练飞行员作为教练每月指导20名飞行员（包括自己在内）进行训练，每名飞行员在完成本月的飞行任务后必须有一个月的带薪休假，然后返回待命可再次投入飞行，已知各项费用平均单价如下表所示（单位：千元）。

第一个月第二个月第三个月第四个月新飞机价格200 195 190 185闲置的熟练飞行员报酬7 6.9 6.8 6.710 9.9 9.8 9.7教练及飞行员报酬和训练费用执行飞行任务的飞行员报9 8.9 9.8 9.7酬休假期的飞行员报酬 5 4.9 4.8 4.7（1）为甲方安排一个总费用最小的飞行计划。

（2）如果每名熟练飞行员作为教练每月指导不超过20名飞行员（包括自己在内）进行训练，相应的模型和安排将会发生怎样的改变？解：（1）设每月初购买飞机数量为d1，d2，d3,d4架，每月闲置飞机数量为y1,y2,y3,y4架，每月教练与新飞行员总数量为a1,a2,a3,a4人，每月闲置熟练飞行员的数量为b1,b2,b3,b4人。

由于每月执行任务的飞行员和休假期的飞行员的数量是固定的，即这部分的花费是固定的，所以在优化目标中可以不必考虑。

数学建模大作业习题答案

数学建模大作业习题答案数学建模大作业习题答案作为一门应用数学课程，数学建模在现代科学研究和工程技术中具有重要的地位和作用。

通过数学建模，我们可以将实际问题转化为数学模型，从而利用数学方法进行分析和求解。

在数学建模的学习过程中，我们经常会遇到一些习题，下面我将为大家提供一些数学建模大作业题目的答案，希望能对大家的学习有所帮助。

1. 题目：某城市的交通拥堵问题解答：针对这个问题，我们可以采用图论的方法进行建模和求解。

首先，我们将城市的道路网络抽象为一个图，图的节点表示交叉口，边表示道路。

然后，我们可以给每条边赋予一个权重，表示道路的通行能力。

接着，我们可以使用最短路径算法，比如Dijkstra算法，来计算从一个交叉口到另一个交叉口的最短路径，从而找到最优的交通路线。

此外，我们还可以使用最小生成树算法，比如Prim算法，来构建一个最小的道路网络，以减少交通拥堵。

2. 题目：某工厂的生产调度问题解答：对于这个问题，我们可以采用线性规划的方法进行建模和求解。

首先，我们可以将工厂的生产任务抽象为一个线性规划模型，其中目标函数表示最大化生产效益，约束条件表示生产能力、物料供应和市场需求等方面的限制。

然后，我们可以使用线性规划求解器，比如Simplex算法或内点法，来求解这个线性规划模型，得到最优的生产调度方案。

此外，我们还可以引入一些启发式算法，比如遗传算法或模拟退火算法，来寻找更好的解决方案。

3. 题目：某股票的价格预测问题解答：对于这个问题，我们可以采用时间序列分析的方法进行建模和求解。

首先，我们可以将股票的价格序列抽象为一个时间序列模型，比如ARIMA模型。

然后，我们可以使用历史数据来拟合这个时间序列模型，并进行参数估计。

接着，我们可以利用这个时间序列模型来预测未来的股票价格。

此外，我们还可以引入其他的预测方法，比如神经网络或支持向量机，来提高预测的准确性。

通过以上的例子，我们可以看到，在数学建模的过程中，我们需要将实际问题抽象为数学模型，然后利用数学方法进行分析和求解。

数学建模作业

数学建模作业————————————————————————————————作者：————————————————————————————————日期：说明：本电子版题目与教材原题不符者以教材为准，教材上没有的做了会适当加分。

教材上有而本电子版题目没有原题的，请同学们自行录入原题。

所有基本题目解答过程均须不少于姜启源先生《数学模型第三版习题参考解答》之答案长度！第1章数学模型引论1.1 在稳定的椅子问题中，如设椅子的四脚连线呈长方形，结论如何？（稳定的椅子问题见姜启源《数学模型》第6页）（小型题目模版）解：模型分析（黑体五号字）：……宋体五号字模型假设与符号说明（黑体五号字）：……宋体五号字模型建立：……宋体五号字模型求解：……宋体五号字程序源代码（如果需要编程）：……宋体五号字程序运行结果（如果有图形或数据）：……宋体五号字模型讨论：……宋体五号字1.2 在商人们安全过河问题中，若商人和随从各四人，怎样才能安全过河呢？一般地，有n 名商人带n 名随从过河，船每次能渡k 人过河，试讨论商人们能安全过河时，n 与k 应满足什么关系。

（商人们安全过河问题见姜启源《数学模型》第7页）1.3 人、狗、鸡、米均要过河，船需要人划，另外至多还能载一物，而当人不在时，狗要吃鸡，鸡要吃米。

问人、狗、鸡、米怎样过河？1.4 有3对阿拉伯夫妻过河，船至多载两人，条件是根据阿拉伯法典，任一女子不能在其丈夫不在的情况下与其他的男子在一起。

问怎样过河？1.5 如果银行存款年利率为5.5%，问如果要求到2010年本利积累为100000元，那么在1990年应在银行存入多少元？而到2000年的本利积累为多少元？1.6 某城市的Logistic 模型为2610251251N N dt dN ⨯-=，如果不考虑该市的流动人口的影响以及非正常死亡。

设该市1990年人口总数为8000000人，试求该市在未来的人口总数。

当∞→t 时发生什么情况。

数学建模历年竞赛试题

目录前言................................................................................................. 错误!未定义书签。

目录........................................................................................................................... - 0 - 一、什么是数学模型............................................................................................... - 3 -2001年B题……公交车调度......................................................................... - 4 - 2001年C题……基金使用计划..................................................................... - 9 - 2002年A题……车灯线光源的优化设计................................................... - 10 - 2002年B题……彩票中的数学................................................................... - 11 - 2003年A题……SARS的传播.................................................................... - 15 - 2003年B题……露天矿生产的车辆安排................................................... - 26 - 2003年D题……抢渡长江........................................................................... - 29 - 2004年C题……饮酒驾车........................................................................... - 32 - 2004年B题……电力市场的输电阻塞管理............................................... - 34 - 电力市场交易规则：............................................................................. - 35 -输电阻塞管理原则：............................................................................. - 36 -表1各机组出力方案（单位：兆瓦，记作MW） ............................ - 39 -表2各线路的潮流值（各方案与表1相对应，单位：MW） ......... - 41 -表3各机组的段容量（单位：MW） ................................................. - 42 -表4各机组的段价（单位：元/兆瓦小时，记作元/MWh）............. - 42 -表5各机组的爬坡速率（单位：MW/分钟） .................................... - 43 -表6各线路的潮流限值（单位：MW）和相对安全裕度 ................. - 43 -2008年B题……高等教育学费标准探讨................................................... - 43 - 2008年D题……NBA赛程的分析与评价 ................................................. - 45 - 2009年A题……制动器试验台的控制方法分析....................................... - 47 - 2009年B题……眼科病床的合理安排....................................................... - 50 - 【附录】2008-07-13到2008-09-11的病人信息 ................................ - 51 - 2009年D题……会议筹备........................................................................... - 77 - 附表1……10家备选宾馆的有关数据................................................. - 78 -附表2……本届会议的代表回执中有关住房要求的信息（单位：人）- 79 -附表3……以往几届会议代表回执和与会情况.................................. - 80 -附图（其中500等数字是两宾馆间距，单位为米）......................... - 81 -二、为什么要学习数学模型................................................................................. - 83 -1、数学模型无处不在，我们的生活、工作、学习都离不开它............... - 83 -例1买房贷款问题................................................................................. - 83 -例2物体冷却过程的数学模型............................................................. - 84 -2、是学好数学用好数学的必经之路........................................................... - 86 -3、是数学教学改革的重要手段和有效路径............................................... - 88 -4、数学建模竞赛所提唱的团队精神是现代大学生必须具备素质........... - 91 -5、数学建模竞赛鼓励学生用跳跃式的、发散式的形象思维方法，这有利于培养学生的创新意识。

数学建模作业4

数学建模作业-----第五章①非线性最小二乘问题（1）最小二乘方法无约束问题为：minz=∑(a+aa aa a−ηa )2aa=1编写lingo程序：sets:quantity/1..15/:x,y;endsetsdata:x=2,5,7,10,14,19,26,31,34,38,45,52,53,60,65;y=54,50,45,37,35,25,20,16,18,13,8,11,8,4,6;enddatamin=sum(quantity:(a+b*exp(x*c)-y)^2);free(a);free(b);free(c);结果：Local optimal solution found.Objective value: 44.78049Extended solver steps: 5Total solver iterations: 68Variable Value Reduced CostA 2.430177 0.000000B 57.33209 0.000000C -0.4460383E-01 0.000000 （2）最小一乘法无约束问题为：min(a,b) z=∑|a+aa aa a−ηa |aa=1编写lingo程序写出相应的LINGO程序如下：sets:quantity/1..15/:x,y;endsetsdata:x=2,5,7,10,14,19,26,31,34,38,45,52,53,60,65;y=54,50,45,37,35,25,20,16,18,13,8,11,8,4,6; enddatamin=sum(quantity:abs(a+b*exp(x*c)-y));free(a);free(b);free(c);运行结果：Local optimal solution found.Objective value: 20.80640Extended solver steps: 2Total solver iterations: 643Variable Value Reduced CostA 3.398267 0.000000B 57.11461 0.000000C -0.4752126E-01 0.000000（3）最大偏差最小的方法编写程序：sets:quantity/1..15/: x,y;endsetsdata:x = 2,5,7,10,14,19,26,31,34,38,45,52,53,60,65;y = 54,50,45,37,35,25,20,16,18,13,8,11,8,4,6;enddatamin=max(quantity: abs(a+b*exp(c*x)-y));free(a); free(b); free(c);结果：Linearization components added:Constraints: 91Variables: 76Integers: 30Local optimal solution found.Objective value: 2.774408Extended solver steps: 38Total solver iterations: 7654Variable Value Reduced CostA 2.885594 0.000000B 55.86246 0.000000C -0.4441314E-01 0.000000（4）画出散点图和曲线图编写matlab程序如下：>> x=[2 5 7 10 14 19 26 31 34 38 45 52 53 60 65];y=[54 50 45 37 35 25 20 16 18 13 8 11 8 4 6];scatter(x,y,'k*');hold onx=0:0.1:100;y=2.43+exp(-0.0446*x).*57.33; plot(x,y, 'g');xlabel('x');ylabel('y'); hold onx=0:0.1:100;y=3.40+exp(-0.048*x).*57.11; plot(x,y, 'r');xlabel('x');ylabel('y'); hold onx=0:0.1:100;y=2.89+exp(-0.044*x).*55.86; plot(x,y, 'b');xlabel('x');ylabel('y'); 得到图形：分析：三条曲线的结果大致相同，但是略有差异。

数学建模训练习题（含代码程序）

5组于金龙王超焦艳彬快速评卷策略摘要本文研究的是快速评卷问题，在保证准确率和公平公正的原则下，使每位评卷人评阅的试卷总数最小，即满足总的工作量最小。

为解决该问题，在考虑系统误差的前提下，本文建立了多目标优化模型和圆桌评卷模型，利用计算机仿真，建立了以下两种方案，并验证了方案的合理性。

对于方案一，采取了截至分数线淘汰制，每一轮我们将试卷尽可能的平均分成8份，根据该轮试卷的期望值设定一个截至分数线，淘汰分数线以下的所有试卷，剩下的试卷带编号进入下一轮。

当最后的试卷数在2W 附近时，停止进行下一轮仿真，将2W 左右份试卷分给每一位评卷老师进行评阅打分，然后各试卷取平均分进行排名，取前三名为最终优胜者，并记录这三份试卷的编号进行对应。

最后我们通过对上述批卷次数进行统计，一组仿真结果如下：总阅卷次数平均阅卷次数准确率 2182797.1%对于方案二，采用了圆桌评卷模型，将所有试卷尽可能平均分成8份（对应8位带有标号的评卷老师），以第一份试卷为例，首先由第一位老师进行评分，淘汰60%，将余下试卷（含试卷标号）交给右手边的第二位老师进行评分，然后将评分与第一位老师的评分取平均值进行排名，淘汰40%，传给右手边的第三位老师进行评分，按照上一回合的排名制继续淘汰，直至该份试卷只剩下一个则不再淘汰。

将这个试卷依次交给右手边未评过此卷的老师，进行平均打分，最后得出此份中的最优试卷分数及标号。

同样方法，得到剩余7份试卷各自的最优试卷份数及标号，最后对所得8份试卷进行排名，取成绩较高者前三名为优胜试卷，并记录这三份试卷的标号，统计评卷总次数，一组仿真结果如下：总阅卷次数平均阅卷次数准确率 21226.598.3%最后对方案进行分析和改进，对于三种变量：试卷数量、评卷人数和优胜者数量，当其中两种变量不变，调整第三变量时，观察各方案准确率的浮动，得出三者变动规律，寻求出最优评卷策略。

关键字关键字：：计算机仿真圆桌模型系统误差多目标优化1.问题重述在确定像数学建模竞赛这种形式比赛的优胜者时，常常要评阅大量的答卷，比如说，有P=100份答卷。

第二篇模糊数学建模习题解答

第六章习题1.解：设~A 表示“乘客满意”，等车时间为论域}0|{≥=t t U ，由分段函数法，可得高峰期：⎩⎨⎧>≤≤-=30303/1)(~t t t t A μ，非高峰期：⎩⎨⎧>≤≤-=50505/1)(~t t t t A μ设~B 表示“公交公司满意”，载客率为论域}0|{>=v v U ，由分段函数法，可得⎪⎩⎪⎨⎧>≤<≤<=2.1/2.12.15.02.1/5.000)(~v v v v v v B μ 注：答案不惟一，不同的方法可有不同的结果，比如利用时间t 的单调减少函数通过曲线拟合的方法建立~A 的隶属度函数，或者利用模糊数学工具箱中的函数。

2.解：由()()x x B A ~~μμ=可解出曲线交点：5.1*=x ，则有：xeA x Rx c2)21(~1--∈-=⎰，(]x exeB A x x x x 22)22(),5.1[)21(5.1,~~--+∞∈--∞-∈⎰⎰+=(]()xexeB A x x x x 22)21(,5.1)22(5.1,~~--+∞∈--∞-∈⎰⎰+=3.解：（1）654321~6.05.04.019.00D u u u u u u +++++=; （2）~~D C ⊆（3）654321~~4.04.04.001.00C A u u u u u u +++++=（4）{}541~0.5,,A u u u =，{}3~0.5B u =，{}632~0.5,,C u u u ={}1~0.7A u =，∅=~0.7B ，{}32~0.7,C u u =4.解：设=~A “高产”，论域U 为各年的产量首先利用MATLAB 软件计算产量数据的均值与标准差，程序如下： q=[4653.5,4304.7,4497.0,4605.7, 5155.1,5458.1,6193.4,7051.5,8235.6,9953.5];mean(q)std(q)[mean(q)- std(q),mean(q)+std(q)]运行结果：均值=6010.8，标准差=1875.8，估计区间],[s x s x +-==[4135，7886.6]然后利用分段函数法建立隶属度函数⎪⎩⎪⎨⎧>≤<--≤=6.788616.7886413541356.7886413541350)(~x x x x x A μ 说明：估计区间也可以是]2,2[s x s x +-或]3,3[s x s x +-等，但要使得区间包含全部样本点，当然区间越短越好。

2023全国大学生数学建模竞赛模拟题

2023全国大学生数学建模竞赛模拟题第一部分：问题描述在2023年全国大学生数学建模竞赛中，我们将考虑以下问题：问题一：某大学计划对校园内的停车管理进行优化。

假设校园内有N个停车位（N为正整数），每个停车位只能停放一辆车。

现在需要设计一个停车系统，使得所有车辆能够尽可能高效地停放在停车位上。

请你们给出一个数学模型，以及相应的优化策略，以满足停车位利用效率最大化的要求。

问题二：某电商公司为了提高货物的配送效率，需要选址一些配送中心，以覆盖尽可能多的用户。

假设已知用户的分布情况和需求量，在这些信息的基础上，请你们设计一个数学模型，并给出选址策略，以最大化用户的满意度，同时尽量减少配送的时间和成本。

第二部分：问题分析与数学模型建立问题一：停车管理优化我们首先定义问题的目标函数，即停车位利用效率的优化目标。

假设停车场内每个停车位的编号为i（i=1,2,...,N），对于每个停车位，我们引入二进制变量x_i，表示该停车位是否被使用，其中x_i=1表示被占用，x_i=0表示空闲。

接着，我们需要确定约束条件。

显然，每个停车位只能被一辆车使用，即∑x_i ≤ 1 （i=1,2,...,N）其中，∑表示求和。

为了使停车位利用效率最大化，我们可以引入一个系数p_i，表示第i个停车位的利用效率，取值范围为[0,1]。

利用效率越高，则p_i越接近1，反之越接近0。

我们可以根据停车位距离出入口的远近、停车位所在区域的拥挤程度等因素来确定p_i的取值。

然后，我们可以构建目标函数：Maximize ∑p_i*x_i （i=1,2,...,N）最后，我们将目标函数和约束条件整合，形成一个数学模型。

问题二：配送中心选址对于问题二，我们可以将用户的需求量作为权重，即需求量越高的用户对配送中心的选择影响越大。

假设有M个可能的配送中心位置（M为正整数），每个位置编号为j（j=1,2,...,M），我们引入二进制变量y_j，表示第j个位置是否选址为配送中心，其中y_j=1表示选址，y_j=0表示不选址。

数学建模基础练习一及参考答案

数学建模基础练习一及参考答案练习1 matlab练习一、矩阵及数组操作1．利用基本矩阵产生3×3和15×8的单位矩阵、全1矩阵、全0矩阵、均匀分布随机矩阵（[-1，1]之间）、正态分布矩阵（均值为1，方差为4）,然后将正态分布矩阵中大于1的元素变为1，将小于1的元素变为0。

2．利用fix及rand函数生成[0，10]上的均匀分布的10×10的整数随机矩阵a，然后统计a中大于等于5的元素个数。

3．在给定的矩阵中删除含有整行内容全为0的行，删除整列内容全为0的列。

4．随机生成10阶的矩阵，要求元素值介于0~1000之间，并统计元素中奇数的个数、素数的个数。

二、绘图5．在同一图形窗口画出下列两条曲线图像，要求改变线型和标记y1=2x+5；y2=x^2-3x+1，并且用legend标注。

6．画出下列函数的曲面及等高线z=sinxcosyexp(-sqrt(x^2+y^2))． 7．在同一个图形中绘制一行三列的子图，分别画出向量x=[1 5 8 10 12 5 3]的三维饼图、柱状图、条形图。

三、程序设计8．编写程序计算（x在[-8，8]，间隔0.5）先新建的，在那上输好，保存，在命令窗口代数；9．用两种方法求数列前15项的和。

10．编写程序产生20个两位随机整数，输出其中小于平均数的偶数。

11．试找出100以内的所有素数。

12．当时，四、数据处理与拟合初步1随机产生由10个两位随机数的行向量A,将A中元素按降序排列为B,再将B重排为A。

14．通过测量得到一组数据t 1 2 3 4 5 6 7 8 9 10 y 842 362 754 368 169 038 034 016 012 005 分别采用y=c1+c2e^(-t)和y=d1+d2te^(-t)进行拟合，并画出散点及两条拟合曲线对比拟合效果。

15．计算下列定积分16．（1）微分方程组当t=0时，x1(0)=1，x2(0)=-0.5，求微分方程t在[0，25]上的解，并画出相空间轨道图像。

数学建模matlab例题参考及练习讲课稿

数学建模m a t l a b例题参考及练习数学实验与数学建模实验报告学院：专业班级：姓名：学号：完成时间：年月日承诺书本人承诺所呈交的数学实验与数学建模作业都是本人通过学习自行进行编程独立完成，所有结果都通过上机验证，无转载或抄袭他人，也未经他人转载或抄袭。

若承诺不实，本人愿意承担一切责任。

承诺人：年月日数学实验学习体会（每个人必须要写字数1200字以上，占总成绩的20%）练习1 一元函数的图形1．画出x y arcsin =的图象.2．画出x y sec =在],0[π之间的图象.3．在同一坐标系中画出x y =，2x y =，3x y =，3x y =，x y =的图象.4．画出3232)1()1()(x x x f ++-=的图象，并根据图象特点指出函数)(x f 的奇偶性.5．画出)2ln(1++=x y 及其反函数的图象.6．画出321+=x y 及其反函数的图象.练习2 函数极限1．计算下列函数的极限. (1)x xx 4cos 12sin 1lim 4-+π→.程序：sym x ;f=(1+sin(2*x))/(1-cos(4*x));limit(f,x,pi/4)运行结果：lx21ans =1(2).程序：sym x ;f=(1+cos(x))^(3*sec(x));limit(f,x,pi/2)运行结果：lx22ans =exp(3)(3)22)2(sin ln lim x xx -ππ→.程序：sym x ;f=log(sin(x))/(pi-2*x)^2;limit(f,x,pi/2)运行结果：lx23ans =x x x sec 3 2 ) cos 1 ( lim + π →-1/8(4)212lim xxex→.程序：sym x;f=x^2*exp(1/x);limit(f,x,0)limit(f,x,0,'right')limit(f,x,0,'left')运行结果：lx24ans =NaNans =Infans =%左极限为零，存在，右极限为无穷大，在x趋近于零时函数没有极限(5))215(lim122xx xx+-∞→.程序：sym x;f=5*x^2/(1-x^2)+2^(1/x);limit(f,x,inf)运行结果：>> lx25ans =(6)x x x x x -+-→32112lim .程序：sym x ;f=(x^2-2*x+1)/(x^3-x);limit(f,x,1)运行结果：>> lx26ans =0 (7)x x x 11lim 20-+→.程序：sym x ;f=(sqrt(1+x^2)-1)/x;limit(f,x,0)运行结果：>> lx27ans =0 (8))3sin(cos 21lim 3π--π→x xx . 程序：sym x ;f=(1-2*cos(x))/sin(x-pi/3);limit(f,x,pi/3)运行结果：>> lx28ans =3^(1/2) (9)tgxx x )1(lim 0+→.程序：sym x ;f=(1/x)^tan(x);limit(f,x,0,'right')运行结果：>> lx29ans =(10)xx arctgx )2(lim π+∞→.程序：sym x ;f=(2/pi*atan(x))^x;limit(f,x,inf,'left')运行结果：>> lx210ans =Inf2．解方程012=-⋅x x .程序：sym x ;X=solve(x*2^x-1)运行结果：>> lx202X =lambertw(0, log(2))/log(2)%方程有两个解3．解方程1sin 3+=x x .程序：sym x ;X=solve(3*sin(x)+1-x)运行结果：>> lx203X =-0.538470451711254993610615326557454．解方程03=++q px x .(p 、q 为实数)程序：X=solve('x^3+p*x+q=0','x')运行结果：X =((p^3/27 + q^2/4)^(1/2) - q/2)^(1/3) - p/(3*((p^3/27 +q^2/4)^(1/2) - q/2)^(1/3))p/(6*((p^3/27 + q^2/4)^(1/2) - q/2)^(1/3)) - ((p^3/27 + q^2/4)^(1/2) - q/2)^(1/3)/2 - (3^(1/2)*i*(p/(3*((p^3/27 + q^2/4)^(1/2) - q/2)^(1/3)) + ((p^3/27 + q^2/4)^(1/2) - q/2)^(1/3)))/2p/(6*((p^3/27 + q^2/4)^(1/2) - q/2)^(1/3)) - ((p^3/27 + q^2/4)^(1/2) - q/2)^(1/3)/2 + (3^(1/2)*i*(p/(3*((p^3/27 + q^2/4)^(1/2) - q/2)^(1/3)) + ((p^3/27 + q^2/4)^(1/2) - q/2)^(1/3)))/2练习 3 导数及偏导数计算1.求下列函数的导数. (1))11)(1(-+=x x y程序：sym x ;f=(sqrt(x)+1)*(1/sqrt(x)-1);diff(f)运行结果：>> lx31ans =(1/x^(1/2) - 1)/(2*x^(1/2)) - (x^(1/2) +1)/(2*x^(3/2))(2)x x x y ln sin =程序：sym x ;f=x*sin(x)*log(x);diff(f)运行结果：>> lx32ans =sin(x) + log(x)*sin(x) + x*cos(x)*log(x)2.求下列参数方程所确定的函数的导数.(1)⎩⎨⎧==t y t x 44程序：sym t ;f1=t^4;f2=4*t;diff(f2)/diff(f1)运行结果：>> lx3211/t^3(2)⎩⎨⎧-=+=arctgt t y t x )1ln(2程序：sym t ;f1=log(1+t^2);f2=t-atan(t);diff(f2)/diff(f1)运行结果：>> lx322ans =-((t^2 + 1)*(1/(t^2 + 1) - 1))/(2*t)3.求下列隐函数的导数. (1)22ln y x x y arctg +=程序：syms x y ;f=atan(y/x)-log(sqrt(x^2+y^2));yx=-diff(f,x)/diff(f,y)运行结果；>> lx331yx =(x/(x^2 + y^2) + y/(x^2*(y^2/x^2 +1)))/(1/(x*(y^2/x^2 + 1)) - y/(x^2 + y^2))(2)x y y x =程序：syms x y ;f=x^y-y^xyx=-diff(f,x)/diff(f,y)运行结果：>> lx332f =x^y - y^x(x^(y - 1)*y - y^x*log(y))/(x*y^(x - 1) - x^y*log(x))4.设x e y x cos =，求)4(y .程序：sym x ;f=exp(x)*sin(x);diff(f,x,4)运行结果：>> lx34ans =(-4)*exp(x)*sin(x)5.验证x e y x sin =满足关系式:022=+'-''y y y程序：sym x ; f=exp(x)*sin(x);y2=diff(f,x,2);y1=diff(f,x,1);y=f;y2-y1*2+2*y=='0'运行结果：>> lx35ans =1%运行结果为1表示y2-y1*2+2*y=='0'成立6.设)ln(y x x u +=，求22x u ∂∂，22y u ∂∂，y x u ∂∂∂2. 程序：syms x y ;f=x*log(x+y);uxx=diff(f,x,2)uyy=diff(f,y,2)f1=diff(f,x);uxy=diff(f1,y)运行结果：>> lx36uxx =2/(x + y) - x/(x + y)^2uyy =-x/(x + y)^2uxy =1/(x + y) - x/(x + y)^27.求下列多元隐函数的偏导数y zx z ∂∂∂∂,.(1)1cos cos cos 222=++z y x程序：syms x y z ;f=(cos(x))^2+(cos(y))^2+(cos(z))^2-1;zx=-diff(f,x)/diff(f,z)zy=-diff(f,y)/diff(f,z)运行结果：>> lx371zx =-(cos(x)*sin(x))/(cos(z)*sin(z))zy =-(cos(y)*sin(y))/(cos(z)*sin(z))(2)xyz e z =程序：syms x y z ;f=exp(z)-x*y*zzx=-diff(f,x)/diff(f,z)zy=-diff(f,y)/diff(f,z)运行结果：>> lx372f =exp(z) - x*y*zzx =(y*z)/(exp(z) - x*y)zy =(x*z)/(exp(z) - x*y)8.证明函数22)()(ln b y a x u -+-=(b a ,为常数)满足拉普拉斯方程: 02222=∂∂+∂∂y u x u （提示：对结果用simplify 化简）练习4 积分计算1.计算下列不定积分. (1)⎰+dx x x 12 (2)⎰+x xdx 2sin 12sin2.计算下列定积分.(1)⎰e xdx x 1ln (2)⎰ππ342sin dxx x3.求⎰+t dx x x x 12)ln (ln 1并用diff 对结果求导.4.求摆线)cos 1(),sin (t a y t t a x -=-=的一拱(π≤≤20t )与x 轴所围成的图形的面积.5.计算二重积分(1)⎰⎰≤++122)(y x dxdy y x (2)⎰⎰≤++x y x dxdy y x 22)(22 6.计算⎰+L ds y x 22 L 为圆周)0(22>=+a ax y x7.计算⎰++-L dy y x dx y x )()(2222，其中L 为抛物线2x y =上从点(0，0)到点(2，4)的一段弧.练习5 matlab 自定义函数与导数应用1.建立函数x x a a x f 3sin 31sin ),(+=，当a 为何值时，该函数在3π=x 处取得极值，它是极大值还是极小值，并求此极值.2.确定下列函数的单调区间.（1）7186223---=x x x y （2）)0(82>+=x x x y3.求下列函数的最大值、最小值.（1）2332xx y -=41≤≤-x （2）312824≤≤-+-=x x x y练习6 matab 矩阵运算与数组运算1．计算（1）⎪⎪⎪⎭⎫ ⎝⎛--521111204321+⎪⎪⎪⎭⎫ ⎝⎛--232002101041221 （2）⎪⎪⎭⎫ ⎝⎛-01301213⎪⎪⎪⎪⎪⎭⎫ ⎝⎛030101020501⎪⎪⎪⎭⎫ ⎝⎛-205101 (3)52422⎪⎪⎭⎫ ⎝⎛- 2．设⎪⎪⎪⎭⎫ ⎝⎛-=243121013A ,⎪⎪⎪⎭⎫ ⎝⎛-=112111201B ,求满足关系B X A =-23的X ．练习7 矩阵与线性方程组１．求下列矩阵的秩．(1)⎪⎪⎪⎭⎫ ⎝⎛-321110021 (2)⎪⎪⎪⎪⎪⎭⎫ ⎝⎛4820322513454947513253947543173125 ２．求下列矩阵的行列式，如可逆，试用不同的方法求其逆矩阵．（1）⎪⎪⎪⎭⎫ ⎝⎛--285421122 （2）⎪⎪⎪⎪⎪⎭⎫ ⎝⎛---6201111121324321 ３．设X ⎪⎪⎪⎭⎫ ⎝⎛-111012111=⎪⎪⎪⎭⎫ ⎝⎛521234311求X ．４．解下列线性方程组．(1)⎪⎪⎩⎪⎪⎨⎧=--+=+-+=+-+=+-+6223312433862344224221432143214321x x x x x x x x x x x x x x x x (2)⎪⎪⎩⎪⎪⎨⎧-=+--=+--=-+-212201432143214321x x x x x x x x x x x x练习8 常微分方程与级数求1-6题微分方程的通解1．1222+='y y y x 2．x y x y dx dy -+= 3．x x x y y +='cos 4．1)2sin cos (='+y y y x 5．x e y y y x 2cos 3=-'+'' 6．x x y y sin 14++=+''求7、8题初值问题的解7．⎪⎩⎪⎨⎧==-++-+=10)2(212222x y dx dy x xy y y xy x8．⎪⎪⎩⎪⎪⎨⎧===++==0000222,02V dt dxx x x a dt dx n dt x d t t9．给出函数x x e x f x x cos 2sin )(+=在点0=x 的7阶taylor 展开式以及在x=1处的 5阶taylor 展开式．10．判别下列级数的敛散性，若收敛求其和． (1)Λ++++7151311 (2)∑∞=+112n n n tg π11．求幂级数∑∞=--22)1(n nn n n x 的和函数． 12．求函数项级数∑∞=-1)2sin )1(n n n n x π的和函数．。

数学建模最优化模型例题

数学建模最优化模型例题好，咱们今天来聊聊数学建模和最优化模型这块儿。

数学建模，这名字听起来就挺高大上的，实际上，咱们日常生活中处处都是它的身影。

想象一下，早上起床，看到窗外阳光明媚，心里琢磨着今天去不去公园，顺便锻炼锻炼。

于是，你心里开始盘算，公园离家有多远，走路要多久，还是骑个单车比较快？这就是在用数学建模，算一算，看看哪个更划算。

再说说最优化模型，这就像是在挑选午饭一样。

你有一大堆选择，米饭、面条、快餐还是外卖，真是眼花缭乱。

你心里想，要是不吃太油腻的，又想吃得饱，还得好吃。

于是开始分析：今天外卖不如自己做，自己做的话，买啥材料比较好，怎么搭配更营养呢？这时候，你的脑子就像一个小计算机，开始进行各种选择。

想想，如果能把所有的选择变成一个数学问题，肯定能算出最优解，嘿，生活简直就像在解题一样，乐趣多多。

再说说商场里打折的那种，真是让人心痒痒的。

假如你打算买新鞋，满心期待。

可是一进商场，各种颜色、各种款式扑面而来，心里顿时就犯了选择困难症。

想要买的那双鞋打折了，可是另外一双颜色也不错，怎么办呢？这时候，最优化模型就可以帮你了。

想一想，你最看重什么，舒适、样式还是价格？用数学的眼光来审视，看看哪双鞋的性价比最高，没准儿就能找到那个最适合自己的了。

有些小伙伴可能会问了，数学建模到底有什么用呢？你知道吗，很多企业在决策的时候都离不开这些模型。

就拿快递公司来说，他们每天都要处理成千上万的包裹，怎么能保证包裹及时送到呢？他们需要用到最优化模型来安排路线，减少运输成本。

想象一下，如果没有这些模型，快递员可能跑了一大圈，最后才发现原来只需要直走就到了。

那可真是得不偿失，没准儿包裹还会晚到，这可就麻烦了。

数学建模的魅力就在于它能把复杂的问题简单化。

我们生活中遇到的各种难题，最终都可以转化为一个个数学问题。

你说这是不是挺神奇的？比如你要规划一次旅行，想去多少个地方，怎么安排最合适，住哪儿能便宜又舒服，这些全都可以用建模来解决。

数学建模上机练习习题及答案教学内容

数学建模上机练习习题及答案教学内容练习1 基础练习⼀、矩阵及数组操作：1．利⽤基本矩阵产⽣3×3和15×8的单位矩阵、全1矩阵、全0矩阵、均匀分布随机矩阵（[-1，1]之间）、正态分布矩阵（均值为1，⽅差为4）。

A=eye(3) B=eye(15,8) C=ones(3) D=ones(15,8) E=zeros(3) F=zeros(15,8) G=(-1+(1-(-1))*rand(3)) H=1+sqrt(4)*randn(5) 2．利⽤fix及rand函数⽣成[0，10]上的均匀分布的10×10的整数随机矩阵a，然后统计a中⼤于等于5的元素个数a=fix(0+(10-0)*rand(10));K=find(a>=5);Num=length(K)或者num=sum(sum(a>=5))num =533．在给定的矩阵中删除含有整⾏内容全为0的⾏，删除整列内容全为0的列。

如已给定矩阵A在给定的矩阵中删除含有整⾏内容全为0的⾏在命令窗⼝中输⼊A(find(sum(abs(A'))==0),:)=[];删除整列内容全为0的列。

A(：，find(sum(abs(A'))==0))=[];⼆、绘图：4．在同⼀图形窗⼝画出下列两条曲线图像： y1=2x+5； y2=x^2-3x+1，并且⽤legend 标注 x=0:0.01:10; y1=2*x+5; y2=x.^2-3*x+1; plot(x,y1,x,y2,'r') legend('y1', 'y2')12345678910-10010203040506070805．画出下列函数的曲⾯及等⾼线： z=x^2+y^2+sin(xy)．在命令窗⼝输⼊： [x,y]=meshgrid(0:0.25:4*pi);z=x.^2+y.^2+sin(x.*y)；contour3(x,y,z); meshc(x,y,z)5101551015100200300400三、程序设计：6．编写程序计算（x在[-3，3]，间隔0.01）建⽴M⽂件d.mx=input('请输⼊x的值:');if x>=-3&x<-1y=(-x.^2-4*x-3)/2;elseif x>=-1&x<1y=-x.^2+1;elseif x>=1&x<=3y=(-x.^2+4*x-3)/2;elsey='error'endy在命令窗⼝输⼊x 的值：7．有⼀列分数序列：求前15项的和。

常见数学建模练习题目及解答

=
3
−
x2
+
x3
+
d
− 4
−
d
+ 4
=
0
x1
+
d
− 5
−
d
+ 5
=
2.4
x2
+
d
− 6
−
d
+ 6
=3
model: sets: level/1..3/:p,z,goal; variable/1..3/:x; s_con_num/1..6/:g,dplus,dminus; s_con(s_con_num,variable):c; obj(level,s_con_num)/1 1,2 2,2 3,2 4,3 5,3 6/:wplus,wminus; endsets data: ctr=?; goal=? ? 0; g=9000 3 3 0 2.4 3; c=500 500 1000 100 -1 1 0 0 -1 1 100 0 1 0; wplus=1 1 1 1 0 0; wminus=0 0 0 0 1 1; enddata min=@sum(level:p*z);
⎨ ⎪⎩
y
2
'
=
−
y2 x
− (1 −
1 4x
2
)
y1
y1=dsolve('x^2*D2y+x*Dy+(x^2-1/4)*y','y(pi/2)=2,Dy(pi/2)=-2/pi',' x'); subplot(1,3,1) ezplot(y1,[pi/2,pi]) f=@(x,y) [y(2);-y(2)/x-(1-0.25./x.^2)*y(1)]; [x,y2]=ode45(f,[pi/2,pi],[2,-2/pi]) subplot(1,3,2) plot(x,y2(:,1)) y3=subs(y1,'x',x); subplot(1,3,3) plot(x,y3)

数学建模作业43508

08年全国大学生数学建模A题 文档

数学建模案例分析-- 模糊数学方法建模5习题八

数学建模答案(完整版)

数学建模题目及答案-数学建模100题

数学建模作业43508

数学建模大作业习题答案

数学建模作业

数学建模历年竞赛试题

数学建模作业4

数学建模训练习题（含代码程序）

第二篇模糊数学建模习题解答

2023全国大学生数学建模竞赛模拟题

数学建模基础练习一及参考答案

数学建模matlab例题参考及练习讲课稿

数学建模最优化模型例题

数学建模上机练习习题及答案教学内容

常见数学建模练习题目及解答

08年全国大学生数学建模A题文档