数学建模案例分析线性代数建模案例20例Word文件下载.docx
《数学建模案例分析线性代数建模案例20例Word文件下载.docx》由会员分享,可在线阅读,更多相关《数学建模案例分析线性代数建模案例20例Word文件下载.docx(53页珍藏版)》请在冰豆网上搜索。
X4100
X3
300
X4300
其增广矩阵
110
1
100
(Ab)=100
初等行变换0
600
(A,b)011
001
由此可得
X
X4
即
X4600.
为了唯一确定未知流量,只要增添X4
统计的值即可.
当X4=350时,确定X1=250,X2;
=250,X3
=50.
若X4=200,贝UX1=100,X2=
:
400,X3=
100<
0.这表明单行线
“③
改为“③④”才合理.
【模型分析】
(1)由(A,b)的行最简形可见,上述方程组中的最后一个方程是多余的.这意味着最后一个方程中的数据“300”可以不用统计.
⑵由X2
可得
200,
x-ix3200
XX3300,这
x4x3300
就是说X1,X2,X3,X4这四个未知量中,任意一个未知量的值统计出来之后都可以确定出其他三个未知量的值•
Matlab实验「题
某城市有下图所示的交通图,每条道路都是单行线,需要调查每条道路每小时的车流量•图中的数字表示该条路段的车流数•如果每个交叉路口进入和离开的车数相等,整个图中进入和离开的车数相等•
150400290
图4某城市单行线车流量
(2)分析哪些流量数据是多余的.
(3)为了唯一确定未知流量,需要增添哪几条道路的流量统计
案例二.配方问题
在化工、医药、日常膳食等方面都经常涉及到配方问题•在不考虑各种成分之间可能发生某些化学反应时,配方问题可以用向量和线性方程组来建模•
模型准备】一种佐料由四种原料A、B、C、D混合而成.这种佐料现有两种规格,这两种规格的佐料中,四种原料的比例分别为2:
3:
1:
1和1:
2:
12现在需要四种原料的比例为4:
7:
5的第三种规格的佐料•问:
第三种规格的佐料能否由前两种规格的佐料按一定比例配制而成?
模型假设】
(1)假设四种原料混合在一起时不发生化学变化•
(2)假设四种原料的比例是按重量计算的•(3)假设前两种规格的佐料分装成袋,比如说第一种规格的佐料每袋净重7克(其中A、B、C、D四种原料分别为2克,3克,1克,1克),第二种规格的佐料每袋净重6克(其中A、B、C、D四种原料分别为1克,2克,1克,2克).
【模型建立】根据已知数据和上述假设,可以进一步假设将x袋第一种规格的佐料与y袋第二种规格的佐料混合在一起,得到的混合物中A、B、C、D四种原料分别为4克,7克,3克,5克,则有以下线性方程组
2xy4,3x2y7,xy3,x2y5.
【模型求解】上述线性方程组的增广矩阵
2
4
3
7
初等行变换
(A,b)=1
5
可见x又因为第一种规格的佐料每袋净重7克,第二种规格的佐料每袋净重6
y2.
克,所以第三种规格的佐料能由前两种规格的佐料按7:
12的比例配制而成.
(1)若令1=(2,3,1,1),2=(1,2,1,1〕,=(4,7,5,3〕,则原问题等价于线性方程组Ax=b是否有解”,也等价于“能否由1,2线性表示”.
(2)若四种原料的比例是按体积计算的,则还要考虑混合前后体积的关系(未必是简单的叠加),因而最好还是先根据具体情况将体积比转换为重量比,然后再按上述方法处理.
(3)上面的模型假设中的第三个假设只是起到简化运算的作用.如果直接设x克第一种规格的佐料与y克第二种规格的佐料混合得第三种规格的佐料,则有下表
表1混合后四种原料的含量
原料佐料规格
A
B
C
D
第一种
_x
_x
第二种
_y
6
第三种
-49(x+y)
19
+y)
19(x+y>
-59(x+y)
因而有如下线性方程组
4“
x
-y
(x
y),
7/
y
—x
5/
y).
【模型检验】把x=7,y=12代入上述方程组(),则各等式都成立.可见模型假设中的第三个假设不影响解的正确性•
Matlab实验题
蛋白质、碳水化合物和脂肪是人体每日必须的三种营养,但过量的脂肪摄入不利于健康.人们可以通过适量的运动来消耗多余的脂肪.设三种食物(脱脂牛奶、大豆面粉、乳清)每100克中蛋白质、碳水化合物和脂肪的含量以及慢跑5分钟消耗蛋白质、碳水化合物和脂肪的量如下表•
表2三种食物的营养成分和慢跑的消耗情况
营养
每100克食物所含营养(克)
慢跑5分钟
每日需要的营
牛奶
大豆面粉
乳清
消耗量(克)
养量(克)
蛋白质
36
51
13
10
33
碳水化合物
52
34
74
20
45
脂肪
15
问怎样安排饮食和运动才能实现每日的营养需求?
案例三.投入产出问题
在研究多个经济部门之间的投入产出关系时,W.Leontief提出了投入产出模型.这为经济学研究提供了强有力的手段•W.Leontief因此获得了1973年的Nobel经济学奖•
【模型准备】某地有一座煤矿,一个发电厂和一条铁路•经成本核算,每生产价值1元钱的煤需消耗0.3元的电;
为了把这1元钱的煤运出去需花费0.2元的运费;
每生产1元的电需0.6元的煤作燃料;
为了运行电厂的辅助设备需消耗本身0.1元的电,还需要花费0.1元的运费;
作为铁路局,每提供1元运费的运输需消耗0.5元的煤,辅助设备要消耗0.1元的电.现煤矿接到外地6万元煤的订货,电厂有10万元电的外地需求,问:
煤矿和电厂各生产多少才能满足需求?
【模型假设】假设不考虑价格变动等其他因素.
【模型建立】设煤矿,电厂,铁路分别产出x元,y元,z元刚好满足需求.则有下表
表3消耗与产出情况
产出(1元)
产出
消耗
订单
煤
电
运
消耗
0.6
0.5
0.6y+0.5z
60000
0.3
0.1
0.3x+0.1y+0.1z
100000
0.2
z
0.2x+0.1y
根据需求,应该有
(0.6y
0.5z)
(0.3x
0.1y
0.1z)100000,
(0.2x
0.1y)
x0.6y
r0.5z
0.3x
0.9y
0.1z100000
0.2x
z0
【模型求解】在Matlab命令窗口输入以下命令
>
A=[1,-0.6,-0.5;
-0.3,0.9,-0.1;
-0.2,-0.1,1];
b=[60000;
100000;
0];
x=A\b
Matlab执行后得
x=
1.0e+005*
1.9966
1.8415
0.5835
可见煤矿要生产1.9966105元的煤,电厂要生产1.8415105元的电恰好满足需求.
x00.60.560000
【模型分析】令x=y,A=0.30.10.1,b=100000,其中x称为总产值列向z0.20.100
量,A称为消耗系数矩阵,b称为最终产品向量,则
00.60.5x0.6y0.5z
Ax=0.30.10.1y=0.3x0.1y0.1z
0.20.10z0.2x0.1y
根据需求,应该有xAx=b,即(EA)x=b.故x=(EA)1b.
某乡镇有甲、乙、丙三个企业.甲企业每生产1元的产品要消耗0.25元乙企业的产品和0.25元丙企业的产品.乙企业每生产1元的产品要消耗0.65元甲企业的产品,0.05元自产的产品和0.05元丙企业的产品.丙企业每生产1元的产品要消耗0.5元甲企业的产品和0.1元乙企业的产品.在一个生产周期内,甲、乙、丙三个企业生产的产品价值分别为100万元,120万元,60万元,同时各自的固定资产折旧分别为20万元,5万元和5万元.
(1)求一个生产周期内这三个企业扣除消耗和折旧后的新创价值.
⑵如果这三个企业接到外来订单分别为50万元,60万元,40万元,那么他们
各生产多少才能满足需求?
案例四.平板的稳态温度分布问题
在热传导的研究中,一个重要的问题是确定一块平板的稳态温度分布•根据…定律,只要测定一块矩形平板四周的温度就可以确定平板上各点的温度•
*<
11530
图8一块平板的温度分布图
【模型准备】如图9所示的平板代表一条金属梁的截面•已知四周8个节点处的温度(单位°
C),求中间4个点处的温度Ti,T2,T3,T4.
80
‘T'
90Tl
*■
T2
■*
60
to.Jh
80T3
60-
T4
5^
rII
50
►
图9一块平板的温度分布图
【模型假设】假设忽略垂直于该截面方向上的热传导,并且每个节点的温度等于与它相邻的四个节点温度的平均值•
【模型建立】根据已知条件和上述假设,
有如下线性方程组
T1(90
T3)
T2(80
T1
T4)
T3(80
T4—(50
【模型求解】将上述线性方程组整理得
4T1
T2T3
190
4T2
T4140
4T3
4T4100
在Matlab命令窗口输入以下命令
A=[4,-1,-1,0;
-1,4,0,-1;
-1,0,4,-1;
0,-1,-1,4];
b=[190;
140;
100];
x=A\b;
x'
ans=
82.916770.833370.833360.4167
可见T1=82.9167,T2=70.8333,T3=70.8333,T4=60.4167.
参考文献
陈怀琛,高淑萍,杨威,工程线性代数,北京:
电子工业出版社,2007.页码:
15-16.
假定下图中的平板代表一条金属梁的截面,并忽略垂直于该截面方向上的热传导.已知平板内部有30个节点,每个节点的温度近似等于与它相邻的四个节点温度的平均值.设4条边界上的温度分别等于每位同学学号的后四位的5倍,例如学号为16308209的同学计算本题时,选择Ti=40,Tu=10,Tr=0,Td=45.
1
Tu
|b—4
►<
■■■■■■■■■■■■■
»
和1
TiT1
0<
T6
<
f■
T26
Tr
Ti
I
T27
i
T7
T5
T10
Td4
图10一块平板的温度分布图
(1)建立可以确定平板内节点温度的线性方程组.
(2)用Matlab软件求解该线性方程组.
(3)用Matlab中的函数mesh绘制三维平板温度分布图.
案例五.CT图像的代数重建问题
X射线透视可以得到3维对象在2维平面上的投影,CT则通过不同角度的X射线得到3维对象的多个2维投影,并以此重建对象内部的3维图像•代数重建方法就是从这些2维投影出发,通过求解超定线性方程组,获得对象内部3维图像的
方法.
图11双层螺旋CT
图12CT图像
这里我们考虑一个更简单的模型,从2维图像的1维投影重建原先的2维图像.一个长方形图像可以用一个横竖均匀划分的离散网格来覆盖,每个网格对应一个像素,它是该网格上各点像素的均值•这样一个图像就可以用一个矩阵表示,其元素就是图像在一点的灰度值(黑白图像).下面我们以33图像为例来说明.
表4消耗与产出情况
33图像
各点的灰度值
水平方向上的叠加值
X1=1
X2=0
X3=0
X1+X2+X3=1
X4=0
X5=0.5
X6=0.5
X4+X5+X6=1
X7=0.5
X8=0
X9=1
X7+X8+X9=1.5
竖直方向上的叠加值
X1+X4+X7
=1.5
X2+X5+X8
=0.5
X3+X6+X9
每个网格中的数字Xi代表其灰度值,范围在[0,1]内.0表示白色,1表示黑色,0.5表示灰色.如果我们不知道网格中的数值,只知道沿竖直方向和水平方向的叠加值,为了确定网格中的灰度值,可以建立线性方程组(含有6个方程,9个未知数)
X5
X6
L
Xg
显然该方程组的解是不唯一的,为了重建图像,必须增加叠加值.如我们增加从右上方到左下方的叠加值,则方程组将增加5个方程
X1=1,
X2+X4=0,
X3+X5+X7=1,
X6+X8=0.5,
X9=1,
和上面的6个方程放在一起构成一个含有11个方程,9个未知数的线性方程组.
【模型准备】设33图像中第一行3个点的灰度值依次为xi,x2,X3,第二行3个点
的灰度值依次为X4,X5,X6,第三行3个点的灰度值依次为X7,X8,X9.沿竖直方向的叠加值依次为1.5,0.5,1.5,沿水平方向的叠加值依次为1,1,1.5,沿右上方到左下方的叠加值依次为1,0,1,0.5,1.确定X1,X2,…,X9的值.
【模型建立】由已知条件可得(含有11个方程,9个未知数的)线性方程组
X2X31
X5X61
X91
A=[1,1,1,0,0,0,0,0,0;
0,0,0,1,1,1,0,0,0;
0,0,0,0,0,0,1,1,1;
1,0,0,1,0,0,1,0,0;
0,1,0,0,1,0,0,1,0;
0,0,1,0,0,1,0,0,1;
1,0,0,0,0,0,0,0,0;
0,1,0,1,0,0,0,0,0;
0,0,1,0,1,0,1,0,0;
0,0,0,0,0,1,0,1,0;
0,0,0,0,0,0,0,0,1];
b=[1;
1;
1.5;
0.5;
0;
1];
x=A\b;
Warning:
Rankdeficient,rank=8tol=4.2305e-015.
1.00000.00000-0.00000.50000.50000.5000-0.00001.0000
可见上述方程组的解不唯一.其中的一个特解为
X1=1,X2=0,X3=0,X4=0,X5=0.5,X6=0.5,X7=0.5,X8=0,X9=1.
【模型分析】上述结果表明,仅有三个方向上的叠加值还不够.可以再增加从左上方到右下方的叠加值.在实际情况下,由于测量误差,上述线性方程组可能是超定的这时可以将超定方程组的近似解作为重建的图像数据.
给定一个33图像的2个方向上的灰度叠加值:
沿左上方到右下方的灰度叠加值依次为0.8,1.2,1.7,0.2,0.3;
沿右上方到左下方的灰度叠加值依次为0.6,0.2,1.6,1.2,0.6.
(1)建立可以确定网格数据的线性方程组,并用Matlab求解.
⑵将网格数据乘以256,再取整,用Matlab绘制该灰度图像.
案例六.平衡结构的梁受力计算
在桥梁、房顶、铁塔等建筑结构中,涉及到各种各样的梁•对这些梁进行受力分析是设计师、工程师经常做的事情.
图14埃菲尔铁塔局部
F面以双杆系统的受力分析为例,说明如何研究梁上各铰接点处的受力情况.
【模型准备】在图15所示的双杆系统中,已知杆1重Gi=200牛顿,长Li=2米与水平方向的夹角为i=/6,杆2重G2=100牛顿,长L2=.2米与水平方向的夹角为2=/4.三个铰接点AB,C所在平面垂直于水平面.求杆1,杆2在铰接
点处所受到的力•
【模型假设】假设两杆都是均匀的•在铰接点处的受力情况如图16所示.
【模型建立】对于杆1:
水平方向受到的合力为零,故N1=N3,
竖直方向受到的合力为零,故N2+N4=G1,
以点A为支点的合力矩为零,故(L1sin1)N3+(L1cos1)N4=(―L1cos1)G1.
N5
图16两杆受力情况
对于杆2类似地有
1N5=N7,N6=N8+G2,(L2Sin2)N7=(L2C0S2)N8+(L2C0S2)G2.
此外还有N3=N7,N4=N8.于是将上述8个等式联立起来得到关于N1,N2,…,N8的线性方程组:
N1N30
N2N4G1
N4N80
G1=200;
L1=2;
theta仁pi/6;
G2=100;
L2=sqrt
(2);
theta2=pi/4;
A=[1,0,-1,0,0,0,0,0;
0,1,0,1,0,0,0,0;
0,0,L1*sin(theta1),L1*cos(theta1),0,0,0,0;
0,0,0,0,1,0,-1,0;
0,0,0,0,0,1,0,-1;
0,0,0,0,0,0,L2*sin(theta2),-L2*cos(theta2);
0,0,1,0,0,0,-1,0;
0,0,0,1,0,0,0,-1];
b=[0;
G1;
0.5*L1*cos(theta1)*G1;
G2;
0.5*L2*cos(theta2)*G2;
Matlab执行后得ans=
95.0962154.903895.096245.096295.0962145.096295.096245.0962
【模型分析】最后的结果没有出现负值,说明图16中假设的各个力的方向与事实一致.如果结果中出现负值,则说明该力的方向与假设的方向相反.
陈怀琛,高淑萍,杨威,工