版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

brqhl_ing

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 262.9
散金: 11
红花: 1
帖子: 115
在线: 19.5小时
虫号: 854693
注册: 2009-09-23
性别: GG
专业: 物理电子学

[交流] 【求助】解常微分方程组已有2人参与

dx(1)=-8/3*x(1)+x(2)*x(3)+0.5*(x(1)-x(2));
dx(2)= -10*x(2)+10*x(3)
dx(3)=-x(1)*x(2)+28*x(2)-x(3)];

dx(4)=-8/3*x(4)+x(5)*x(6)
dx(5)=-10*x(5)+10*x(6)
dx(6)=-x(3)*x(4)+28*x(4)-x(6)
这两个其实就是两个系统的耦合如何去解这个方程组？

回复此楼

» 猜你喜欢

孩子确诊有中度注意力缺陷已经有12人回复
2025冷门绝学什么时候出结果已经有3人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有4人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复
AI论文写作工具：是科研加速器还是学术作弊器？已经有3人回复
2026博士申请-功能高分子，水凝胶方向已经有6人回复
论文投稿，期刊推荐已经有4人回复
硕士和导师闹得不愉快已经有13人回复
请问2026国家基金面上项目会启动申2停1吗已经有5人回复
同一篇文章，用不同账号投稿对编辑决定是否送审有没有影响？已经有3人回复

» 本主题相关价值贴推荐，对您同样有帮助:

求助四波混频耦合波方程数值求解已经有5人回复
请问如何求解二元一阶微分方程组已经有24人回复
高金求助matlab解微分方程组已经有12人回复
【高奖励】消灭5月份0回复已经有20人回复
【求助】请教非齐次常微分方程组的解析解法已经有4人回复
【求助】如何求解下列偏微分方程已经有13人回复
【求助】急请微分方程高手看下这个方程组！已经有14人回复
【求助】偏微分方程和常微分方程组成的方程组怎么求，大家帮帮忙啊已经有3人回复
【求助】matlab怎么求解偏微分方程组啊，先谢谢了已经有13人回复
【求助】如何用Runge-Kutta迭代求解二阶常微分方程组【已解决】已经有9人回复

1楼2010-06-04 17:06:38

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

brqhl_ing

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 262.9
散金: 11
红花: 1
帖子: 115
在线: 19.5小时
虫号: 854693
注册: 2009-09-23
性别: GG
专业: 物理电子学

引用回帖:

Originally posted by brqhl_ing at 2010-06-04 17:06:38:
dx(1)=-8/3*x(1)+x(2)*x(3)+0.5*(x(1)-x(2));
dx(2)= -10*x(2)+10*x(3)
dx(3)=-x(1)*x(2)+28*x(2)-x(3)];

dx(4)=-8/3*x(4)+x(5)*x(6)
dx(5)=-10*x(5)+10*x(6)
dx(6)=-x(3)*x(4)+28*x(4)-x(6)
这两个其实就 ...

PS：初值可以任意给定

赞一下

回复此楼

2楼2010-06-04 17:07:41

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

change0618

铁杆木虫 (著名写手)

方丈大师

应助: 44 (小学生)
金币: 17724.5
红花: 17
帖子: 2413
在线: 546.7小时
虫号: 496517
注册: 2008-01-19
专业: 化学反应工程

★ ★
robert2020(金币+2):多谢提供解决方案，欢迎常来仿真模拟版交流！ 2010-06-04 23:18:29
zzuwangshilei(金币+2):帮助他人解决问题 2010-06-05 10:46:11
adu886886(金币-2):金币加重复了，不好意思 2010-06-05 18:34:47

function test
clear
clc
x0=eps*ones(1,3);
[t,x]=ode45(@OdeFun,[0:0.02:8],x0);
plot(t,x(:,1),'r.-',t,x(:,2),'b.-',t,x(:,3),'k.-')

function dx = OdeFun(t,x)
dx1=-8/3*x(1)+x(2)*x(3)+0.5*(x(1)-x(2));
dx2= -10*x(2)+10*x(3);
dx3=-x(1)*x(2)+28*x(2)-x(3);
dx=[dx1;dx2;dx3];

赞一下(3人)

回复此楼

3楼2010-06-04 23:09:48

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

brqhl_ing

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 262.9
散金: 11
红花: 1
帖子: 115
在线: 19.5小时
虫号: 854693
注册: 2009-09-23
性别: GG
专业: 物理电子学

zzuwangshilei:鼓励继续讨论 2010-06-05 10:46:36

引用回帖:

Originally posted by change0618 at 2010-06-04 23:09:48:
function test
clear
clc
x0=eps*ones(1,3);
[t,x]=ode45(@OdeFun,[0:0.02:8],x0);
plot(t,x(:,1),'r.-',t,x(:,2),'b.-',t,x(:,3),'k.-')

function dx = OdeFun(t,x)
dx1=-8/3*x(1)+x(2)*x(3)+0.5*(x(1) ...

谢谢，你的帮助！我现在想对第一个方程和第四个方程进行希尔伯特变换，来定义相位该如何进行？是不是把解出的x y z 的值按照他的时间序列代入到第一个和第四个方程把所选时间序列的dx1 和dx4
的值求出来在进行变换？？

赞一下(1人)

回复此楼

4楼2010-06-05 09:46:02

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