应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » 帮我看一下这个程序啊为什么输不出来，菜鸟求问
81/1返回列表
查看: 499  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

Abla

新虫 (初入文坛)

[求助] 帮我看一下这个程序啊为什么输不出来，菜鸟求问

\[Epsilon] =
  1 - (2*Pi*12.8*10^9*2*Pi*12.8*10^9 -
    2*Pi*10.3*10^9*2*Pi*10.3*10^9)/(
   67.2*10^9*67.2*10^9 - 2*Pi*10.3*10^9*2*Pi*10.3*10^9 +
    I*0.02*10^9*67.2*10^9);
\[Mu] = 1 - (
   2*Pi*10.95*10^9*2*Pi*10.95*10^9 -
    2*Pi*10.05*10^9*2*Pi*10.05*10^9)/(
   67.2*10^9*67.2*10^9 - 2*Pi*10.05*10^9*2*Pi*10.05*10^9 +
    I*0.02*10^9*67.2*10^9);
\[Chi] = Sqrt[\[Epsilon]]/Sqrt[\[Mu]];
w = 67.2*10^9;
c = 3.0*10^8;
d = 0.01;
\[Lambda] = (2 Pi*c)/w;
q[\[Theta]_] =
  Sqrt[\[Epsilon]]/Sqrt[\[Mu]]*Sqrt[
   1 - (Sin[\[Theta]]*Sin[\[Theta]])/(\[Epsilon]*\[Mu])];
k[\[Theta]_] = (67.2*10^9)/(3*10^8)*Sqrt[\[Epsilon]]*Sqrt[\[Mu]]*Sqrt[
   1 - (Sin[\[Theta]]*Sin[\[Theta]])/(\[Epsilon]*\[Mu])];

r[\[Theta]_] := (-I*Cos[\[Theta]]*Cos[\[Theta]]*
    Sin[k[\[Theta]]*d]/q[\[Theta]] +
   I*q[\[Theta]]*Sin[k[\[Theta]]*d])/(
  2*Cos[\[Theta]]*Cos[k[\[Theta]]*d] -
   I*Cos[\[Theta]]*Cos[\[Theta]]*Sin[k[\[Theta]]*d]/q[\[Theta]] -
   I*q[\[Theta]]*Sin[k[\[Theta]]*d]);
t[\[Theta]_] := (2 Cos[\[Theta]])/(
  2*Cos[\[Theta]]*Cos[k[\[Theta]]*d] -
   I*Cos[\[Theta]]*Cos[\[Theta]]*Sin[k[\[Theta]]*d]/q[\[Theta]] -
   I*q[\[Theta]]*Sin[k[\[Theta]]*d]);
c[\[Theta]_] = Re[r[\[Theta]]];
d[\[Theta]_] = Im[r[\[Theta]]];
a[\[Theta]_] := D[Im[r[\[Theta]]], \[Theta]];
a1[\[Theta]_] = a[\[Theta]];
b[\[Theta]_] := D[Re[r[\[Theta]]], \[Theta]];
b1[\[Theta]] = b[\[Theta]];
c1[\[Theta]] = Re[t[\[Theta]]];
d1[\[Theta]] = Im[t[\[Theta]]];
a2[\[Theta]_] := D[Im[t[\[Theta]]], \[Theta]];
a3[\[Theta]_] = a2[\[Theta]];
b2[\[Theta]_] := D[Re[t[\[Theta]]], \[Theta]];
b3[\[Theta]_] = b2[\[Theta]];
\[CapitalDelta]r[\[Theta]_] = -(\[Lambda]/(2 Pi))*1/(
  c[\[Theta]]*c[\[Theta]] +
   d[\[Theta]]*d[\[Theta]])*(c[\[Theta]]*a1[\[Theta]] -
    d[\[Theta]]*b1[\[Theta]])
\[CapitalDelta]t[\[Theta]_] = -(\[Lambda]/(2 Pi))*1/(
  c[\[Theta]]*c[\[Theta]] +
   d[\[Theta]]*d[\[Theta]])*(c1[\[Theta]]*a3[\[Theta]] -
    d1[\[Theta]]*b3[\[Theta]])
Plot[\[CapitalDelta]r[\[Theta]], {\[Theta], 0, Pi/2}]
Plot[\[CapitalDelta]t[\[Theta]], {\[Theta], 0, Pi/2}]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-03-02 11:46:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Abla

新虫 (初入文坛)
为什么输不出来图形啊
赞一下 回复此楼
2楼2013-03-02 11:54:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Abla

