应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 求助关于非线性最小二乘法进行参数识别
181/2返回列表12下一页
查看: 1837  |  回复: 17
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

涡振怎么振？

金虫 (小有名气)

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

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

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

@月只蓝

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


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

» 本帖@通知

@月只蓝

» 猜你喜欢
云总会开，月终会明。
1楼 2015-08-07 14:48:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

涡振怎么振？

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

涡振怎么振？

金虫 (小有名气)
@feixiaolin
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的回帖
相关版块跳转 我要订阅楼主 涡振怎么振？ 的主题更新
181/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 中青基了要发朋友圈吗？ +5 349506619 2026-08-28 5/250 2026-08-31 08:22 by huihoudubo
[基金申请] 投票:  有多少人是今天查系统知道结果的？ +17 爱看书的可乐 2026-08-26 19/950 2026-08-31 08:09 by xiangy672
[基金申请] 能否申诉？ +4 echo8914667 2026-08-30 5/250 2026-08-30 21:11 by echo8914667
[基金申请] 国自然面上复盘~欢迎讨论 (金币+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
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[基金申请] 我不理解！ +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 国合现在查不到了吗？ +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
信息提示
关闭
请填处理意见
关闭