function悲伤的鸭子/p clear all
clc
k0=[0. 50.50.50.50.50.50.5 ] %参数初始值
lb=[00000]； %参数下限
UB=[INFINFINFINF]； %参数上限
使用%函数fmincon (进行参数估计
使用frintf (函数fmincon ) )估计的参数值为：(n ) )
(frintf((TK1=%.4f(n )，k (1) ) ) ) ) ) ) ) ) )
frintf((TK2=%.4f(n )，k (2) ) ) ) ) ) ) ) ) )。
(frintf((TK3=%.4f(n )，k (3) ) ) ) ) ) ) ) ) )
frintf((TK4=%.4f(n )，k ) )
(frintf((TK5=%.4f(n )，k (5) ) ) ) ) ) ) ) ) )
(fr intf (thesumofthesquaresis : %.1en (n )，fval ) ) ) ) ) ) ) )。
k_fmincon=k；
------------- -请参阅
功能f=obj func4fmin con (k，x0，yexp ) )。
[ tx ]=ode 45 (关键点eqs，tspan，x0，[]，k )；
y 65:1 )=x 65:1 )；
y65:2:4(x65:4:6 )；
f=sum ((y 653360，1 )-yexp ) :1 ) )2) sum ) ) y ) :2 )-yexp ) :2 ).^2)。
sum(y ) :3 )-yexp ) :3 ) )2) sum ) ) y ) :4 )-yexp ) :4 ).^2)；
------------- -请参阅
functiondxdt=kineticeqs(t，x，k ) )。
q=8.75k(5；
dxdt=.
[k(5)-q*x(1)1)-k (1) x (1) * x (2)-k (4) * x (1) x (6) * sqrt (0.9 ) ]
(7.0-q*x(2)-k )1) x )1) x )2) k )2) *x )3) )
(1.75-q*x(3)-k )2(x ) x (3) )
(-q*x(4)2) k )1) x )1) x )2)-k )3) x )4) *x )5) )
(-q*x(5)3) k )2) x )2) x )3) k )3) x )4) x )5) )
(-q*x(6)2*k )3) x )4) x )5)-k )4) x )1) x )6) sqrt ) 0.9 ) )
(-q*x(7)2) k )4) x )1) x )6) sqrt ) 0.9 ) )
]；