| 查看: 226 | 回复: 0 | ||
lyjade2014金虫 (小有名气)
|
[求助]
matlab菜鸟一枚，请走过路过的进来指导一下，跪谢！！！
|
|
公式为： 1=x•(1+10^(9.27-pH)+[Ni]•10^2.8+[Mn]•10^1.0)+x^2•(2•[Ni]•10^4.04+2•[Mn]•10^1.54) +x^3•(3•[Ni]•10^5.77+3•[Mn]•10^1.7)+x^4•(4•[Ni]•10^6.96+4•[Mn]•10^1.3) +x^5•(5•[Ni]•10^2.71)+x^6•(6•[Ni]•10^2.74) 其中[Ni][Co][Mn]是上一步计算出来的，可以看做是已知的，而pH是个变量，我希望给出一个pH就能得到一个x（是正实数解）。以下是我编写的程序，能不能帮我看看我写的对不对呢？跪谢！！求回复！！ %计算出x a1=1+10^(9.27-pH)+[Ni]*10^2.8+[Mn]*10; a2=2*[Ni]*10^4.04+2*[Mn]*10^1.54; a3=3*[Ni]*10^5.77+3*[Mn]*10^1.70; a4=4*[Ni]*10^6.96+4*[Mn]*10^1.30; a5=5*[Ni]*10^2.71; a6=6*[Ni]*10^2.74; p1=[a6 a5 a4 a3 a2 a1 -1]; x=roots(p1); % f=sym('x*a1+x^2*a2+x^3*a3+x^4*a4+x^5*a5+x^6*a6-1'); % solve('f','x'); for k=1:length(x) if x(k)>0&&isreal(x(k)) index=k; else index=[]; end end if ~isempty(index) x=x(index); end |
回复此楼» 猜你喜欢
谈谈两天一夜的“延安行”
已经有15人回复
-
EST投稿状态问题
已经有6人回复
-
职称评审没过，求安慰
已经有15人回复
-
垃圾破二本职称评审标准
已经有11人回复
-
投稿Elsevier的Neoplasia杂志，到最后选publishing options时页面空白，不能完成投稿
已经有16人回复
-
毕业后当辅导员了，天天各种学生超烦
已经有4人回复
-
聘U V热熔胶研究人员
已经有10人回复
-
求助文献
已经有3人回复
-
投稿返修后收到这样的回复，还有希望吗
已经有8人回复
-
三无产品还有机会吗
已经有6人回复