 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1141 | 回復(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ù)此樓» 猜你喜歡
346求調(diào)劑[0856]
已經(jīng)有7人回復(fù)
-
一志愿山大07化學(xué) 332分 四六級(jí)已過 本科山東雙非 求調(diào)劑！
已經(jīng)有3人回復(fù)
-
310求調(diào)劑
已經(jīng)有3人回復(fù)
-
機(jī)械專碩299求調(diào)劑至材料
已經(jīng)有4人回復(fù)
-
070300化學(xué)319求調(diào)劑
已經(jīng)有7人回復(fù)
-
08工科 320總分 求調(diào)劑
已經(jīng)有6人回復(fù)
-
一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑
已經(jīng)有12人回復(fù)
-
307求調(diào)劑
已經(jīng)有3人回復(fù)
-
265求調(diào)劑
已經(jīng)有3人回復(fù)
-
301求調(diào)劑
已經(jīng)有10人回復(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ā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿 西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。 +3 | 晨昏線與星海 2026-03-18 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境專碩，總分308求調(diào)劑 +3 | 墨墨漠 2026-03-18 | 3/150 |
|
|
|
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +8 | 小材化本科 2026-03-18 | 8/400 |
|
|
|
[考研] 288求調(diào)劑 +16 | 于海海海海 2026-03-19 | 16/800 |
|
|
|
[考研] 317求調(diào)劑 +5 | 申子申申 2026-03-19 | 9/450 |
|
|
|
[考研] 一志愿 西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。 +4 | 晨昏線與星海 2026-03-19 | 4/200 |
|
|
|
[考研]
 材料與化工 322求調(diào)劑
+4
|
然11 2026-03-19 | 4/200 |
|
|
|
[考研] 藥學(xué)383 求調(diào)劑 +3 | 藥學(xué)chy 2026-03-15 | 5/250 |
|
|
|
[考研] 265求調(diào)劑 +12 | 梁梁校校 2026-03-19 | 13/650 |
|
|
|
[考研] 一志愿華中農(nóng)業(yè)071010，總分320求調(diào)劑 +3 | 困困困困坤坤 2026-03-20 | 3/150 |
|
|
|
[考研] 一志愿西安交通大學(xué) 學(xué)碩 354求調(diào)劑 +3 | 我想要讀研究生 2026-03-20 | 3/150 |
|
|
|
[考研] 一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑 +11 | Ymlll 2026-03-18 | 15/750 |
|
|
|
[考研] 0856調(diào)劑，是學(xué)校就去 +8 | sllhht 2026-03-19 | 9/450 |
|
|
|
[考研] 招收調(diào)劑碩士 +4 | lidianxing 2026-03-19 | 12/600 |
|
|
|
[考博] 招收博士1-2人 +3 | QGZDSYS 2026-03-18 | 3/150 |
|
|
|
[論文投稿]
申請(qǐng)回稿延期一個(gè)月，編輯同意了。但系統(tǒng)上的時(shí)間沒變，給編輯又寫郵件了，沒回復(fù)
 10+3
|
wangf9518 2026-03-17 | 4/200 |
|
|
|
[考研] 梁成偉老師課題組歡迎你的加入 +9 | 一鴨鴨喲 2026-03-14 | 11/550 |
|
|
|
[考研] 286求調(diào)劑 +6 | lemonzzn 2026-03-16 | 10/500 |
|
|
|
[考研] 312求調(diào)劑 +8 | 陌宸希 2026-03-16 | 9/450 |
|
|
|
[考研] 0856專碩279求調(diào)劑 +5 | 加油加油！? 2026-03-15 | 5/250 |
|