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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 求助关于非线性最小二乘法进行参数识别
181/2返回列表12下一页
查看: 1686  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 336求调剂 +5 rmc8866 2026-03-21 5/250 2026-03-21 17:24 by 学员8dgXkO
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +6 嘉年新程 2026-03-15 6/300 2026-03-21 17:07 by Dream007008
[考研] 求调剂 +3 .m.. 2026-03-21 4/200 2026-03-21 16:25 by barlinike
[考研] 求调剂 +4 十三加油 2026-03-21 4/200 2026-03-21 16:20 by barlinike
[考研] 286分人工智能专业请求调剂愿意跨考！ +4 lemonzzn 2026-03-17 7/350 2026-03-21 16:09 by lemonzzn
[考研] 070300化学319求调剂 +7 锦鲤0909 2026-03-17 7/350 2026-03-21 03:46 by JourneyLucky
[考研] 08工科 320总分 求调剂 +6 梨花珞晚风 2026-03-17 6/300 2026-03-21 03:40 by JourneyLucky
[考研] 材料工程（专）一志愿985 初试335求调剂 +3 hiloiy 2026-03-17 4/200 2026-03-21 03:04 by JourneyLucky
[考研] 274求调剂 +10 S.H1 2026-03-18 10/500 2026-03-20 23:51 by JourneyLucky
[考研] 22408 344分 求调剂 一志愿 华电计算机技术 +4 solanXXX 2026-03-20 4/200 2026-03-20 23:49 by alg094825
[考研] 南京大学化学376求调剂 +3 hisfailed 2026-03-19 6/300 2026-03-20 23:43 by hisfailed
[考研] 317求调剂 +5 申子申申 2026-03-19 9/450 2026-03-20 22:26 by JourneyLucky
[考研] 350求调剂 +5 weudhdk 2026-03-19 5/250 2026-03-20 22:04 by luoyongfeng
[考研] 材料学硕297已过四六级求调剂推荐 +11 adaie 2026-03-19 11/550 2026-03-20 21:30 by laoshidan
[考研] 295材料求调剂，一志愿武汉理工085601专硕 +5 Charlieyq 2026-03-19 5/250 2026-03-20 20:35 by JourneyLucky
[考研] 081700化工学硕调剂 +3 【1】 2026-03-16 3/150 2026-03-19 23:40 by edmund7
[考研] 328求调剂，英语六级551，有科研经历 +4 生物工程调剂 2026-03-16 12/600 2026-03-19 11:10 by 生物工程调剂
[考研] 0703化学336分求调剂 +6 zbzihdhd 2026-03-15 7/350 2026-03-18 09:53 by zhukairuo
[考研] 有没有道铁/土木的想调剂南林，给自己招师弟中～ +3 TqlXswl 2026-03-16 7/350 2026-03-17 15:23 by TqlXswl
[考研] 11408 一志愿西电，277分求调剂 +3 zhouzhen654 2026-03-16 3/150 2026-03-17 07:03 by laoshidan
信息提示
关闭
请填处理意见
关闭