我正在尝试绘制下面的等式。这是Kraemer和Lagenbach-Belz公式,如果有人熟悉的话。参数是 0< t< 1 0< cA< 4 0< cS< 4
我想我只是没有正确地增加数组。也许我需要将Wq存储为数组?或发表声明?我对MATLAB不是很熟悉。我也试过没有if语句而没有运气。矩阵错误。
t = 0:.1:1;
cA= 0:.1:4;
cS= 0:.1:4;
p = t;
if(cA <= 1)
Wq = ((t*p)/(1-p))*((cA+cS)/2)*exp(-(2*(1-p))/(3*p))*((1-cA)^2/(cA+cS));
else
Wq = ((t*p)/(1-p))*((cA+cS)/2)*exp((-(1-p))*((cA-1)/(cA+4*cS)));
end
plot(Wq,t)