 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1143 | 回復(fù): 4 | ||
沒有實(shí)力銀蟲 (小有名氣)
|
[求助]
初學(xué)mathematica，琢磨了一天實(shí)在不行了，求教高手！
|
|
先把我寫的內(nèi)容直接貼出來，問題在后面： y1 = 40; y2 = 20; b = 0.18; Subscript[c, 1] = 10^-3; Subscript[c, 2] = 10^-4; Subscript[d, 1] = Subscript[d, 2] = 40; e = 1; g = 0.72; h = 10^-14; (*前面是已知的常數(shù)*) f[x_, y1_, y2_] = (y1/b)/( 1 + Subscript[c, 1]/x + Subscript[c, 1]*Subscript[d, 1]*e/(b*x)) + (y2/b)/( 1 + Subscript[c, 2]/x + Subscript[c, 2]*Subscript[d, 2]*e/(b*x)) == g + h/x - x; t = x /. Solve[f[x, y1, y2]][[4]]; pH = -Log10[t] (*這里得出我要的結(jié)果*) Ai1 = (y1/b)/( 1 + Subscript[c, 1]/t + Subscript[c, 1]*Subscript[d, 1]*e/(b*t)); Ai2 = (y2/b)/( 1 + Subscript[c, 2]/t + Subscript[c, 2]*Subscript[d, 2]*e/(b*t)); y1 = t*Ai1/Subscript[c, 1] + Ai1; y2 = t*Ai2/Subscript[c, 2] + Ai2; (*這里重新得到新的y1和y2*) f[x_, y1_, y2_] = (y1/b)/( 1 + Subscript[c, 1]/x + Subscript[c, 1]*Subscript[d, 1]*e/(b*x)) + (y2/b)/( 1 + Subscript[c, 2]/x + Subscript[c, 2]*Subscript[d, 2]*e/(b*x)) == g + h/x - x; t = x /. Solve[f[x, y1, y2]][[4]]; pH = -Log10[t] (*這里重復(fù)了一次，繼續(xù)得到我要的結(jié)果*) 前面是直接復(fù)制出來的，不知道高手們能不能用，我再傳一個(gè)圖片。 我的問題是，我需要重復(fù)計(jì)算這個(gè)pH值，重復(fù)兩百次，然后畫出pH值隨重復(fù)次數(shù)的變化曲線圖，目前方程可能還有點(diǎn)問題，得到的曲線圖不一定好看。關(guān)鍵問題是現(xiàn)在我一不會(huì)寫循環(huán)語句，二不會(huì)寫作圖的命令。初學(xué)mathematica，請(qǐng)勿見笑，麻煩高手指教！如果能幫上忙僅有的三十余枚金幣全部奉上。多謝幫助！ 截圖mathematica [ Last edited by 沒有實(shí)力 on 2012-5-24 at 10:56 ] |
回復(fù)此樓» 猜你喜歡
278求調(diào)劑
已經(jīng)有7人回復(fù)
-
26考研一志愿中國(guó)石油大學(xué)(華東)305分求調(diào)劑
已經(jīng)有5人回復(fù)
-
280求調(diào)劑
已經(jīng)有9人回復(fù)
-
279求調(diào)劑
已經(jīng)有5人回復(fù)
-
286求調(diào)劑
已經(jīng)有5人回復(fù)
-
材料與化工（0856）304求 B區(qū) 調(diào)劑
已經(jīng)有3人回復(fù)
-
268求調(diào)劑
已經(jīng)有8人回復(fù)
-
265求調(diào)劑
已經(jīng)有14人回復(fù)
-
22 350 本科985求調(diào)劑，求老登收留
已經(jīng)有3人回復(fù)
-
【考研調(diào)劑】化學(xué)專業(yè) 281分，一志愿四川大學(xué)，誠(chéng)心求調(diào)劑
已經(jīng)有9人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 求高手分析PS的TG-DTA圖，本人菜鳥，望指點(diǎn)
已經(jīng)有3人回復(fù)
- 求教各位高手！
已經(jīng)有3人回復(fù)
- 公派成功，但是接下來的問題也很迷茫，求教大家！
已經(jīng)有11人回復(fù)
- 菜鳥求教MS復(fù)制粘貼問題
已經(jīng)有4人回復(fù)
- 求教mathematica高手
已經(jīng)有5人回復(fù)
- SDS-PAGE凝膠電泳求教
已經(jīng)有14人回復(fù)
- 加溶劑后計(jì)算，不收斂，請(qǐng)高手指點(diǎn)
已經(jīng)有5人回復(fù)
- gaussian 09 做好輸入文件后總是關(guān)聯(lián)不上，總說找不到文件，請(qǐng)高手指點(diǎn)！
已經(jīng)有16人回復(fù)
- 用氣相色譜測(cè)反式脂肪酸結(jié)果什么峰都沒有，急求高手解答啊！
已經(jīng)有7人回復(fù)
- 【求助】請(qǐng)教FLUNT DPM 高手
已經(jīng)有3人回復(fù)
- 求助高手關(guān)于文章作圖，附件中的圖是怎么做出來的
已經(jīng)有9人回復(fù)
- 【求助】求教高手指導(dǎo)柵格數(shù)據(jù)投影
已經(jīng)有8人回復(fù)
wangww2011
木蟲 (著名寫手)
- 程序強(qiáng)帖: 13
- 應(yīng)助: 11 (小學(xué)生)
- 金幣: 4023.1
- 散金: 2709
- 紅花: 18
- 沙發(fā): 1
- 帖子: 1915
- 在線: 1537.1小時(shí)
- 蟲號(hào): 772953
- 注冊(cè): 2009-05-17
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
2樓2012-05-24 19:56:09
mshwangg
至尊木蟲 (正式寫手)
- 程序強(qiáng)帖: 5
- 應(yīng)助: 206 (大學(xué)生)
- 金幣: 10702.8
- 散金: 100
- 紅花: 19
- 帖子: 597
- 在線: 195.4小時(shí)
- 蟲號(hào): 576702
- 注冊(cè): 2008-06-21
- 專業(yè): 物理學(xué)I
3樓2012-05-24 22:36:11
xiaoqing8569
鐵桿木蟲 (著名寫手)
奧林匹亞光學(xué)院院長(zhǎng)
- 應(yīng)助: 37 (小學(xué)生)
- 金幣: 7414.2
- 散金: 2240
- 紅花: 14
- 帖子: 1822
- 在線: 696.6小時(shí)
- 蟲號(hào): 458097
- 注冊(cè): 2007-11-13
- 性別: GG
- 專業(yè): 應(yīng)用光學(xué)
【答案】應(yīng)助回帖
感謝參與，應(yīng)助指數(shù) +1
|
y1 = 40; y2 = 20; b = 0.18; Subscript[c, 1] = 10^-3; Subscript[c, 2] = 10^-4; Subscript[d, 1] = Subscript[d, 2] = 40; e = 1; g = 0.72; h = 10^-14; (*前面是已知的常數(shù)*) DD = (yy1/b)/(1 + Subscript[c, 1]/x + Subscript[c, 1]*Subscript[d, 1]*e/(b*x)) + (yy2/b)/(1 + Subscript[c, 2]/x + Subscript[c, 2]*Subscript[d, 2]*e/(b*x)) == g + h/x - x; phData = {}; For[i = 1, i <= 200, i++, t = x /. Solve[Evaluate[DD /. {yy1 -> y1, yy2 -> y2}]][[4]]; pH = -Log[10, t]; phData = Append[phData, pH]; Ai1 = (y1/b)/(1 + Subscript[c, 1]/t + Subscript[c, 1]*Subscript[d, 1]*e/(b*t)); Ai2 = (y2/b)/(1 + Subscript[c, 2]/t + Subscript[c, 2]*Subscript[d, 2]*e/(b*t)); y1 = t*Ai1/Subscript[c, 1] + Ai1; y2 = t*Ai2/Subscript[c, 2] + Ai2] phData; ListLinePlot[phData, PlotStyle -> Directive[Red, Thickness[.0035]], Joined -> True, PlotRange -> All, AspectRatio -> .5] Export["D:\\pH.gif", %, ImageSize -> {900, 650}] |
4樓2012-05-24 22:47:37
 |
