 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1164 | 回復(fù): 4 | |||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | |||
沒有實(shí)力銀蟲 (小有名氣)
|
[求助]
初學(xué)mathematica，琢磨了一天實(shí)在不行了，求教高手！
|
||
|
先把我寫的內(nèi)容直接貼出來(lá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ù)制出來(lái)的，不知道高手們能不能用，我再傳一個(gè)圖片。 我的問題是，我需要重復(fù)計(jì)算這個(gè)pH值，重復(fù)兩百次，然后畫出pH值隨重復(fù)次數(shù)的變化曲線圖，目前方程可能還有點(diǎn)問題，得到的曲線圖不一定好看。關(guān)鍵問題是現(xiàn)在我一不會(huì)寫循環(huán)語(yǔ)句，二不會(huì)寫作圖的命令。初學(xué)mathematica，請(qǐng)勿見笑，麻煩高手指教！如果能幫上忙僅有的三十余枚金幣全部奉上。多謝幫助！ 截圖mathematica [ Last edited by 沒有實(shí)力 on 2012-5-24 at 10:56 ] |
回復(fù)此樓» 猜你喜歡
0703化學(xué)321分求調(diào)劑
已經(jīng)有4人回復(fù)
-
一志愿北京化工大學(xué)材料與化工（085600）296求調(diào)劑
已經(jīng)有19人回復(fù)
-
298求調(diào)劑
已經(jīng)有5人回復(fù)
-
346求調(diào)劑 一志愿070303有機(jī)化學(xué)
已經(jīng)有6人回復(fù)
-
288資源與環(huán)境專碩求調(diào)劑，不限專業(yè)，有學(xué)上就行
已經(jīng)有3人回復(fù)
-
324求調(diào)劑
已經(jīng)有11人回復(fù)
-
307求調(diào)劑
已經(jīng)有13人回復(fù)
-
071010 323 分求調(diào)劑
已經(jīng)有3人回復(fù)
-
材料調(diào)劑
已經(jīng)有9人回復(fù)
-
材料專碩調(diào)劑
已經(jīng)有10人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 求高手分析PS的TG-DTA圖，本人菜鳥，望指點(diǎn)
已經(jīng)有3人回復(fù)
- 求教各位高手！
已經(jīng)有3人回復(fù)
- 公派成功，但是接下來(lái)的問題也很迷茫，求教大家！
已經(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)不上，總說(shuō)找不到文件，請(qǐng)高手指點(diǎn)！
已經(jīng)有16人回復(fù)
- 用氣相色譜測(cè)反式脂肪酸結(jié)果什么峰都沒有，急求高手解答啊！
已經(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ù)
xiaoqing8569
鐵桿木蟲 (著名寫手)
奧林匹亞光學(xué)院院長(zhǎng)
- 應(yīng)助: 37 (小學(xué)生)
- 金幣: 7419.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
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
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 284求調(diào)劑 +14 | junqihahaha 2026-03-26 | 15/750 |
|
|---|---|---|---|---|
|
|
[考研] 322求調(diào)劑 +8 | 宋明欣 2026-03-27 | 8/400 |
|
|
[考研] 085701求調(diào)劑初試286分 +5 | secret0328 2026-03-28 | 5/250 |
|
|
|
[基金申請(qǐng)] 面上5B能上會(huì)嗎？ +7 | redcom 2026-03-29 | 7/350 |
|
|
|
[考研] 283求調(diào)劑（080500） +14 | A child 2026-03-27 | 14/700 |
|
|
|
[考研] 296求調(diào)劑 +10 | 彼岸t 2026-03-29 | 10/500 |
|
|
|
[考研] 一志愿哈爾濱工業(yè)大學(xué)材料與化工方向336分 +11 | 辰沐5211314 2026-03-26 | 11/550 |
|
|
|
[考研] 322求調(diào)劑：一志愿湖南大學(xué) 材料與化工（085600），已過(guò)六級(jí)。 +4 | XX小鄧 2026-03-29 | 4/200 |
|
|
|
[考研] 一志愿鄭州大學(xué)，080500學(xué)碩，總分317分求調(diào)劑 +8 | 舉個(gè)栗子oi 2026-03-24 | 9/450 |
|
|
|
[考研] 一志愿華北電力大學(xué)能動(dòng)專碩，293，求調(diào)劑 +3 | 15537177284 2026-03-23 | 5/250 |
|
|
|
[考研] 085701環(huán)境工程，267求調(diào)劑 +16 | minht 2026-03-26 | 16/800 |
|
|
|
[考研] 081200-314 +3 | LILIQQ 2026-03-27 | 4/200 |
|
|
|
[考研] 085701環(huán)境工程求調(diào)劑 +9 | 多久上課 2026-03-27 | 9/450 |
|
|
|
[考研] 275求調(diào)劑 +10 | jjjjjjjjjjl 2026-03-27 | 10/500 |
|
|
|
[考研] 285求調(diào)劑 +4 | AZMK 2026-03-27 | 7/350 |
|
|
|
[考研] 一志愿南師大0703化學(xué) 275求調(diào)劑 +4 | Ripcord上岸 2026-03-27 | 4/200 |
|
|
|
[考研] 298調(diào)劑 +3 | jiyingjie123 2026-03-27 | 3/150 |
|
|
|
[論文投稿] Journal of Mechanical Science and Technology +3 | Russ_ss 2026-03-25 | 5/250 |
|
|
|
[考研] 求b區(qū)院校調(diào)劑 +4 | 周56 2026-03-24 | 5/250 |
|
|
|
[考研]
 085404電子信息284分求調(diào)劑
+4
|
13659058978 2026-03-24 | 4/200 |
|