24小时热门版块排行榜    

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

e531

新虫 (初入文坛)

我编写程序应该是有问题的,因为sqrt函数是对实数开根号,对于复数开根号我也不知道怎么解决。
4楼2018-06-25 12:17:01
已阅   回复此楼   关注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的回帖

ju5200

木虫 (正式写手)

那就用MATLAB啊  从头开始学一下
5楼2018-07-28 14:37:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 上海交大招收材料化学方向科研助理/“申请考核”博士(请勿回复帖子或站内投条) +3 灵梦and紫 2024-06-12 4/200 2024-06-15 20:58 by 1822836277
[基金申请] Nature 11日发文,中国著名学者们称造假迫不得已 +5 babu2015 2024-06-14 5/250 2024-06-15 20:40 by lc231001
[找工作] 江西双非一本和四川双一流高校如何选择? 5+7 寒山敲钟 2024-06-12 22/1100 2024-06-15 20:00 by LittleBush
[基金申请] 博后基金,博管会会提前知道消息吗? +3 yuyiang 2024-06-13 3/150 2024-06-15 19:52 by Lynn212
[硕博家园] 博士毕业高校和就业的相关问题 +4 SCITOPPP 2024-06-14 6/300 2024-06-15 18:54 by SCITOPPP
[有机交流] 车间生产,真空度很高,温度很高,但减压蒸馏速度很慢。 10+12 召唤鬼泣lL 2024-06-13 35/1750 2024-06-15 16:14 by 老钱秀
[基金申请] 2024国社科通讯评审 +9 qsd10086 2024-06-13 14/700 2024-06-15 15:51 by thesuna
[教师之家] 饶议:什么制度能保障大学普通教师不用为领导拎包,不用看领导脸色 +8 zju2000 2024-06-12 14/700 2024-06-15 13:59 by chemhua
[考博] 2025考博 +4 自强不息a?a 2024-06-15 6/300 2024-06-15 12:36 by 自强不息a?a
[基金申请] 博后基金,以往的结果点不开,怎么回事呢?最后一次机会了,两次都没中前面。 +7 kyukitu 2024-06-14 13/650 2024-06-15 06:46 by 我是王小帅
[基金申请] E12面上申请 +4 汉风之遗 2024-06-13 4/200 2024-06-14 15:28 by 天外飞去来
[硕博家园] 关于硕博连读的一些疑问? +4 Lwenter 2024-06-14 4/200 2024-06-14 14:32 by ou0551
[有机交流] ππ堆积会发生在有机溶剂中吗 5+3 zibuyu0420 2024-06-13 4/200 2024-06-14 14:17 by 小肉干
[基金申请] 国自然基金公布的时候基金号有吗 +8 潇洒怡惜 2024-06-13 11/550 2024-06-14 11:24 by JRfei
[基金申请] 化学会年会改成一年开一次了? +5 babu2015 2024-06-12 13/650 2024-06-14 10:46 by babu2015
[教师之家] 请各位老师谈谈身边非升即走的例子 +10 Lucas1121 2024-06-09 25/1250 2024-06-13 21:52 by Lucas1121
[找工作] 成都产品质量检测研究院 200+3 鲸鱼663 2024-06-11 8/400 2024-06-13 08:10 by 加纳居士
[硕博家园] 申博 +3 悦悦小小鱼 2024-06-12 3/150 2024-06-12 15:11 by chen5805
[论文投稿] water research状态咨询 5+3 Flyyawa 2024-06-10 6/300 2024-06-11 09:45 by bobvan
[教师之家] 公办双非,学生论文升学就业都不行。一本大学的论文升学就业没问题吧? +5 河西夜郎 2024-06-09 5/250 2024-06-10 17:59 by yyallen2003
信息提示
请填处理意见