24小时热门版块排行榜    

CyRhmU.jpeg
查看: 3703  |  回复: 19
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

ylzhang508

金虫 (小有名气)

[交流] 【求助】复超越方程的解法已有11人参与

各位大侠,我在解金属波导时解出来的色散关系是复数的超越方程,不知道怎么才能求解,请各位不吝赐教!方程如下:

这里,d是已知的数值,比如50nm;eps_d是已知的介质的介电常数,为实数,比如为1.7。eps_m是金属的介电常数,是电磁场频率的复值函数
kd和km满足
整个方程,对每一个频率有一组(beta,w)的解。因为eps_m(w)是复值的,所以beta也是复值函数,。我想知道怎么解这个方程啊?在文献里看到说是用的Nelder-Mead单纯形求最小值方法解的,请大家给我提一些建议吧。
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

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

ylzhang508

金虫 (小有名气)

引用回帖:
Originally posted by 小木虫 at 2010-05-26 21:41:55:
eps_m=eps_m(w)具体是什么,另外给我一个beta的特解(当然也要给出特解时的参数值)

多谢版主。eps_m是光波段的金属Ag的介电常数

其中,eps_inf=3.7,w_p=1.38e16(Hz),gamma=2.73e13(Hz),因此对于光波波长为400-1200 (nm)的光来说,w=2*pi*f=2*pi*c/lambda=6.28*[0.25~0.75]e15[Hz],对于给定的角频率来说,这不过是一个复数。

另外,我前面给出的公式有点问题,根号下那两个eps_d,eps_m后面要乘以(w/c)^2,其中c是光速,见下面


前面的公式d=50 (nm)最好用国际单位制的米表示d=5e-8.

文献中的曲线在这幅图里面,灰色的曲线就是我想画的,a--->反对称,s--->对称

这幅图的横坐标是beta的实部,纵坐标实际是角频率,作者用能量单位表示的,把频率乘以普朗克常数就是能量值。

[ Last edited by ylzhang508 on 2010-5-26 at 22:07 ]
10楼2010-05-26 22:01:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 20 个回答

kuangfg1987

至尊木虫 (职业作家)


ylzhang508(金币+1):还是谢谢你 2010-05-21 21:00:53
tongling921(金币+1):确实 2010-05-26 20:57:58
用matlab解,不过我也不太清楚
2楼2010-05-21 17:34:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ylzhang508

金虫 (小有名气)

谢谢你啊。自己顶一下,请大家多多指点
3楼2010-05-24 16:43:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tower598

铁杆木虫 (职业作家)

上网达人

ylzhang508(金币+2):3q 2010-05-26 19:57:43
fortran的二分法应该能解出来
4楼2010-05-26 19:24:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复(可上传附件)
信息提示
请填处理意见