| 查看: 4316 | 回复: 19 | |||
[交流]
【求助】复超越方程的解法 已有11人参与
|
|||
各位大侠，我在解金属波导时解出来的色散关系是复数的超越方程，不知道怎么才能求解，请各位不吝赐教！方程如下： 这里,d是已知的数值，比如50nm;eps_d是已知的介质的介电常数，为实数,比如为1.7。eps_m是金属的介电常数，是电磁场频率的复值函数  kd和km满足  整个方程，对每一个频率有一组(beta,w)的解。因为eps_m(w)是复值的，所以beta也是复值函数,  。我想知道怎么解这个方程啊？在文献里看到说是用的Nelder-Mead单纯形求最小值方法解的，请大家给我提一些建议吧。 |
回复此楼» 猜你喜欢
欢迎发来filecode的Mz6后的代码验证其规律
已经有123人回复
-
今天基金会出结果吗？20260819
已经有13人回复
-
时间戳变了，能看出什么问题？
已经有16人回复
-
filecode，4个jtjc了
已经有13人回复
-
为什么余额宝的年化利率越来越低？主要原因有哪些？
已经有5人回复
-
今天放榜没戏了吧
已经有5人回复
-
今天放榜吗？
已经有15人回复
-
filecode=后面第一个是大写字母
已经有12人回复
-
2027广东省杰青
已经有10人回复
-
朋友圈看到的
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 请问我求解的数值解合理吗？谁能帮我验证一下？
已经有7人回复
- 求助：超越方程的解法
已经有13人回复
- 求助：如何求含贝塞尔函数的超越方程的复数解
已经有3人回复
- 求不等式方程的解法
已经有11人回复
- 高数关于不显含x和y的微分方程求解？
已经有8人回复
- 偏微分方程组的求解
已经有3人回复
- matlab如何求解多元多次方程，包括绘图
已经有7人回复
- 高斯方程积分编程matlab或C语言
已经有15人回复
- 用matlab求解方程
已经有6人回复
- 2阶偏微分方程的解法
已经有11人回复
- fortran循环求助
已经有3人回复
- 用fortran程序遗传算法解非线性方程组
已经有7人回复
- 这个方程组怎么解啊
已经有14人回复
- 求三角函数的最大值···
已经有7人回复
- 如何求解下列方程：xlnx=a(x-b)
已经有14人回复
- 各位朋友，如何用matlab编写这个方程，并求解？
已经有15人回复
- 函数求导后，不会解，求助大侠！
已经有5人回复
- 【求助】帮忙求解一个一元四次方程
已经有9人回复
- 求各位前辈帮忙！matlab解复数方程，出现图中的警告，怎来修改呢？谢谢了啊
已经有5人回复
- MATLAB解方程出现??? Subscripted assignment dimension mismatch.
已经有3人回复
3楼2010-05-24 16:43:17
小木虫:)
荣誉版主 (著名写手)
- 应助: 0 (幼儿园)
- 贵宾: 2.542
- 金币: 9239.2
- 散金: 130
- 红花: 22
- 帖子: 1211
- 在线: 343.9小时
- 虫号: 800712
- 注册: 2009-06-29
- 性别: GG
- 专业: 基础物理学
- 管辖: 物理
7楼2010-05-26 21:16:28
kuangfg1987
至尊木虫 (职业作家)
- 应助: 0 (幼儿园)
- 金币: 21318.2
- 散金: 30
- 红花: 25
- 帖子: 3028
- 在线: 405小时
- 虫号: 891745
- 注册: 2009-11-02
- 专业: 原子和分子物理
2楼2010-05-21 17:34:45
4楼2010-05-26 19:24:38
小木虫:)
荣誉版主 (著名写手)
- 应助: 0 (幼儿园)
- 贵宾: 2.542
- 金币: 9239.2
- 散金: 130
- 红花: 22
- 帖子: 1211
- 在线: 343.9小时
- 虫号: 800712
- 注册: 2009-06-29
- 性别: GG
- 专业: 基础物理学
- 管辖: 物理
5楼2010-05-26 19:40:59
6楼2010-05-26 19:57:09
|
哇，谢谢您，不过我现在没装mathematica， 还不能验证。不过您这个写法好像只解一个单beta的方程，实际上我的方程里面eps_m是一个复值的函数eps_m=eps_m(w)，因此这里有两个变量频率w和beta。我希望得到的解是一条曲线，横坐标是Re(beta)，纵坐标是w，这里w是频率。在文献里说是用matlab里面的helder mead最优化方法求解的，我怀疑是用fminsearch这个命令，但是fminsearch只能处理输入参数是实数的情况，所以...... [ Last edited by ylzhang508 on 2010-5-26 at 21:41 ] |
8楼2010-05-26 21:28:10
小木虫:)
荣誉版主 (著名写手)
- 应助: 0 (幼儿园)
- 贵宾: 2.542
- 金币: 9239.2
- 散金: 130
- 红花: 22
- 帖子: 1211
- 在线: 343.9小时
- 虫号: 800712
- 注册: 2009-06-29
- 性别: GG
- 专业: 基础物理学
- 管辖: 物理
血饮玄冥9楼2010-05-26 21:41:55
|
多谢版主。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
