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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 求助关于非线性最小二乘法进行参数识别
181/2返回列表12下一页
查看: 1693  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 本人女孩 +7 吼吼， 2026-04-10 9/450 2026-04-11 14:45 by ACS Nano——
[考研] 280求调剂 +6 兮兮夜夜 2026-04-09 9/450 2026-04-11 12:16 by zhq0425
[考研] 0854调剂 +11 长弓傲 2026-04-09 12/600 2026-04-11 11:16 by zhq0425
[考研] 0854调剂 +8 950824he@ 2026-04-09 8/400 2026-04-11 10:11 by zhq0425
[考研] 求调剂 +10 翩翩一书生 2026-04-09 10/500 2026-04-11 09:32 by 益达核桃露
[考研] 284求调剂 +19 梵@@ 2026-04-06 21/1050 2026-04-10 21:12 by zhouxiaoyu
[考研] 求调剂288 +6 ioodiiij 2026-04-10 8/400 2026-04-10 21:07 by zhouxiaoyu
[硕博家园] 0856材料化工求调剂，一志愿211，初试成绩349 +5 江淮北月 2026-04-05 5/250 2026-04-10 16:26 by 高维春
[考研] 考研二轮调剂 +8 故人?? 2026-04-09 8/400 2026-04-10 09:44 by 青梅duoduo
[考研] 08600生物与医药-327 +10 18755400796 2026-04-05 10/500 2026-04-10 08:14 by kangsm
[硕博家园] 有没有学校材料专业收跨调(一志愿085410) +5 momo(上岸版) 2026-04-06 8/400 2026-04-09 15:07 by only周
[考研] 化学工程与技术专业一志愿哈工程 291分B区 国家级大创负责人 有一作论文 +13 Emmy~ 2026-04-09 13/650 2026-04-09 14:47 by only周
[考研] 286求调剂 +19 Faune 2026-04-08 20/1000 2026-04-09 08:36 by 哦哦123
[考研] 085404，334分，求调剂 +5 sunjie8888 2026-04-08 8/400 2026-04-09 07:26 by sunjie8888
[考研] 372分材料与化工（085600）英二数二求调剂 +4 蓝笺片 2026-04-06 4/200 2026-04-07 12:30 by dongzh2009
[考研] 297分083200求助 +9 aekx 2026-04-05 9/450 2026-04-06 20:57 by flysky1234
[考研] 362求调剂一志愿中国石油大学 +4 我要考大 2026-04-06 6/300 2026-04-06 14:11 by 无际的草原
[考研] 331求调剂 +8 于征yz 2026-04-05 8/400 2026-04-06 00:54 by fmesaito
[考研] 一志愿北交大材料工程总分358求调剂 +6 cs0106 2026-04-05 6/300 2026-04-05 16:34 by imissbao
[考研] 295求调剂 +4 A你好研究生 2026-04-04 5/250 2026-04-04 22:46 by yu221
信息提示
关闭
请填处理意见
关闭