24小时热门版块排行榜    

查看: 2329  |  回复: 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的回帖

ju5200

木虫 (正式写手)

那就用MATLAB啊  从头开始学一下
5楼2018-07-28 14:37:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答

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的回帖
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 硕博巨婴,也许才刚刚开始 +30 SNaiL1995 2024-06-12 65/3250 2024-06-16 02:10 by 青古
[文学芳草园] 累并快乐着 +13 MYHLD521 2024-06-14 13/650 2024-06-15 22:59 by zeolitess
[考博] 上海交大招收材料化学方向科研助理/“申请考核”博士(请勿回复帖子或站内投条) +3 灵梦and紫 2024-06-12 4/200 2024-06-15 20:58 by 1822836277
[基金申请] BO4的YQ答辩通知发布了吗? +6 博学笃行 2024-06-11 6/300 2024-06-15 16:04 by 悲催科研狗
[基金申请] 2024国社科通讯评审 +9 qsd10086 2024-06-13 14/700 2024-06-15 15:51 by thesuna
[找工作] 杭电、天津科技、青农和宁波工程学院如何选? +9 味道很好啊 2024-06-13 15/750 2024-06-15 09:59 by yulan19992002
[论文投稿] 投了一篇4区的SCI,审稿人一个拒稿,一个小修,编辑给了大修。 +9 安稳22123 2024-06-13 10/500 2024-06-14 23:45 by jurkat.1640
[食品] 食品博士导师 +6 小李醒yy 2024-06-11 9/450 2024-06-14 23:37 by 小李醒yy
[基金申请] 面上基金有一个(两个)C是不是就没戏了.... 5+3 zzzm116 2024-06-13 17/850 2024-06-14 22:26 by lzt8076
[考研] 物理化学一对一辅导 +3 林大diao 2024-06-12 5/250 2024-06-14 20:57 by 林大diao
[基金申请] 博士后基金需要结题吗? +8 zhouchuck 2024-06-13 8/400 2024-06-14 17:27 by liuyupu132
[基金申请] 工材E口JQ有消息了吗 +4 babyduck 2024-06-11 4/200 2024-06-14 17:23 by firepick
[基金申请] 国自然基金公布的时候基金号有吗 +8 潇洒怡惜 2024-06-13 11/550 2024-06-14 11:24 by JRfei
[论文投稿] 投稿kbs被拒后,系统提供的推荐的其他期刊有用吗? +3 chenrui2015 2024-06-10 3/150 2024-06-14 10:30 by Tr.sjx1997
[论文投稿] 最近写了一篇控制优化领域的文章,可以投哪里啊?有没有水一些的期刊推荐 +7 香瓜木香 2024-06-12 13/650 2024-06-14 07:05 by 香瓜木香
[论文投稿] 文章proof要求使用机构的邮箱 5+3 不可不信缘 2024-06-11 11/550 2024-06-14 07:00 by 3001160025
[基金申请] 连续两年医学口青年项目初审体会 +11 进击的荣耀 2024-06-09 18/900 2024-06-13 17:27 by 进击的荣耀
[考博] 博导选择 +3 bing85977 2024-06-12 3/150 2024-06-13 15:34 by 我是邱尧
[有机交流] 原料反应完了,怎么知道是产物还是中间体 +6 小胡在努力 2024-06-11 8/400 2024-06-13 13:33 by 091602
[教师之家] 公办双非,学生论文升学就业都不行。一本大学的论文升学就业没问题吧? +5 河西夜郎 2024-06-09 5/250 2024-06-10 17:59 by yyallen2003
信息提示
请填处理意见