应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » 新人求助解MATLAB超越方程，方程中含有复数
51/1返回列表
查看: 3078  |  回复: 9
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

小雨滴001

新虫 (初入文坛)

[求助] 新人求助解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]

新人求助解MATLAB超越方程，方程中含有复数
N16Y7`CLCLKS{PN)53U@6JX.png
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2015-05-08 14:54:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xzczd

木虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★
小雨滴001: 金币+5, ★★★很有帮助, 谢谢您的指教，对我非常有帮助 2015-05-09 09:47:26
引用回帖:
6楼: Originally posted by 小雨滴001 at 2015-05-08 17:38:18
嗯嗯 我看到了 其实是我下面给出的公式里打错了 程序里是对的 谢谢您哦...

啥，程序里才是对的？那你确定你这方程没错？：
CODE:
With[{N1 = x + I y},
Plot3D[Abs[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], {x, -20, 20}, {y, -20, 20}(*,
  PlotRange\[Rule]{0,1}*), PlotPoints -> 100]]

新人求助解MATLAB超越方程，方程中含有复数-1
赞一下 回复此楼
小木虫Mathematica版块已毁（当然原本也不咋的），建议大家前往百度贴吧或Stackexchange。
8楼2015-05-08 17:59:03
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答

小雨滴001

新虫 (初入文坛)
忘记说清楚，n是可以取从0开始的任意整数的
赞一下 回复此楼
2楼2015-05-08 14:56:51
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

eidie

木虫之王 (文坛精英)

文献杰出贡献

【答案】应助回帖

感谢参与，应助指数 +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
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xzczd

木虫 (小有名气)

【答案】应助回帖

感谢参与，应助指数 +1
别的先不说，你的代码里，第二个ArcTan中某处的n2写成了n3你发现了没有？
赞一下 回复此楼
小木虫Mathematica版块已毁（当然原本也不咋的），建议大家前往百度贴吧或Stackexchange。
4楼2015-05-08 17:23:19
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