应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 用MATLAB求解四元一次方程组
11/1返回列表
查看: 3311  |  回复: 0
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

小犀牛来了

铁虫 (初入文坛)

[求助] 用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)
回复此楼

» 猜你喜欢
勇往直前
1楼 2015-04-20 21:25:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

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

科研从小木虫开始，人人为我，我为人人
相关版块跳转 我要订阅楼主 小犀牛来了 的主题更新
11/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