| 查看: 305 | 回复: 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 |
回复此楼» 猜你喜欢
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有23人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有15人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有6人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有6人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有6人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有7人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有7人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有10人回复
50