24小时热门版块排行榜    

查看: 1863  |  回复: 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

铁杆木虫 (职业作家)

【答案】应助回帖

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的回帖
查看全部 7 个回答

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的回帖

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 人气不行了 +9 fansofjerry 2026-08-21 9/450 2026-08-24 17:48 by zhaozhengxin
[基金申请] 没有任何消息-是不是就凉了 +6 图啦图啦 2026-08-24 6/300 2026-08-24 17:44 by 小小红帆
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +5 scalable 2026-08-24 5/250 2026-08-24 17:03 by spiderzhuzhi
[基金申请] 2026国自然函评费到账 +17 羊腰板 2026-08-21 19/950 2026-08-24 16:52 by iaeyuan
[基金申请] 能否退出参与的面上项目解除限项 +19 koalala 2026-08-24 21/1050 2026-08-24 16:48 by Wei_ren
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 范进中举一文的中心思想 +6 炎黄贵胄 2026-08-22 7/350 2026-08-24 11:58 by 6543yes
[基金申请] 让我中一个面上吧! +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[基金申请] 明天应该可查了!? +4 chengyan1220 2026-08-23 4/200 2026-08-24 08:47 by gloomy6159
[教师之家] 跳槽后在研项目怎么办? +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 张春生
[基金申请] 93BebMhtakh前后11位开头都是大写 +7 且听虎啸 2026-08-17 8/400 2026-08-22 21:55 by 医学老男孩
[基金申请] 只有每年这种时候来逛逛小木虫 +24 yaoyewhu2008 2026-08-20 26/1300 2026-08-22 17:43 by kammury
[基金申请] 今天基金会出结果吗?20260819 +16 kkkl_v 2026-08-19 17/850 2026-08-22 16:12 by 阿布Abu
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 重要消息,中午系统在维护 +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 家与远方
信息提示
请填处理意见