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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » 金属狭缝内的有效折射率的方程求解
北京石油化工学院2026年研究生招生接收调剂公告
283/3返回列表上一页123
查看: 1973  |  回复: 27
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

月只蓝

主管区长 (职业作家)

【答案】应助回帖

引用回帖:
19楼: Originally posted by 想要_回家 at 2013-07-15 14:53:34
no是什么？

不是no，是n0，在代码的第9行，你手动修改，然后run即可。n0是迭代初值，就是预先对解的大概数值有个猜测。
我大概知道，为什么n出现负值了：因为n解出来虚数，在解方程时，a+bi与-a-bi 的平方是相等的，所以若给n加个绝对值，就不会出现图像波动的情况，同时也不需要对初值那么严格，得到类似17楼的结果。
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
21楼2013-07-15 14:59:49
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

引用回帖:
20楼: Originally posted by 想要_回家 at 2013-07-15 14:57:04
怎么我改成1.0就不知道是什么东西了，这个n0要怎么个改法？

对于解非线性方程而言，初值的选取非常关键，根据你的方程，若果n0=1的话，方程两端立刻陷入0=0的死循环，如果不人为干预，是解不出方程的。n0取大于1即可，稍微验证发现，n的取值其实随意性很大的，n取2到35，同时对n取绝对值，出来的结果都是满意的
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
22楼2013-07-15 15:03:40
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

想要_回家

铜虫 (小有名气)
但是我看了一下具体数据好像有波动的地方他们的绝对值是差不多的，
赞一下 回复此楼
。。。。。。。
23楼2013-07-15 15:05:49
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)
引用回帖:
23楼: Originally posted by 想要_回家 at 2013-07-15 15:05:49
但是我看了一下具体数据好像有波动的地方他们的绝对值是差不多的，

是的，加上绝对值后，就会严格地按照单调递减的趋势来走图了。
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
24楼2013-07-15 15:07:00
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

想要_回家

铜虫 (小有名气)
谢了啊，金币你应该收到了吧
赞一下 回复此楼
。。。。。。。
25楼2013-07-15 15:08:47
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

tgvictor99

木虫 (知名作家)

迷糊领域荣誉领主
楼主，我来领金币了
程序
k0 = 6.082e6;
kesim = -42.7616+(-1.308i);
kesid = 1;
ww = 1e-9;
syms n
pn=[];
%重复调用sovle函数，会有点慢。
for i = 1:1:100
    i %看循环的圈数
    w = ww*i;
    f = -tanh(w*k0/2*sqrt(n*n-kesid))  -kesid*sqrt(n*n-kesim)/kesim/sqrt(n*n-kesid);
    f = vpa(f);
    pn = [pn,solve(f)];
