24小时热门版块排行榜    

查看: 1858  |  回复: 6

luke13vivid

铜虫 (初入文坛)

[求助] 想请教一个拟合得一阶德拜方程参数的问题 已有1人参与

我的测量对象为生物实验里常用的PBS溶液。现在需要知道它的弛豫频率、电导率、静电场介电常数和高频率的介电常数这4个参数。
然后用安捷伦的网络分析仪对这个溶液进行了测量,分别得到了1GHz~6GHz频率范围内的实部和虚部值。
现在的问题是不知道如何对数据进行拟合,以此来得到提到的4个参数。感觉拟合的时候,肯定需要同时满足实部和虚部的值。使用matlab的话,我使用的是cftool这个命令的工具箱,但是这里面没办法实现2组数据的同时拟合,仅对虚部拟合的话,R方为负数;使用1stopt这个工具,可以满足二者的同时拟合,但是每次拟合得到的结果中,有2个参数的值总是不稳定,偏差很大,基本不能用。
所以现在想来请教各位,你们一般如何处理呢?
附件为数据。谢谢各位!
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : PBS_ask.xlsx
  • 2015-09-17 12:36:24, 13.9 K

» 猜你喜欢

笑策无非
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
拟合公式给出来啊,文本格式
2楼2015-09-18 16:21:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luke13vivid

铜虫 (初入文坛)

引用回帖:
2楼: Originally posted by dingd at 2015-09-18 16:21:56
拟合公式给出来啊,文本格式

y1: = a1/(1+x^2/(a3^2));
y2: = a4/(8.85e-12*x*2*pi) + (a2-a1)*x*2*pi/(2*pi*a3+2*pi*x^2/a3);
a1是高频介电常数,a2是静电场介电常数,a3是频率,a4是电导率。
另,我标题写了是一阶德拜方程的表达式。
笑策无非
3楼2015-09-18 23:00:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

1stOpt结果:

均方差(RMSE):0.190588164869755
残差平方和(SSE):7.26476971768416
相关系数(R): 0.994346014509003
相关系数之平方(R^2): 0.988723996569938
确定系数(DC): 0.98865337083679
F统计(F-Statistic): 1755.12267972599

参数                  最佳估算
--------------------        -------------
a1        64.5035113000666
a3        -20286617999.6817
a4        1.55496534593737
a2        6.32823281735981
4楼2015-09-22 12:09:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luke13vivid

铜虫 (初入文坛)

引用回帖:
4楼: Originally posted by dingd at 2015-09-22 12:09:12
1stOpt结果:

均方差(RMSE):0.190588164869755
残差平方和(SSE):7.26476971768416
相关系数(R): 0.994346014509003
相关系数之平方(R^2): 0.988723996569938
确定系数(DC): 0.98865337083679
F统计(F-Stati ...

谢谢啦,不过a3是频率,拟合出负数的话。。。
笑策无非
5楼2015-09-22 18:42:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
luke13vivid: 金币+10, 有帮助 2015-09-24 19:52:59
实际上有两组解,不为负的条件应该提前说清:

均方差(RMSE):0.190588164869755
残差平方和(SSE):7.26476971768418
相关系数(R): 0.994346014498193
相关系数之平方(R^2): 0.988723996548442
确定系数(DC): 0.988653370853386
F统计(F-Statistic): 1755.12268446855

参数                  最佳估算
--------------------        -------------
a1        64.5035112881546
a3        20286618063.223
a4        1.55496534595833
a2        122.67878996248
6楼2015-09-22 21:36:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luke13vivid

铜虫 (初入文坛)

引用回帖:
6楼: Originally posted by dingd at 2015-09-22 21:36:29
实际上有两组解,不为负的条件应该提前说清:

均方差(RMSE):0.190588164869755
残差平方和(SSE):7.26476971768418
相关系数(R): 0.994346014498193
相关系数之平方(R^2): 0.988723996548442
确定系数(DC): 0 ...

哦可,谢谢啦
笑策无非
7楼2015-09-24 19:52:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 luke13vivid 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 范进中举一文的中心思想 +5 炎黄贵胄 2026-08-22 6/300 2026-08-24 06:34 by sincosx
[基金申请] 2026国自然函评费到账 +16 羊腰板 2026-08-21 17/850 2026-08-23 23:02 by tianxiaochun
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +4 process2012 2026-08-23 5/250 2026-08-23 19:58 by jurkat.1640
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 今天放榜吗? +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[基金申请] 只有每年这种时候来逛逛小木虫 +24 yaoyewhu2008 2026-08-20 26/1300 2026-08-22 17:43 by kammury
[基金申请] filecode,4个jtjc了 +13 ziyangfang 2026-08-19 16/800 2026-08-22 17:08 by WH3796
[基金申请] 人气不行了 +8 fansofjerry 2026-08-21 8/400 2026-08-22 16:30 by zyqchem
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +16 医学老男孩 2026-08-20 20/1000 2026-08-21 21:13 by Ldrop2023
[基金申请] 科研孤儿太难了 +17 我4大白菜 2026-08-20 18/900 2026-08-21 20:57 by zhangev
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 我面上完蛋了 +7 且听虎啸 2026-08-20 8/400 2026-08-21 12:31 by 酷酷墨镜
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[论文投稿] 投稿咨询 +5 wwm09 2026-08-17 7/350 2026-08-21 10:11 by 期刊论文帮手
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 估计是周四 +3 archvillain 2026-08-18 3/150 2026-08-21 01:48 by jnhyjjm
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
请填处理意见