24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1135  |  回复: 7
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

brqhl_ing

银虫 (小有名气)

[交流] 【求助】解常微分方程组已有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)
这两个其实就是两个系统的耦合  如何去解这个方程组?
回复此楼

» 猜你喜欢

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

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

captainmeng

木虫 (小有名气)

这里高手真多啊!!!!
8楼2013-11-22 00:19:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答

brqhl_ing

银虫 (小有名气)

引用回帖:
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

铁杆木虫 (著名写手)

方丈大师

★ ★
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楼2010-06-04 23:09:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

brqhl_ing

银虫 (小有名气)

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
的值求出来 在进行变换??
4楼2010-06-05 09:46:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复(可上传附件)
信息提示
请填处理意见