| 查看: 1532 | 回复: 16 | ||||
[交流]
【求助】帮忙作一个图~
|
||||
 相关的mathematica 公式已经给写好，就是不清楚怎么做，请高手帮忙，给出源代码和结果(最好用7.0的版本)~ C1 = 120; C2 = 100; C3= 60; Aa = - (Cos[alpha] Cos[v] + Sin[alpha] Sin[v])/Tan [beta]; Dir = Cos[beta]^2 (Cos[alpha]^2 Cos[v]^2 + Sin[alpha]^2 Sin[v]^2)/(1 + Aa^2) + (Sin[beta] Aa)^2/(1 + Aa^2); p = (Cos[alpha] Cos[beta])^4 + (Sin[alpha] Cos[beta])^4 + Sin[beta]^4; You = 1/(1/3/(C1 + 2*C2) - (1 - 3*p)/3/(C1 - C2) + (1 - p)/2/C3); rad = ((C1 + 2*C2)*(C1 - C2 - 2*C3)*Dir + 2*C2*C3)*You/(2*(C1 + 2*C2)*(C1 - C2)*C3); [ Last edited by onexf on 2011-3-23 at 19:50 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 大开眼界 |
» 猜你喜欢
本科新能源科学与工程，一志愿华理能动285求调剂
已经有5人回复
-
一志愿华北电力大学能动专硕，293，求调剂
已经有5人回复
-
调剂310
已经有7人回复
-
食品工程专硕一志愿中海洋309求调剂
已经有6人回复
-
275求调剂
已经有14人回复
-
283求调剂
已经有3人回复
-
085602 化工专硕 338分 求调剂
已经有12人回复
-
312，生物学求调剂
已经有3人回复
-
化学0703 调剂 306分 一志愿211
已经有4人回复
-
311（085601）求调剂
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求帮忙作图，急！！！
已经有16人回复
- 液相色谱图主峰本来很好，一个月后突然出现肩峰，求大侠帮忙，附图。
已经有49人回复
- 请帮忙分析两张色谱图
已经有4人回复
- 实验得到一张质谱图，麻烦各位大哥帮忙解析一下这，含有结构式。
已经有12人回复
- 已有核磁图，但有一个峰不知是什么，求高手帮忙，谢谢！
已经有17人回复
- 帮忙用XCrySDen出个图吧.
已经有4人回复
- 大神来帮忙解释一下，上下两个图是什么意思
已经有10人回复
- 帮忙看一下下面这个谱图是不是我要的？为什么11ppm左右有一个峰？
已经有8人回复
- 求高手帮解一个图谱，谢谢
已经有7人回复
- 帮忙分析一个HR-TEM图！
已经有11人回复
- 以前没做过，求帮忙看下DMA图
已经有8人回复
- 帮忙分析一个XRD图
已经有4人回复
- 急急急，哪位大神帮忙分析一下透镜图~~看不懂~汗一个先
已经有20人回复
- 编辑来信，让提供一个图作为TOC，这是什么意思？
已经有8人回复
- 做一个核磁图谱 ，有提供原始数据帮忙做成图！
已经有36人回复
- 请问XRD的局部放大的图是怎么做出来的？
已经有4人回复
- 用CASTEP做出来的一个能带图，大家帮忙看看它的导电性
已经有21人回复
- 【请教】帮忙分析一张N2气吸脱附的图，不闭合
已经有9人回复
- LaTex使用问题，帮忙画两个图标
已经有3人回复
- 【求助】请高手帮忙分析一下这张固化DSC曲线图，非常感谢
已经有18人回复
- 【求助】谁能帮忙分析一张多肽二级质谱图，必有重谢！
已经有17人回复
» 抢金币啦！回帖就可以得到:
-
【招聘】温州医科大学/附属医院/瓯江实验室李校堃院士团队诚聘博士后和科研人员
+1/190
-
三峡大学材料与化工学院环境功能材料课题组招生硕士调剂生
+1/160
-
山东建筑大学 （市政与环境工程学院）硕士研究生招生啦 静电纺丝纤维与盐水分离应用
+1/88
-
青岛科技大学生物工程学院招调剂硕士研究生
+1/71
-
招聘︱国家纳米科学中心丁宝全杰青团队诚聘博士后（特别研究助理）
+2/42
-
欢迎08或部分07开头的有缘人报考
+1/30
-
双一流大学湘潭大学“化工过程模拟与强化”国家地方联合工程研究中心招收各类博士生
+1/16
-
食品学硕考研调剂
+1/16
-
烟台大学先进陶瓷方向 招收材料专业调剂研究生
+1/11
-
西华大学材料学院表面科学与工程技术科研团队2026年招收研究生
+1/11
-
韩国庆熙大学 Yu Taekyung 教授招2026年博士研究生
+1/10
-
诚挚招生（硕士+博士）
+1/9
-
材料科学与工程学院招生
+1/8
-
武汉理工大学船海与能源动力工程学院 刘鹏飞院士团队专任教师招聘启事
+1/5
-
三峡大学材化学院磷硅精细化工课题组招收化学类、材料类、材料与化工类硕士
+1/5
-
中国科学技术大学纳米科学技术学院张振教授/江雷院士诚招2026级博士
+1/3
-
26年环境博士申请
+1/3
-
哈工大韩飞课题组招收3D微纳加工、光学工程、纳米化学等方向博士生及博后（长期有效）
+1/3
-
上海师范大学有机化学专业研究生招生
+1/2
-
2026南京林业大学双一流建设高校——化学工程学院招收化学/材料相关专业硕士生
+1/2
★ ★ ★ ★ ★ ★
onexf(金币+1):谢谢参与
cenwanglai(金币+5): 谢谢，呵呵，很有耐心。等lz给你bb哈！ 2011-03-23 23:52:26
onexf(金币+30): 多谢，麻烦继续给力！ 2011-03-24 19:53:57
ben_ladeng: 2011-09-05 20:55:54
onexf(金币+1):谢谢参与
cenwanglai(金币+5): 谢谢，呵呵，很有耐心。等lz给你bb哈！ 2011-03-23 23:52:26
onexf(金币+30): 多谢，麻烦继续给力！ 2011-03-24 19:53:57
ben_ladeng: 2011-09-05 20:55:54
|
源程序如下: C1 = 120; C2 = 100; C3 = 60; Aa = -(Cos[alpha] Cos[v] + Sin[alpha] Sin[v])/Tan[beta]; Dir = Cos[ beta]^2 (Cos[alpha]^2 Cos[v]^2 + Sin[alpha]^2 Sin[v]^2)/(1 + Aa^2) + (Sin[beta] Aa)^2/(1 + Aa^2); p = (Cos[alpha] Cos[beta])^4 + (Sin[alpha] Cos[beta])^4 + Sin[beta]^4; You = 1/(1/3/(C1 + 2*C2) - (1 - 3*p)/3/(C1 - C2) + (1 - p)/2/C3); rad = ((C1 + 2*C2)*(C1 - C2 - 2*C3)*Dir + 2*C2*C3)* You/(2*(C1 + 2*C2)*(C1 - C2)*C3); Plot3D[NMaximize[{rad, 0 < v < 2 \[Pi]}, v][[1]], {alpha, 0, 2 \[Pi]}, {beta, -\[Pi]/2, \[Pi]/2}] Plot3D[NMinimize[{rad, 0 < v < 2 \[Pi]}, v][[1]], {alpha, 0, 2 \[Pi]}, {beta, -\[Pi]/2, \[Pi]/2}] Quit[] 计算结果图:  就是算得很慢 |
2楼2011-03-23 23:45:40
★
onexf(金币+1):谢谢参与
onexf(金币+1):谢谢参与
|
本帖内容被屏蔽 |
4楼2011-03-24 08:21:38
6楼2011-03-24 08:47:20
8楼2011-03-24 19:52:58
那位虫子教教我？最近比较忙，脑筋有点不够用了
9楼2011-03-29 22:13:09
★ ★ ★
onexf(金币+1):谢谢参与
微尘、梦想(金币+2): 谢谢回复！ 2011-05-01 18:25:08
onexf(金币+30): xx! 2011-05-02 13:13:59
onexf(金币+1):谢谢参与
微尘、梦想(金币+2): 谢谢回复！ 2011-05-01 18:25:08
onexf(金币+30): xx! 2011-05-02 13:13:59
|
I have writen the following code, and the singluar points are disppeared. Try again, Plz! C11 = 120; C12 = 100; C44 = 60; n1 = Cos[a] Cos; n2 = Sin[a] Cos; n3 = Sin; m1 = Cos[v] Cos[w]; m2 = Sin[v] Cos[w]; m3 = Sin[w]; AA = - (Cos[a] Cos[v] + Sin[a] Sin[v])/Tan; Dir = Cos[ b]^2 (Cos[a]^2 Cos[v]^2 + Sin[a]^2 Sin[v]^2)/(1 + AA^2) + (Sin AA)^2/(1 + AA^2); p = n1^4 + n2^4 + n3^4; you = 1/(1/3/(C11 + 2*C12) - (1 - 3*p)/3/(C11 - C12) + (1 - p)/2/C44); ((C11 + 2*C12)*(C11 - C12 - 2*C44)*Dir + 2*C12*C44)/(2*(C11 + 2*C12)*(C11 - C12)*C44)*you; radi = ((C11 + 2*C12)*(C11 - C12 - 2*C44)*Dir + 2*C12*C44)* you/(2*(C11 + 2*C12)*(C11 - C12)*C44); ParametricPlot3D[{Abs[ NMaximize[{radi, 0 + 0.01 < v < \[Pi] - 0.01}, v][[1]]]*n1, Abs[NMaximize[{radi, 0 + 0.01 < v < \[Pi] - 0.01}, v][[1]]]*n2, Abs[NMaximize[{radi, 0 < v < \[Pi]}, v][[1]]]*n3}, {a, 0, 2 \[Pi]}, {b, -\[Pi]/2 + 0.01, \[Pi]/2}, PlotPoints -> 40]  |
10楼2011-05-01 13:27:47
都这么厉害啊，11楼2011-06-06 16:50:11
12楼2011-06-09 15:36:51
13楼2011-09-05 11:08:19
16楼2011-09-06 03:05:50
17楼2012-02-07 11:03:42
简单回复
2011-03-24 00:00
onexf(金币+1):谢谢参与
牛
k-li5楼
2011-03-24 08:41
onexf(金币+1):谢谢参与
2011-03-24 08:47
lxzhq14楼
2011-09-05 12:30
onexf(金币+1):谢谢参与
八戒的小弟15楼
2011-09-05 17:23
onexf(金币+1):谢谢参与
