| 查看: 1534 | 回复: 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 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 大开眼界 |
» 猜你喜欢
食品工程专硕一志愿中海洋309求调剂
已经有8人回复
-
求调剂
已经有7人回复
-
311（085601）求调剂
已经有4人回复
-
305求调剂
已经有4人回复
-
一志愿北化085600材料专硕275|有文章专利｜求调剂
已经有7人回复
-
一志愿华理，数一英一285求A区调剂
已经有12人回复
-
289求调剂
已经有11人回复
-
生物学学硕，一志愿湖南大学，初试成绩338
已经有6人回复
-
321求调剂
已经有7人回复
-
343求调剂
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求帮忙作图，急！！！
已经有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人回复
» 抢金币啦！回帖就可以得到:
-
基金祈福
+3/1679
-
【招聘】温州医科大学/附属医院/瓯江实验室李校堃院士团队诚聘博士后和科研人员
+1/194
-
耐热氧化耐热黄变效果好的聚醚多元醇制备的聚氨酯都有哪些潜在应用领域？
+2/148
-
多通道数据采集/线性功率放大器租赁｜自动化控制实验
+1/95
-
焊接专业应届毕业生招聘
+1/83
-
艾里卡特 (Alicat) - 燃料电池测试系统与质量流量和压力仪表
+3/49
-
济南大学水利与环境学院赵艳侠简介
+1/42
-
天津城建大学招收物理、化学或者生物调剂生一名，前往天津大学联合培养
+1/38
-
浙江省人民医院博士后岗位招聘
+1/37
-
【青岛大学】2026年生物与医药申请考核制博士生招生（含少数民族骨干人才）
+1/21
-
齐齐哈尔大学化学与化学工程学院招收专硕调剂
+1/12
-
烟台大学先进陶瓷方向 招收材料专业调剂研究生
+1/11
-
中国科学院上海硅酸盐研究所招聘科研助理（纳米材料与生物应用方向）
+1/7
-
南京林业大学化学工程学院柏惺峰课题组诚聘科研助理、水杉英才E类/F类青年教师或博后
+1/6
-
【博士招生】广东工业大学国家优青课题组招收2026年环境方向博士生
+1/5
-
欢迎各位同学申请参加澳门大学健康科学学院的2026年夏令营，报名截止4月12号
+1/4
-
福建师范大学国家级人才团队招收2026年博士研究生
+1/4
-
招收2026级研究生
+1/4
-
26年招收07 08学术型/专业型硕士研究生
+1/3
-
苏州高校招收08方向调剂
+1/3
★ ★ ★ ★ ★ ★
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):谢谢参与
