跪求各位亲请教非自治系统的LYapunov指数程序相关问题

mumianhua · 发表于 2014-6-17 10:06

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

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

x

各位，亲们，本人做非自治系统的LYapunov指数计算，遇到问题，状态化后八个方程，其中还有非线性油膜力以及其他非线性项（正余弦函数），在计算雅克比矩阵时，即对方程各个变量分别进行求导时，计算结果过于复杂，不知道是自己的方法有问题，还是别的，恳请大家帮助。I need your hellp

mumianhua · 发表于 2014-6-17 21:12

希望能够得到大家的帮助谢谢 I need your help

mumianhua · 发表于 2014-6-17 21:14

各位，亲们，本人做非自治系统的LYapunov指数计算，遇到问题，状态化后八个方程，其中还有非线性油膜力以及其他非线性项（正余弦函数），在计算雅克比矩阵时，即对方程各个变量分别进行求导时，计算结果过于复杂，不知道是自己的方法有问题，还是别的，恳请大家帮助。I need your hellp

mumianhua · 发表于 2014-6-18 11:22

亲怎么没人回应呢你们很忙吗？？？？？？？

mumianhua · 发表于 2014-6-22 09:31

不用大家帮忙了我已经找到方法了结果还不错想要呈现者@@@@@我

wuzhongshanhe · 发表于 2014-6-25 07:54

我看未必！

wystar · 发表于 2014-6-30 15:58

我现在遇到了同样问题，4个状态方程，分母中还含有未知量，在计算雅克比矩阵时，也是相当的复杂。求帮助啊。

mumianhua · 发表于 2014-6-30 17:16

wystar 发表于 2014-6-30 15:58
我现在遇到了同样问题，4个状态方程，分母中还含有未知量，在计算雅克比矩阵时，也是相当的复杂。求帮助啊 ...

x其实计算Jac很容易的你不要对方程中的各项求导，只需要利用matlab编程即可以求出Jac 你试试

wystar · 发表于 2014-6-30 17:32

x其实计算Jac很容易的你不要对方程中的各项求导，只需要利用matlab编程即可以求出Jac 你试试

我尝试过对方程中的各项求导,也尝试过用matlab的符号计算进行编程自动求导，Jac求出来之后计算，求李雅普诺夫指数存在问题，无法求解结果。采用方法是Wolf 等提出的基于Gram_Schmidt正交化求解方法。下面是错误提示。
Warning: Failure at t=1.386101e-001.  Unable to meet integration tolerances without reducing the step size below the smallest value
allowed (4.440892e-016) at time t.
> In ode15s at 753
  In lyapunov at 113
  In run_lyap at 4
Elapsed time is 3.801703 seconds.

mumianhua · 发表于 2014-6-30 20:19

wystar 发表于 2014-6-30 17:32
我尝试过对方程中的各项求导,也尝试过用matlab的符号计算进行编程自动求导，Jac求出来之后计算，求李雅普 ...

你试试利用matlab编程计算Jac矩阵呀很简单你试试先

wystar · 发表于 2014-6-30 20:32

你试试利用matlab编程计算Jac矩阵呀很简单你试试先

已经试过了的啊，可惜求不出来啊。这个是我做的间隙机构的相轨迹图，不知道是否正确。D:\

wystar · 发表于 2014-6-30 20:35

D:\ode45相图.JPG

wystar · 发表于 2014-6-30 20:41

晕了，图片一直没有上传上去

mumianhua · 发表于 2014-7-1 08:59

wystar 发表于 2014-6-30 20:32
已经试过了的啊，可惜求不出来啊。这个是我做的间隙机构的相轨迹图，不知道是否正确。

那你贴图吧我来看看贴程序看看

mumianhua · 发表于 2014-7-1 08:59

wystar 发表于 2014-6-30 20:35

？？？？？？？？

账号		自动登录	找回密码
密码			我要加入

[稳定性与分岔] 跪求各位亲请教非自治系统的LYapunov指数程序相关问题

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

本帖被以下淘专辑推荐:

点评

评分

评分

浏览过的版块

[稳定性与分岔] 跪求 各位亲 请教非自治系统的LYapunov指数程序相关问题

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

本帖被以下淘专辑推荐:

点评

评分

评分

浏览过的版块

[稳定性与分岔] 跪求各位亲请教非自治系统的LYapunov指数程序相关问题