24小时热门版块排行榜    

查看: 2210  |  回复: 8

pianfeng1986

铁虫 (初入文坛)

[求助] 一个关于数值拟合求常系数的问题,求思路

最近遇到一个数值拟合求常系数的问题,绞尽脑汁实在搞不定,故来此请教下论坛的前辈。
     具体是这样,实际实验中测试得到一组随频率 f 变化的导纳数值 y ,这样可以得到电导曲线(实部)和电纳曲线(虚部)如下:


     
      我根据理论计算,导纳y=g(f)这个函数中包含有 a,b 两个常系数,我的目的就是把这两系数 a,b 求出来。本来考虑用matlab进行符号运算,从y的解析式找规律,但是算出来的 y 解析式实在是太复杂,大概有四五十项,里面包含a^2, b^2,.......a^5, b^5, a^2*b^2, a^2*b^2/(a+b)......各种多次元形式都有,而且还包含两个cos(h(a,b)),sin(h(a,b)的三角函数项。 基本很难找规律。 所以想请教下各位前辈或学长,这两条曲线拟合我该怎么入手? 希望能提供一个有效的思路!谢谢!
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

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

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
pianfeng1986: 金币+5, ★有帮助 2013-04-25 21:33:47
csgt0: 金币+1, 谢谢 2013-04-26 15:57:25
应该属于复数拟合吧,如果模型公式已确定直接用数值拟合就行了。
2楼2013-04-25 17:35:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pianfeng1986

铁虫 (初入文坛)

引用回帖:
2楼: Originally posted by dingd at 2013-04-25 17:35:31
应该属于复数拟合吧,如果模型公式已确定直接用数值拟合就行了。

谢谢回复,模型公式的项数太多了,不同次元的项数有四五十个
这种情况下如果非要拟合,有什么好方法吗?
3楼2013-04-25 20:32:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★
pianfeng1986: 金币+5, ★★★很有帮助 2013-04-30 22:48:10
与模型公式的项数没什么关系吧,只不过计算耗时多些而已。

[ 发自手机版 http://muchong.com/3g ]
4楼2013-04-25 21:14:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pianfeng1986

铁虫 (初入文坛)

引用回帖:
4楼: Originally posted by dingd at 2013-04-25 21:14:53
与模型公式的项数没什么关系吧,只不过计算耗时多些而已。

厚着脸皮再问下~~ 具体该怎么做呢,我以前只会用matlab中的cftool做过拟合,用这种函数处理的话拟合式的项数是有限制的。
有没有什么计算方法或函数可以进行这种很复杂的拟合?
5楼2013-04-25 21:32:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wlxydgq

铁虫 (正式写手)

★ ★ ★ ★ ★
pianfeng1986(csgt0代发): 金币+5, lz代发 2013-05-03 11:40:10
继续关注中,我拟合的时候,也需要拟合类似第二张图的图像,楼主,咱们常交流呀。
6楼2013-04-25 23:22:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wlxydgq

铁虫 (正式写手)

★
csgt0: 金币+1, 谢谢 2013-05-02 09:12:27
你可以用1stopt试试看,这个软件上手容易,但我感觉计算速度慢些!
7楼2013-04-25 23:24:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

★ ★ ★ ★ ★
pianfeng1986(csgt0代发): 金币+5, lz代发 2013-05-03 11:39:58
引用回帖:
5楼: Originally posted by pianfeng1986 at 2013-04-25 21:32:48
厚着脸皮再问下~~ 具体该怎么做呢,我以前只会用matlab中的cftool做过拟合,用这种函数处理的话拟合式的项数是有限制的。
有没有什么计算方法或函数可以进行这种很复杂的拟合?...

楼上有同学说用1stopt,是非常好的选择。
cftool有项数限制的话,可以试着自己写个程序试试。
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
8楼2013-04-26 08:51:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pianfeng1986

铁虫 (初入文坛)

引用回帖:
5楼: Originally posted by pianfeng1986 at 2013-04-25 21:32:48
厚着脸皮再问下~~ 具体该怎么做呢,我以前只会用matlab中的cftool做过拟合,用这种函数处理的话拟合式的项数是有限制的。
有没有什么计算方法或函数可以进行这种很复杂的拟合?...

原来是同道中人~
9楼2013-04-26 10:51:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 pianfeng1986 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 5sJFDKp2kwuI 2026-11-06 10/500 2026-11-09 02:59 by l48hFd8Nwae2
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 5sJFDKp2kwuI 2026-11-06 11/550 2026-11-09 02:48 by l48hFd8Nwae2
[考研] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 5sJFDKp2kwuI 2026-11-06 9/450 2026-11-09 02:41 by l48hFd8Nwae2
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 5sJFDKp2kwuI 2026-11-06 8/400 2026-11-08 23:01 by l48hFd8Nwae2
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 5sJFDKp2kwuI 2026-11-06 9/450 2026-11-08 22:46 by l48hFd8Nwae2
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +10 ncTVBM6AOJHF 2026-11-04 22/1100 2026-11-06 17:18 by 5sJFDKp2kwuI
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 ELBZ6WwEaRj6 2026-11-05 5/250 2026-11-06 16:43 by 5sJFDKp2kwuI
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 ELBZ6WwEaRj6 2026-11-04 14/700 2026-11-06 16:32 by 5sJFDKp2kwuI
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 ELBZ6WwEaRj6 2026-11-04 10/500 2026-11-06 16:23 by 5sJFDKp2kwuI
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 16:12 by 5sJFDKp2kwuI
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 ttH3f1oA6ace 2026-11-04 11/550 2026-11-06 16:10 by 5sJFDKp2kwuI
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +7 ncTVBM6AOJHF 2026-11-04 17/850 2026-11-06 13:47 by BkA9Bi5cMhkD
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +8 ncTVBM6AOJHF 2026-11-04 18/900 2026-11-06 13:36 by BkA9Bi5cMhkD
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 ELBZ6WwEaRj6 2026-11-05 5/250 2026-11-06 09:01 by 2lRNluKrJpwL
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ELBZ6WwEaRj6 2026-11-04 10/500 2026-11-06 05:32 by qL1cU46YIdK0
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-06 04:47 by qL1cU46YIdK0
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 ELBZ6WwEaRj6 2026-11-04 8/400 2026-11-06 02:24 by qL1cU46YIdK0
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ELBZ6WwEaRj6 2026-11-04 9/450 2026-11-06 02:15 by qL1cU46YIdK0
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 7/350 2026-11-06 02:13 by qL1cU46YIdK0
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +6 ncTVBM6AOJHF 2026-11-04 10/500 2026-11-06 01:51 by qL1cU46YIdK0
信息提示
请填处理意见