声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1175|回复: 6

[综合讨论] 请教!非线性微分方程组,给定末时刻条件

[复制链接]
发表于 2008-10-27 21:55 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
请教各位高手们,这样的方程组用MATLAB怎么解,ODE解的都是给定初始状态的解,但是这个是给定末时刻的状态

方程组

方程组

,对t求导,

要求解的变量

要求解的变量
末时刻状态为

末时刻状态

末时刻状态

谢谢了。。。

[ 本帖最后由 ChaChing 于 2009-3-7 12:15 编辑 ]
回复
分享到:

使用道具 举报

发表于 2008-10-28 16:06 | 显示全部楼层
终值问题,只需要把[t0,tf]中的t0设置成终止时间,tf设置成开始时间就可以了。

  1. function dx=zhongzhiode(t,x)
  2. dx=[2*x(2)^2-2;
  3. -x(1)+2*x(2)*x(3)-1;
  4. -2*x(2)+2*x(3)^2-4];
复制代码


  1. [t,x]=ode45(@zhongzhiode,[3,0],[1;0;2]);plot(t,x)
复制代码
终值问题.jpg

评分

1

查看全部评分

 楼主| 发表于 2008-10-28 17:17 | 显示全部楼层
非常感谢sogooda的回复!原来可以这样用的。。。谢谢!
发表于 2009-3-7 09:45 | 显示全部楼层
这个问题用simulink怎么解决呢?
发表于 2009-3-8 08:52 | 显示全部楼层

回复 地板 Hyowinner 的帖子

到simulink版问问看。我对simulink基本一窍不通。
发表于 2009-3-9 14:38 | 显示全部楼层

回复 沙发 sogooda 的帖子

为何我改成初值问题击出问题?
[t,x]=ode45(@zhongzhiode,[0,3],[1;0;2]);plot(t,x)
发表于 2009-3-10 08:45 | 显示全部楼层

回复 6楼 ChaChing 的帖子

可能和初值有关,也许在那个初始值时,方程有些特殊吧。我把方程的初值稍微改了一下,就可以了。
  1. [t,x]=ode45(@zhongzhiode,[0,3],[2.5;1;1.7]);plot(t,x)
复制代码
untitled.jpg

评分

1

查看全部评分

您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-11 23:12 , Processed in 0.068138 second(s), 24 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表