| 查看: 456 | 回复: 7 | ||
[求助]
求帮我看看这个程序哪里出了问题啊为什么输出的都是1这条直线
|
|
nh = 2.22 nl = 1.41 c = 3*10^8 dh = (692*10^-9)/(4*nh) dl = (692*10^-9)/(4*nl) ps = Cos[Pi/6] a = Sqrt[1 - (1/(2*2.22))^2] b = Sqrt[1 - (1/(2*1.41))^2] p1 = a/2.22 p2 = b/1.41 k1[w_] = w/(3*10^8)*2.22*a k2[w_] = w/(3*10^8)*1.41*b MatrixForm[n1[w_] = {{Cos[k1[w]*dh], I*1/p1*Sin[k1[w]*dh]}, {I*p1*Sin[k1[w]*dh], Cos[k1[w]*dh]}}] MatrixForm[n2[w_] = {{Cos[k2[w]*dl], I*1/p2*Sin[k2[w]*dl]}, {I*p2*Sin[k2[w]*dl], Cos[k2[w]*dl]}}] n[w_] = n1[w].n2[w].n1[w].n2[w].n1[w].n2[w].n1[w].n2[w].n1[w].n2[ w].n1[w]; x11[w_] = Part[n[w], 1, 1]; x12[w_] = Part[n[w], 1, 2]; x21[w_] = Part[n[w], 2, 1]; x22[w_] = Part[n[w], 2, 2]; t[w_] = (2*Sqrt[3]/2)/(Sqrt[3]/2*x22[w] + Sqrt[3]/2*x11[w] - 3/4*x12[w] - x2[w]); t1[w_] = Re[t[w]] t2[w_] = Im[t[w]] st[w_] = Sqrt[t1[w]^2 + t2[w]^2]; T[w_] = st[w]^2; Plot[T[w], {w, 0, 2}, AxesOrigin -> {0, 0}] 这个程序哪里有错误也请多多指正啊，本人初学者表示很弱 [ Last edited by Abla on 2013-1-9 at 14:08 ] |
回复此楼» 猜你喜欢
投稿精细化工
已经有6人回复
-
博士读完未来一定会好吗
已经有36人回复
-
之前让一硕士生水了7个发明专利，现在这7个获批发明专利的维护费可从哪儿支出哈？
已经有10人回复
-
博士申请都是内定的吗？
已经有9人回复
-
心脉受损
已经有8人回复
-
读博
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
一个gjf文件，总是pending,各位帮我看看是什么问题？
已经有16人回复- PCR好多次了，一直没结果，麻烦各位帮我看看是哪出问题了？
已经有15人回复
- 求助matlab一程序运行的问题，帮忙看看哪里不对
已经有4人回复
- 麻烦塔菲尔曲线高手帮我看看这几个图，说明什么问题，或者能从中得出什么信息？
已经有12人回复
- 我编的Simpson积分法fortran程序给不出结果，大侠们看看哪里出了问题？
已经有4人回复
- 请有经验的虫友帮我看看SCI期刊IEEE TRANSACTIONS ON PLASMA SCIENCE的费用支付问题
已经有4人回复
- 帮忙看一下这个席夫碱的合成，问题出在哪里？多谢了！
已经有7人回复
- 请教高手，帮我看下贴中图片的水体是什么问题，谢谢
已经有19人回复
- 【求助】大家帮忙给看看，到底是哪里出了问题？谢谢了！！！
已经有20人回复
- 【问题求助】说说我的心里状态，大家帮我看看是不是该去看心理医生了？
已经有28人回复
- 【求助】求UDF高手帮忙看看我的问题
已经有18人回复
- 【求助】谁能帮我看看这张豚鼠皮肤病理切片,能说明什么问题吗?帮忙解读一下吧
已经有4人回复
- 【求助】大家帮我看看这护肤品成分有没有问题
已经有7人回复
- 大家帮我看看这3个Offer选择问题
已经有55人回复
- 哪位同学帮忙试试这个网页里文字能不能复制，看我电脑有没有问题
已经有6人回复
- 【问题求助】帮我看看这种情况传不传染？
已经有11人回复
- 【求助】能不能帮我看看这份配方有什么问题？
已经有23人回复
- 【讨论】大伙帮满看看这个TG-DTA基线文件的问题
已经有12人回复
2楼2013-01-09 14:20:36
mshwangg
至尊木虫 (正式写手)
- 程序强帖: 5
- 应助: 206 (大学生)
- 金币: 10702.8
- 散金: 100
- 红花: 19
- 帖子: 597
- 在线: 195.4小时
- 虫号: 576702
- 注册: 2008-06-21
- 专业: 物理学I
【答案】应助回帖
★ ★ ★
感谢参与，应助指数 +1
jjdg: 金币+1, 感谢参与 2013-01-10 10:38:19
Abla: 金币+2 2013-01-10 14:29:57
感谢参与，应助指数 +1
jjdg: 金币+1, 感谢参与 2013-01-10 10:38:19
Abla: 金币+2 2013-01-10 14:29:57
|
从现在的程序上看，这一句 t[w_] = (2*Sqrt[3]/2)/(Sqrt[3]/2*x22[w] + Sqrt[3]/2*x11[w] - 3/4*x12[w] - x2[w]); 中存在未定义的x2，会导致画不出图，检查一下是否有错误。 将x2改成x21或者x22之后，plot可以画图是一条直线。此时说明程序语法上没有问题，为什么是一条直线需要你从程序逻辑上和众多参数和函数定义上考虑是否有问题。 这样的问题非本领域的人几乎无法帮你了，因为涉及到专业的知识未必都知道。 |
3楼2013-01-10 06:34:00
walk1997
金虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 4676.2
- 红花: 22
- 帖子: 1066
- 在线: 798.1小时
- 虫号: 416039
- 注册: 2007-06-29
- 性别: GG
- 专业: 粒子物理学和场论
★ ★
jjdg: 金币+2, 辛苦了 2013-01-10 10:38:44
jjdg: 金币+2, 辛苦了 2013-01-10 10:38:44
|
这个程序里面 你仔细看下 你的输入 n1[w] n2[w] n[w] 在小w情况下 差不多就是个单位矩阵 你最后算出来的自然差不多是个直线 把画的范围改下就好 不过 估计你输入的参数可能有点问题 －-------------------------------------------------------------------------- Clear["Global`*"] nh = 2.22 nl = 1.41 c = 3*10^8 dh = (692*10^-9)/(4*nh) dl = (692*10^-9)/(4*nl) ps = Cos[Pi/6] a = Sqrt[1 - (1/(2*2.22))^2] b = Sqrt[1 - (1/(2*1.41))^2] p1 = a/2.22 p2 = b/1.41 k1[w_] := w/(3*10^8)*2.22*a k2[w_] := w/(3*10^8)*1.41*b n1[w_] := {{Cos[k1[w]*dh], I*1/p1*Sin[k1[w]*dh]}, {I*p1*Sin[k1[w]*dh], Cos[k1[w]*dh]}} n2[w_] := {{Cos[k2[w]*dl], I*1/p2*Sin[k2[w]*dl]}, {I*p2*Sin[k2[w]*dl], Cos[k2[w]*dl]}} n[w_] := n1[w].n2[w].n1[w].n2[w].n1[w].n2[w].n1[w].n2[w].n1[w].n2[ w].n1[w]; x11[w_] := Part[n[w], 1, 1]; x12[w_] := Part[n[w], 1, 2]; x21[w_] := Part[n[w], 2, 1]; x22[w_] := Part[n[w], 2, 2]; t[w_] := (2*Sqrt[3]/2)/(Sqrt[3]/2*x22[w] + Sqrt[3]/2*x11[w] - 3/4*x12[w] - x21[w]); t1[w_] := Re[t[w]] t2[w_] := Im[t[w]] st[w_] := Sqrt[t1[w]^2 + t2[w]^2]; T[w_] := st[w]^2; Plot[T[w], {w, 1*10^16, 1*10^16 + 10^15}] －---------------------------------------------------------------------------- |
4楼2013-01-10 09:22:51
walk1997
金虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 4676.2
- 红花: 22
- 帖子: 1066
- 在线: 798.1小时
- 虫号: 416039
- 注册: 2007-06-29
- 性别: GG
- 专业: 粒子物理学和场论
5楼2013-01-10 09:23:54
6楼2013-01-10 09:33:15
7楼2013-01-10 09:39:30
我把参数调大了好几个数量级输出来了，怎么木有办法给你赏钱哦O(∩_∩)O哈哈哈~8楼2013-01-10 09:44:39