24小时热门版块排行榜    

查看: 1466  |  回复: 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 的主题更新
信息提示
请填处理意见