3-数学分支中相关数学模型.ppt

第三篇 数学分支中的相关数学模型 §1 高等数学相关模型 1.1卫星轨道长度 1.2射击命中概率 1.3人口增长率 §2 线性代数相关模型 2.1投入产出综合平衡分析 2.2输电网络 §3 概率统计相关模型 3.1合金强度与碳含量 3.2年龄与运动能力 3.3商品销售量与价格 §1 高等数学相关模型 问题 1.1 卫星轨道长度 人造地球卫星轨道可视为平面上的椭圆. 分析 卫星轨道椭圆的参数方程 椭圆长度 分别是长、短半轴 椭圆积分 无法解析计算 近地点距地球表面439km. 远地点距地球表面2384km. 地球半径6371km. 求该卫星的轨道长度. 输出 MATLAB程序 function y=x5(t) a=8755;b=6810; y=sqrt(a^2*sin(t).^2+b^2*cos(t).^2); clear t=0:pi/10:pi/2; y1=x5(t); L1=4*trapz(t,y1) L2=4*quad(x5,0,pi/2,1e-6) L1=4.908996526785276e+004 L2=4.908996531830460e+004 输出 求解 梯形公式 辛普森公式 评注 问题 1.2 射击命中概率 射击目标为正椭圆形区域,弹着点与中心有随机偏差. 分析 设目标中心x=0,y=0, 无法解析计算 弹着点围绕中心成二维正态分布,偏差在X、Y方向独立. 求炮弹落在椭圆形区域内的概率. 则弹着点(x,y)概率密度函数 炮弹命中椭圆形区域的概率 椭圆在X方向半轴长120m,Y方向半轴长80m. 设弹着点偏差的均方差在X和Y方向均为100m. 求解：蒙特卡罗方法 作变换 以100(m)为1单位，则 MATLAB程序 输出 clear, n=100000; a=1.2;b=0.8;m=0;z=0; for i=1:n x=rand(1,2); y=0; if x(1)^2+x(2)^2=1 y=exp(- 0.5*(a^2* P=0.3752, m=78552 x(1)^2+b^2*x(2)^2)); z=z+y;m=m+1; end end p=(4*a*b/2/pi) *(z /n),m 评注 问题1 1.3 人口增长率 20世纪美国人口数据(106 ), 计算各年份人口增长率. 记时刻t人口为x(t)，则人口相对增长率为 分析 记1900年为k=0 求解：数值微分三点公式 年增长率 2.20 1.66 1.46 1.02 1.04 1.58 1.49 1.16 1.05 1.04 评注 计算程序 问题2 已知某地区20世纪70年代的人口增长率,且1970年人口为210（百万）， 试估计1980年的人口. 记时刻t人口为x(t)，则人口增长满足微分方程 分析 记1970年为k=0 求解 评注 1980年该地区人口为230.2（百万） 数值积分梯形公式 计算程序 clear x=[76 92 106.5 123.2 131.7 150.7 179.3 204 226.5 251.4]; r(1)=(-3*x(1)+4*x(2)-x(3))/20/x(1); for i=2:9 r(i)=(x(i+1)-x(i-1))/20/x(i); end r(10)=(x(8)-4*x(9)+3*x(10))/20/x(10); r=100*r 数值微分计算程序：people_model.m 数值积分计算程序：people_model.m clear,x0=210;t=1970:2:1980; r=[0.87 0.85 0.89 0.91 0.95 1.10]; s=trapz(t,r); x1980=x0*exp(s/100) 为算出瑞士的国土面积，首先对瑞士地图作如下测量：以由西向东方向为x轴，由南到北方向为y轴，选择方便的原点，并将从最西边界点到最东边界点在x轴上的区间适当地划分为若干段，在每个分点的y方向测出南边界点和北边界点的坐标，得到数据（mm）见下页. 习题： 国土面积问题 根据地图比例,18mm相当于40km，试由测量数据计算瑞士国土的近似面积，与它的精确值41288km比较. 计算程序（数据见下页,省略） clear,x=[7.0 … 158];y1=[44 … 68];y2=[44 … 68]; s=trapz(x,y2)-trapz(x,y1); s=s*40^2/18^2, error=s-412