24小时热门版块排行榜    

查看: 3434  |  回复: 0

小犀牛来了

铁虫 (初入文坛)

[求助] 用MATLAB求解四元一次方程组

方程组为
-m*n*k*p_1/(1-g)+m*n*k*p_2/(1-g)-b*s-h*m*n*k=0
-2*p_1/(1-g)-2*p_1/g+2*p_2/(1-g)-m*n*k*s/(1-g)-h/g-n*k/(1-g)=0
2*p_1/(1-g)-2*p_2/(1-g)+m*n*k*s/(1-g)+h+1+n*k/(1-g)=0
-p_1/g+p_2-m*n*k*s-n*k=0
求p_1,p_2,s,h的值
我自己的用MATLAB求解过程及结果如下,但不知道哪儿出了问题,我求得的p_1、p_2、s、h值竟然分别是我看到的一篇论文上的h、p_1、p_2、s的值 (方程组都是上面这个)
syms p_1 p_2 s h m n g k b;
>> [p_1,p_2,s,h]=solve('-m*n*k*p_1/(1-g)+m*n*k*p_2/(1-g)-b*s-h*m*n*k=0','-2*p_1/(1-g)-2*p_1/g+2*p_2/(1-g)-m*n*k*s/(1-g)-h/g-n*k/(1-g)=0','2*p_1/(1-g)-2*p_2/(1-g)+m*n*k*s/(1-g)+h+1+n*k/(1-g)=0','-p_1/g+p_2-m*n*k*s-n*k=0','p_1,p_2,s,h')

p_1 =

n*k*g*(2*b+m^2*n*k)/(-2*b*g+2*b-m^2*n^2*k^2+m^2*n^2*k^2*g)

p_2 =

g*(-b+b*n*k+m^2*n^2*k^2)/(-2*b+m^2*n^2*k^2)

s =

-(1+n*k)*b/(-2*b+m^2*n^2*k^2)

h =

-(1+n*k)*m*n*k/(-2*b+m^2*n^2*k^2)
回复此楼
勇往直前
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

找到一些相关的精华帖子,希望有用哦~

科研从小木虫开始,人人为我,我为人人
相关版块跳转 我要订阅楼主 小犀牛来了 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 师弟论文见刊大半年才想起来申请专利,还能抢救一下吗? +5 13108017953 2026-09-21 5/250 2026-09-27 16:46 by Leogzhya
[教师之家] 某top大学教授说“能够在市场中兑现的能力才是真能力”无比同意! +3 zju2000 2026-09-26 3/150 2026-09-27 16:26 by lala1008
[硕博家园] 上海工程技术大学 激光智能制造课题组 2027级博士研究生招生 +5 oxidpl 2026-09-20 7/350 2026-09-27 10:43 by Equinoxhua
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 ViqFlLSxHDGt 2026-09-26 3/150 2026-09-27 10:26 by 4IW0sJvtEqX8
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 ViqFlLSxHDGt 2026-09-26 3/150 2026-09-27 10:21 by 4IW0sJvtEqX8
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 ViqFlLSxHDGt 2026-09-26 3/150 2026-09-27 10:15 by 4IW0sJvtEqX8
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 ViqFlLSxHDGt 2026-09-26 3/150 2026-09-27 10:03 by 4IW0sJvtEqX8
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 GxESXJmptOtQ 2026-09-25 4/200 2026-09-27 06:44 by XgCC8uTcwILl
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 GxESXJmptOtQ 2026-09-25 4/200 2026-09-27 06:42 by XgCC8uTcwILl
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 GxESXJmptOtQ 2026-09-25 4/200 2026-09-27 06:39 by XgCC8uTcwILl
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 GxESXJmptOtQ 2026-09-25 3/150 2026-09-27 06:32 by XgCC8uTcwILl
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 vQZoDrrm7mUF 2026-09-25 3/150 2026-09-27 06:08 by XgCC8uTcwILl
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 vQZoDrrm7mUF 2026-09-25 3/150 2026-09-27 05:56 by XgCC8uTcwILl
[论文投稿] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 vQZoDrrm7mUF 2026-09-25 4/200 2026-09-27 05:56 by XgCC8uTcwILl
[找工作] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +5 vQZoDrrm7mUF 2026-09-25 5/250 2026-09-27 05:53 by XgCC8uTcwILl
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 vQZoDrrm7mUF 2026-09-25 4/200 2026-09-27 05:50 by XgCC8uTcwILl
[考研] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 vQZoDrrm7mUF 2026-09-25 5/250 2026-09-27 05:44 by XgCC8uTcwILl
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 vQZoDrrm7mUF 2026-09-25 4/200 2026-09-27 05:44 by XgCC8uTcwILl
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 vQZoDrrm7mUF 2026-09-25 3/150 2026-09-27 05:44 by XgCC8uTcwILl
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 vQZoDrrm7mUF 2026-09-25 6/300 2026-09-27 05:42 by XgCC8uTcwILl
信息提示
请填处理意见