end
pn = vpa(pn',6)%%显示结果,6位数字
plot(1:max(size(pn)),pn);
计算得到的结果，就是对应于1~100nm， n的值。（都是复数哦）
ppn =

     0.180014*i + 9.45122
   - 0.081308*i - 5.90488
  - 0.0538191*i - 4.64296
  - 0.0411873*i - 3.96265
  - 0.0339041*i - 3.52568
  - 0.0291296*i - 3.21643
    0.0257322*i + 2.98365
    0.0231743*i + 2.80079
      229.545*i - 0.85346
  - 0.0195433*i - 2.52953
  - 0.0181963*i - 2.42544
  - 0.0170574*i - 2.33601
  - 0.0160788*i - 2.25821
    0.0152269*i + 2.18979
    0.0144769*i + 2.12907
  - 0.0138103*i - 2.07477
    - 121.548*i - 0.45136
    0.0126739*i + 1.98156
     54.3757*i + 0.401496
    0.0117372*i + 1.90428
  - 0.0113269*i - 1.87034
  - 0.0109486*i - 1.83902
  - 0.0105986*i - 1.81002
     43.0431*i + 0.316457
   0.00997055*i + 1.75799
   - 39.7302*i - 0.291513
     38.2403*i - 0.280406
     110.691*i + 0.273958
   - 71.2479*i - 0.263855
     0.254983 - 68.8568*i
     133.289*i - 0.247581
   0.00830834*i + 1.62072
     31.2965*i + 0.227711
      30.375*i + 0.220702
- 0.00777043*i - 1.57661
     28.6857*i + 0.207825
   0.00745278*i + 1.55066
     27.1616*i - 0.196263
  - 0.0071626*i - 1.52705
   0.00702666*i + 1.51602
   0.00689628*i + 1.50545
     24.5711*i - 0.176353
     23.9987*i - 0.171934
     23.4523*i - 0.167708
     0.163662 - 22.9302*i
     0.166962 - 202.117*i
- 0.00621305*i - 1.45044
    - 21.504*i - 0.152499
    - 63.2466*i - 0.15598
     20.6325*i - 0.145754
   - 60.7657*i - 0.149796
   - 59.5968*i - 0.146881
     253.399*i + 0.144938
     95.6546*i + 0.141968
     56.3451*i + 0.138769
     110.678*i - 0.136979
     108.736*i - 0.134567
      53.4298*i + 0.13149
     52.5159*i - 0.129227
     51.6482*i + 0.127039
      16.9099*i + 0.11616
     16.6286*i - 0.113943
     16.3712*i + 0.111808
     16.1144*i + 0.109725
     15.8582*i - 0.107689
      15.617*i - 0.105719
       15.3829*i - 0.1038
     15.1558*i - 0.101932
     14.9351*i - 0.100112
    14.7276*i + 0.0983479
    14.5192*i + 0.0966174
    14.3165*i + 0.0949289
    14.1195*i + 0.0932805
    13.9277*i + 0.0916704
     13.741*i + 0.0900972
    13.5593*i + 0.0885593
    13.3822*i + 0.0870551
    13.2036*i - 0.0855722
    13.0355*i - 0.0841315
    12.8715*i - 0.0827207
    12.7117*i - 0.0813385
     12.5557*i - 0.079984
    12.4092*i + 0.0786679
    12.2605*i + 0.0773654
    12.1153*i + 0.0760873
    11.9734*i + 0.0748328
    11.8293*i - 0.0735882
    11.6939*i - 0.0723778
     11.567*i + 0.0712012
     11.4375*i + 0.070032
    11.3108*i + 0.0688821
    11.1817*i - 0.0677375
    11.0605*i - 0.0666239
    10.9418*i - 0.0655275
  - 0.00355334*i - 1.2437
   - 0.003523*i - 1.24143
  - 0.00349321*i - 1.2392
- 0.00346396*i - 1.23702
    20.8487*i + 0.0736971
     0.0728772 - 20.635*i
图已经画好，无奈网速太慢，传不了图，所以给数据了。程序运行自动能得到图。
赞一下 回复此楼
每个人的人生都是时代的产物。每个人的血统、命运和梦想组成了现在这个时代。
26楼2013-07-15 16:46:34
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

tgvictor99

木虫 (知名作家)

迷糊领域荣誉领主
晕，这网速，有人已经领了金币了，才 显示出来。晕啊。
赞一下 回复此楼
每个人的人生都是时代的产物。每个人的血统、命运和梦想组成了现在这个时代。
27楼2013-07-15 16:49:01
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

tgvictor99

木虫 (知名作家)

迷糊领域荣誉领主
引用回帖:
13楼: Originally posted by 月只蓝 at 2013-07-15 14:06:40
我用的比较传统的直接解方程的办法，算出不同w对应的n，然后作图，作图的时候，n的虚数部分去掉。
可见看见在w=0和50nm处有跳跃。程序如下：

function sn_20130715
clear all;clc
global  w k0 ed em
k0=6.0 ...

这个是matlab专业户么，用的这么溜了:tiger06::tiger06::tiger06::tiger06::tiger06::tiger06::tiger06::tiger06::tiger06::tiger06::tiger06::tiger06::tiger06::tiger06::tiger06::tiger06::tiger06::tiger06::tiger06::tiger06::tiger06::tiger06::tiger06::tiger06::tiger06::tiger06::tiger06::tiger06::tiger06::tiger06::tiger06::tiger06::tiger06::tiger06::tiger05::tiger05:
赞一下 回复此楼
每个人的人生都是时代的产物。每个人的血统、命运和梦想组成了现在这个时代。
28楼2013-07-15 16:53:19
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 想要_回家 的主题更新
283/3返回列表上一页123
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂 +3 wos666 2026-04-03 3/150 2026-04-04 05:16 by gswylq
[考研] 305求调剂 +3 77Qi 2026-04-03 3/150 2026-04-03 23:01 by qzxyhcsy
[考研] 294求调剂 +6 Grey_Ey 2026-04-03 6/300 2026-04-03 20:46 by 欣喜777
[考研] 一志愿北京工业大学，324分求调剂 +8 零八# 2026-03-28 8/400 2026-04-03 14:15 by 百灵童888
[考研] 321求调剂 +17 y-yh 2026-04-01 20/1000 2026-04-03 12:57 by y-yh
[考研] 309求调剂 +14 呆菇不是戴夫 2026-04-02 14/700 2026-04-03 09:42 by 蓝云思雨
[考研] 一志愿厦门大学材料工程专硕354找调剂！！！ +8 贝呗钡钡 2026-03-30 8/400 2026-04-03 09:41 by hypershenger
[考研] 一志愿南开大学0710生物学359求调剂 +6 兔兔兔111223314 2026-03-29 8/400 2026-04-02 22:37 by louise0220
[考研] 考研调剂 +3 李木子0120 2026-04-02 5/250 2026-04-02 21:45 by dongzh2009
[考研] 0856材料与化工调剂，339 +14 10213207 2026-03-31 14/700 2026-04-02 21:01 by 1104338198
[考研] 279求调剂 +6 学而思兮知 2026-04-01 6/300 2026-04-02 09:16 by vgtyfty
[考研] 298求调剂 +4 什么是胖头鱼 2026-03-30 6/300 2026-04-01 22:06 by 客尔美德
[考研] 材料与化工（0856）304求B区调剂 +8 邱gl 2026-03-30 16/800 2026-04-01 17:58 by 邱gl
[考研] 求调剂，一志愿北林食品与营养095500，301分，已过六级，有科研经历 +4 快乐储蓄罐 2026-03-31 4/200 2026-04-01 09:26 by JourneyLucky
[考研] 339求调剂 +5 zjjkt 2026-03-31 5/250 2026-04-01 09:18 by JourneyLucky
[考研] 本2一志愿C9-333分，材料科学与工程，求调剂 +9 升升不降 2026-03-31 9/450 2026-03-31 18:01 by 无际的草原
[考研] 考研调剂求助 +7 13287130938 2026-03-31 7/350 2026-03-31 16:39 by 690616278
[考研] 材料工程专硕求调剂 +10 hyl3153942 2026-03-29 10/500 2026-03-31 16:31 by hypershenger
[考研] 08工科求调剂286 +5 tgs_001 2026-03-28 5/250 2026-03-31 08:18 by 一只好果子?
[考研] 11408总分309，一志愿东南大学求调剂，不挑专业 +5 天赋带到THU 2026-03-29 6/300 2026-03-30 20:49 by dick_runner
信息提示
关闭
请填处理意见
关闭