本帖最后由 sangzhenyu 于 2020-10-30 13:57 编辑
function G=MiuCal(X)
%参数:Miu:生长速率;MiuMax:生长速率最大值;Ks:底物饱和常数;
% Ki:底物抑制常数;S:极限底物浓度;
%[t,y,te,ye,ie] = ode23s(odefun,tspan,y0,options)
%ODE求解微分过程,然而在参数模拟中涉及的是最后的积分状态,不需要用到
%i=1,2,3,4
CPi=X(1);ni=X(2);a=X(3);b=X(4);c=X(5);
P=[ 1.13 10.481 42.06 79.635 96.7567 113.793 141.438 192.6075 231.086];
S=[28.2534 25.987 24.518 25.979 24.9225 25.5405 24.3915 29.4365 28.931];
MiuN=zeros(1,9);
for i=1:length(S)
MiuN(i)=a*S(i)/(S(i)+b+S(i)^2/c);
for j=1:length(P)
MiuN(j)=MiuN(i)*(1-(P(j)/CPi))^ni;
end
end
% 需要计算Ks、Ki、MiuMax、CPi、ni+
% 规划法求最小参数:实验次数、实验变量数、实验时间点
% minG=求和(((预测值-实际值)/实际值