![用MATLAB求解微分方程.ppt_第1页](http://file.renrendoc.com/FileRoot1/2019-11/13/0a7d9605-cb87-4c7d-8b7d-6fd57ca745c4/0a7d9605-cb87-4c7d-8b7d-6fd57ca745c41.gif)
![用MATLAB求解微分方程.ppt_第2页](http://file.renrendoc.com/FileRoot1/2019-11/13/0a7d9605-cb87-4c7d-8b7d-6fd57ca745c4/0a7d9605-cb87-4c7d-8b7d-6fd57ca745c42.gif)
![用MATLAB求解微分方程.ppt_第3页](http://file.renrendoc.com/FileRoot1/2019-11/13/0a7d9605-cb87-4c7d-8b7d-6fd57ca745c4/0a7d9605-cb87-4c7d-8b7d-6fd57ca745c43.gif)
![用MATLAB求解微分方程.ppt_第4页](http://file.renrendoc.com/FileRoot1/2019-11/13/0a7d9605-cb87-4c7d-8b7d-6fd57ca745c4/0a7d9605-cb87-4c7d-8b7d-6fd57ca745c44.gif)
![用MATLAB求解微分方程.ppt_第5页](http://file.renrendoc.com/FileRoot1/2019-11/13/0a7d9605-cb87-4c7d-8b7d-6fd57ca745c4/0a7d9605-cb87-4c7d-8b7d-6fd57ca745c45.gif)
已阅读5页,还剩8页未读, 继续免费阅读
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1 微分方程的解析解 求微分方程 组 的解析解命令 dsolve 方程1 方程2 方程n 初始条件 自变量 结果 u tan t c 用MATLAB求解微分方程 解输入命令 dsolve Du 1 u 2 t 解输入命令 y dsolve D2y 4 Dy 29 y 0 y 0 0 Dy 0 15 x 结果为 y 3e 2xsin 5x 解输入命令 x y z dsolve Dx 2 x 3 y 3 z Dy 4 x 5 y 3 z Dz 4 x 4 y 2 z t x simple x 将x化简y simple y z simple z 结果为 x c1 c2 c3 c2e 3t c3e 3t e2ty c1e 4t c2e 4t c2e 3t c3e 3t c1 c2 c3 e2tz c1e 4t c2e 4t c1 c2 c3 e2t 2 用Matlab求常微分方程的数值解 t x solver f ts x0 options 1 在解n个未知函数的方程组时 x0和x均为n维向量 m 文件中的待解方程组应以x的分量形式写成 2 使用Matlab软件求数值解时 高阶微分方程必须等价地变换成一阶微分方程组 注意 解 令y1 x y2 y1 1 建立m 文件vdp1000 m如下 functiondy vdp1000 t y dy zeros 2 1 dy 1 y 2 dy 2 1000 1 y 1 2 y 2 y 1 2 取t0 0 tf 3000 输入命令 T Y ode15s vdp1000 03000 20 plot T Y 1 3 结果如图 解1 建立m 文件rigid m如下 functiondy rigid t y dy zeros 3 1 dy 1 y 2 y 3 dy 2 y 1 y 3 dy 3 0 51 y 1 y 2 2 取t0 0 tf 12 输入命令 T Y ode45 rigid 012 011 plot T Y 1 T Y 2 T Y 3 3 结果如图 图中 y1的图形为实线 y2的图形为 线 y3的图形为 线 导弹追踪问题 设位于坐标原点的甲舰向位于x轴上点A 1 0 处的乙舰发射导弹 导弹头始终对准乙舰 如果乙舰以最大的速度v0 是常数 沿平行于y轴的直线行驶 导弹的速度是5v0 求导弹运行的曲线方程 又乙舰行驶多远时 导弹将它击中 解法一 解析法 由 1 2 消去t整理得模型 解法二 数值解 1 建立m 文件eq1 mfunctiondy eq1 x y dy zeros 2 1 dy 1 y 2 dy 2 1 5 sqrt 1 y 1 2 1 x 2 取x0 0 xf 0 9999 建立主程序ff6 m如下 x0 0 xf 0 9999 x y ode15s eq1 x0 xf 00 plot x y 1 b holdony 0 0 01 2 plot 1 y b 结论 导弹大致在 1 0 2 处击中乙舰 令y1 y y2 y1 将方程 3 化为一阶微分方程组 解法三 建立参数方程求数值解 设时刻t乙舰的坐标为 X t Y t 导弹的坐标为 x t y t 3 因乙舰以速度v0沿直线x 1运动 设v0 1 则w 5 X 1 Y t 4 解导弹运动轨迹的参数方程 建立m 文件eq2 m如下 functiondy eq2 t y dy zeros 2 1 dy 1 5 1 y 1 sqrt 1 y 1 2 t y 2 2 dy 2 5 t y 2 sqrt 1 y 1 2 t y 2 2 取t0 0 tf 2 建
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 一年级下数学教案-期末复习(2)-苏教版
- 五年级上册英语-Module 9 Unit 2 I feel happy.外研版(三起)
- 习作:《奇妙的想象》第一课时(教案)部编版三年级语文下册
- 2024车辆赔偿协议书
- 面积 复习课(教案)-三年级下册数学人教版
- 乐器采购合同范本(2024版)
- 2024苏泊尔用户服务协议
- 家长会教师讲话稿
- 分析新能源电动汽车充电设施的建设与运营模式
- 高中毕业典礼致辞
- 分子生物学实验的常见问题与解决方案
- 2023年中国电信校园招聘笔试参考题库附带答案详解
- 水沸腾时气泡的变化
- pa800plus生物制药分析系统试剂盒说明和应用报告
- 加油站油气回收治理设施日常管理制度及检查维护台账(模板)
- 我国民事诉讼中调解制度研析的开题报告
- 钢铁公司安全生产事故隐患排查、治理管理制度
- 西学中中医基础理论试题含答案含四个季度
- 2023-2024学年广西壮族自治区柳州市小学语文四年级期末深度自测试卷
- 国家开放大学电大《监督学》论述题题库及答案(试卷号2528)
- GB/T 2550-2016气体焊接设备焊接、切割和类似作业用橡胶软管
评论
0/150
提交评论