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

登录注册

应《网络安全法》要求，自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)
这两个其实就是两个系统的耦合如何去解这个方程组？

回复此楼

» 猜你喜欢

职称评审没过，求安慰已经有49人回复
26申博自荐已经有3人回复
A期刊撤稿已经有4人回复
垃圾破二本职称评审标准已经有17人回复
投稿Elsevier的Neoplasia杂志，到最后选publishing options时页面空白，不能完成投稿已经有22人回复
EST投稿状态问题已经有7人回复
毕业后当辅导员了，天天各种学生超烦已经有4人回复
三无产品还有机会吗已经有6人回复

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

求助四波混频耦合波方程数值求解已经有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的回帖