5樓2012-05-25 12:03:46
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 沒有實(shí)力 的主題更新
| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 278求調(diào)劑 +7 | 煙火先于春 2026-03-17 | 7/350 |
|
|---|---|---|---|---|
|
|
[考研] 302求調(diào)劑 +11 | 呼呼呼。。。。 2026-03-17 | 11/550 |
|
|
[考研] 二本跨考鄭大材料306英一數(shù)二 +3 | z1z2z3879 2026-03-17 | 3/150 |
|
|
|
[考研] 一志愿華南師大 070300（化學(xué)）304分求調(diào)劑 +3 | 0703武芊慧雪304 2026-03-18 | 3/150 |
|
|
|
[考研] 330求調(diào)劑 +4 | 小材化本科 2026-03-18 | 4/200 |
|
|
|
[考研] 085600材料與化工 +8 | 安全上岸！ 2026-03-16 | 8/400 |
|
|
|
[考研] 一志愿中南化學(xué)（0703）總分337求調(diào)劑 +8 | niko- 2026-03-19 | 9/450 |
|
|
|
[考研] 北科281學(xué)碩材料求調(diào)劑 +5 | tcxiaoxx 2026-03-20 | 5/250 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 ，六級(jí)已過，有科研經(jīng)歷 +13 | 曦熙兮 2026-03-15 | 13/650 |
|
|
|
[考研] 廣西大學(xué)家禽遺傳育種課題組2026年碩士招生（接收計(jì)算機(jī)專業(yè)調(diào)劑） +3 | 123阿標(biāo) 2026-03-17 | 3/150 |
|
|
|
[考研]
 320求調(diào)劑0856
+3
|
不想起名字112 2026-03-19 | 3/150 |
|
|
|
[考研] 085601材料工程專碩求調(diào)劑 +10 | 慕寒mio 2026-03-16 | 10/500 |
|
|
|
[考研] 085600材料與化工調(diào)劑 324分 +10 | llllkkkhh 2026-03-18 | 12/600 |
|
|
|
[考研] 085600材料與化工求調(diào)劑 +6 | 緒幸與子 2026-03-17 | 6/300 |
|
|
|
[考研] 材料工程專碩調(diào)劑 +5 | 204818@lcx 2026-03-17 | 6/300 |
|
|
|
[考研] 收復(fù)試調(diào)劑生 +4 | 雨后秋荷 2026-03-18 | 4/200 |
|
|
|
[考研]
0854，計(jì)算機(jī)類招收調(diào)劑
+3
|
胡辣湯放糖 2026-03-15 | 6/300 |
|
|
|
[考研] 考研化學(xué)學(xué)碩調(diào)劑，一志愿985 +4 | 張vvvv 2026-03-15 | 6/300 |
|
|
|
[考研] 有沒有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3 | TqlXswl 2026-03-16 | 7/350 |
|
|
|
[考研] 一志愿211 0703方向310分求調(diào)劑 +3 | 努力奮斗112 2026-03-15 | 3/150 |
|