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

辽宁工程技术大学上机实验报告 成绩 实验名称 Matlab 求解方程的数值解和解析解 院系 专业 班级 姓名 学号 日期 简述本次实验目的： 实验 1、熟悉 MATLAB软件环境； 目的 2、熟悉 MATLAB的常用运算符； 3、了解 MATLAB的一些常用函数； 你为本次实验做了哪些准备： 实验 提前熟悉线性代数中的方程求解相关运算； 准备 提前熟悉 Matlab 中的方程求解相关的命令； 实验 本次共有 4 个练习，完成 4 个。 进度 本次实验的收获、体会、经验、问题和教训： 通过本次实验我发现，在 Matlab 中一些算法会变得很简单，有时候并不需要我 们去了解具体的程序内部的算法，只要我们学会如何熟练运用 Matlab 软件就好。 学会如何运用 Matlab 中的算法会对我们研究一些问题带来很大的方便，解决问 题会变得很方便，免去了一些手动难以解决的问题。 实验 总结 教师 评语 用 MATLAB求解质点振动方程 振动是日常生活和工程技术中常见的一种运动形式。利用常系数线性微分方程的理论来 讨论有关自由振动和强迫振动的相关问题。 利用 MA TLAB 数学软件大致可分四类情况： （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 ）的通解