初值未知且需要传递中间变量的微分方程参数拟合

最新推荐文章于 2021-10-02 09:27:47 发布
原创 最新推荐文章于 2021-10-02 09:27:47 发布 · 526 阅读 · 0
标签
#经验分享
本文分享了在使用1stopt进行7参数拟合过程中遇到的问题，即如何处理一个包含微分方程和代数方程的系统，其中涉及一个需要传递的中间变量lamb和未知初值lambv。通过gsl的Lu扩展数学库，实现了传递中间变量的功能，并探讨了利用时间插值提高求解精度的方法。同时，将lambv的初始值作为额外的拟合变量。最终得到了拟合结果并展示了相关图形。

问题来源：1stopt 7参数拟合求助 - 计算模拟 - 程序代码 - 小木虫论坛-学术科研互动平台 (muchong.com)

根据两个式子（一个微分方程，一个代数方程）拟合7个未知参数fac1~fac7。微分方程中lamb为需要传递的中间变量，另外，初值lambv未知。

lambv'=((1.0/3/fac7)*(fac1*((lamb/lambv')^fac2-(lambv'/lamb)^(0.5*fac2))+fac3*((lamb/lambv')^fac4-(lambv'/lamb)^(0.5*fac4))+fac5*((lamb/lambv')^fac6-(lambv'/lamb)^(0.5*fac6))));

p=(pa+fac1*((lamb/lambv)^fac2/lamb-(lambv/lamb)^(fac2*0.5)/lamb)+fac3*((lamb/lambv)^fac4/lamb-(lambv/lamb)^(fac4*0.5)/lamb)+fac5*((lamb/lambv)^fac6/lamb-(lambv/lamb)^(fac6*0.5)/lamb));

实验数据：

t

lamb

pa

p

0

0.959926667

-0.218715994

-0.060966222

0.25

0.963063333

-0.200840559

-0.046836278

0.5

0.96612

-0.183553472

-0.035042256

0.75

0.969016667

-0.167289626

-0.024438988

1

0.97171

-0.152268899

-0.01706765

1.25

0.97409

-0.139075651

-0.010939712

1
最低0.47元/天 解锁文章
matlab微分方程未知参数求解拟合,使用matlab最小二乘法拟合求解常微分方程组未知参数...
weixin_34791959的博客
03-19 2117
function k1k2k32format longclear allclctspan = [0 6 24 44 68 72 74 92 104 116]';%%这是时间yexp= [3.111,3.639,3.887,4.289,4.658,5.531,6.218,6.979,7.111,7.114]';%%%这是菌落总数x0 = [3.111];k0 = [0.1345 4 37 8...
matlab微分方程未知参数求解拟合,怎么在MATLAB 拟合常微分方程中的系数，有5个未知量，急......
weixin_34183908的博客
03-19 720
参考下：均方差(RMSE): 0.244595927796274残差平方和(SSR): 21.4181261062382相关系数(R): 0.754273774767991相关系数之平方(R^2): 0.568928927302755修正R平方(Adj. R^2): 0.572359944199415确定系数(DC): 0.536839307201338F统计(F-Statistic): 160....
缺少部分参数的微分方程拟合之一
Lu的专栏
07-31 452
问题来源：如何拟合微分方程组的参数？ – MATLAB中文论坛 (ilovematlab.cn) 微分方程组如下： dx/dt=a*x-b*x*y dy/dt=-c*y+d*x*y 数据如下： t x(t) y(t) 11 45.79 41.40 12 53.03 38.90 13 64.05 36.78 14 75.40 36.04 15 90.36 33.78 16 107.14 35.40 17 127.79 34.68 18 150.77 36.61 19 179.65 37.71 20 ..
美赛整理之带参数的常微分方程拟合问题研究
赛亚茂的博客
10-03 7099
带参数的常微分方程拟合问题研究一.问题的背景：二.提出一个较为简单，但是很有代表性的一个问题：三.求解的基本原理：四.求解的基本算法：1.利用matlab遗传算法求解出其最优误差函数值：2.差分后用LINGO求解： 一.问题的背景： 一般我们平常学习到的拟合问题大部分都是已知了函数表达式后，对于函数表达式里的系数进行的拟合问题。大部分的教材上其实并没有见讲解该如何对不能求解出相关表达式的函数关系的参数拟合求解。而且LINGO也是没有现成的微分符号表达2阶微分。因此，我通过在网上查阅部分资料提出了自己的一些小
隐函数式微分方程参数拟合
wlfc的专栏
08-09 1506
问题来源：1stopt 7参数拟合求助 - 计算模拟 - 程序代码 - 小木虫论坛-学术科研互动平台 (muchong.com) 根据两个式子（一个隐函数式微分方程，一个代数方程）拟合7个未知参数fac1~fac7。式中lambv'=d(lambv)/d(lamb)，lambv为中间变量。lambv'=((1.0/3/fac7)*(fac1*((lamb/lambv')^fac2-(lambv'/lamb)^(0.5*fac2))+fac3*((lamb/lambv')^fac4-(lambv'/lamb
微分方程参数拟合-过拟合现象及处理
Lu的专栏
07-31 2260
原问题地址：matlab 分式拟合,matlab 微分方程组的参数拟合_催眠神兔的博客-CSDN博客 微分方程式： x'=dx/dt=a*0.0321*(b-x)-d*x-dy/dt, y'= dy/dt=0.25*p1*exp(-p1*t)*x , 四个待求参数：a、b、d、p1 t、x、y数据见下面： //0 0 0//这是初值 0,0,0, 0.1,0.486966799,0.048018378, 0.167,1.6657,0.05823, 0.2,0.860306078,0...
Lu 微分方程参数优化（拟合）
Lu的专栏
07-30 3006
欢迎访问Lu程序设计 Lu 微分方程参数优化（拟合） 例子1数学模型： dy/dt=k*y*z+0.095*b*z dz/dt=-b*z-0.222*z 实验数据（ti; yi）： ti yi0.012.3291015630.683.8512927831.14.500939361.636.7491722472.079.1120628722.679.6916573663.0911.16928...
matlab 微分方程组参数拟合,Matlab 微分方程组参数拟合求助！
weixin_42303078的博客
03-16 2661
//若长时间没有结果，在任意可接受输入的窗口，按Ctrl+Alt+q退出Forcal运行。!using["XSLSF","sys","io"]; //使用命名空间XSLSF、sys和io//输出一维数组，该函数看不懂不要紧。将%14.6e改为%25.16e可提高输出精度i: OutVector(p:k,i)= k=FCDLen(p),printf{"\r\n"},i=0,(i//函数定义，用于计算...
带参数传递的微分方程拟合
wlfc的专栏
08-10 1152
问题来源：求大佬用1stopt 跑下程序 我运行总是出错 不知道哪里错了 - 仿真模拟 - 其他 - 小木虫论坛-学术科研互动平台 (muchong.com) 微分方程：y'=p1*y+p2*x 拟合参数：p1,p2 数据： t y x 0 1.25664E-07 0 0.005 0.114139113 0.156918191 0.01 0.224758442 0.31286893 0...
matlab 微分方程组参数拟合,如何拟合微分方程组的参数？
weixin_31889919的博客
03-16 2823
用5阶龙格-库塔法，结果如下，常微分方程算法: 五阶龙格-库塔法(Fith Order Runge-Kutta Method)优化算法: 通用全局优化法(UGO1)计算结束原因: 达到收敛判定标准计算用时(时:分:秒:微秒): 00:00:02:672均方差(RMSE): 3.01974154808009残差平方和(SSE): 1623.1533450618相关系数(R): 0.999906678...
matlab 微分方程组参数拟合,拟合常微分方程 (ODE)
weixin_39622217的博客
03-16 3821
洛仑兹方程组：定义和数值解洛仑兹方程组是常微分方程组(请参阅洛仑兹方程组)。对于实数常量 σ,ρ,β，方程组为dxdt=σ(y-x)dydt=x(ρ-z)-ydzdt=xy-βz.一个敏感系统的洛仑兹参数值是 σ=10,β=8/3,ρ=28。从 [x(0),y(0),z(0)] = [10,20,10] 开始启动系统，查看系统从时间 0 到 100 的演变。sigma = 10;beta = 8/...
matlab 微分方程组参数拟合,matlab拟合微分方程组中的参数
weixin_30616309的博客
03-16 2954
CODE:function k1k2k3format longclear allclctspan = 0:12:144;x0 = [0.350 99.8];k0 = [0.5 8 1.5 2 -0.5 0.1 0.5 5 -0.01]; %k的初值，最需要调节lb = [1 1 1 1 1 1 1 1 1 ]*-inf;ub = [1 1 1 1 1 1 1 1 1]*inf;d...
用 GSL 解常微分方程初值问题
Ivan 的专栏
05-12 5825
用 GSL 解常微分方程初值问题 GNU Scientific Library (GSL) 是一个用于科学计算的 C 语言类库。这里介绍如何用 GSL 来进行常微分方程初值问题的计算。 n 维的常微分方程组一般可以描述为： \[ \frac{{dy_i (t)}}{{dx}} = f_i (t,y_1 (t), \ldots y_n (t)) \] 其中 i = 1,\ldots,n
Matlab拟合解微分方程,求大神指点matlab用拟合的方式解延迟微分方程组参数
weixin_35141129的博客
03-20 410
该楼层疑似违规已被系统折叠隐藏此楼查看此楼dy(1)=-k*y(1)*y(2);dy(2)=Z(3,1)-a*y(2)-q*y(2);dy(3)=k*y(1)*y(2)-Z(3,1);dy(4)=a*y(2);dy(5)=q*y(2);t=12345678910111213141516171819202122232425262728293031323334353637383940y2=08563...
C/C++类库大全（附github连接）
热门推荐
风轻云淡的博客
12-04 1万+
C++ 资源大全中文版 浏览中发现一篇总结的很全的c/c++类库，在此做记录，以备不时之需   github地址（内有各个库的源码连接）  https://github.com/jobbole/awesome-cpp-cn   本文源自  http://www.sohu.com/a/111825907_464041，   标准库     C++标准库，包括了STL容器，算法和函数等...
GSL中的常微分方程
yxmlhc的专栏
10-02 599
常微分方程 本章描述求解常微分方程初值问题的函数。本库提供了各种低级方法，如Runge-Kutta和Burirsch-Stoer例程，以及用于自适应步长控制的高级组件。用户可以组合这些组件来实现所需的解决方案，并完全访问任何中间步骤。驱动程序对象可以用作高级包装器，以便于使用低级函数。 这些函数声明在头文件gsl_odeiv2.h中。这是版本1.15中的一个新接口，所有函数都使用前缀gsl_odeiv2。推荐使用gsl_odeiv.h中定义的以前的gsl_odeiv开头的实现。为了向后兼...
微分方程模型——偏微分方程
hanhan的博客
07-14 9524
1. 简介 微分方程：描述自然界中存在的物理现象和普遍规律。 常微分方程（ODE） 偏微分方程（PDE） 偏微分方程理论： 物理/工程问题————翻译（建模）/物理工程规律————》数学问题（PDE） 物理/工程问题————求解/数学理论————》数学结果 物理/工程问题————分析————》数学公式/物理意义 偏微分方程的基本概念： 定义：未知函数及其偏导数所满足的方程；F(x, u(x), Du, D2u,…, Dnu) = 0; 阶数：偏微分方程中偏导数的最高阶数，有n阶就为n。 线性偏微分
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值