24小时热门版块排行榜    

查看: 2313  |  回复: 4
【悬赏金币】回答本帖问题,作者e531将赠送您 5 个金币

e531

新虫 (初入文坛)

[求助] 复函数拟合 介电函数Lorentz模型拟合已有1人参与

万能的虫友们,你们好!
我现在急需拟合一个材料的介电函数的Lorentz模型,模型如图(1)。文献中只有材料的反射率-波数图形,已知反射率和介电函数的关系如图(2)。我将其转换为反射率-频率关系,用getdata提取数据点,如下。但是用1stopt拟合的参数量级和文献中的相差过大,并且计算值没变化。
我希望虫友们的帮助:
(1)帮我指出所编函数是否有问题;
(2)帮忙拟合一下介电函数模型参数。
真的非常感谢你/你们。
注:
所编代码
Title   拟合ALAS;
ComplexStr = i;
Parameter A,W1,W2,V1,V2;
Variable x,y;
Function y=abs((sqrt(A*(W1^2-x^2+i*V1*x)/(W2^2-x^2+i*V2*x))-1)/(sqrt(A*(W1^2-x^2+i*V1*x)/(W2^2-x^2+i*V2*x))+1));
Data;
6.5940679E+13        0.0765394
6.6140106E+13        0.0779875
6.6339631E+13        0.0800653
6.6539001E+13        0.0811535
6.6738385E+13        0.0823317
6.6941348E+13        0.0840498
6.7140860E+13        0.0860377
6.7340286E+13        0.0874857
6.7539854E+13        0.0898335
6.7739393E+13        0.0920013
6.7938891E+13        0.0938992
6.8138487E+13        0.0964269
6.8337984E+13        0.0983248
6.8540990E+13        0.1003129
6.8740501E+13        0.1023007
6.8939928E+13        0.1037488
6.9139524E+13        0.1062765
6.9339318E+13        0.1100638
6.9542805E+13        0.1151108
6.9739699E+13        0.1226766
6.9940766E+13        0.1345611
7.0138721E+13        0.1488746
7.0341020E+13        0.1685864
7.0541408E+13        0.1983745
7.0741344E+13        0.2030615
7.0940629E+13        0.1813877
7.1138697E+13        0.1519765
7.1342114E+13        0.1343514
7.1539573E+13        0.1232937
7.1737924E+13        0.1179042
7.1940096E+13        0.1145840
7.2142536E+13        0.1129733
7.2341552E+13        0.1118122
7.2540667E+13        0.1112810
7.2739895E+13        0.1114695
7.2939152E+13        0.1118379
7.3141889E+13        0.1121165
7.3341004E+13        0.1115853
7.3540317E+13        0.1123136
7.3739631E+13        0.1130419
7.3942453E+13        0.1138604
7.4141653E+13        0.1138689
7.4340768E+13        0.1133377
7.4539939E+13        0.1131663
7.4739111E+13        0.1129949
7.4941834E+13        0.1131836
7.5141020E+13        0.1131021
7.5340220E+13        0.1131107
7.5539420E+13        0.1131193
7.5742129E+13        0.1132180
7.5941329E+13        0.1132265
7.6140458E+13        0.1127852
7.6339715E+13        0.1131537
7.6542367E+13        0.1128925
7.6740846E+13        0.1083126
7.6939508E+13        0.1049023
7.7139854E+13        0.1121984
7.7339549E+13        0.1153559
7.7542202E+13        0.1150947
7.7741501E+13        0.1157330
7.7940743E+13        0.1160115
7.8139915E+13        0.1158401
7.8339186E+13        0.1162985
7.8538428E+13        0.1165770
7.8737656E+13        0.1167655
7.8936941E+13        0.1173139
7.9139678E+13        0.1175925
7.9338977E+13        0.1182308
7.9538305E+13        0.1190491
7.9737604E+13        0.1196875
7.9940327E+13        0.1198761
8.0139513E+13        0.1197947
8.0338840E+13        0.1206130
8.0538097E+13        0.1209814
8.0740820E+13        0.1211701
8.0940062E+13        0.1214485
8.1139262E+13        0.1214571
8.1338420E+13        0.1211958
8.1537662E+13        0.1214742
8.1740315E+13        0.1212130
8.1939430E+13        0.1206818
8.2138516E+13        0.1199706
8.2337561E+13        0.1189895
8.2540086E+13        0.1179186
8.2739229E+13        0.1175673
8.2938429E+13        0.1175758
8.3137629E+13        0.1175844
8.3340494E+13        0.1186727
8.3539906E+13        0.1200308
8.3739318E+13        0.1213889
8.3938801E+13        0.1231969
8.4138171E+13        0.1242851
8.4341021E+13        0.1252834
8.4540420E+13        0.1265516
8.4739704E+13        0.1270999

