24小时热门版块排行榜    

查看: 525  |  回复: 2

309065816

木虫 (正式写手)

[求助] 采用solve函数求解非线性方程却无法运行,求指导!急急急

%求定态解和Jacobi矩阵,分析方程稳定性
[x,y,z]=solve('D*(c-x)-k1*x*(1-y-z)+k01*y-k4*x*z+k5*z=0','k1*x*(1-y-z)-k01*y-k02*z-k3*y*z=0','k2*y-k02*z+k3*y*z-k4*x*z-k5*z=0')
%求动力学方程组的零解的解析解
x1=x;y1=y;z1=z;
syms D c k1 k01 k2 k02 k3 k4 k5 V x y z ze zel;
ze=jacobian([D*(c-x)-k1*x*(1-y-z)+k01*y-k4*x*z+k5*z;...
k1*x*(1-y-z)-k01*y-k02*z-k3*y*z;k2*y-k02*z+k3*y*z-k4*x*z-k5*z],[x y z]);
%求动力学方程组的Jacobi矩阵
zel=subs(ze,{x,y,z},{x1(2),y1(2),z1(2)})
%将非平凡零解代入Jacobi矩阵
非线性方程组如图所示,运算时出现警告:Warning: Explicit solution could not be found.
这是怎么回事呢?

采用solve函数求解非线性方程却无法运行,求指导!急急急
S@UT`4J5~[%PHU@(PY0W(16.jpg
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cooooldog

铁杆木虫 (著名写手)

ส็็็

你是用符号方法求解的;
但是这个规模的问题,不存在"显式解";或至少用当前matlab的solve函数的符号解法不能给出"显式解"

也就是,解的解析形式是无法这样得到,或这个问题不存在解析解.
数值计算的稳定性是不能这样分析的; 只能通过前向或后向扰动分析的方法从数值理论的角度分析; 很多问题都不能有符号解的.
ส็็็็็็็็็็็็็็็็็็็็
2楼2014-07-28 12:28:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

309065816

木虫 (正式写手)

引用回帖:
2楼: Originally posted by cooooldog at 2014-07-28 12:28:19
你是用符号方法求解的;
但是这个规模的问题,不存在"显式解";或至少用当前matlab的solve函数的符号解法不能给出"显式解"

也就是,解的解析形式是无法这样得到,或这个问题不存在解析解.
数值 ...

嗯嗯,谢谢。用solve不能得到显示解,所以老师提出要手算并进行化简,运用求根公式进行计算。
稳定性的话,是原先文献里有用过这样的方法,看特征值是否是虚数什么的进行判断。但是文献里的方程比较简单,求出来的特征值也很简单,可以很容易的进行判断。但是上述方程过于复杂,目前准备采用其他方法进行判断。
3楼2014-07-29 12:22:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 309065816 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿中国海洋大学,生物学,301分,求调剂 +5 1孙悟空 2026-03-17 5/250 2026-03-19 18:03 by zcl123
[考研] 一志愿福大288有机化学,求调剂 +3 小木虫200408204 2026-03-18 3/150 2026-03-19 13:31 by houyaoxu
[考研] 287求调剂 +3 晨昏线与星海 2026-03-19 4/200 2026-03-19 12:32 by peike
[考研] 0703化学调剂,求各位老师收留 +10 秋有木北 2026-03-14 10/500 2026-03-19 05:52 by anny19840123
[考研] 330求调剂 +3 小材化本科 2026-03-18 3/150 2026-03-18 21:55 by 无懈可击111
[考研] 一志愿武理材料305分求调剂 +5 想上岸的鲤鱼 2026-03-18 6/300 2026-03-18 17:53 by 无际的草原
[考研] 295求调剂 +3 一志愿京区211 2026-03-18 5/250 2026-03-18 17:03 by zhaoqian0518
[考研] 085601专硕,总分342求调剂,地区不限 +5 share_joy 2026-03-16 5/250 2026-03-18 14:48 by haxia
[考研] 工科材料085601 279求调剂 +6 困于星晨 2026-03-17 6/300 2026-03-18 10:21 by kkcoco25
[考研] 301求调剂 +9 yy要上岸呀 2026-03-17 9/450 2026-03-18 08:58 by 无际的草原
[考研] 277调剂 +5 自由煎饼果子 2026-03-16 6/300 2026-03-17 19:26 by 李leezz
[考研] 326求调剂 +5 上岸的小葡 2026-03-15 6/300 2026-03-17 17:26 by ruiyingmiao
[考研] 26考研求调剂 +6 丶宏Sir 2026-03-13 6/300 2026-03-17 16:13 by 醉在风里
[考博] 26申博 +4 八6八68 2026-03-16 4/200 2026-03-17 13:00 by 轻松不少随
[考研] 275求调剂 +4 太阳花天天开心 2026-03-16 4/200 2026-03-17 10:53 by 功夫疯狂
[考研] 考研调剂 +3 淇ya_~ 2026-03-17 5/250 2026-03-17 09:25 by Winj1e
[考研] 278求调剂 +3 Yy7400 2026-03-13 3/150 2026-03-17 08:24 by laoshidan
[考研] 304求调剂 +3 曼殊2266 2026-03-14 3/150 2026-03-16 16:39 by houyaoxu
[考研] 327求调剂 +6 拾光任染 2026-03-15 11/550 2026-03-15 22:47 by 拾光任染
[考研] 中科大材料专硕319求调剂 +3 孟鑫材料 2026-03-13 3/150 2026-03-14 18:10 by houyaoxu
信息提示
请填处理意见