新虫 (初入文坛)
直接粘贴到mathmatica上好像没有这么乱
赞一下 回复此楼
3楼2013-03-02 12:01:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

walk1997

金虫 (著名写手)
程序里面c,d各出现了2次定义
1次是做为变量 另外又是做为函数出现
把第二个c,d改一下名字吧 另外 很多的定义 cc[\theta]=xxx[\theta]也改下
赞一下 回复此楼
4楼2013-03-06 13:35:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Abla

新虫 (初入文坛)
= 1 - (2*Pi*12.8*10^9*2*Pi*12.8*10^9 -
    2*Pi*10.3*10^9*2*Pi*10.3*10^9)/(
   67.2*10^9*67.2*10^9 - 2*Pi*10.3*10^9*2*Pi*10.3*10^9 +
    I*0.02*10^9*67.2*10^9);
\[Mu] = 1 - (
   2*Pi*10.95*10^9*2*Pi*10.95*10^9 -
    2*Pi*10.05*10^9*2*Pi*10.05*10^9)/(
   67.2*10^9*67.2*10^9 - 2*Pi*10.05*10^9*2*Pi*10.05*10^9 +
    I*0.02*10^9*67.2*10^9);
\[Chi] = Sqrt[\[Epsilon]]/Sqrt[\[Mu]];

\[Lambda] = (2 Pi*3*10^8)/(67.2*10^9);
q[\[Theta]_] =
  Sqrt[\[Epsilon]]/Sqrt[\[Mu]]*Sqrt[
   1 - (Sin[\[Theta]]*Sin[\[Theta]])/(\[Epsilon]*\[Mu])];
k[\[Theta]_] = (67.2*10^9)/(3*10^8)*Sqrt[\[Epsilon]]*Sqrt[\[Mu]]*Sqrt[
   1 - (Sin[\[Theta]]*Sin[\[Theta]])/(\[Epsilon]*\[Mu])];

r[\[Theta]_] := (-I*Cos[\[Theta]]*Cos[\[Theta]]*
    Sin[k[\[Theta]]*0.01]/q[\[Theta]] +
   I*q[\[Theta]]*Sin[k[\[Theta]]*0.01])/(
  2*Cos[\[Theta]]*Cos[k[\[Theta]]*d] -
   I*Cos[\[Theta]]*Cos[\[Theta]]*Sin[k[\[Theta]]*d]/q[\[Theta]] -
   I*q[\[Theta]]*Sin[k[\[Theta]]*d]);
t[\[Theta]_] := (2 Cos[\[Theta]])/(
  2*Cos[\[Theta]]*Cos[k[\[Theta]]*0.01] -
   I*Cos[\[Theta]]*Cos[\[Theta]]*Sin[k[\[Theta]]*0.01]/q[\[Theta]] -
   I*q[\[Theta]]*Sin[k[\[Theta]]*0.01]);
c[\[Theta]_] = Re[r[\[Theta]]];
d[\[Theta]_] = Im[r[\[Theta]]];
a[\[Theta]_] := D[Im[r[\[Theta]]], \[Theta]];
a1[\[Theta]_] = a[\[Theta]];
b[\[Theta]_] := D[Re[r[\[Theta]]], \[Theta]];
b1[\[Theta]] = b[\[Theta]];
c1[\[Theta]] = Re[t[\[Theta]]];
d1[\[Theta]] = Im[t[\[Theta]]];
a2[\[Theta]_] := D[Im[t[\[Theta]]], \[Theta]];
a3[\[Theta]_] = a2[\[Theta]];
b2[\[Theta]_] := D[Re[t[\[Theta]]], \[Theta]];
b3[\[Theta]_] = b2[\[Theta]];
\[CapitalDelta]r[\[Theta]_] = -(\[Lambda]/(2 Pi))*1/(
   c[\[Theta]]*c[\[Theta]] +
    d[\[Theta]]*d[\[Theta]])*(c[\[Theta]]*a1[\[Theta]] -
     d[\[Theta]]*b1[\[Theta]]);
\[CapitalDelta]t[\[Theta]_] = -(\[Lambda]/(2 Pi))*1/(
   c[\[Theta]]*c[\[Theta]] +
    d[\[Theta]]*d[\[Theta]])*(c1[\[Theta]]*a3[\[Theta]] -
     d1[\[Theta]]*b3[\[Theta]]);
