| 查看: 3644 | 回复: 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 |
回复此楼» 猜你喜欢
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有4人回复
-
某top大学教授说“能够在市场中兑现的能力才是真能力”无比同意！
已经有13人回复
-
同一个分子，一条来自文献，一条来自AI——不告诉你答案，你会选哪条？
已经有7人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有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人回复
6楼2015-05-08 17:38:18
2楼2015-05-08 14:56:51
eidie
木虫之王 (文坛精英)
- 应助: 15 (小学生)
- 贵宾: 1.016
- 金币: 71294.4
- 散金: 8785
- 红花: 240
- 沙发: 1
- 帖子: 13888
- 在线: 1583.2小时
- 虫号: 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
4楼2015-05-08 17:23:19
