数学建模章绍辉版作业 SANY标准化小组 #QS8QHH-HHGX8Q8-GNHHJ8-HHMHGN#
第四章作业
第二题：
针对严重的交通情况，国家质量监督检验检疫局发布的国家标准，车辆驾驶人员血液中的酒精含量大于或等于20mg/100ml,小于80mg/100ml 为饮酒驾车，血液中的酒精含量大于或等于80mg/100ml 的为醉酒驾车。
下面分别考虑大李在很短时间内和较长时间内（如2个小时）喝了三瓶啤酒，多长时间内驾车就会违反新的国家标准。
1、 问题假设
大李在短时间内喝下三瓶啤酒后，酒精先从吸收室（肠胃）吸收进中心室（血液和体液），然后从中心室向体外排除，忽略喝酒的时间，根据生理学知识，假设
（1） 吸收室在初始时刻t=0时，酒精量立即为032
D ;在任意时刻，酒精从吸收室吸收进中心室的速率（吸收室在单位时间内酒精含量的减少量）与吸收室的酒精含量成正比，比例系数为1k ；
（2） 中心室的容积V 保持不变；在初始时刻t=0时，中心室的酒精含量为0；在任意时
刻，酒精从中心室向体外排除的速率（中心室在单位时间内酒精含量的减少量）与中心室的酒精含量成正比，比例系数为2k ；
（3） 在大李适度饮酒没有酒精中毒的前提下，假设1k 和2k 都是常量，与饮酒量无关。
2、 符号说明
酒精量是指纯酒精的质量，单位是毫克；
酒精含量是指纯酒精的浓度，单位是毫克/百毫升；
~t 时刻（小时）；
()~x t 在时刻t 吸收室（肠胃）内的酒精量（毫克）；
0~D 两瓶酒的酒精量（毫克）；
(t)~c 在时刻t 吸收室（血液和体液）的酒精含量（毫克/百毫升）；
2()~c t 在时刻t 中心室（血液和体液）的酒精含量（毫克/百毫升）；
~V 中心室的容积（百毫升）；
1~k 酒精从吸收室吸收进中心室的速率系数（假设其为常数）；
2~k 酒精从中心室向体外排除的速率系数（假设其为常数）；
3~k 在短时间喝下三瓶酒的假设下是指短时间喝下的三瓶酒的酒精总量除以中心室体积，即03/2D V ；而在较长时间内（2小时内）喝下三瓶酒的假设下就特指03/4D V .
3、 模型建立和求解
（1） 酒是在很短时间内喝的：
记喝酒时刻为0t =（小时），设(0)0c =，可用()2113212
()k t k t k k c t e e k k --=--来计算血液中的酒精含量，此时12k k 、为假设中所示的常数，而033155.792D k V ??== ???
. 下面用MATLAB 程序画图展示血液中酒精含量随时间变化并且利用fzero 函数和fminbnd 函数来得到饮酒驾车醉酒驾车对应的时间段，以及血液中酒精含量最高的时刻。
MATLAB 程序如下：
k1=;k2=;k3=;
c=@(t)(k1.*k3)./(k1-k2).*(exp(-k2.*t)-exp(-k1.*t));
f=@(t)c(t)-20;
g=@(t)c(t)-80;
h=@(t)-c(t);
t1(1)=fzero(f,1);t1(2)=fzero(f,12),
t2(1)=fzero(g,1);t2(2)=fzero(g,12)
[t3,c3]=fminbnd(h,0,24)
fplot(c,[0,20],'k')
hold on
plot([0,20],[20,20],'k',[0,20],[80,80],'k')
hold off
xlabel('时刻t （小时），从开始喝酒算起')
ylabel('血液中的酒精含量（mg/100ml ）')
title('短时间喝下三瓶酒时，血液中酒精含量随时间的变化过程')
gtext(',20)')
gtext(',20)')
gtext(',80)')
gtext(',80)')
gtext(',')
运行结果如下：
t1 =
t2 =
t3 =
c3 =
所绘图形如下：
结果分析：
所以，当[0.06891,0.38052)(4.1125,11.589]t ∈?时，20()80c t ≤<，属饮酒驾车。当
[0.38052,4.1125]t ∈时，属醉酒驾驶；当 1.307t =时，血液中的酒精含量最高为毫克/百毫升。
（2） 酒是在2小时内喝的：
可假设三瓶啤酒是在2小时内匀速喝的. 同样记喝酒时刻为0t =（小时），设(0)0c =，则吸收室的酒精量1()x t 满足分段的初值问题
解得
于是中心室内的酒精含量2()c t 满足分段的初值问题
解得
其中
0334D k V
=，3412k k k k =-，1452k k k k =，362= k k k ，1273(1)k k k e =- 12228456k k k k e k e k --=-+，7912
k k k k =-，22122()1089k k k k k e k e -=+ 因为1 2.0079k =，20.1855k =以及0/103.86D V =，所以
377.896k =，442.743k =，5462.66k =，6419.92k =
74243.1k =，8101.43k =，92328.3k =，10207.82k =
下面用MATLAB 程序画图展示血液中酒精含量随时间变化并且利用fzero 函数和fminbnd 函数来得到饮酒驾车醉酒驾车对应的时间段，以及血液中酒精含量最高的时刻。
MATLAB 程序如下：
k1=;k2=;k3=;
k4=;k5=;k6=;k9=;k10=;
c1=@(t)(k4.* exp(-k1.*t)-k5.*exp(-k2.*t)+k6).*(t>=0&t<=2)+...
( k10.* exp(-k2.*t)-k9.*exp(-k1.*t)).*(t>2);
f1=@(t)c1(t)-20;
g1=@(t)c1(t)-80;
h1=@(t)-c1(t);
t1(1)=fzero(f1,1);t1(2)=fzero(f1,12),
t2(1)=fzero(g1,1);t2(2)=fzero(g1,12),
[t3,c3]=fminbnd(h1,0,20)
fplot(c1,[0,20],'k')
hold on
plot([0,20],[20,20],'k',[0,20],[80,80],'k')
hold off
xlabel('时刻t （小时），从开始喝酒算起')
ylabel('血液中的酒精含量（mg/100ml ）')
title('短时间喝下三瓶酒时，血液中酒精含量随时间的变化过程')
gtext(',20)')
gtext(',20)')
gtext(',80)')
gtext(',80)')
gtext(',')
运行结果如下：
t1 =
t2 =
t3 =
c3 =
所绘图形如下：
结果分析：
所以，当[0.62321,1.6366)(5.1412,12.62]
≤<，属饮酒驾车。当
c t
t∈?时，20()80
t=时，血液中的酒精含量最高，为毫克/百毫升.
t∈时，属醉酒驾驶；当 2.6328
[1.6366,5.1412]
下面用图形比较两种不同假设下血液中酒精含量的变化过程。
MATLAB程序如下：
k1=;k2=;k3=;
k4=;k5=;k6=;k9=;k10=;
c=@(t)(k1.*k3)./(k1-k2).*(exp(-k2.*t)-exp(-k1.*t));
c1=@(t)(k4.* exp(-k1.*t)-k5.*exp(-k2.*t)+k6).*(t>=0&t<=2)+...
( k10.* exp(-k2.*t)-k9.*exp(-k1.*t)).*(t>2);
plot(0::20,c(0::20),'--k',0::20,c1(0::20),'k',2,c1(2),'.k')
xlabel('时刻t（小时），从开始喝酒算起')
ylabel('血液中的酒精含量（mg/100ml）')
title('短时间喝下三瓶酒时，血液中酒精含量随时间的变化过程')
legend('很短时间内喝三瓶啤酒','两小时内匀速喝下三瓶啤酒','函数的分段点')
所绘图形如下：
第四题：
研究将鹿群放入草场后，草和鹿两个种群的相互作用，草的生长服从Logistic规律，年固有增长率，最大密度为3000个密度单位，在草最茂盛时，每只鹿每年吃掉个密度单位的草，若没有草，鹿群的年死亡率高达,而在草最茂盛的时候草对鹿的死亡的补偿率为。
1、建立差分方程组模型，比较将100只鹿放入密度为1000和密度为3000的两种草场的情况下，草和鹿两个种群的数量演变过程。
（1）符号说明：
x~
第k年草场的密度单位
k
~
y第k年草场上鹿的数量
k
r~草场上草的年固有增长率
a~由于捕食导致的草的密度单位减少的速度大小
d~如果没有草，鹿群的年死亡率
b~草对鹿群死亡的补偿率
N~草的最大密度单位
（2）模型的建立与求解：
基于以上假设，由于草的生长服从Logistic模型，建立差分方程组模型如下所示：
令1k k x x x +==，1k k y y y +==，与上述方程组联立得到平衡点为0(0,0)P 、1(,0)P N 、2()(
,)dN rN b d P b ab
- 以下用MATLAB 实现差分方程组模型。
MATLAB 程序如下：
n=20;r=;a=;b=;d=;
N=3000;t=1:n+1;
x1(1)=1000;y1(1)=100;
for k=1:n
??? x1(k+1)=x1(k)+r*x1(k)*(1-x1(k)/N)-a*x1(k)*y1(k)/N
??? y1(k+1)=y1(k)-d*y1(k)+b*x1(k)*y1(k)/N
end
subplot(2,2,1),plot(t,x1,'k^',t,y1,'kv'),axis([-1,21,0,3000]), xlabel('第k 年'),ylabel('数量'),gtext('草密度'),gtext('鹿数量'), title('草和鹿随时间的演变');
subplot(2,2,2),plot(x1,y1,'ko'),axis([1000,3000,0,1000]), xlabel('草密度'),ylabel('鹿数量'),title('相平面图');
x2(1)=3000;y2(1)=100;
for k=1:n
??? x2(k+1)=x2(k)+r*x2(k)*(1-x2(k)/N)-a*x2(k)*y2(k)/N
??? y2(k+1)=y2(k)-d*y2(k)+b*x2(k)*y2(k)/N
end
subplot(2,2,3),plot(t,x2,'k^',t,y2,'kv'),axis([-1,21,0,3000]), xlabel('第k 年'),ylabel('数量'),gtext('草密度'),gtext('鹿数量'), title('草和鹿随时间的演变');
subplot(2,2,4),plot(x2,y2,'ko'),axis([1000,3000,0,1000]), xlabel('草密度'),ylabel('鹿数量'),title('相平面图');
运行结果如下：
x1 =
Columns 1 through 8
1000 1480
Columns 9 through 16
1550
Columns 17 through 21
y1 =
Columns 1 through 8
100 60
Columns 9 through 16
Columns 17 through 21
x2 =
Columns 1 through 8
3000 2840 2570 2149
Columns 9 through 16
1698 1907
Columns 17 through 21
y2 =
Columns 1 through 8
100 160
Columns 9 through 16
Columns 17 through 21
所绘图像如下： 由图像中可以看出2()(,)dN rN b d P b ab
-为渐进稳定的平衡点。 2、建立常微分方程组模型，重做以上问题。
记草和鹿在第t 年的数量分别记为()x x t =和()y y t =，其余的符号假设与（1）相同，建立常微分方程组模型为： 令(1)0x axy rx N N --=且0k bxy dy N
-+=，得到常微分方程组的临界点为0(0,0)P 、1(,0)P N 、2()(,)dN rN b d P b ab
-，以下运用数值解，方向场和特征线等技巧，用matlab 绘制常微分方程组的解曲线图，并加以分析。
为了消去方程中的参数t ，将两式相除，得到：()dy Ndy bxy dx rx N x axy -+=-- 当将100头鹿放入密度为1000与密度为3000的草场中，初始值分别为(0)1000x =，(0)100y =与(0)3000x =，(0)100y =.
利用matlab 实现上述过程的源程序与运行结果如下：
函数m 文件：
function dx=fun(t,x)
dx=zeros(2,1);
dx(1)=x(1)**(1-x(1)/3000)*x(1)*x(2)/3000;
dx(2)=*x(2)+*x(1)*x(2)/3000;
主程序：
[t1,x1]=ode45('fun',0:20,[1000 100]);
[t1,x1]
subplot(2,2,1)
plot(t1,x1(:,1),'.-k',t1,x1(:,2),'.-k'),grid on
axis([0 20 0 3000])
title('草和鹿随时间的演变（x_0=1000）')
gtext('草密度'),gtext('鹿数量')
xlabel('第k 年'),ylabel('数量')
subplot(2,2,2)
plot(x1(:,1),x1(:,2),'.-k'),grid on
axis([1000 3000 0 800])
title('相平面图（x_0=1000）')
xlabel('草密度'),ylabel('鹿数量')
[t2,x2]=ode45('fun',0:20,[3000 100]);
[t2,x2]
subplot(2,2,3)
plot(t2,x2(:,1),'.-k',t2,x2(:,2),'.-k'),grid on
axis([0 20 0 3000])
title('草和鹿随时间的演变（x_0=3000）')
gtext('草密度'),gtext('鹿数量')
xlabel('第k 年'),ylabel('数量')
subplot(2,2,4)
plot(x2(:,1),x2(:,2),'.-k'),grid on
axis([1000 3000 0 800])
title('相平面图（x_0=3000）')
xlabel('草密度'),ylabel('鹿数量')
所绘图像如下： 通过图像发现该常微分方程组的平衡点2()(,)dN rN b d P b ab
-是渐进稳定的，而0(0,0)P 、1(,0)P N 不是渐进稳定的平衡点。
3差分方程相当于取步长1h =时用欧拉方法计算微分方程组的数值解。
数学建模部分课后习题解答 中国地质大学 能源学院 华文静 1.在稳定的椅子问题中，如设椅子的四脚连线呈长方形，结论如何？ 解： 模型假设 （1） 椅子四条腿一样长，椅脚与地面接触处视为一点，四脚的连线呈长方形 （2） 地面高度是连续变化的，沿任何方向都不会出现间断（没有像台阶那样的情况）， 即从数学角度来看，地面是连续曲面。这个假设相当于给出了椅子能放稳的必要条件 （3） 椅子在任何位置至少有三只脚同时着地。为了保证这一点，要求对于椅脚的间 距和椅腿的长度而言，地面是相对平坦的。因为在地面上椅脚间距和椅腿长度的尺寸大小相当的范围内，如果出现深沟或凸峰（即使是连续变化的），此时三只脚是无法同时着地的。 模型建立 在上述假设下，解决问题的关键在于选择合适的变量，把椅子四只脚同时着地表示出来。首先，引入合适的变量来表示椅子位置的挪动。生活经验告诉我们，要把椅子通过挪动放稳，通常有拖动或转动椅子两种办法，也就是数学上所说的平移与旋转变换。然而，平移椅子后问题的条件没有发生本质变化，所以用平移的办法是不能解决问题的。于是可尝试将椅子就地旋转，并试图在旋转过程中找到一种椅子能放稳的情形。 注意到椅脚连线呈长方形，长方形是中心对称图形，绕它的对称中心旋转180度后，椅子仍在原地。把长方形绕它的对称中心旋转，这可以表示椅子位置的改变。于是，旋转角度θ这一变量就表示了椅子的位置。为此，在平面上建立直角坐标系来解决问题。 设椅脚连线为长方形ABCD,以对角线AC 所在的直线为x 轴，对称中心O 为原点，建立平面直角坐标系。椅子绕O 点沿逆时针方向旋转角度θ后，长方形ABCD 转至A1B1C1D1的位置，这样就可以用旋转角）0（πθθ≤≤表示出椅子绕点O 旋转θ后的位置。 其次，把椅脚是否着地用数学形式表示出来。当椅脚与地面的竖直距离为零时，椅脚就着地了，而当这个距离大于零时，椅脚不着地。由于椅子在不同的位置是θ的函数，因此，椅脚与地面的竖直距离也是θ的函数。 由于椅子有四只脚，因而椅脚与地面的竖直距离有四个，它们都是θ的函数，而由假设（3）可知，椅子在任何位置至少有三只脚同时着地，即这四个函数对于任意的θ，其函数值至少有三个同时为0。因此，只需引入两个距离函数即可。考虑到长方形ABCD 是对称中心图形，绕其对称中心O 沿逆时针方向旋转180度后，长方形位置不变，但A,C 和B,D 对换了。因此，记A ，B 两脚与地面竖直距离之和为)(θf ，C,D 两脚之和为 )(θg ，其中[]πθ，0∈，使得)()(00θθg f =成立。 模型求解 如果0)0()0(== g f ，那么结论成立。
兰州交通大学 数学建模大作业 学院：机电工程学院 班级：车辆093 学号：200903812 姓名：刘键学号：200903813 姓名：杨海斌学号：200903814 姓名：彭福泰学号：200903815 姓名：程二永学号：200903816 姓名：屈辉
高速公路问题 1 实验案例 (2) 1.1 高速公路问题(简化) (2) 1.1.1 问题分析 (3) 1.1.2 变量说明 (3) 1.1.3 模型假设 (3) 1.1.4 模型建立 (3) 1.1.5 模型求解 (4) 1.1.6 求解模型的程序 (4) 1实验案例 1.1 高速公路问题(简化) A城和B城之间准备建一条高速公路，B城位于A城正南20公里和正东30公里交汇处，它们之间有东西走向连绵起伏的山脉。公路造价与地形特点有关，图4.2.4给出了整个地区的大致地貌情况，显示可分为三条沿东西方向的地形带。 你的任务是建立一个数学模型，在给定三种地形上每公里的建造费用的情况下，确定最便宜的路线。图中直线AB显然是路径最短的，但不一定最便宜。而路径ARSB过山地的路段最短，但是否是最好的路径呢？ A B 图8.2 高速公路修建地段
1.1.1 问题分析 在建设高速公路时，总是希望建造费用最小。如果要建造的起点、终点在同一地貌 中，那么最佳路线则是两点间连接的线段，这样费用则最省。因此本问题是一个典型的最优化问题，以建造费用最小为目标，需要做出的决策则是确定在各个地貌交界处的汇合点。 1.1.2 变量说明 i x ：在第i 个汇合点上的横坐标（以左下角为直角坐标原点），i ＝1，2，…，4；x 5＝30（指目的地B 点的横坐标） x=[x 1，x 2，x 3，x 4]T l i ：第i 段南北方向的长度（i ＝1，2， (5) S i ：在第i 段上地所建公路的长度（i ＝1，2， (5) 由问题分析可知， () ()() () 2 542552 432442 322332212 222 1211x x l S x x l S x x l S x x l S x l S -+=-+=-+=-+=+= C 1：平原每公里的造价（单位：万元/公里） C 2：高地每公里的造价（单位：万元/公里） C 3：高山每公里的造价（单位：万元/公里） 1.1.3 模型假设 1、 假设在相同地貌中修建高速公路，建造费用与公路长度成正比； 2、 假设在相同地貌中修建高速公路在一条直线上。在理论上，可以使得建造费用最少， 当然实际中一般达不到。 1.1.4 模型建立 在A 城与B 城之间建造一条高速公路的问题可以转化为下面的非线性规划模型。优化目标是在A 城与B 城之间建造高速公路的费用。 () 4,3,2,1300. .)(min 5142332211=≤≤++++=i x t s S C S C S C S C S C x f i
·综 述· 辐射剂量数学模型在医学影像学的应用及研究进展＊ 刘 潇综述，曾勇明△审校 （重庆医科大学附属第一医院放射科 ４０００１６） 关键词：辐射剂量；医学影像学；数学模型；蒙特卡洛；仿真人体体模 ｄｏｉ：１０．３９６９／ｊ．ｉｓｓｎ．１６７１－８３４８．２０１３．１４．０３３文献标识码：Ａ文章编号：１６７１－８３４８（２０１３）１４－１６５０－０３ 随着医学的不断发展，现代医学影像技术越来越多的应用于临床实践中，尤其是在ＣＴ、ＤＳＡ的临床应用呈逐年上升趋势，辐射剂量问题已引起全世界的关注。有效实施辐射剂量检测是保证医学影像学检查合理使用的基本要求。当前，临床上主要采用影像设备的剂量测试工具来获得辐射剂量数据，并评估患者的辐射剂量，但不能前瞻性的评价和预估某一放射学检查时的辐射水平。近年来，数学模型开始应用于医学影像学领域，对研究辐射剂量的科学实验带来便利。本文就辐射剂量数学模型的临床应用及进展综述如下。 １ 辐射剂量数学模型 为了使描述更具科学性、逻辑性、客观性和可重复性，人们采用一种普遍认为比较严格的语言来描述各种现象，这种语言就是数学。使用数学语言描述的事物就称为数学模型。数学建模是一种数学的思考方法，是运用数学的语言和方法，通过抽象、简化建立能近似刻画并“解决”实际问题的一种强有力的数学手段。数学建模是用数学语言描述实际现象的过程，不但包括外在形态，内在机制的描述，也包括预测，试验和解释实际现象等内容［１］。 以蒙特卡洛（Ｍｏｎｔｅ Ｃａｒｌｏ）为代表的数学模拟方法是一种通过设定随机过程，反复生成时间序列，计算参数估计量和统计量，进而研究其分布特征的方法。具体的，当系统中各个单元的可靠性特征量已知，但系统的可靠性过于复杂，难以建立可靠性预计的精确数学模型或模型太复杂而不便应用时，可用随机模拟法近似计算出系统可靠性的预计值；随着模拟次数的增多，其预计精度也逐渐增高。由于涉及时间序列的反复生成，蒙特卡洛模拟法是以高容量和高速度的计算机为前提条件的，因此，只是在近些年才得到广泛推广。蒙特卡洛模拟方法的原理是当问题或对象本身具有概率特征时，可以用计算机模拟的方法产生抽样结果，根据抽样计算统计量或者参数的值；随着模拟次数的增多，可以通过对各次统计量或参数的估计值求平均的方法得到稳定结论［２］。 在医学影像学中，基于蒙特卡洛模拟技术开发的软件的临床应用近年来有较大发展，如Ｉｍｐａｃｔ ＭＣ软件包（ＶＡＭＰ Ｇｍ－ｂＨ，Ｅｒｌａｎｇｅｎ，Ｇｅｒｍａｎｙ）功能独特，目前科研中，提供快速的三维剂量分布计算，该软件可以适用于多种任务，包括普通放射学、ＣＴ、Ｃ型臂（基于平板探测器）ＣＴ等。在科研中成功的剂量分布计算已经在３０多个专业领域的国际刊物也有极好的反馈［３］。还有一些通用的软件工具常在实验研究中应用，如用来模拟辐射ＣＴ剂量沉积的基于蒙特卡洛的软件ＭＣＮＰＸ ｅｘ－ｔｅｎｄｅｄ ｖ２．６，在洛杉矶洛斯阿拉莫斯国家实验室执行模拟［４］。国内应用较广的免费软件，如Ｇｅａｎｔ４［５－６］或ＭＣＮＰ ＥＧＳ４［７］，这些软件可执行辐射剂量估算。器官剂量估算软件ＰＣＸＭＣ（ＳＴＵＫ，Ｆｉｎｌａｎｄ）是基于蒙特卡洛计算方法，用于估算人体器官所受吸收剂量（ａｂｓｏｒｂｅｄｄｏｓｅ，ＡＤ）和全身有效剂量（ｅｆｆｅｃｔｉｖｅｄｏｓｅ，ＥＤ）的常用计算软件［８］。 ２ 数学模型在ＣＴ检查中的应用 在ＣＴ检查中减少辐射剂量是医学影像研究的热点问题，常用于评价ＣＴ检查的ＥＤ通过剂量长度乘积（ｄｏｓｅ ｌｅｎｇｔｈｐｒｏｄｕｃｔ，ＤＬＰ）乘以权重因子获得［９］，但与利用仿真体模检测辐射剂量的方法比较，其值不够准确［１０］。应用数学模型软件，模拟患者的辐射剂量，可避免不必要的重复照射。通过在软件中加入ＣＴ的扫描参数及患者的性别、体质量指数（ＢＭＩ）、心率等因素，因而更具个性化。辐射剂量数学模型在ＣＴ的应用已越来越受到重视。 有研究采用数学模型评估冠状动脉造影患者接受的辐射剂量，模型模拟固定管电流下ＥＤ，与常规心电门控管电流自动调制技术接受的剂量相比较。可以得到心电门控管电流自动调制技术（预设１００ｍＡｓ）的ＥＤ为（７．１±２．１）ｍＳｖ，而模拟固定管电流（１００ｍＡｓ）下肺组织的ＥＤ为（１２．５±５．３）ｍＳｖ；并证明应用心电门控管电流自动调制技术后辐射剂量减少了５２％［１１］。 Ｉｍｐａｃｔ ＭＣ软件生成的三维剂量分布是其特点，可涉及到器官剂量的估算和计算患者个体风险的ＥＤ水平。在对每个采集的参数和重建的容积数据的基础上，进行了蒙特卡洛模拟，以计算每个像素的沉积与光子相互作用方面的剂量。它可模拟现代ＣＴ系统的所有参数，比如蝶形过滤器、管电流调制、双源ＣＴ设置和动态Ｚ轴准直等。Ｉｍｐａｃｔ ＭＣ软件的可视仿真体模（ＮＶＩＤＩＡ ＧＰＵ）功能，模拟一个高精度的ＣＴ检查环境，因此Ｉｍｐａｃｔ ＭＣ是最快最全面的蒙特卡洛模拟软件包之一。为了确保最好的结果，Ｉｍｐａｃｔ ＭＣ已在三个不同的ＣＴ系统（西门子、ＧＥ、飞利浦公司产品）验证［１２］。 ＭＣＮＰＸ ｅｘｔｅｎｄｅｄ ｖ２．６软件能模拟以１ｋｅＶ的低能量辐射剂量为基准的剂量，这种软件可使用１２０ｋＶｐ、３００ｍＡ的条件下模拟全身ＣＴ扫描。针对普通患者，扫描范围可扩大，从头顶的底部到耻骨随意调节。利用蒙特卡洛技术模拟的人体数学模型，以现场调查（与临床应用相适应）所得的ＣＴ技术参数和几何条件为输入参数，从理论上估算了成人ＣＴ冠状动脉检查所接受各器官组织的吸收剂量［４］。一些免费软件（如Ｇｅａｎｔ４）缺乏灵活性，难以适应ＣＴ扫描技术的复杂多变，这些原因促进了开发以蒙特卡洛技术为仿真基础的应用于放射诊断的软件，尤其是与ＣＴ检查相关的应用软件［１３－１７］。 ３ 数学模型在介入治疗的应用 介入治疗是临床、医学与工程技术紧密结合，相互依存而发展起来的前沿学科，它具有微创、简便、安全等优点，为过去 ０ ５ ６ １重庆医学２０１３年５月第４２卷第１４期 ＊基金项目：重庆市卫生局科研基金资助项目（２０１０－２－０５５）。 作者简介：刘潇（１９８１～），技师，在读硕士研究生，主要从事医学影像技术研究。 △ 通讯作者，Ｔｅｌ：１３６０８３３８４８８；Ｅ－ｍａｉｌ：ｚｅｎｇ－ｙｍ＠ｖｉｐ．ｓｉｎａ．ｃｏｍ。
第一部分课后习题 1.学校共1000名学生，235人住在A宿舍，333人住在B宿舍，432人住在C宿舍。学 生们要组织一个10人的委员会，试用下列办法分配各宿舍的委员数： （1）按比例分配取整数的名额后，剩下的名额按惯例分给小数部分较大者。 （2）2.1节中的Q值方法。 （3）d’Hondt方法：将A，B，C各宿舍的人数用正整数n=1，2，3，…相除，其商数如下表： 将所得商数从大到小取前10个（10为席位数），在数字下标以横线，表中A，B，C行有横线的数分别为2，3，5，这就是3个宿舍分配的席位。你能解释这种方法的道理吗。 如果委员会从10人增至15人，用以上3种方法再分配名额。将3种方法两次分配的结果列表比较。 （4）你能提出其他的方法吗。用你的方法分配上面的名额。 2.在超市购物时你注意到大包装商品比小包装商品便宜这种现象了吗。比如洁银牙膏50g 装的每支1.50元，120g装的3.00元，二者单位重量的价格比是1.2：1。试用比例方法构造模型解释这个现象。 （1）分析商品价格C与商品重量w的关系。价格由生产成本、包装成本和其他成本等决定，这些成本中有的与重量w成正比，有的与表面积成正比，还有与w无关的因素。 （2）给出单位重量价格c与w的关系，画出它的简图，说明w越大c越小，但是随着w的增加c减少的程度变小。解释实际意义是什么。 3.一垂钓俱乐部鼓励垂钓者将调上的鱼放生，打算按照放生的鱼的重量给予奖励，俱乐部 只准备了一把软尺用于测量，请你设计按照测量的长度估计鱼的重量的方法。假定鱼池中只有一种鲈鱼，并且得到8条鱼的如下数据（胸围指鱼身的最大周长）： 先用机理分析建立模型，再用数据确定参数 4.用宽w的布条缠绕直径d的圆形管道，要求布条不重叠，问布条与管道轴线的夹角 应 多大（如图）。若知道管道长度，需用多长布条（可考虑两端的影响）。如果管道是其他形状呢。
郑重声明： 本作业仅供参考，可能会有错误，请自己甄别。 应用运筹学作业 6.某工厂生产A，B，C，D四种产品，加工这些产品一般需要经刨、磨、钻、镗四道工序，每种产品在各工序加工时所需设备台时如表1-18所示，设每月工作25天，每天工作8小时，且该厂有刨床、磨床、钻床、镗床各一台。问：如何安排生产，才能使月利润最大？又如A，B，C，D四种产品，每月最大的销售量分别为300件、350件、200件和400件，则该问题的线性规划问题又该如何？ 1234 四种产品的数量，则得目标函数： Max=(200?150)x1+(130?100)x2+(150?120)x3+(230?200)x4 =50x1+30x2+30x3+30x4 生产四种产品所用时间： (0.3+0.9+0.7+0.4)x1+(0.5+0.5+0.5+0.5)x2+(0.2+0.7+0.4+ 0.8)x3+(0.4+0.8+0.6+0.7)x4≤25×8 即：2.3x1+2.0x2+2.1x3+2.5x4≤200 又产品数量不可能为负，所以：x i≥0(i=1,2,3,4) 综上，该问题的线性规划模型如下： Max Z=50x1+30x2+30x3+30x4 S.T.{2.3x1+2.0x2+2.1x3+2.5x4≤200 x i≥0(i=1,2,3,4) 下求解目标函数的最优解： max=50*x1+30*x2+30*x3+30*x4; 2.3*x1+2.0*x2+2.1*x3+2.5*x4<200; Global optimal solution found. Objective value: 4347.826 Total solver iterations: 0 Variable Value Reduced Cost X1 86.95652 0.000000 X2 0.000000 13.47826 X3 0.000000 15.65217
习题2作业讲评 1. 继续考虑 2.2节的“汽车刹车距离”案例，请问“两秒准则”和“一车长度准则”一样吗？“两秒准则”是否足够安全？对于安全车距，你有没有更好的建议？（“两秒准则”，即后车司机从前车经过某一标志开始，默数2秒之后到达同一标志，而不管车速如何. 刹车距离与车速的经验公式 20.750.082678d v v =+，速度单位为m/s ，距离单位为m ） 解答 （1）“两秒准则”表明前后车距与车速成正比例关系. 引入以下符号： D ～ 前后车距（m ）；v ～ 车速（m/s ）； 于是“两秒准则”的数学模型为22D K v v ==. 与“一车长度准则”相比是否一样，依赖于一车长度的选取. 比较2 0.750.082678d v v =+与2D v =，得： ()0.082678 1.25d D v v -=- 所以当15.12 m/s v <（约合54.43 km/h ）时，有d
v=(20:5:80).*0.44704; d2=[18,25,36,47,64,82,105,132,162,196,237,283,334 22,31,45,58,80,103,131,165,202,245,295,353,418 20,28,40.5,52.5,72,92.5,118,148.5,182,220.5,266,318,376]; d2=0.3048.*d2; k1=0.75; k2=0.082678; K2=2; d1=[v;v;v].*k1; d=d1+d2; plot([0,40],[0,K2*40],'k') hold on 51015 2025 303540 车速v （m/s ） 距离（m ） 图1
西南大学2014年春《数学建模》作业及答案(已整理) 第一次作业 1：[填空题] 名词解释: 1．原型 2．模型 3．数学模型 4．机理分析 5．测试分析 6．理想方法 7．计算机模拟 8．蛛网模型 9．群体决策 10．直觉 11．灵感 12．想象力 13．洞察力 14．类比法 15．思维模型 16．符号模型 17．直观模型 18．物理模型19．2倍周期收敛20．灵敏度分析21．TSP问题22．随机存储策略23．随机模型24．概率模型25．混合整数规划26．灰色预测 参考答案： 1．原型：原型指人们在现实世界里关心、研究或者从事生产、管理的实际对象。2．模型：指为某个特定目的将原形的某一部分信息简缩、提炼而构造的原型替代物。3．数学模型：是由数字、字母或其它数字符号组成的，描述现实对象数量规律的数学公式、图形或算法。4．机理分析：根据对客观事物特性的认识，找出反映内部机理的数量规律，建立的模型常有明显的物理意义或现实意义。5．测试分析：将研究对象看作一个"黑箱”系统，通过对系统输入、输出数据的测量和统计分析，按照一定的准则找出与数据拟合得最好的模型。6．理想方法：是从观察和经验中通过想象和逻辑思维，把对象简化、纯化，使其升华到理状态，以其更本质地揭示对象的固有规律。7．计算机模拟：根据实际系统或过程的特性，按照一定的数学规律用计算机程序语言模拟实际运行情况，并依据大量模拟结构对系统或过程进行定量分析。8．蛛网模型：用需求曲线和供应曲线分析市场经济稳定性的图示法在经济学中称为蛛网模型。9．群体决策：根据若干人对某些对象的决策结果，综合出这个群体的决策结果的过程称为群体决策。10．直觉：直觉是人们对新事物本质的极敏锐的领悟、理解或推断。11．灵感：灵感是指在人有意识或下意识思考过程中迸发出来的猜测、思路或判断。12．想象力：指人们在原有知识基础上，将新感知的形象与记忆中的形象相互比较、重新组合、加工、处理，创造出新形象，是一种形象思维活动。13．洞察力：指人们在充分占有资料的基础上，经过初步分析能迅速抓住主要矛盾，舍弃次要因素，简化问题的层次，对可以用那些方法解决面临的问题，以及不同方法的优劣作出判断。14．类比法：类比法注意到研究对象与以熟悉的另一对象具有某些共性，比较二者相似之处以获得对研究对象的新认识。15．思维模型：指人们对原形的反复认识，将获取的知识以经验的形式直接储存于人脑中，从而可以根据思维或直觉作出相应的决策。16．符号模型：是在一定约束条件或假设下借助于专门的符号、线条等，按一定形式组合起来描述原型。17．直观模型：指那些供展览用的实物模型以及玩具、照片等，通常是把原型的尺寸按比例缩小或放大，主要追求外观上的逼真。18．物理模型：主要指科技工作者为一定的目的根据相似原理构造的模型，它不仅可以显示原型的外形或某些特征，而且可以用来进行模拟实验，间接地研究原型的某些规律。19．2倍周期收敛：在离散模型中，如果一个数列存在两个收敛子列就称为2倍周期收敛。20．灵敏度分析：系数的每个变化都会改变线性规划问题，随之也会影响原来求得的最优解。为制定一个应付各种偶然情况的全能方法，必须研究以求得的最优解是怎样随输入系数的变化而变化的。这叫灵敏性分析。21．TSP问题：在加权图中寻求最佳推销员回路的问题可以转化为在一个完备加权图中寻求最佳哈密顿圈的问题，称为TSP问题。22．随机存储策略：商店在订购货物时采用的一种简单的策略，是制定一个下界s和一个上界S，当周末存货不小于s时就不定货；当存货少于s 时就订货，且定货量使得下周初的存量达到S，这种策略称为随机存储策略。23．随机模型：如果随机因素对研究对象的影响必须考虑，就应该建立随机性的数学模型，简称为随机模型。24．概
数学建模 第一章 线性规划 §1 线性规划 在人们的生产实践中，经常会遇到如何利用现有资源来安排生产，以取得最大经济效益的问题。此类问题构成了运筹学的一个重要分支—数学规划，而线性规划(Linear Programming 简记LP)则是数学规划的一个重要分支。自从1947年G. B. Dantzig 提出求解线性规划的单纯形方法以来，线性规划在理论上趋向成熟，在实用中日益广泛与深入。特别是在计算机能处理成千上万个约束条件和决策变量的线性规划问题之后，线性规划的适用领域更为广泛了，已成为现代管理中经常采用的基本方法之一。 1.1 线性规划的实例与定义 例1 某机床厂生产甲、乙两种机床，每台销售后的利润分别为4000元与3000元。生产甲机床需用B A 、机器加工，加工时间分别为每台2小时和1小时；生产乙机床需用C B A 、、三种机器加工，加工时间为每台各一小时。若每天可用于加工的机器时数分别为A 机器10小时、B 机器8小时和C 机器7小时，问该厂应生产甲、乙机床各几台，才能使总利润最大？ 上述问题的数学模型：设该厂生产1x 台甲机床和2x 乙机床时总利润最大，则21,x x 应满足 （目标函数） 2134m ax x x z += (1) s.t. （ 约 束 条 件 ） ?????? ?≥≤≤+≤+0 ,781022122 121x x x x x x x （2） 这里变量21,x x 称之为决策变量，（1）式被称为问题的目标函数，（2）中的几个不等式是问题的约束条件，记为s.t.(即subject to)。
上述即为一规划问题数学模型的三个要素。由于上面的目标函数及约束条件均为线性函数，故被称为线性规划问题。 总之，线性规划问题是在一组线性约束条件的限制下，求一线性目标函数最大或最小的问题。 在解决实际问题时，把问题归结成一个线性规划数学模型是很重要的一步，但往往也是困难的一步，模型建立得是否恰当，直接影响到求解。而选取适当的决策变量，是我们建立有效模型的关键之一。 1.2 线性规划的Matlab 标准形式 线性规划的目标函数可以是求最大值，也可以是求最小值，约束条件的不等号可以是小于号也可以是大于号。为了避免这种形式多样性带来的不便，Matlab 中规定线性规划的标准形式为 b Ax x c x T ≤ that such min 其中c 和x 为n 维列向量，b 为m 维列向量，A 为n m ?矩阵。 例如线性规划 b Ax x c x T ≥ that such max 的Matlab 标准型为 b Ax x c x T -≤-- that such min 1.3 线性规划问题的解的概念 一般线性规划问题的标准型为 ∑==n j j j x c z 1min (3) ∑==≤n j i j ij m i b x a 1,,2,1 s.t.Λ (4) 可行解 满足约束条件（4）的解),,,(21n x x x x Λ=，称为线性规划问题的可行解，而使目标函数（3）达到最小值的可行解叫最优解。
数学建模章绍辉版作业集团标准化工作小组 #Q8QGGQT-GX8G08Q8-GNQGJ8-MHHGN#
第四章作业 第二题： 针对严重的交通情况，国家质量监督检验检疫局发布的国家标准，车辆驾驶人员血液中的酒精含量大于或等于20mg/100ml,小于80mg/100ml 为饮酒驾车，血液中的酒精含量大于或等于80mg/100ml 的为醉酒驾车。 下面分别考虑大李在很短时间内和较长时间内（如2个小时）喝了三瓶啤酒，多长时间内驾车就会违反新的国家标准。 1、 问题假设 大李在短时间内喝下三瓶啤酒后，酒精先从吸收室（肠胃）吸收进中心室（血液和体液），然后从中心室向体外排除，忽略喝酒的时间，根据生理学知识，假设 （1） 吸收室在初始时刻t=0时，酒精量立即为 32 D ;在任意时刻，酒精从吸收室吸收进中心室的速率（吸收室在单位时间内酒精含量的减少量）与吸收室的酒精含量成正比，比例系数为1k ； （2） 中心室的容积V 保持不变；在初始时刻t=0时，中心室的酒精含量为0； 在任意时刻，酒精从中心室向体外排除的速率（中心室在单位时间内酒精含量的减少量）与中心室的酒精含量成正比，比例系数为2k ； （3） 在大李适度饮酒没有酒精中毒的前提下，假设1k 和2k 都是常量，与饮酒量 无关。 2、 符号说明 酒精量是指纯酒精的质量，单位是毫克； 酒精含量是指纯酒精的浓度，单位是毫克/百毫升； ~t 时刻（小时）； ()~x t 在时刻t 吸收室（肠胃）内的酒精量（毫克）；
0~D 两瓶酒的酒精量（毫克）； (t)~c 在时刻t 吸收室（血液和体液）的酒精含量（毫克/百毫升）； 2()~c t 在时刻t 中心室（血液和体液）的酒精含量（毫克/百毫升）； ~V 中心室的容积（百毫升）； 1~k 酒精从吸收室吸收进中心室的速率系数（假设其为常数）； 2~k 酒精从中心室向体外排除的速率系数（假设其为常数）； 3~k 在短时间喝下三瓶酒的假设下是指短时间喝下的三瓶酒的酒精总量除以中心室体 积，即03/2D V ；而在较长时间内（2小时内）喝下三瓶酒的假设下就特指03/4D V . 3、 模型建立和求解 （1） 酒是在很短时间内喝的： 记喝酒时刻为0t =（小时），设(0)0c =，可用()2113 212 ()k t k t k k c t e e k k --= --来计算血液中的酒精含量，此时12k k 、为假设中所示的常数，而033155.792D k V ?? == ??? . 下面用MATLAB 程序画图展示血液中酒精含量随时间变化并且利用fzero 函数和fminbnd 函数来得到饮酒驾车醉酒驾车对应的时间段，以及血液中酒精含量最高的时刻。 MATLAB 程序如下： k1=;k2=;k3=; c=@(t)(k1.*k3)./(k1-k2).*(exp(-k2.*t)-exp(-k1.*t)); f=@(t)c(t)-20; g=@(t)c(t)-80; h=@(t)-c(t); t1(1)=fzero(f,1);t1(2)=fzero(f,12), t2(1)=fzero(g,1);t2(2)=fzero(g,12) [t3,c3]=fminbnd(h,0,24) fplot(c,[0,20],'k') hold on plot([0,20],[20,20],'k',[0,20],[80,80],'k') hold off
09级数模试题 1. 把四只脚的连线呈长方形的椅子往不平的地面上一放，通常只有三只脚着地，放不稳，然后稍微挪动几次，就可以使四只脚同时着地，放稳了。试作合理的假设并建立数学模型说明这个现象。（15分） 解：对于此题，如果不用任何假设很难证明，结果很 可能是否定的。 因此对这个问题我们假设： （1）地面为连续曲面 （2）长方形桌的四条腿长度相同 （3）相对于地面的弯曲程度而言，方桌的腿是足够长的 （4）方桌的腿只要有一点接触地面就算着地。 那么，总可以让桌子的三条腿是同时接触到地面。 现在，我们来证明：如果上述假设 条件成立，那么答案是肯定的。以长方 桌的中心为坐标原点作直角坐标系如图 所示，方桌的四条腿分别在A、B、C、D 处，A、、D的初始位置在与x轴平行，再 假设有一条在x轴上的线,则也与A、B，C、D平行。当方桌绕中心0旋转时，对角线与x轴的夹角记为θ。 容易看出，当四条腿尚未全部着地时，腿到地面的距离是不确定的。为消除这一不确定性，令() fθ为A、B离地距离之和，
()g θ为C 、D 离地距离之和，它们的值由θ唯一确定。由假设（1）， ()f θ,()g θ均为θ的连续函数。又由假设（3） ，三条腿总能同时着地， 故()f θ()g θ=0必成立（?θ）。不妨设(0)0f =(0)0g >（若(0)g 也为0，则初始时刻已四条腿着地，不必再旋转），于是问题归结为： 已知()f θ,()g θ均为θ的连续函数，(0)0f =,(0)0g >且对任意θ有00()()0f g θθ=，求证存在某一0θ，使00()()0f g θθ=。 证明：当θ=π时，与互换位置，故()0f π>，()0g π=。作()()()h f g θθθ=-，显然，()h θ也是θ的连续函数，(0)(0)(0)0h f g =-<而()()()0h f g πππ=->，由连续函数的取零值定理，存在0θ，00θπ<<，使得0()0h θ=，即00()()f g θθ=。又由于00()()0f g θθ=，故必有00()()0f g θθ==，证毕。 2.学校共1000名学生，235人住在A 宿舍，333人住在B 宿舍，432人住在C 宿舍。学生 们要组织一个10人的委员会，试用合理的方法分配各宿舍的委员数。（15分） 解：按各宿舍人数占总人数的比列分配各宿舍的委员数。设：A 宿舍的委员数为x 人，B 宿舍的委员数为y 人，C 宿舍的委员数为z 人。计算出人数小数点后面的小数部分最大的整数进1，其余取整数部分。 则 10； 10=235/1000；
数学建模第二章作业答案章绍辉
习题2作业讲评 1. 继续考虑 2.2节的“汽车刹车距离”案例，请问“两秒准则”和“一车长度准则”一样吗？“两秒准则”是否足够安全？对于安全车距，你有没有更好的建议？（“两秒准则”，即后车司机从前车经过某一标志开始，默数2秒之后到达同一标志，而不管车速如何. 刹车距离与车速的经验公式 20.750.082678d v v =+，速度单位为m/s ，距离单位为m ） 解答 （1）“两秒准则”表明前后车距与车速成正比例关系. 引入以下符号： D ～ 前后车距（m ）；v ～ 车速（m/s ）； 于是“两秒准则”的数学模型为22D K v v ==. 与“一车长度准则”相比是否一样，依赖于一车长度的选取. 比较2 0.750.082678d v v =+与2D v =，得： ()0.082678 1.25d D v v -=- 所以当15.12 m/s v <（约合54.43 km/h ）时，有d
v=(20:5:80).*0.44704; d2=[18,25,36,47,64,82,105,132,162,196,237,283,334 22,31,45,58,80,103,131,165,202,245,295,353,418 20,28,40.5,52.5,72,92.5,118,148.5,182,220.5,266,318,376]; d2=0.3048.*d2; k1=0.75; k2=0.082678; K2=2; d1=[v;v;v].*k1; d=d1+d2; plot([0,40],[0,K2*40],'k') hold on plot(0:40,polyval([k2,k1,0],0:40),':k') plot([v;v;v],d,'ok','MarkerSize',2) title('比较刹车距离实测数据、理论值和两秒准则') legend('两秒准则','刹车距离理论值',... '刹车距离的最小值、平均值和最大值',2) xlabel('车速v （m/s ）') ylabel('距离（m ）') hold off 51015 2025 303540 020406080100120 140160180比较刹车距离实测数据、理论值和两秒准则 车速v （m/s ） 距离（m ） 两秒准则 刹车距离理论值 刹车距离的最小值、平均值和最大值 图1
1.你要在雨中从一处沿直线走到另一处，雨速是常数，方向不变。 你是否走得越快，淋雨量越少呢？ 2.假设在一所大学中，一位普通教授以每天一本的速度开始从图书 馆借出书。再设图书馆平均一周收回借出书的1/10，若在充分长的时间内，一位普通教授大约借出多少年本书？ 3.一人早上6：00从山脚A上山，晚18：00到山顶B；第二天，早 6：00从B下山，晚18：00到A。问是否有一个时刻t,这两天都在这一时刻到达同一地点？ 4.如何将一个不规则的蛋糕I平均分成两部分？ 5.兄妹二人沿某街分别在离家3公里与2公里处同向散步回家，家 中的狗一直在二人之间来回奔跑。已知哥哥的速度为3公里/小时，妹妹的速度为2公里/小时，狗的速度为5公里/小时。分析半小时后，狗在何处？ 6.甲乙两人约定中午12：00至13：00在市中心某地见面，并事先 约定先到者在那等待10分钟，若另一个人十分钟内没有到达，先到者将离去。用图解法计算，甲乙两人见面的可能性有多大？ 7.设有n个人参加某一宴会，已知没有人认识所有的人，证明：至 少存在两人他们认识的人一样多。 8.一角度为60度的圆锥形漏斗装着10 端小孔的 面积为0.5 9.假设在一个刹车交叉口，所有车辆都是由东驶上一个1/100的斜
坡，计算这种情 下的刹车距离。如果汽车由西驶来，刹车距离又是多少？ 10. 水管或煤气管经常需要从外部包扎以便对管道起保护作用。包扎时用很长的带子缠绕在管道外部。为了节省材料，如何进行包扎才能使带子全部包住管道而且带子也没有发生重叠。 :顶=1：a:b ，选坐v>0,而设语雨速 L( 1q -+v x ),v≤x Q(v)= L( v x -q +1),v>x 2.解：由于教授每天借一本书，即一周借七本书，而图书馆平均每周
数学建模期末大作业论文 题目：A题美好的一天 组长：何曦（2014112739） 组员：李颖（2014112747）张楚良（2014112740） 班级：交通工程三班 指导老师：陈崇双
美好的一天 摘要 关键字：Dijkstra算法多目标规划有向赋权图 MATLAB SPSS
1 问题的重述 Hello！大家好，我是没头脑,住在西南宇宙大学巨偏远的新校区（节点22）。明天我一个外地同学来找我玩，TA叫不高兴，是个镁铝\帅锅，期待ing。我想陪TA在城里转转，当然是去些不怎么花钱的地方啦~~。目前想到的有林湾步行街(节点76)、郫郫公园（节点91），大川博物院（节点72）。交通嘛，只坐公交车好了,反正公交比较发达，你能想出来的路线都有车啊。另外，进城顺便办两件事，去老校区财务处一趟（节点50），还要去新东方（节点34）找我们宿舍老三，他抽奖中了两张电影票，我要霸占过来明晚吃了饭跟TA一起看。电影院嘛，TASHIWODE电影院（节点54）不错，比较便宜哈。我攒了很久的钱，订了明晚开心面馆（节点63）的烛光晚餐，额哈哈，为了TA，破费一下也是可以的哈。哦，对了，老三说了，他明天一整天都上课，只有中午休息的时候能接见我给我票。 我主要是想请教一下各位大神： 1）明天我应该怎么安排路线才能够让花在坐车上的时间最少？ 2）考虑到可能堵车啊，TA比较没耐心啊，因为TA叫不高兴嘛。尤其是堵车啊，等车啊，这种事，万一影响了气氛就悲剧了。我感觉路口越密的地方越容易堵，如果考虑这个，又应该怎么安排路线呢？ 3）我们城比较挫啊，连地图也没有，Z老师搞地图测绘的，他有地图，跟他要他不给，只给了我一个破表格（见附件，一个文件有两页啊），说“你自己画吧”。帮我画一张地图吧，最好能标明我们要去的那几个地方和比较省时的路线啊，拜托了~ 2 问题的分析 2.1 对问题一的分析 问题一要求安排路线使得坐车花费的时间最少。 对于问题一，假设公交车的速度维持不变，要使花费的时间最少，则将问题转化为对最短路径的求解。求解最短路径使用Dijkstra算法很容易进行求解，在运用MATLAB编程，得到最优的一条路径，则这条路径所对应的时间即为最少用时。 2.2 对问题二的分析 问题二要求在考虑堵车的情况下，路口越密越容易发生拥堵，安排路线是乘车时间最短。 对于问题二，在问题的基础上增加了附加因素，即公交车的速度会因道路的密集程度而发生改变，从而问题一建立的基本Dijkstra算法对于问题二就不再适用了，因此对问题一的基本Dijkstra算法进行改进，并结合蚁群算法的机理与特点，运用MATLAB求解出最短路径，保证了花费时间的最少性。 2.3 对问题三的分析 问题三要求根据提供的附件，画出一张地图，标明要去的那几个地方和比较省时的路线。 对于问题三，在问题一和问题二的基础上，根据求解的结果，运用SPSS软件画出地图。
1．（1） n=101; x1=linspace(-1,1,n); x2=linspace(-2,2,n); y1=[sqrt(1-x1.^2);-sqrt(1-x1.^2)]; y2=[sqrt(4-x2.^2);-sqrt(4-x2.^2);sqrt(1-(x2.^2)/4);-sqrt(1-(x2.^2)/4)]; plot(x1,y1) … hold on; plot(x2,y2) title('椭圆x^2/4+y^2=1的内切圆和外切圆') axis equal -2.5 -2-1.5-1-0.500.51 1.52 2.5 -2-1.5-1-0.500.511.5 2椭圆x 2/4+y 2=1的内切圆和外切圆 （2） x1=linspace(-2,2,101); / x2=linspace(-2,8); axis equal plot(exp(x1),x1,x1,exp(x1),x2,x2) title('指数函数y=exp(x)和对数函数y=ln(x)关于y=x 对称')
-2 -1 1 2 3 4 5 6 7 8 -2-101234567 8指数函数y=exp(x)和对数函数y=ln(x)关于y=x 对称 （3） hold on — q=input('请输入一个正整数q;') for i=1:q for j=1:i if rem(j,i) plot(j/i,1/i) end end end @
0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 00.050.10.150.20.250.30.350.40.45 0.5 3．代码如下： n=input('请输入实验次数n=') k=0; for i=1:n 。 x=ceil(rand*6)+ceil(rand*6); if x ==3|x==11 k=k+1; elseif x~=2&x~=7&x~=12 y= ceil(rand*6)+ceil(rand*6); while y~=x&y~=7 y=ceil(rand*6)+ceil(rand*6); end if y==7 ； k=k+1; end end end
数学建模作业 姓名：叶勃 学号： 班级：024121
一：层次分析法 1、 分别用和法、根法、特征根法编程求判断矩阵 1261/2141/61/41A ????=?????? 11/2433 217551/4 1/711/21/31/31/52111/31/5 3 1 1A ????????=? ?????? ? 的特征根和特征向量 （1）冪法求该矩阵的特征根和特征向量 程序为： #include X[i]=0; for(j=0;j 数学建模模拟试题及答案 一、填空题（每题5分，共20分） 1. 若,, x z z y ∝∝则y 与x 的函数关系是. 2. 在超级市场的收银台有两条队伍可选择，队1有1m 个顾客，每人都买了1n 件商品，队2有2m 个顾客，每人都买了2n 件商品，假设每个人付款需p 秒，而扫描每件商品需t 秒，则加入较快队1的条件是 . 3. 马尔萨斯与罗捷斯蒂克两个人口增长模型的主要区别是假设了 4. 在研究猪的身长与体重关系时，我们通过与已知其相关性质的的弹性梁作 的方法建立了模型. 二、分析判断题（每小题15分，满分30分） 1. 要为一所大学编制全校性选修课程表，有哪些因素应予以考虑？试至少列出5种. 2. 一起交通事故发生3个小时后，警方测得司机血液中酒精的含量是 ),m l /m g (100/56 又过两个小时，含量降为),m l /m g (100/40试判断，当事故发生时，司 机是否违反了酒精含量的规定（不超过80/100)m l /m g (. （提示：不妨设开始时刻为)(,0t C t =表示t 时刻血液中酒精的浓度，则依平衡原理，在时间间隔],[t t t ?+内酒精浓度的改变量为 t t kC t C t t C ??=??+)()()( 其中0>k 为比例常数，负号则表示了浓度随时间的推移是递减的.） 三、计算题（每题25分，满分50分） 1. 一个毛纺厂使用羊毛、兔毛和某种纤维生产甲、乙两种混纺毛料，生产一个单位产品甲需要的三种原料依次为3、2、8个单位，产值为580元；生产一个单位产品乙需要的三种原料依次为2、3、5个单位，产值为680元，三种原料在计划期内的供给量依次为90、30和80单位.试建立线性规划模型以求一个生产方案，使得总产值达到最大，并由此回答： （1） 最优生产方案是否具有可选择余地？若有请至少给出两个，否则说明理由. （2） 原材料的利用情况.数学建模模拟试题及答案.pdf
相关主题
文本预览