Plot[\[CapitalDelta]r[\[Theta]], {\[Theta], 0, Pi/2}]
Plot[\[CapitalDelta]t[\[Theta]], {\[Theta], 0, Pi/2}]
我把常量直接带进去了 不存在那个问题 还是输不出结果
赞一下 回复此楼
5楼2013-03-19 09:38:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Abla

新虫 (初入文坛)
引用回帖:
4楼: Originally posted by walk1997 at 2013-03-06 13:35:59
程序里面c,d各出现了2次定义
1次是做为变量 另外又是做为函数出现
把第二个c,d改一下名字吧 另外 很多的定义 cc=xxx也改下

我把常量已经带进去了 不存在那个问题 还是输不出结果
赞一下 回复此楼
6楼2013-03-19 09:38:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

walk1997

金虫 (著名写手)
如上面所提的 你第个代码  d 还是有问题 出现过d 和d[\theta]

简单看下输出 \[CapitalDelta]r[\[Theta]]/{\[Theta]->1.0}  结果就能清楚看到 为何没图形
赞一下 回复此楼
7楼2013-03-20 10:37:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Abla

新虫 (初入文坛)
引用回帖:
7楼: Originally posted by walk1997 at 2013-03-20 10:37:06
如上面所提的 你第个代码  d 还是有问题 出现过d 和d

简单看下输出 \r/{\->1.0}  结果就能清楚看到 为何没图形

d=0.01
回复此楼
8楼2013-04-12 14:51:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Abla 的主题更新
81/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +6 3dfhjxgsh7 2026-02-22 9/450 2026-02-23 07:49 by w4l55oybr1
[考博] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +4 khieu8v8m0 2026-02-22 4/200 2026-02-23 06:46 by jsjzfl
[论文投稿] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 w89i99eaeh 2026-02-22 4/200 2026-02-23 06:36 by w4l55oybr1
[博后之家] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +4 khieu8v8m0 2026-02-22 5/250 2026-02-23 06:34 by w4l55oybr1
[公派出国] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 khieu8v8m0 2026-02-22 5/250 2026-02-23 06:29 by w4l55oybr1
[考研] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +4 usprnugpzw 2026-02-21 10/500 2026-02-23 04:58 by 5jlh3qtdvx
[论文投稿] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 khieu8v8m0 2026-02-22 6/300 2026-02-23 02:08 by 5jlh3qtdvx
[考博] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +5 3dfhjxgsh7 2026-02-22 6/300 2026-02-23 02:04 by 5jlh3qtdvx
[教师之家] 版面费该交吗 +7 苹果在哪里 2026-02-22 8/400 2026-02-22 22:37 by otani
[基金申请] 基金正文30页指的是报告正文还是整个申请书 +5 successhe 2026-02-16 6/300 2026-02-22 21:38 by 山西悬空寺空悬�
[基金申请] 面上可以超过30页吧？ +4 阿拉贡aragon 2026-02-22 4/200 2026-02-22 21:22 by 山西悬空寺空悬�
[教师之家] 为什么中国大学教授们水了那么多所谓的顶会顶刊，但还是做不出宇树机器人？ +5 欢乐颂叶蓁 2026-02-21 5/250 2026-02-22 21:15 by 山西悬空寺空悬�
[论文投稿] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +4 usprnugpzw 2026-02-21 6/300 2026-02-22 19:48 by w89i99eaeh
[考研] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 3dfhjxgsh7 2026-02-22 4/200 2026-02-22 16:52 by khieu8v8m0
[找工作] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 usprnugpzw 2026-02-22 3/150 2026-02-22 16:37 by khieu8v8m0
[公派出国] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 usprnugpzw 2026-02-21 4/200 2026-02-22 16:27 by khieu8v8m0
[基金申请] “人文社科而论，许多学术研究还没有达到民国时期的水平” +4 苏东坡二世 2026-02-18 5/250 2026-02-22 16:07 by liangep1573
[基金申请] 什么是人一生最重要的？ +4 瞬息宇宙 2026-02-21 4/200 2026-02-22 11:44 by huagongfeihu
[基金申请] 今年春晚有几个节目很不错，点赞！ +11 瞬息宇宙 2026-02-16 12/600 2026-02-21 21:14 by lq493392203
[基金申请] 体制内长辈说体制内绝大部分一辈子在底层，如同你们一样大部分普通教师忙且收入低 +9 瞬息宇宙 2026-02-20 12/600 2026-02-21 10:39 by 欢乐颂叶蓁
信息提示
关闭
请填处理意见
关闭