| 查看: 420 | 回复: 1 | ||
[求助]
python
|
| 考虑一个三层对称平板波导,包层材料为二氧化硅,芯层材料为硅。计算波导TE基膜的模场宽度(E场的半高全宽)随波导厚度变化的趋势曲线。 |
» 猜你喜欢
售SCI一区文章,我:8 O5 51O 54,科目齐全
已经有4人回复
售SCI一区文章,我:8 O5 51O 54,科目齐全
已经有4人回复
售SCI一区文章,我:8 O5 51O 54,科目齐全
已经有5人回复
售SCI一区文章,我:8 O5 51O 54,科目齐全
已经有4人回复
售SCI一区文章,我:8 O5 51O 54,科目齐全
已经有6人回复
售SCI一区文章,我:8 O5 51O 54,科目齐全
已经有6人回复
售SCI一区文章,我:8 O5 51O 54,科目齐全
已经有6人回复
售SCI一区文章,我:8 O5 51O 54,科目齐全
已经有7人回复
售SCI一区文章,我:8 O5 51O 54,科目齐全
已经有7人回复
售SCI一区文章,我:8 O5 51O 54,科目齐全
已经有7人回复
» 本主题相关商家推荐: (我也要在这里推广)
|
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() 这是我写的,希望指点。。。。。 |
2楼2016-11-26 22:42:11













回复此楼
10