| 查看: 567 | 回复: 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. 这是怎么回事呢？  S@UT`4J5~[%PHU@(PY0W(16.jpg |
回复此楼» 猜你喜欢
FileCode能看出啥？
已经有23人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有36人回复
-
关于Filecode分析方法
已经有9人回复
-
2019年青年基金涵评意见，大家看看几个A，几个B?
已经有10人回复
-
不应该看fileCode
已经有6人回复
-
好奇怪的filecode
已经有6人回复
-
长年满屏的广告，版主太不责任了。基金也等的急
已经有4人回复
-
filecode
已经有15人回复
-
职称评审，求友友推荐见刊最快的期刊
已经有4人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有10人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求fortran计算方法指导 14元非线性方程组 之前8元的使用bfs法能解14元解不了呀求指导
已经有4人回复
- matlab 求解非线性方程组未果
已经有10人回复
- 求解这类方程 用什么软件编程较好····
已经有8人回复
- 求教matlab傅里叶变换求解
已经有5人回复
- 二阶非线性微分方程求解
已经有14人回复
- 求助！matlab用fsolve函数求解非线性方程组的问题！
已经有19人回复
- 毕设，数学相关专业，求指导！
已经有18人回复
- 【求助】求解隐函数
已经有6人回复
cooooldog
铁杆木虫 (著名写手)
ส็็็
- 应助: 237 (大学生)
- 金币: 6101.9
- 散金: 1114
- 红花: 39
- 帖子: 1380
- 在线: 553.8小时
- 虫号: 506699
- 注册: 2008-02-18
- 专业: 模式识别
2楼2014-07-28 12:28:19
309065816
木虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 2494.2
- 散金: 180
- 红花: 3
- 帖子: 516
- 在线: 142.3小时
- 虫号: 1227451
- 注册: 2011-03-09
- 性别: MM
- 专业: 电化学
3楼2014-07-29 12:22:24
