24小时热门版块排行榜    

查看: 1836  |  回复: 17

涡振怎么振?

金虫 (小有名气)

[求助] 求助关于非线性最小二乘法进行参数识别 已有1人参与

图1中是要求解的八元一次方程组,里面的有八个待解参数(就是小括号里出ti的八个参数,ti是时间),其他参数都是已知的(yi是位移时程,yi上一点是速度时程,D是高度,U是速度,减数fvi(ti)是力的时程),被减数是一个有八个待解参数的函数,函数形式如图2,本质就是想通过已知的位移时程、速度时程、力的时程来求得未知参数,从而拟合得到被减数这一函数。

如何通过已知的位移时程、速度时程、力的时程,求得这八个参数,进而拟合得到图2中的函数图像?

@月只蓝

求助关于非线性最小二乘法进行参数识别
图1.jpg


求助关于非线性最小二乘法进行参数识别-1
图2.jpg
回复此楼

» 本帖@通知

» 猜你喜欢

云总会开,月终会明。
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

涡振怎么振?

金虫 (小有名气)

其实也就是要对减数fvi(ti)这个力的时程曲线进行拟合,只是拟合的函数形式已经给定,就是图2中的公式,这个公式里除了八个待解参数,其他都是已知的。
云总会开,月终会明。
2楼2015-08-07 15:01:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

涡振怎么振?

金虫 (小有名气)

3楼2015-08-07 15:02:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

涡振怎么振?

金虫 (小有名气)

下面是数据,含有时间、位移时程、速度时程和力的时程。
U=9.1,D=0.175,密度rou为51.6.

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : 数据.txt
  • 2015-08-07 15:32:37, 1.24 M
云总会开,月终会明。
4楼2015-08-07 15:34:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖


感谢参与,应助指数 +1
ben_ladeng: 金币+1, 辛苦了 2015-08-11 08:16:14
就拟合而言,首先要判断方程形式是代数式还是常微分,然后再考虑方程是线性还是非线性。从图2来看,问题涉及的方程是非线性的常微分方程,其中最难解决的问题是y的导数y'本身又是隐式的。
隐式的常微分方程的拟合问题,我从未解决过,而且你的数据量达到三万多行,这个问题我难以解决。
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
5楼2015-08-07 15:54:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

涡振怎么振?

金虫 (小有名气)

引用回帖:
5楼: Originally posted by 月只蓝 at 2015-08-07 15:54:00
就拟合而言,首先要判断方程形式是代数式还是常微分,然后再考虑方程是线性还是非线性。从图2来看,问题涉及的方程是非线性的常微分方程,其中最难解决的问题是y的导数y'本身又是隐式的。
隐式的常微分方程的拟合问 ...

这一个方程组是令误差函数为0,然后由误差函数中各个参数求偏导得到的。

y’是速度时程,是已知的。
求助关于非线性最小二乘法进行参数识别-2
误差函数.jpg

云总会开,月终会明。
6楼2015-08-07 16:06:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

涡振怎么振?

金虫 (小有名气)

引用回帖:
4楼: Originally posted by 涡振怎么振? at 2015-08-07 15:34:29
下面是数据,含有时间、位移时程、速度时程和力的时程。
U=9.1,D=0.175,密度rou为51.6.

密度写错了,是50.61
云总会开,月终会明。
7楼2015-08-07 16:10:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

也许用高版本的1stOpt软件可以拟合隐式的常微分方程,但我手头的版本有4的参数限制,你的问题涉及8个参数。如果有高版本1stopt这个也许能高效解决。

另外,MATLAB编程应该能做,我只能提供一点思路:一般显式的常微分方程拟合,用 ode45函数和 lsqnonlin函数联用可以求解,其中ode45是求解显式常微分方程的,MATLAB有提供求解隐式常微分方程的函数ode15i,因此你的问题用ode15i 函数和 lsqnonlin函数联用可以一试,不过关键还是你的数据量太大,程序代码解决之后,计算量又是一个问题。可以把按一定间距提取数据,比如每隔10个取一个点,这样数据量缩减为1/10。
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
8楼2015-08-07 16:11:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

涡振怎么振?

金虫 (小有名气)

引用回帖:
8楼: Originally posted by 月只蓝 at 2015-08-07 16:11:01
也许用高版本的1stOpt软件可以拟合隐式的常微分方程,但我手头的版本有4的参数限制,你的问题涉及8个参数。如果有高版本1stopt这个也许能高效解决。

另外,MATLAB编程应该能做,我只能提供一点思路:一般显式的常 ...

开始尝试用origin没有搞出来,matlab还在初学阶段。。。我先试试吧

谢谢!
云总会开,月终会明。
9楼2015-08-07 16:21:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

引用回帖:
6楼: Originally posted by 涡振怎么振? at 2015-08-07 16:06:09
这一个方程组是令误差函数为0,然后由误差函数中各个参数求偏导得到的。

y’是速度时程,是已知的。

误差函数.jpg
...

y'是一个已知的数 还是在你给出的数据中 的一列数?
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
10楼2015-08-07 16:49:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 涡振怎么振? 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 国社科又开始会评了,不知道这次命运如何 +5 雨打竹帘 2026-08-30 9/450 2026-08-30 22:23 by 余韵清
[基金申请] 为什么到现在没收到通知? +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[基金申请] 面上意见出来了 +9 黄鸟于飞Chao 2026-08-29 18/900 2026-08-30 16:47 by 黄鸟于飞Chao
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 中青基了要发朋友圈吗? +4 349506619 2026-08-28 4/200 2026-08-29 22:41 by alongwaytogo
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 基金未中,这种答复是模板吗? +5 zhaosm1982 2026-08-27 6/300 2026-08-27 16:00 by lfy8008
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
信息提示
请填处理意见