| 查看: 1159 | 回复: 4 | |||
[求助]
初学mathematica，琢磨了一天实在不行了，求教高手！
|
|
先把我写的内容直接贴出来，问题在后面： 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; (*前面是已知的常数*) 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] (*这里得出我要的结果*) 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] (*这里重复了一次，继续得到我要的结果*) 前面是直接复制出来的，不知道高手们能不能用，我再传一个图片。 我的问题是，我需要重复计算这个pH值，重复两百次，然后画出pH值随重复次数的变化曲线图，目前方程可能还有点问题，得到的曲线图不一定好看。关键问题是现在我一不会写循环语句，二不会写作图的命令。初学mathematica，请勿见笑，麻烦高手指教！如果能帮上忙仅有的三十余枚金币全部奉上。多谢帮助！  截图mathematica [ Last edited by 没有实力 on 2012-5-24 at 10:56 ] |
回复此楼» 猜你喜欢
调剂310
已经有11人回复
-
300求调剂，材料科学英一数二
已经有9人回复
-
296求调剂
已经有3人回复
-
一志愿南昌大学324求调剂
已经有3人回复
-
299求调剂
已经有9人回复
-
22408 359分调剂
已经有6人回复
-
085600 材料与化工 329分求调剂
已经有12人回复
-
一志愿北京化工大学材料与化工（085600）296求调剂
已经有14人回复
-
085601材料工程找调剂
已经有11人回复
-
总分293求调剂
已经有9人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求高手分析PS的TG-DTA图，本人菜鸟，望指点
已经有3人回复
- 求教各位高手！
已经有3人回复
- 公派成功，但是接下来的问题也很迷茫，求教大家！
已经有11人回复
- 菜鸟求教MS复制粘贴问题
已经有4人回复
- 求教mathematica高手
已经有5人回复
- SDS-PAGE凝胶电泳求教
已经有14人回复
- 加溶剂后计算，不收敛，请高手指点
已经有5人回复
- gaussian 09 做好输入文件后总是关联不上，总说找不到文件，请高手指点！
已经有16人回复
- 用气相色谱测反式脂肪酸结果什么峰都没有，急求高手解答啊！
已经有7人回复
- 【求助】请教FLUNT DPM 高手
已经有3人回复
- 求助高手关于文章作图，附件中的图是怎么做出来的
已经有9人回复
- 【求助】求教高手指导栅格数据投影
已经有8人回复
wangww2011
木虫 (著名写手)
- 程序强帖: 13
- 应助: 11 (小学生)
- 金币: 4023.1
- 散金: 2709
- 红花: 18
- 沙发: 1
- 帖子: 1915
- 在线: 1537.1小时
- 虫号: 772953
- 注册: 2009-05-17
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
2楼2012-05-24 19:56:09
mshwangg
至尊木虫 (正式写手)
- 程序强帖: 5
- 应助: 206 (大学生)
- 金币: 10702.8
- 散金: 100
- 红花: 19
- 帖子: 597
- 在线: 195.4小时
- 虫号: 576702
- 注册: 2008-06-21
- 专业: 物理学I
3楼2012-05-24 22:36:11
xiaoqing8569
铁杆木虫 (著名写手)
奥林匹亚光学院院长
- 应助: 37 (小学生)
- 金币: 7414.2
- 散金: 2240
- 红花: 14
- 帖子: 1822
- 在线: 696.6小时
- 虫号: 458097
- 注册: 2007-11-13
- 性别: GG
- 专业: 应用光学
【答案】应助回帖
感谢参与，应助指数 +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; (*前面是已知的常数*) 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
20