| 查看: 1710 | 回复: 11 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
wanyshenger金虫 (正式写手)
|
[求助]
请教大神：matlab怎么求解有符号变量的方程组 已有2人参与
|
|
|
MATLAB刚刚入门，不太会用，老板让我计算绝热火焰温度，现在遇到的问题是燃烧生成的CO2,H2O会分解，但化学平衡的方程组不会解，简单来说就是方程组里有符号变量如何求解 例如： x*[c+0.5*(x+y)]^0.5/(a-x)=K1; y*[c+0.5*(x+y)]^0.5/(b-y)=K2; 其中a,b,c,K1,K2 均为过程中赋值的变量，求解x，y，请问这种情况应该如何用MATLAB求解。请各位大神不吝指教，感谢感谢！ |
回复此楼» 猜你喜欢
综述论文作为代表作会不会影响评审专家的印象分？
已经有12人回复
-
2019年青年基金涵评意见，大家看看几个A，几个B?
已经有4人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有29人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有9人回复
-
是这周出结果还是下周出结果？
已经有3人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有9人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有4人回复
-
帮忙看看fileCode
已经有13人回复
-
为什么网上很多人说本周 12号出结果
已经有7人回复
-
读博的好处
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 救助关于几率分布、图形拟合的问题
已经有1人回复
断剑残雪
铁杆木虫 (职业作家)
- 应助: 25 (小学生)
- 贵宾: 0.057
- 金币: 8517.4
- 散金: 9307
- 红花: 191
- 沙发: 5
- 帖子: 3566
- 在线: 1232.8小时
- 虫号: 2187840
- 注册: 2012-12-15
- 专业: 代数学
3楼2015-07-30 08:15:49
断剑残雪
铁杆木虫 (职业作家)
- 应助: 25 (小学生)
- 贵宾: 0.057
- 金币: 8517.4
- 散金: 9307
- 红花: 191
- 沙发: 5
- 帖子: 3566
- 在线: 1232.8小时
- 虫号: 2187840
- 注册: 2012-12-15
- 专业: 代数学
2楼2015-07-30 08:14:25
断剑残雪
铁杆木虫 (职业作家)
- 应助: 25 (小学生)
- 贵宾: 0.057
- 金币: 8517.4
- 散金: 9307
- 红花: 191
- 沙发: 5
- 帖子: 3566
- 在线: 1232.8小时
- 虫号: 2187840
- 注册: 2012-12-15
- 专业: 代数学
4楼2015-07-30 08:17:42
断剑残雪
铁杆木虫 (职业作家)
- 应助: 25 (小学生)
- 贵宾: 0.057
- 金币: 8517.4
- 散金: 9307
- 红花: 191
- 沙发: 5
- 帖子: 3566
- 在线: 1232.8小时
- 虫号: 2187840
- 注册: 2012-12-15
- 专业: 代数学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★
wanyshenger: 金币+8, ★★★很有帮助 2015-07-30 13:00:56
wanyshenger: 金币+8, ★★★很有帮助 2015-07-30 13:00:56
|
最终版： syms z a b c K1 K2 x y; a=1 ;b=1 ;c=1 ;K1=1 ;K2=1; %比如已知它们的值，要求x，y。那么把前面的注释去掉再运行即可，它们的值也可以自己设定。 [x,y]=solve(x*sqrt(c+0.5*(x+y))/(a-x)==K1 , y*sqrt(c+0.5*(x+y))/(b-y)==K2,x,y); %给定方程组，并返回x，y的符号表达式; z1 = solve(z^4 + z^3*(K1 + K2) - z^2*(c - K1*K2) - z*((K1*a)/2 + (K2*b)/2 + K1*c + K2*c) - K1*K2*c - (K1*K2*b)/2 - (K1*K2*a)/2==0, z) x=subs(x,'z1',z1) %返回x的值 y=subs(y,'z1',z1) %返回y的值 |
5楼2015-07-30 08:34:25
