Question

我写了一个简单的代码来解决稳态发酵问题。由于变量太多，它看起来很多。

function y = f (x)

mumax = 1.10; %h-1
Ksx = 1.32; %gl^-1
Pix=1.39; %gl^-1
Pmx=49.9; %gl^-1
Kix=304; %gl^-1
Pis=47.1; %gl^-1
Pms=95.5; %gl^-1
Kis=140; %gl^-1
qsmax=3.42; %gg^-1h^-1
Kss=2.05; %gl^-1
alp=0.39;
qpmax=3.02; %gg^-1h^-1
Ksp=2.05;  %gl^-1
Kip=140;  %gl^-1
Pip=47.1;  %gl^-1
Pmp=95.5;  %gl^-1
F=240;
S=40;
V=120;
D=0.5;
mu= (mumax*x(2)) / (Ksx+x(2));


  y=[x(1)*(-D+mu*(1-(x(3)-Pix)/(Pmx-Pix))*(Kix/(Kix+x(2))));
  D*(S-x(2))-(qsmax*(x(2)/(Kss+x(2)))*(1-((x(3)-Pis)/(Pms-Pis)))*(Kis/(Kis+x(2))))*x(1);
  -x(3)*D+x(1)*(-D+mu*(1-(x(3)-Pix)/(Pmx-Pix))*(Kix/(Kix+x(2))))*alp+qpmax*(x(2)/(Ksp+x(2)))*(1-(x(3)-Pip)/(Pmp-Pip))*x(1)*(Kip/(Kip+x(2)));]; 

endfunction


[x, fval,  info] = fsolve (@f, [2; 10; 30])

我将D定义为0.5但实际上我需要在0和1之间的区间内使用D的解，然后绘制所有x（1），x（2），x（3）和D. 我试过像

这样的东西

表示i = 0：0.1：1 d = num2str（I）

但它不起作用也许我把它放错了？最好的方法是将所有内容保存到一个矩阵中以便轻松绘制。

Answer 1

这个解决方法来自MATLAB帮助fsolve处理多个输入参数。将D定义为参数：

function y = fffff (x,D)

然后你可以创建一个匿名函数来撬掉fsolve的参数。尝试像

这样的东西

ivals=linspace(0,1);
for j=ivals;D=j;
  [x, fval,  info] = fsolve (@(x) fffff(x,D), [2; 10; 30]);
  Q=[Q x];
end

如果你想测试D的很多值（或者如果你正在进行时间演化）那么你就不会以这种方式构建Q，但是像这样使用fsolve应该可以完成工作。

求解变量变量的方程

1 个答案: