24小时热门版块排行榜

返回列表

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

[交流] 国庆快乐

回复此楼

» 猜你喜欢

售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急已经有3人回复
浙江师范大学是怎么坑我的已经有7人回复
求教各位大神：2026教育部人文社科青年基金项目何时公示呀？已经有6人回复
面上有专家说实验设备不是我们单位的已经有8人回复
科研求助已经有4人回复
博士去二本高校当辅导员，如何调整心态？已经有11人回复
江苏省自然基金什么时候出结果已经有6人回复
闲聊已经有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()