用龙格库塔算法解二阶常微分方程,利用c++编程 另z=y'原来的方程就可以化简成以z'和y,x的方程和y'=z 带入两个初始条件,就可以进行迭代了用matlab编程实现四阶龙格库塔解二元二阶微分方程组 求解二阶微分方程,初始条件还需要给出y1'(0)和y2'(0)。这里暂时按照0处理。function zd530003514a=0.1;b=0.1;Y0=[b-1;0;b;0];解方程[t,Y]=ode45(@ode,[0 10],Y0);y1=Y(:,1);y2=Y(:,3);绘图subplot 211plot(t,y1);subplot 212plot(t,y2);微分方程定义function dY=ode(t,Y)L1=5;L2=0.01;a0=2;b0=2;c0=2;y1=Y(1);y2=Y(3);dY=Y(2);(a0*y2+b0*y2^2+c0*y2^3)-L1^2*L2*y1-L1^2*y1;Y(4);(a0*y2+b0*y2^2+c0*y2^3)-L1^2*L2*y1;用四阶龙格-库塔求微分方程 fun=(s,f)10./(s.^4+8*s.^3+36*s.^2+40*s+10);s0=0;send=1;ds=0.001;s=s0:ds:send;f=0;初值i=1;for ss=s(1:end-1)k1=fun(ss,f(i));k2=fun(ss+ds/2,f(i)+ds/2*k1);k3=fun(ss+ds/2,f(i)+ds/2*k2);k4=fun(ss+ds,f(i)+ds*k3);f(i+1)=f(i)+ds/6*(k1+2*k2+2*k3+k4);i=i+1;endplot(s,f);如何用四阶龙格库塔法ode45求取常微分方程的参数 我通过查文献知道这里们要用到四阶龙格库塔法进行ode45函数法拟合,即求取同时满足。可以给定初值的!数据可以变成这样t=[0 7 32 57 82 107 187 307 467 607];x=[0.9748 0.0105 0.0108 0.0039;0.8088 0.1032 0.0143 0.0737;0.2481 0.1258 0.0941 0.5103;0.0771 0.1264 0.1310 0.6655;0.0307 0.0931 0.1395 0.7367;0.0000 0.0714 0.1155 0.8131;0.0000 0.0599 0.0986 0.8415;0.0000 0.0588 0.0914 0.8498;0.0000 0.0506 0.0715 0.8798];x=[0.9919 0.0043 0.0029 0.0008];跪求四阶龙格库塔公式求 常微分方程的方法 不要代码 只要具体的推倒或者是 解出的方程组 http://baike.baidu.com/link?url=dTC8pz6PJXWpZdsua4rFN_panquCZBoHgwqPNwdvVqgZL_aAEXaCy7vebNX-Sk_vOX9qU5NZUcUtmbnD8hTsA3EjhlRPwwk-rPmjGGne0mNuaYTurUhqP90-jM7DMiVXfBYS1GBXwLd4HVmbBri7UkalIJIcXS9318VED0fqDXsc02gB_qfKXnpwaxfKcyUt就有推导啊指出楼主的一个认识错误 龙格库塔法仅仅是数值求解 并不能说是解出了方程组 数值解再怎么说也只是解析解的近似描述(个人认为 仅供参考)龙格-库塔方法求解三阶常微分方程 第一步:将高阶常微分方程转换成常微分方程组,func(t,x)第二步:调用runge_kutta(@func,y0,h,a,b)例如:二阶常微分方程func。mfunction z=func(t,y)z=[y(2);(1-y(1)^2)*y(2)-y(1)];main。mclear all;close all;clcy0=[0.25;0];h=0.1;a=0;b=20;[t1 y1]=runge_kutta(@rhs_7,y0,h,a,b)
随机阅读
- 乙氧基甲叉氰乙酸乙蜡 反流性食管炎需要吃什么药
- 站着气割怎样走枪 请问哪位气割工师傅,气割枪氧气阀门坏了关不严,一直跑氧气,这时怎样才能点着火?
- 少年的心儿永不老什么意思 歌词里有一句“阿里山的姑娘美如水啊”,这首歌叫什么名字,谁唱的?
- 兴安盟乌兰浩特集散点能送到跃进马场吗? 乌兰浩特到跃进马场找人
- wta 皇冠明珠 WTA皇冠明珠赛的比赛积分
- 盛龙路口离龙园远不 赤湾西龙岗坪山沙湖
- 武汉电动车上牌的标准是什么? 汉阳王派电动车电话
- 福建省福州市连江县东岱镇龙山村龙兴南路邮编是什么? 连江东岱镇东水村
- 莫道桑榆晚为霞尚满天启示 莫道桑榆晚 为霞尚满天是什么意思
- 去商丘旅游有什么值得一尝的特色美食? 郑州纬四路家和米线
- 桂林市农田灌溉试验中心站 赣粤运河和湘桂运河可行性如何?
- 谁能告诉我东莞东城哪家海鲜酒楼既好吃价格又合理呢? 东莞食为先海鲜大酒楼
- 最近听到一个词,叫 房价对家庭收入比
- 儿童节我们通常干什么用英文来表达 儿童节人们通常干什么 用英语
- 龙之崛起沙漠建不了水井 《皇帝-龙之崛起》关于圣地的问题如何解决?
- 洗车污水净化的原理 洗车污水处理方法
- 求张艺谋的《我的父亲母亲》的迅雷下载地址?请原创回答 他母亲的房子 迅雷
- 首长好?近到乡下:听说湟源县波航乡纳隆村把2014年的贫困户底保金轮流着发放,不管贫富,其中有位三组的百头牛养殖户杨先生还在领底保金,这些乡村干部是无做为还是党的扶贫政策当儿戏? 百里杜鹃普底乡永丰村九牛组李光后
- 系噶嘻嘻一觉觉啦啦来次够抖音音乐 前几天在抖音听了一首节奏感很强的纯音乐,歌词就是“噔噔 噔 噔噔噔” 记住这歌词,最后三个一个节拍
- 四海一家综合缴费软件,一站通便民综合缴费系统.是不是真的象他们说的那样里面的功能都全还是就是个骗人