数学建模实验

1、数学建模实验数学建模课程实验报告专题实验7班级数财系班学号3姓名李丛文实验题目常微分方程数值解实验目的1掌握用matl求微分方程初值问题数值解的方法；.通过实例学习微分方程模型解决简化的实际问题；3.了解欧拉方法和龙格库塔方法的基本思想。实验内容(包括分析过程、方法、和代码,结果）. 用欧拉方法和龙格库塔方法求下列微分方程初值问题的数值解,画出解的图形，对结果进行分析比较解;m文件funci f=f（x，)f=+;程序；c;clar;0;b=; %求解区间x1,yr=45(f,b,1); %调用龙格库塔求解函数求解数值解;% 以下利用eue方法求解y(1）=；n0;=（b-a)/n;x=a：:

2、b;r =:n y(i1)y（i)+h*f(x(i)，y();dfgure（1）plot(x1,y_r,*,x，y,b+,x,*exp()x2,-);%数值解与真解图itl(数值解与真解图);eend(r,euer,真解);xlal(x);ylae(y);igure(2)ot(x1,b(y_r(3*exp（1)-2*x1-2),-);%龙格库塔方法的误差ttle(龙格库塔方法的误差）lbel(）;ylae(error)；igure(3)plo(x,abs(-(3*ex(x)-2x-)）,r)%euler方法的误差it（euer方法的误差)xlae(x);ylabe(eror);4. 单摆运动是

3、一个我们熟悉的物理模型，可以看作工程技术中一些振动问题的简化,图中一根长l的(无弹性的)细线，一端固定,另一端悬挂一质量为m的小球，在重力作用下小球处于竖直的平衡位置，使小球偏离平衡位置一根小的角度,然后让它无初速度的放开,小球就会沿圆弧摆动,在不考虑空气阻力的情况下建立关于时间的微分方程，设=25cm，在等于和两种情况下求方程的数值解,并与近似解比较。解;实验原理与数学模型:在小球摆动过程中的任一位置,小球所受重力沿运动轨迹方向的分力为-mgi(负号表示力的方向与的正方向相反）,利用牛顿第二定律即得微分方程描述单摆运动规律的微分方程（)是2阶微分方程,无解析解，但可用mtlb或其它软件编程求

4、其数值解，但都需要先将它化成方程组的形式。令 则微分方程 (1) 化为 初始条件转化为 在前面的两式中，g=98，l.25,x10为1o=0.74（弧度）及30o=0536(弧度)两种情况. 周期 根据上原理，可以建立模型,用matlab编程求解，过程如下,程序 以单摆的两个周期来计算作图。文件；function dx=abi(t,)9.8；=0.;dxx()；-g/lsin（x(1);程序；一.当时0:05:；a0=.174;x0=a0,;t,=ode23(dnbai,s,x）;y=a0*cos(sqr(4).*t）；，x(：，）,splot(1,2,2）,pot (t,x(:,）,k)，t

5、le(摆角10度数值解)subplt（，2,1）,lot（t,y，b*）,lot（t,，-*）,titl(摆角10度近似解）ans = 0.7 0145 .0500 0661 0.168 0.10 .141 0.1 1500 0.13 0.1017 0200 0.0550 0.052 .250 0.0015 -.018 0.3000 -.022 0.056 0350 -0.1009 -0.1046 0.000 -01397 -.1429 0.50 -.169 -0.69 05000 -0.174 -0.145 0.50 .1667 -01647 .6000 -1430 -0.186 .650

6、-.4 -0.098 0.7000 -.57 -0.91 07 -00044 0.04 0.800 0094 0.94 0.500 983 175 .90 01377 0.1449 0.500637 01691.0000 .1738 0174二. 当时,=0:0.5：2;a0=0.523;x0=a，0;t,=oe2（danai,s,x0)；y=a*cos(qrt(40).*t)；t,x(:,）,suplot(，2,2),plt (t，x(:,1）,-k*),titl(摆角30度数值解）supot(1，2,1), plt（t,y,b*),plot(，y，*)，title(摆角3度近似解)ans

7、= 0 .5236 0.5236 0050 0.4993 0.9 0.100 0283 04223 0.150 0.310 0.3051 0.2000 0.17 .177 0.2500 0.0 -0.005 0.00 .141 -0.168 0350 -0.289 0.319 0.4000 -.407 .42 0.4500 -0.4874 -0500 0.5000 -05217 -05 0.5500 -0.05 -0.942 0600 -0.4460 -0158 0.00 -0.2 0.963 07000 0.206 -147 07500 -.00 0016 .800 10 .1782 0.8

8、500 0202 .325 0.0 0.385 0.448 0.9500 4736 0039 1.000 0.5177 0522 1.050 05136 0.45 1.1000 0.4617 0.4092 1.1500 .366 .3 .0 .265 .169 .200 083 .0271 .30 -0.70 -0.188 3500 -.33 033 1400 -0.1 -0.4407 .500 -0.480 -08 1.5000 0.5116 -0.526 .500 0.5177-4 1.600 -0.475 0.4023 1.650 -0389 -0.2782 00 -0.2654 -0.124 .7500 -1164 007 .800 0.0437 0.1984 1.500 0.995 033 1.9000 0.331 065 1.950 0405 0.5094 2000 0036 5218从数据