 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1165 | 回復(fù): 4 | ||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||
沒(méi)有實(shí)力銀蟲(chóng) (小有名氣)
|
[求助]
初學(xué)mathematica，琢磨了一天實(shí)在不行了，求教高手！
|
|
|
先把我寫(xiě)的內(nèi)容直接貼出來(lái)，問(wèn)題在后面： 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ù)制出來(lái)的，不知道高手們能不能用，我再傳一個(gè)圖片。 我的問(wèn)題是，我需要重復(fù)計(jì)算這個(gè)pH值，重復(fù)兩百次，然后畫(huà)出pH值隨重復(fù)次數(shù)的變化曲線圖，目前方程可能還有點(diǎn)問(wèn)題，得到的曲線圖不一定好看。關(guān)鍵問(wèn)題是現(xiàn)在我一不會(huì)寫(xiě)循環(huán)語(yǔ)句，二不會(huì)寫(xiě)作圖的命令。初學(xué)mathematica，請(qǐng)勿見(jiàn)笑，麻煩高手指教！如果能幫上忙僅有的三十余枚金幣全部奉上。多謝幫助！ 截圖mathematica [ Last edited by 沒(méi)有實(shí)力 on 2012-5-24 at 10:56 ] |
回復(fù)此樓» 猜你喜歡
295材料工程專(zhuān)碩求調(diào)劑
已經(jīng)有8人回復(fù)
-
求調(diào)劑
已經(jīng)有8人回復(fù)
-
293求調(diào)劑
已經(jīng)有5人回復(fù)
-
279求調(diào)劑
已經(jīng)有8人回復(fù)
-
328求調(diào)劑
已經(jīng)有8人回復(fù)
-
322求調(diào)劑：一志愿湖南大學(xué) 材料與化工（085600），已過(guò)六級(jí)。
已經(jīng)有9人回復(fù)
-
307求調(diào)劑
已經(jīng)有17人回復(fù)
-
286求調(diào)劑
已經(jīng)有5人回復(fù)
-
085600 295分求調(diào)劑
已經(jīng)有6人回復(fù)
-
287求調(diào)劑
已經(jīng)有13人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 求高手分析PS的TG-DTA圖，本人菜鳥(niǎo)，望指點(diǎn)
已經(jīng)有3人回復(fù)
- 求教各位高手！
已經(jīng)有3人回復(fù)
- 公派成功，但是接下來(lái)的問(wèn)題也很迷茫，求教大家！
已經(jīng)有11人回復(fù)
- 菜鳥(niǎo)求教MS復(fù)制粘貼問(wèn)題
已經(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)不上，總說(shuō)找不到文件，請(qǐng)高手指點(diǎn)！
已經(jīng)有16人回復(fù)
- 用氣相色譜測(cè)反式脂肪酸結(jié)果什么峰都沒(méi)有，急求高手解答��！
已經(jīng)有7人回復(fù)
- 【求助】請(qǐng)教FLUNT DPM 高手
已經(jīng)有3人回復(fù)
- 求助高手關(guān)于文章作圖，附件中的圖是怎么做出來(lái)的
已經(jīng)有9人回復(fù)
- 【求助】求教高手指導(dǎo)柵格數(shù)據(jù)投影
已經(jīng)有8人回復(fù)
mshwangg
至尊木蟲(chóng) (正式寫(xiě)手)
- 程序強(qiáng)帖: 5
- 應(yīng)助: 206 (大學(xué)生)
- 金幣: 10702.8
- 散金: 100
- 紅花: 19
- 帖子: 597
- 在線: 195.4小時(shí)
- 蟲(chóng)號(hào): 576702
- 注冊(cè): 2008-06-21
- 專(zhuān)業(yè): 物理學(xué)I
3樓2012-05-24 22:36:11
wangww2011
木蟲(chóng) (著名寫(xiě)手)
- 程序強(qiáng)帖: 13
- 應(yīng)助: 11 (小學(xué)生)
- 金幣: 4023.1
- 散金: 2709
- 紅花: 18
- 沙發(fā): 1
- 帖子: 1915
- 在線: 1537.1小時(shí)
- 蟲(chóng)號(hào): 772953
- 注冊(cè): 2009-05-17
- 性別: GG
- 專(zhuān)業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
2樓2012-05-24 19:56:09
xiaoqing8569
鐵桿木蟲(chóng) (著名寫(xiě)手)
奧林匹亞光學(xué)院院長(zhǎng)
- 應(yīng)助: 37 (小學(xué)生)
- 金幣: 7419.2
- 散金: 2240
- 紅花: 14
- 帖子: 1822
- 在線: 696.6小時(shí)
- 蟲(chóng)號(hào): 458097
- 注冊(cè): 2007-11-13
- 性別: GG
- 專(zhuān)業(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
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 300求調(diào)劑，材料科學(xué)英一數(shù)二 +18 | leaflight 2026-03-24 | 18/900 |
|
|---|---|---|---|---|
|
[考研] 求調(diào)劑 +7 | 研研，接電話 2026-03-24 | 8/400 |
|
|
|
[考研] 311求調(diào)劑 +10 | lin0039 2026-03-26 | 10/500 |
|
|
|
[考研] 327求調(diào)劑 +4 | 小卡不卡. 2026-03-29 | 4/200 |
|
|
|
[考研] 327求調(diào)劑 +6 | 汲亦昊 2026-03-29 | 6/300 |
|
|
|
[考研] 349求調(diào)劑 +6 | 李木子啊哈哈 2026-03-25 | 6/300 |
|
|
|
[考研] 330分求調(diào)劑 +5 | qzenlc 2026-03-29 | 5/250 |
|
|
|
[考研] 食品工程專(zhuān)碩一志愿中海洋309求調(diào)劑 +4 | 小張zxy張 2026-03-26 | 8/400 |
|
|
|
[考研] 一志愿北化085600材料專(zhuān)碩275|有文章專(zhuān)利｜求調(diào)劑 +7 | Micky11223 2026-03-25 | 7/350 |
|
|
|
[考研] 283求調(diào)劑 +3 | A child 2026-03-28 | 3/150 |
|
|
|
[考研] 085701環(huán)境工程，267求調(diào)劑 +16 | minht 2026-03-26 | 16/800 |
|
|
|
[考研] 286求調(diào)劑 +12 | PolarBear11 2026-03-26 | 12/600 |
|
|
|
[材料工程] 一志愿C9材料與化工專(zhuān)業(yè)總分300求調(diào)劑 +8 | 曼111 2026-03-24 | 9/450 |
|
|
|
[考研]
 一志愿 西北大學(xué) 總分282 英語(yǔ)一62 求調(diào)劑
+7
|
18419759900 2026-03-25 | 8/400 |
|
|
|
[考研] 314求調(diào)劑 +3 | 溪云珂 2026-03-26 | 3/150 |
|
|
|
[考研]
材料調(diào)劑
 5+4
|
想要一壺桃花水 2026-03-25 | 10/500 |
|
|
|
[考研] 中國(guó)科學(xué)院深圳先進(jìn)技術(shù)研究院-光纖傳感課題組招生-中國(guó)科學(xué)院大學(xué)、深圳理工大學(xué)聯(lián)培 +5 | YangTyu1 2026-03-26 | 5/250 |
|
|
|
[考研] 一志愿北化315 求調(diào)劑 +3 | akrrain 2026-03-24 | 3/150 |
|
|
|
[考研] 材料專(zhuān)碩找調(diào)劑 +5 | 哈哈哈吼吼吼哈 2026-03-23 | 5/250 |
|
|
|
[考研] 材料專(zhuān)碩331求調(diào)劑 +4 | 鮮當(dāng)牛 2026-03-24 | 4/200 |
|