| 1 | 1/1 | 返回列表 |
| 查看: 238 | 回復: 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(是正實數(shù)解)。以下是我編寫的程序,能不能幫我看看我寫的對不對呢?跪謝。∏蠡貜停! %計算出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 |
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0703化學調(diào)劑 ,六級已過,有科研經(jīng)歷 +12 | 曦熙兮 2026-03-15 | 12/600 |
|
|---|---|---|---|---|
|
[考研] 【考研調(diào)劑】化學專業(yè) 281分,一志愿四川大學,誠心求調(diào)劑 +5 | 吃吃吃才有意義 2026-03-19 | 5/250 |
|
|
[考研] 266求調(diào)劑 +5 | 陽陽哇塞 2026-03-14 | 10/500 |
|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-18 | 3/150 |
|
|
[考研] 346求調(diào)劑[0856] +3 | WayneLim327 2026-03-16 | 6/300 |
|
|
[考研] 0703化學調(diào)劑 +8 | 妮妮ninicgb 2026-03-15 | 12/600 |
|
|
[考研] 材料專業(yè)求調(diào)劑 +5 | hanamiko 2026-03-18 | 5/250 |
|
|
[考研] 085601專碩,總分342求調(diào)劑,地區(qū)不限 +5 | share_joy 2026-03-16 | 5/250 |
|
|
[考研] 280求調(diào)劑 +6 | 咕嚕曉曉 2026-03-18 | 7/350 |
|
|
[考研] 工科材料085601 279求調(diào)劑 +6 | 困于星晨 2026-03-17 | 6/300 |
|
|
[考研] 334求調(diào)劑 +3 | 志存高遠意在機?/a> 2026-03-16 | 3/150 |
|
|
[考博] 26博士申請 +3 | 1042136743 2026-03-17 | 3/150 |
|
|
[考研] 一志愿蘇州大學材料工程(085601)專碩有科研經(jīng)歷三項國獎兩個實用型專利一項省級立項 +6 | 大火山小火山 2026-03-16 | 8/400 |
|
|
[考研] 211本,11408一志愿中科院277分,曾在中科院自動化所實習 +6 | Losir 2026-03-12 | 7/350 |
|
|
[考研] 275求調(diào)劑 +4 | 太陽花天天開心 2026-03-16 | 4/200 |
|
|
[考研]
|
zhouzhen654 2026-03-16 | 3/150 |
|
|
[考研] 0703一志愿211 285分求調(diào)劑 +5 | ly3471z 2026-03-13 | 5/250 |
|
|
[考研] 材料與化工 323 英一+數(shù)二+物化,一志愿:哈工大 本人本科雙一流 +4 | 自由的_飛翔 2026-03-13 | 5/250 |
|
|
[基金申請] 現(xiàn)在如何回避去年的某一個專家,不知道名字 +3 | zk200107 2026-03-12 | 6/300 |
|
|
[碩博家園] 085600 260分求調(diào)劑 +3 | 天空還下雨么 2026-03-13 | 5/250 |
|