24小时热门版块排行榜    

查看: 2319  |  回复: 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个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 来用亲身体会,一起说说年轻老师的辛苦 +20 zylfront 2024-05-17 26/1300 2024-05-23 07:53 by luwangba
[论文投稿] word转成pdf之后公式里面的字体变了,正文字体没变。 +9 1255037206 2024-05-20 11/550 2024-05-23 05:54 by tjushede
[考博] 大扩招 +12 version211 2024-05-21 18/900 2024-05-22 23:02 by 村安全顾问
[论文投稿] 求能源类,故障诊断,机器学习结合的EI期刊。 36+3 wj0318 2024-05-19 5/250 2024-05-22 22:05 by 知识产权服务
[论文投稿] 期刊投稿进度 +6 jianhuang9 2024-05-21 8/400 2024-05-22 18:29 by sakuraai
[考博] 换导师 +16 是柠檬呀! 2024-05-18 29/1450 2024-05-22 16:29 by oooooo?o
[基金申请] 国自然的面上项目,5个审稿人,5个B能上会吗? 4+10 lancet0903 2024-05-20 32/1600 2024-05-22 15:45 by 潇湘之迷
[论文投稿] 投稿journal of energy chemistry +3 小可爱嘿呀 2024-05-21 3/150 2024-05-22 11:11 by xs74101122
[考博] 化工三篇sci只想有书读,急寻博导! +3 8139173 2024-05-21 5/250 2024-05-22 10:56 by 安塔瓦拉多
[基金申请] 国自然等 80+4 胖虎 2024-05-21 12/600 2024-05-22 09:47 by nono2009
[硕博家园] 耐高温垫片求购 +8 Sexyflea 2024-05-16 11/550 2024-05-22 08:12 by tan151646
[基金申请] 肿瘤口 +5 hml2091 2024-05-20 5/250 2024-05-21 22:48 by Liu131415
[论文投稿] 论文一审意见回来后发现实验程序编错了论证分析部分可能要大改 5+4 hshhenb 2024-05-20 5/250 2024-05-21 11:03 by bnullh
[论文投稿] SENSORS AND ACTUATORS B-CHEMICAL 15+4 荣小撇 2024-05-20 4/200 2024-05-21 10:07 by ca0yan9
[硕博家园] 又见奇葩学生,可咋办 +25 appleapple2 2024-05-17 30/1500 2024-05-21 09:13 by qyws
[考博] 双非博士还是985科研助理,然后再读博 +6 lxdatj123 2024-05-18 13/650 2024-05-21 08:19 by lxdatj123
[基金申请] 去年申请基金的评审意见是ChatGPT在国内是禁止的,研究方案中有使用ChatGPT不合理 +4 瞬息宇宙 2024-05-19 4/200 2024-05-21 00:37 by dxcharlary
[基金申请] 基金评审 +4 阿呆不呆 2024-05-20 4/200 2024-05-20 15:10 by 一路向东
[论文投稿] 推荐转投( transfer pending)是否有用? 50+3 lily5289 2024-05-17 7/350 2024-05-19 15:11 by wanghuawei
[硕博家园] 五氯化铌怎么溶解啊 +3 南南枝枝 2024-05-17 5/250 2024-05-17 11:37 by ad_fish
信息提示
请填处理意见