| 查看: 1545 | 回复: 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 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 大开眼界 |
» 猜你喜欢
299求调剂
已经有8人回复
-
一志愿北京理工大学本科211材料工程294求调剂
已经有6人回复
-
300求调剂，材料科学英一数二
已经有8人回复
-
招收生物学/细胞生物学调剂
已经有5人回复
-
070305高分子化学与物理 304分求调剂
已经有7人回复
-
289求调剂
已经有13人回复
-
一志愿哈尔滨工业大学材料与化工方向336分
已经有9人回复
-
081200-11408-276学硕求调剂
已经有6人回复
-
调剂求院校招收
已经有5人回复
-
调剂310
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求帮忙作图，急！！！
已经有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/268
-
天津科技大学-（新晋院士团队）先进纤维与纸基功能材料团队B+专业
+1/93
-
焊接专业应届毕业生招聘
+1/83
-
青岛农业大学海洋学院招调剂硕士研究生
+1/73
-
青岛大学固态电池研究团队招收材料、化学、物理等方向博士研究生
+2/46
-
【全额奖学金+英国博士学位】国科大杭高院Nathan课题组招收2026级联培博士生1名！
+1/44
-
南通大学国家级人才团队招收卫生检验、医学检验相关方向调剂研究生
+2/42
-
西华大学材料学院表面科学与工程技术科研团队2026年招收研究生
+1/33
-
南京林业大学-国家级青年人才团队 招2026级博士、硕士（合成化学、植物化学方向）
+1/20
-
赣南师范大学智能制造与未来能源学院2026年硕士调剂信息
+1/16
-
烟台大学精准材料高等研究院26年材料方向研究生招生
+1/11
-
26年申博自荐-计算机视觉
+1/10
-
食品调剂-总分291-有中核论文导师一作、本人二作-一志愿南昌大学学硕-介绍详细
+1/8
-
国家双一流高校-国家级青年人才课题组博士招生
+1/7
-
诚聘博士科研人员（有机合成化学方向）
+1/7
-
同济大学环境学院 肖倩特聘研究员课题组 招聘硕士/博士（长期有效）
+1/6
-
上海交通大学陈倩栎课题组招收2026年申请-考核博士生1名（能源材料方向）
+1/5
-
5科大招收材料专业硕士
+1/3
-
香港中文大学 生物医学学院/组织工程与再生医学研究所 招聘博士后/研究助理
+1/2
-
北京石油化工学院康举、张慧杰老师团招2026年研究生复试调剂考生，欢迎联系
+1/1
★ ★ ★
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
★ ★ ★ ★ ★ ★
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
8楼2011-03-24 19:52:58
简单回复
2011-03-24 00:00
onexf(金币+1):谢谢参与
牛
