| 查看: 1540 | 回复: 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 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 大开眼界 |
» 猜你喜欢
0703化学调剂，求导师收
已经有10人回复
-
各位老师好，我的一志愿为北京科技大学085601材料专硕
已经有5人回复
-
一志愿太原理工安全工程300分，求调剂
已经有6人回复
-
0856求调剂
已经有4人回复
-
11408软件工程求调剂
已经有3人回复
-
一志愿211 初试270分 求调剂
已经有8人回复
-
300求调剂，材料科学英一数二
已经有6人回复
-
学硕274求调剂
已经有9人回复
-
本科双非材料，跨考一志愿华电085801电气，283求调剂，任何专业都可以
已经有7人回复
-
317分 一志愿南理工材料工程 本科湖工大 求调剂
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求帮忙作图，急！！！
已经有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/288
-
【招聘】温州医科大学/附属医院/瓯江实验室李校堃院士团队诚聘博士后和科研人员
+1/190
-
武汉工程大学化环学院程老师课题组招收化学类学硕、化学工程与环境工程专硕调剂
+2/74
-
欢迎报考赣南师范大学智能制造与未来能源学院
+2/64
-
首都师范大学招收材料与化工专硕调剂生--中科院半导体所共培
+1/40
-
中山大学生医工学院课题组招硕士生（纳米材料/高分子/分析传感/荧光免疫检测背景）
+2/38
-
韩国庆熙大学 Yu Taekyung 教授招2026年博士研究生
+1/10
-
中国地震局地质研究所+地球物理专业+接收调剂
+1/10
-
齐齐哈尔大学化学与化学工程学院招收专硕调剂
+1/10
-
西安工程大学环境与化学工程学院能源环境催化与先进材料课题组招生
+1/9
-
2026年中国科学院城市环境所与福建农林大学联合培养硕士研究生项目介绍
+1/6
-
墨尔本大学伯明翰大学联合全奖国际博士生
+1/5
-
中国科学院大学纳米科学与技术学院-罗聃课题组诚招硕士研究生
+1/5
-
计算机(11408)考研361分求调剂
+1/4
-
26博士申请
+1/4
-
杭州电子科技大学 磁电功能材料创新团队 招收硕士调剂生
+1/4
-
环境微生物与污染控制研究所招收2026硕士研究生
+1/3
-
【MP期刊】青年编委招募开启！加入我们，一起推动材料科学前沿
+1/2
-
【博士招生】上海交通大学机动学院-招2026年秋“申请-考核制”博士生1名-3D打印方向
+1/1
-
北京石油化工学院康举、张慧杰老师团招2026年研究生复试调剂考生，欢迎联系
+1/1
★ ★ ★ ★ ★ ★
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):谢谢参与
