24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

[交流] 国庆快乐

回复此楼

» 猜你喜欢

参与限项已经有3人回复
假如你的研究生提出不合理要求已经有7人回复
实验室接单子已经有4人回复
全日制（定向）博士已经有4人回复
对氯苯硼酸纯化已经有3人回复
求助:我三月中下旬出站，青基依托单位怎么办？已经有12人回复
不自信的我已经有12人回复
所感已经有4人回复
要不要辞职读博？已经有7人回复
北核录用已经有3人回复

» 抢金币啦！回帖就可以得到:

查看全部散金贴

1楼 2016-10-01 10:08:17

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

蓝青Ching

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 10.5
帖子: 10
在线: 3小时
虫号: 4631710

★
小木虫: 金币+0.5, 给个红包，谢谢回帖

老师，请问能不能帮我解决一下python的问题？
考虑一个三层对称平板波导，包层材料为二氧化硅，芯层材料为硅。计算波导TE基膜的模场宽度（E场的半高全宽）随波导厚度变化的趋势曲线。
这是我写的，希望指点。。。。。

from camfr import *
from matplotlib.pyplot import *
from numpy import *

set_lambda(1.55)
set_N(1)
set_polarisation(TE)

sio2=Material(1.45)
silicon=Material(3.46)
air=Material(1.0)

w2=arange(0.01,6,0.8)
transmission=[]

for tw in w2:
slab=Slab(sio2(4.995)+silicon(0.01)+sio2(4.995))
slab=Slab(sio2(5-w2/2)+silicon(w2)+sio2(5-w2/2)

x=arrange(0,tw,0.8)

field1=[slab.mode(0).field((coord(tx,0,0).E1() for tx in x]
field2=[slab.mode(0).field((coord(tx,0,0)).E1()*(slab.n(Coord(tx,0,0)))**2 for tx in x]
index=[slab.n((coord(tx,0,0)) for tx in x]

M=max(field)
y=M/2
[m,n]=find(abs(M-y)<0.00001)
Q=m-n

plot(x,real(Q))
plot(x,absolute(index))
show()