数学建模实验报告。.doc

专业班级 实验报告 姓 名 学 号 实验日期、时间 同组学生姓名 课程名称 数学建模与数学实验 指导老师 实验地点 实验名称 微分方程与线性方程实验 实验类型 其他 实验成绩 内 容： 一、实验目的和要求 二、实验内容和原理 三、主要仪器设备 四、操作方法和实验步骤 五、实验数据记录和处理 六、实验结果与分析 七、讨论、心得  一、实验目的和要求： 1.掌握求解微分方程的方法； 2.熟练掌握矩阵的基本运算； 3.熟练掌握线性方程组的求解。 二、实验内容和原理 1.根据实际问题建立线性方程、微分方程模型。 2.利用Matlab求解方程模型。 三、主要仪器设备 计算机， matlab 四、操作方法和实验步骤 1.求解下列方程组 （1） （2）的解。 （3） 解：(1)作增广矩阵A，在命令窗口输入命令： A=[2 5 -8 8;4 3 -9 9;2 3 -5 7;1 8 -7 12]; rref(A) （2）在命令窗口输入命令： A=[2 -1 3;2 1 1;4 1 2]; rref(A) (3) 作增广矩阵A，在命令窗口输入命令： A=[1 1 -3 -1 1;3 -1 -3 4 4;1 5 -9 -8 0]; rref(A) 2. 求下列微分方程的解析解 （1） （2） （3） 解：（1）在命令窗口输入命令： dsolve(Dy=a*y+b,x) （2）在命令窗口输入命令： y=dsolve(D2y=sin(2*x)-y,y(0)=0,Dy(0)=1,x) （3）在命令窗口输入命令： dsolve((1+x^2)*D2y=2*x*Dy,y(0)=1,Dy(0)=3,x) 3. 求下列微分方程的特解 （1）在初始条件下的特解． （2）在初始条件下的特解. (3) 解：（1）在命令窗口输入命令： y=dsolve(x*Dy+y-exp(x)=0,y(1)=2*exp(1),x) （2）在命令窗口输入命令： y=dsolve((x^2-1)*Dy+2*x*y-cos(x)=0,y(0)=0,x) （3）在命令窗口输入命令： [x,y]=dsolve(Dx=3*x+4*y,Dy=5*x-7*y,x(0)=0,y(0)=1, t); x=simple(x) y=simple(y) 4. 假设每年森林被火灾吞食面积的速度由以下方程决定： （单位：百万亩/年） 这里t年是从现在开始计算的时间．那么在未来3年里，将有多少森林面积被火灾破坏？ （1）模型建立 解：根据题意，每年森林被火灾吞食面积速度由方程决定，现在森林被吞食的面积为0百万亩，设被吞食的面积为y百万亩，即： 其中初始条件为：y(0)=0 (2)模型求解 在命令窗口输入命令： dsolve(Dy=0.25*exp(0.05*t),y(0)=0,t) 得出含未知数的表达式： ans =5*exp(1/20*t)-5 将t=3带入这个表达式: 5*exp(1/20*3)-5 ans = 0.8092 5. 统计学家对当地一家教养所的统计分析得出:,其中x（年）是从当前开始计算的时间．若目前监狱里收容了2000人,现有700间监狱,每间最多可容纳6人． (1)请预测10年后将有多少人被收容? (2)十年后,该监狱能容纳下收容人员吗?若不能容纳,则至少还需要修多少间房间? （1）模型建立 解：根据题意，监狱里收容人员将以每年的速度增加，目前监狱里已经收容了2000人，设监狱里收容了y人： 即： (2)模型求解 解：①在命令窗口输入命令： dsolve(Dy=280*exp(0.2*x),y(0)=2000,x) ans = 1400*exp(1/5*x)+600 然后将x=10带入ans，得： 1400*exp(1/5*10)+600 ans = 1.0945e+004 ②十年后，不能容纳下收容人员，至少还需要修房间： 五、实验数据记录和处理 1、（1）ans =