MATLAB求解方程解析解和数值解.docx

辽宁工程技术大学上机实验报告 成绩 成绩 实验名称 Matlab 求解方程的数值解和解析解 院系 专业 班级 姓名 学号 日期 实验 简述本次实验目的： 1、熟悉MATLAB 软件环境； 目的 2、熟悉MATLAB 的常用运算符； 3、了解MATLAB 的一些常用函数； 你为本次实验做了哪些准备： 实验 提前熟悉线性代数中的方程求解相关运算； 准备 提前熟悉Matlab 中的方程求解相关的命令； 实验 本次共有 4 个练习，完成 4 个。进度 本次实验的收获、体会、经验、问题和教训： 通过本次实验我发现，在 Matlab 中一些算法会变得很简单，有时候并不需要我们去了解具体的程序内部的算法，只要我们学会如何熟练运用Matlab 软件就好。学会如何运用 Matlab 中的算法会对我们研究一些问题带来很大的方便，解决问题会变得很方便，免去了一些手动难以解决的问题。 实验总结 教师评语 用MATLAB 求解质点振动方程 振动是日常生活和工程技术中常见的一种运动形式。利用常系数线性微分方程的理论来讨论有关自由振动和强迫振动的相关问题。利用 MATLAB 数学软件大致可分四类情况：（1）无阻尼自由振动情况；（2）有阻尼自由振动；（3）无阻尼强迫振动；（4）有阻尼强迫振动 求其数值解和解析解； MATLAB 软件求解微分方程解析解的命令“dsolve()” 求通解的命令格式：（’微分方程’,’自变量’） 注：微分方程在输入时，一阶导数y’应输入 Dy，y’’应输入 D2y 等，D 应大写。 1，无阻尼自由振动情况： 常见的数学摆的无阻尼微小振动方程代码如下： t=0:pi/50:2*pi; y=2*sin(3*t+2); plot(t,y,b) 2 1.5 1 0.5 0 -0.5 -1 -1.5 -2 0 1 2 3 4 5 6 7 2，有阻尼自由振动 由无阻尼振动的通解可以看出，无阻尼振动是按照正弦规律运动的，摆动似乎可以无限期的 进行下去，但事实上，空气从在阻力，在运动时，我们必须把空气阻力考虑在内，所以我们得到有阻尼摆动方程为： 记 u/m=2n,g/l=w^2,这里 n,w 是正常数，所以： y=dsolve(D2y+2*n*Dy+w^2*y=0,t); （4.43） 解得: y = C3*exp(-t*(n + ((n + w)*(n - w))^(1/2))) + C2*exp(-t*(n - ((n + w)*(n - w))^(1/2))) （1）小阻尼情形:nw 时，方程（4.43）的通解为： y=exp(-n*t)*(c1*cos(w1*t)+c2*sin(w1*t)) 和前面无阻尼的情形一样，可以把上式的通解改写为一下形式： y=A*exp(-n*t)*sin(w1*t+Q), (4.45) 这里的 A，Q 为任意常数。用 matlab 操作得到： t=0:0.1:10; y=3*exp(-0.1*t).*sin(5*t+4); plot(t,y,k-) 如图： 3 2 1 0 -1 -2 -3 0 1 2 3 4 5 6 7 8 9 10 由（4.45）可见，摆动的运动不是周期的，振动的幅度随着时间的增加而不断减小。 （2）大阻尼情形：你w 时；r2r10; 方程（4.43）的通解： Y=c1*exp(r1*t)+c2*exp(r2*t), 这里的 c1,c2 为任意常数； （3）临界的情形：即，n=w 的情形，方程（4.43）的通解解为 y=exp(-n*t)*(c1+c2*t), 这里的，c1，c2 为任意常数； 由 MAYLAB 绘制图像得： t=0:0.1:100; yy=exp(-0.2*t).*(0.5-0.2*t); plot(t,yy,r-) 0.5 0.4 0.3 0.2 0.1 0 -0.1 -0.2 -0.3  0 10 20 30 40 50 60 70 80 90 100 有阻尼自由振动 无阻尼自由振动和有阻尼自由振动都属于自由振动，它对应于一个二阶常数齐次线性微分方程。当一个振动系统还经常受到一个外力作用时，这种振动称为强迫振动。 ? ＝Ａｓｉｎ（? ｔ＋? ）＋Ｈ／（? ２－ｐ２）ｓｉｎ（ｐｔ） 取Ａ＝２； ＝５；ｐ＝３； t=0:pi/50:2*pi; y=2*sin(5*t+2)+1/9*sin(3*t); plot(t,y,k) 2.5 2 1.5 1 0.5 0 -0.5 -1 -1.5 -2 -2.5 0 1 2 3 4 5 6 7 有阻尼强迫振动 这时摆动的运动方程（1.11）变为： D2y+2*n*Dy+w.^2*y=H*sin(p*t), （4.52） 根据实际情况，我们只讨论小阻尼的情况，即 nw 的情形，这时的（4.52）对应的齐次线性方程的通解为