| 查看: 3072 | 回复: 9 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
新人求助解MATLAB超越方程，方程中含有复数已有2人参与
|
|||
mathematica和matlab,我都试着写了方程，是用迭代初值的方法，可是结果总是随着初值变化，我也试了先画图，但是也无法画出来，耽误了好几天了，真的好烦躁 。1stopt这个软件我也没得。。。。一共就5个金币了，能回答的大神帮帮忙吧！拜托了！这是的mathematica的程序 n1 = 1.445; n2 = 3.445; n3 = 0.1453 + I*11.3587; n4 = 1; h = 100*10^-9; lamada = 1550*10^-9; k0 = 2*Pi/lamada; b = 70*10^-9; n = 1; N0 = 0.1453 + I*11.3587; FindRoot[n*Pi + ArcTan[Sqrt[N1^2 - n2^2]/Sqrt[n1^2 - N1^2]] + ArcTan[Sqrt[ N1^2 - n3^2]*(1 + (Sqrt[N1^2 - n4^2] - Sqrt[N1^2 - n3^2])* Exp[-2*k0*Sqrt[N1^2 - n3^2]*h]/(Sqrt[N1^2 - n4^2] + Sqrt[N1^2 - n3^2]))/(1 - (Sqrt[N1^2 - n4^2] - Sqrt[N1^2 - n3^2])* Exp[-2*k0*Sqrt[N1^2 - n3^2]*h]/(Sqrt[N1^2 - n4^2] + Sqrt[N1^2 - n3^2]))/Sqrt[n1^2 - N1^2]] - k0*Sqrt[(n1^2 - N1^2)]*b == 0, {N1, N0}, WorkingPrecision -> 24]  N16Y7`CLCLKS{PN)53U@6JX.png |
回复此楼» 猜你喜欢
请问哪里可以有青B申请的本子可以借鉴一下。
已经有4人回复
-
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗
已经有6人回复
-
孩子确诊有中度注意力缺陷
已经有14人回复
-
三甲基碘化亚砜的氧化反应
已经有4人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
请问有评职称，把科研教学业绩算分排序的高校吗
已经有5人回复
-
2025冷门绝学什么时候出结果
已经有3人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
-
AI论文写作工具：是科研加速器还是学术作弊器？
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
方程变换形式时matlab求出的解不同
已经有3人回复- matlab对含贝塞尔函数的超越方程求根
已经有5人回复
- 求助：超越方程的解法
已经有13人回复
- 求助：如何求含贝塞尔函数的超越方程的复数解
已经有3人回复
- matlab如何求解多元多次方程，包括绘图
已经有7人回复
- matlab解方程组求助攻
已经有4人回复
- matlab解一元三次方程，得到的都是复数根。
已经有7人回复
- 用matlab求解非线性方程组说无解，一定是方程组本身无解，还是有可能程序有问题呢？
已经有11人回复
- 求助，用matlab解微分方程组，希望高手能给予指点
已经有17人回复
- matlab解偏微分方程求助
已经有12人回复
- 各位朋友，如何用matlab编写这个方程，并求解？
已经有15人回复
- 求助 MATLAB解方程组-fslove
已经有7人回复
- 求助matlab解非线性方程。求代码
已经有5人回复
- 【求助】帮忙求解一个一元四次方程
已经有9人回复
- 求各位前辈帮忙！matlab解复数方程，出现图中的警告，怎来修改呢？谢谢了啊
已经有5人回复
- 用Matlab求解方程
已经有6人回复
- MATLAB解方程出现??? Subscripted assignment dimension mismatch.
已经有3人回复
- 【求助】matlab 求解微分方程中的未知参数
已经有20人回复
- 【求助】matlab 解矩阵方程 重金悬赏解决方案
已经有10人回复
- 【求助】求助用matlab解偏微分方程
已经有15人回复
- 【求助】复超越方程的解法
已经有19人回复
- 【求助】求教matlab解非线性方程组
已经有9人回复
eidie
木虫之王 (文坛精英)
- 应助: 15 (小学生)
- 贵宾: 1.016
- 金币: 71318.4
- 散金: 8785
- 红花: 240
- 沙发: 1
- 帖子: 13887
- 在线: 1583.1小时
- 虫号: 2259480
- 注册: 2013-01-24
- 专业: 应用光学
【答案】应助回帖
感谢参与，应助指数 +1
|
patpat~ 不要哭。 %%%%%%%%%%%%%%%%%%% %不知道对不对，妹妹先试试看 clear; close; clc; n1 = 1.445; n2 = 3.445; n3 = 0.1453 + i*11.3587; n4 = 1; h = 100*10^-9; lamada = 1550*10^-9; k0 = 2*pi/lamada; b = 70*10^-9; n = 1; N0 = 0.1453 + i*11.3587; syms N1 delta3=(1+((N1^2-n4^2)^(1/2)-(N1^2-n3^2)^(1/2))/((N1^2-n4^2)^(1/2)+(N1^2-n3^2)^(1/2))*exp(-2*k0*(N1^2-n3^2)^(1/2)*h)) ... /(1-((N1^2-n4^2)^(1/2)-(N1^2-n3^2)^(1/2))/((N1^2-n4^2)^(1/2)+(N1^2-n3^2)^(1/2))*exp(-2*k0*(N1^2-n3^2)^(1/2)*h)); f=-k0*(n1^2-N1^2)^(1/2)+n*pi+atan((N1^2-n2^2)^(1/2)/(n1^2-N1^2)^(1/2))-atan((N1^2-n2^2)^(1/2)/(n1^2-N1^2)^(1/2)*delta3); ezplot(f); solve(f,'N1') |
3楼2015-05-08 16:40:29
2楼2015-05-08 14:56:51
4楼2015-05-08 17:23:19
5楼2015-05-08 17:36:50