运行结果,如下图。
文献中的参数如下。

复函数拟合   介电函数Lorentz模型拟合
介电函数和反射率.png


复函数拟合   介电函数Lorentz模型拟合-1
运行结果.jpg


复函数拟合   介电函数Lorentz模型拟合-2
文献参数.png
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
盗版1stOpt 1.5版的无法进行复数拟合计算的吧,何来正确结果。
2楼2018-06-24 22:10:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

e531

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by dingd at 2018-06-24 22:10:31
盗版1stOpt 1.5版的无法进行复数拟合计算的吧,何来正确结果。

请问你能不能帮我拟合一下啊?还有你知不知道怎么提取图片中同一横坐标对应的实部曲线和虚部曲线的数据点?真的非常感谢啊!
3楼2018-06-25 12:09:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

e531

新虫 (初入文坛)

我编写程序应该是有问题的,因为sqrt函数是对实数开根号,对于复数开根号我也不知道怎么解决。
4楼2018-06-25 12:17:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ju5200

木虫 (正式写手)

那就用MATLAB啊  从头开始学一下
5楼2018-07-28 14:37:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 e531 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 科研不是打打杀杀,科研就是人情世故 +9 杞天大圣 2024-05-08 9/450 2024-05-11 07:59 by hyjcs
[基金申请] 数理口函评了吗?多大比例呀 +5 Axvdvbfs 2024-05-09 17/850 2024-05-11 06:52 by Higraduate
[基金申请] 啥也不想干 +5 lr小机灵鬼 2024-05-10 5/250 2024-05-10 23:21 by tjmamhtk
[基金申请] 是这样的? +16 xiangsu121 2024-05-10 18/900 2024-05-10 22:10 by cuicyge
[硕博家园] 盲审 +13 暮色恋伊人 2024-05-07 15/750 2024-05-10 21:28 by 暮色恋伊人
[基金申请] 面上项目代表作必须有通讯作者文章吗? +17 zlxnnx111 2024-05-09 26/1300 2024-05-10 18:58 by 星火12
[基金申请] E05送审了么? +4 wly78528 2024-05-08 4/200 2024-05-10 16:51 by 路遥还有谁
[电化学] 硫化锂正极TEM测试如何制样? +3 发仔子呀 2024-05-07 3/150 2024-05-10 13:14 by 姜尚真。
[基金申请] AI派发是根据关键词、基金题目、研究者研究方向匹配? 180(金币+5) +3 flllgend 2024-05-04 15/750 2024-05-09 18:24 by crona酱
[硕博家园] 矿大本-中科大硕(推免)-24年科研助理求职-方向不限 +7 Kelaizhang 2024-05-08 13/650 2024-05-09 16:56 by 安德里斯
[电化学] 水系电池正极制作 +3 shangfujun 2024-05-09 3/150 2024-05-09 12:03 by jacky601
[有机交流] 求助酰胺的合成 +4 洛悠灬 2024-05-08 5/250 2024-05-09 11:29 by 小肉干
[基金申请] 博士后面上基金 +3 xlwr 2024-05-07 3/150 2024-05-09 09:38 by msjy
[访问学者] NUS Bala教授招CSC联培 +4 山高水远来日方 2024-05-07 5/250 2024-05-08 20:33 by lucas1553
[有机交流] 低沸点物质浓缩 100+3 桃桃PXS 2024-05-04 8/400 2024-05-08 19:41 by 谢大奇
[基金申请] 想问下,F口信息的函评是不是还没开始? +3 wowhow 2024-05-07 5/250 2024-05-08 12:31 by cuicyge
[论文投稿] 请问如何进行查重? 6+6 fc5878 2024-05-06 8/400 2024-05-08 11:36 by TE225
[论文投稿] Optics Express一直卡在Final Decision +4 MissPro 2024-05-05 11/550 2024-05-08 08:21 by MissPro
[高分子] 求推荐研究U V热熔胶的高手 +4 千年一虫 2024-05-04 4/200 2024-05-07 09:16 by Sunrisepay
[论文投稿] 数据共用 +4 Zzd-zzd 2024-05-06 4/200 2024-05-06 15:23 by topedit
信息提示
请填处理意见