matlab常微分方程2次初值问题,MATLAB求解二阶常微分方程初值问题
主程序如下:
>> global m1 m2 k1 c1 c2 v0
>> m1=4080;m2=12.5e3;k1=7.35e6;c1=17.56e3;c2=0.225e3;v0=3.2/60
>> [t,y]=ode45(@doublemass,[0 5],[-0.0154;-0.0049;0;0])
其中的doublemass的函数M文件定义如下:
function dyy=doublemass(t,y)
global m1 m2 k1 c1 c2 v0
dyy=[y(3);(-4.48e6*v0*t/(6-y(1))-(k1+4.48e6/(6-y(1)))*y(1)-(c1+c2)*y(2)+c2*y(4)+4.48e6*y(2)/(6-y(1)))/m1;y(4);(4.48e6*v0*t/(6-y(1))-m2*10+c2*y(3)+4.48e6*y(1)/(6-y(1))-c2*y(4)-4.48e6*y(2)/(6-y(1)))/m2]
但计算出的结果却不对,前面两个初值不是前面设置好的-0.0154;-0.0049,而是 -0.0000 -0.0000,如下图:
y =
1.0e+009 *
-0.0000 -0.0000 0 0
-0.0000 -0.0000 -0.0000 -0.0000
-0.0000 -0.0000 -0.0000 -0.0000
-0.0000 -0.0000 -0.0000 -0.0000
-0.0000 -0.0000 -0.0000 -0.0000
-0.0000 -0.0000 -0.0000 -0.0000
-0.0000 -0.0000 -0.0000 -0.0000
-0.0000 -0.0000 -0.0000 -0.0000
-0.0000 -0.0000 -0.0000 -0.0000
-0.0000 0.0000 -0.0000 -0.0000
-0.0000 0.0000 -0.0000 -0.0000
-0.0000 0.0000 -0.0000 -0.0000
-0.0000 0.0000 -0.0000 -0.0000
-0.0000 0.0000 -0.0000 -0.0000
-0.0000 0.0000 -0.0000 -0.0000
-0.0000 0.0000 -0.0000 -0.0000
-0.0000 0.0000 -0.0000 -0.0000
-0.0000 0.0000 -0.0000 -0.0000
-0.0000 0.0000 -0.0000 -0.0000
-0.0000 0.0000 -0.0000 -0.0000
-0.0000 0.0000 -0.0000 -0.0000
-0.0000 0.0000 -0.0000 -0.0000
-0.0000 0.0000 -0.0000 -0.0000
-0.0000 0.0000 -0.0000 -0.0000
-0.0000 0.0000 -0.0000 -0.0000
-0.0000 0.0000 -0.0000 -0.0000
-0.0000 0.0000 -0.0000 -0.0000
-0.0000 0.0000 -0.0000 -0.0000
-0.0000 0.0000 -0.0000 -0.0000
-0.0000 0.0000 -0.0000 -0.0000
-0.0000 0.0000 -0.0000 -0.0000
-0.0000 0.0000 -0.0000 -0.0000
-0.0000 0.0000 -0.0000 -0.0000
-0.0000 0.0000 -0.0000 -0.0000
-0.0000 0.0000 -0.0000 -0.0000
-0.0000 0.0000 -0.0000 -0.0000
-0.0000 0.0000 -0.0000 -0.0000
-0.0000 0.0000 -0.0000 -0.0000
-0.0000 0.0000 -0.0000 -0.0000
-0.0000 0.0000 -0.0000 -0.0000
-0.0000 0.0000 -0.0000 -0.0000
-0.0000 0.0000 -0.0000 -0.0000
-0.0000 0.0000 -0.0000 -0.0000
-0.0000 0.0000 -0.0000 -0.0000
-0.0000 0.0000 -0.0000 -0.0000
-0.0000 0.0000 -0.0000 -0.0000
-0.0000 0.0000 -0.0000 -0.0000
-0.0000 0.0000 -0.0000 -0.0000
-0.0000 0.0000 -0.0000 -0.0000
-0.0000 0.0000 -0.0000 -0.0000
-0.0000 0.0000 -0.0000 -0.0000
-0.0000 0.0000 -0.0000 -0.0000
-0.0000 0.0000 -0.0000 -0.0000
-0.0000 0.0000 -0.0000 -0.0000
-0.0000 0.0000 -0.0000 -0.0000
-0.0000 0.0000 -0.0000 -0.0000
-0.0000 0.0000 -0.0000 -0.0000
-0.0000 0.0000 -0.0000 -0.0000
-0.0000 0.0000 -0.0000 -0.0000
-0.0000 0.0000 -0.0000 -0.0000
-0.0000 0.0000 -0.0000 -0.0000
-0.0000 0.0001 -0.0000 -0.0000
-0.0000 0.0001 -0.0000 -0.0000
-0.0000 0.0001 -0.0000 -0.0000
-0.0000 0.0001 -0.0000 -0.0000
-0.0000 0.0002 -0.0000 -0.0001
-0.0000 0.0002 -0.0000 -0.0001
-0.0000 0.0002 -0.0000 -0.0001
-0.0000 0.0003 -0.0000 -0.0001
-0.0000 0.0003 -0.0000 -0.0001
-0.0000 0.0004 -0.0000 -0.0001
-0.0000 0.0005 -0.0000 -0.0002
-0.0000 0.0005 -0.0000 -0.0002
-0.0000 0.0006 -0.0000 -0.0002
-0.0000 0.0006 -0.0000 -0.0002
-0.0000 0.0007 -0.0000 -0.0002
-0.0000 0.0007 -0.0000 -0.0002
-0.0000 0.0007 -0.0000 -0.0003
-0.0000 0.0008 -0.0000 -0.0003
-0.0000 0.0008 -0.0000 -0.0003
-0.0000 0.0008 -0.0000 -0.0003
-0.0000 0.0008 -0.0000 -0.0003
-0.0000 0.0008 -0.0000 -0.0003
-0.0000 0.0008 -0.0000 -0.0003
-0.0000 0.0008 -0.0000 -0.0003
-0.0000 0.0008 -0.0000 -0.0003
-0.0000 0.0008 -0.0000 -0.0003
-0.0000 0.0008 -0.0000 -0.0003
-0.0000 0.0008 -0.0000 -0.0003
-0.0000 0.0008 -0.0000 -0.0003
-0.0000 0.0008 -0.0000 -0.0003
-0.0000 0.0008 -0.0000 -0.0003
-0.0000 0.0008 -0.0000 -0.0003
-0.0000 0.0010 -0.0001 -0.0003
-0.0000 0.0013 -0.0001 -0.0003
-0.0000 0.0017 -0.0001 -0.0004
-0.0000 0.0022 -0.0001 -0.0004
-0.0000 0.0030 -0.0001 -0.0004
-0.0000 0.0041 -0.0001 -0.0004
-0.0000 0.0054 -0.0001 -0.0004
-0.0000 0.0069 -0.0002 -0.0004
-0.0000 0.0094 -0.0002 -0.0004
-0.0001 0.0123 -0.0002 -0.0004
-0.0001 0.0157 -0.0002 -0.0004
-0.0001 0.0197 -0.0003 -0.0004
-0.0001 0.0261 -0.0003 -0.0004
-0.0001 0.0335 -0.0003 -0.0004
-0.0002 0.0420 -0.0003 -0.0004
-0.0002 0.0516 -0.0004 -0.0004
-0.0002 0.0678 -0.0004 -0.0004
-0.0003 0.0864 -0.0005 -0.0004
-0.0004 0.1076 -0.0005 -0.0004
-0.0004 0.1314 -0.0006 -0.0005
-0.0005 0.1585 -0.0007 -0.0005
-0.0006 0.1885 -0.0007 -0.0005
-0.0007 0.2214 -0.0008 -0.0005
-0.0008 0.2575 -0.0008 -0.0005
-0.0009 0.2967 -0.0009 -0.0005
-0.0010 0.3392 -0.0010 -0.0006
-0.0011 0.3851 -0.0010 -0.0006
-0.0013 0.4345 -0.0011 -0.0006
-0.0014 0.4875 -0.0012 -0.0006
-0.0016 0.5441 -0.0013 -0.0006
-0.0017 0.6045 -0.0013 -0.0006
-0.0019 0.6688 -0.0014 -0.0007
-0.0021 0.7372 -0.0015 -0.0007
-0.0023 0.8096 -0.0016 -0.0007
-0.0025 0.8863 -0.0017 -0.0007
-0.0027 0.9673 -0.0018 -0.0007
-0.0029 1.0528 -0.0019 -0.0007
-0.0032 1.1428 -0.0020 -0.0008
-0.0034 1.2376 -0.0021 -0.0008
-0.0037 1.3371 -0.0021 -0.0008
-0.0040 1.4416 -0.0023 -0.0008
-0.0043 1.5512 -0.0024 -0.0008
-0.0046 1.6659 -0.0025 -0.0009
-0.0049 1.7859 -0.0026 -0.0009
-0.0051 1.8795 -0.0027 -0.0009
-0.0054 1.9762 -0.0027 -0.0009
-0.0056 2.0761 -0.0028 -0.0009
-0.0059 2.1792 -0.0029 -0.0009
不知是何原因?
matlab常微分方程2次初值问题,MATLAB求解二阶常微分方程初值问题相关推荐
- matlab编写求解二阶常微分方程,求一个复杂的二阶常微分方程的数值解,不会写代码= =...
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼 如下是我写的计算步骤,稍后解释:(不想看可以跳过~) %main %input some parameters m1=input('please inpu ...
- python解常微分方程龙格库_求解二阶常微分方程的RungeKutta四阶方法
我试着做一个简谐振子的例子,它将用龙格-库塔四阶法求解.要求解的二阶常微分方程(ODE)和初始条件为: y''+y=0 y(0)=0和y'(0)=1/pi 范围在0到1之间,共有100步.我用u作为辅 ...
- python求解四阶微分方程_用Python求解二阶常微分方程组的RungeKutta四阶解
本文试图用龙格库塔四阶法数值求解两个常微分方程组. 初始系统: 要解决的系统: 我有非常奇怪的解图... 我有: 正确的图形: 我在我的龙格库塔找不到麻烦.请帮帮我.在 我的代码在这里:dt = 0. ...
- 简单典型二阶系统_MIT—微分方程与线性代数笔记2.1 二阶常微分方程
§2.1 二阶常微分方程 2.1 Second Order Equations 微分方程 MIT公开课<微分方程和线性代数> 2.1 二阶微分方程v.youku.com 二阶常微分方程的 ...
- 一阶欧拉近似matlab,MATLAB改进欧拉法与四阶龙格-库塔求解一阶常微分方程.doc
MATLAB改进欧拉法与四阶龙格-库塔求解一阶常微分方程 姓名:樊元君 学号:2012200902 日期:2012.11.06 一.实验目的 掌握改进欧拉法与四阶龙格-库塔求解一阶常微分方程的初值问题 ...
- 用四阶RungeKutta方程解二阶常微分方程,并计算船舶在规则波中的横摇角(附Matlab代码)
前几天接到师姐分派的任务,让我求解一艘船模的横摇角的时间历程曲线,为后期的减摇控制准备. 1 首先冷静分析一下,原方程如下: 我们要求解的就是theta角和时间t之间的关系曲线,这是一道典型的二阶常微 ...
- matlab求二阶微分方程的通解,二阶常微分方程matlab的数值解和解析解分析总报告.ppt...
方程求解程序清单 a=-1,b=2,c=-1; w=1; m=2; n=1; h = 0.02; t=0:h:30; s1=dsolve('a*D2y+b*Dy+c*y=sin(w*t)','y(0) ...
- matlab微分方程组边值,matlab求解常微分方程边值问题的方法
matlab求解常微分方程边值问题的方法 Matlab 求解常微分方程边值问题的方法:bvp4c 函数 常微分方程的边值问题,即 boundary value problems ,简称 BVP 问题, ...
- 二阶边值问题的数值解matlab,《二阶常微分方程边值问题的数值解法》-毕业论文.doc...
w 摘 要 本文主要研究二阶常微分方程边值问题的数值解法.对线性边值问题,我们总结了两类常用的数值方法,即打靶法和有限差分方法,对每种方法都列出了详细的计算步骤和Matlab程序代码,通过具体的算例对 ...
- 四阶龙格库塔方程(Rungekutta)解二阶常微分方程组并计算船舶在迎浪下的纵摇埀荡耦合运动方程-附Matlab代码
今年年初的时候给师姐做了DDPG算法的船舶减横摇控制算法,师姐还有想法要让我把纵摇-埀荡两个自由度的减摇也做出来,这个任务归我了.实际上不管是多少个自由度的减摇,其实都需要解运动方程,当初做单自由度横 ...
最新文章
- 世界上有哪些代码量很少,但很牛逼很经典的算法或项目案例?
- Facebook曝至今最严重安全漏洞,超过5000万用户受影响
- LintCode-56.两数之和
- 《数据安全法》今日实施,中国信通院联合百度等企业发起“数据安全推进计划”
- 数据科学竞赛-人脸表情识别
- 诺奖得主:中国机制促成抗疫成功经济复苏
- SpringCloud Ribbon中的7种负载均衡策略!
- css 动态生成圆形区域内扇形个数_CSS实用技巧总结
- Ubuntu开启允许root用户远程登录
- 虚拟机安装centos7后无法联网
- android背景颜色动态修改,Android自定义TextView带圆角及背景颜色(动态改变圆角背景颜色)...
- 《A Sub-Pixel Edge Detector: an Implementation of the Canny/Devernay Algorithm》
- 第三课 protel学习系列——protel功能拓展
- chrome 内核的浏览器pdf 字体显示问题
- 爬取天涯帖子(单个帖)
- 电脑版微信打开内置浏览器
- 原神换服报错(二级地址解析失败),换服失败,无反应
- 01、如何用DAP仿真器下载程序
- Trajectory Similarity Join in Spatial Networks
- 微信分享解决wx not defined