Question

我的代码中出现此错误：

使用odearguments时出错（第91行）（。C1 *（（V1 P_OPEN）+ V2）（CA_ER-C））@（T，C） - （（V3 *（C. ^ 2））/（C。 ^ 2 +（K3 ^ 2）））必须返回列向量。

但是在MATLAB文档中有一个例子Example 3，其中向量作为输入给出，但它可以正常工作。为什么我的代码出错？
这是我的代码：

Ca_ER   = 10e-6; 
c0  = 2e-6;   
c1  = .185; 
v1  = 6; 
v2  = .11;
v3  = .09e6;
v4  = 1.2;
k3  = .1e-6;
a1 = 400e6;
a2 = 0.2e6;
a3 = 400e6;
a4 = 0.2e6;
a5 = 20e6;
b2 = .21;
d1 = 0.13e-6;
d2 = b2/a2;
d3 = 943.4e-9;
d4 = d1*d2/d3; 
d5 = 82.34e-9; 
IP= .5e-6; 
Ca=.001e-6:.01e-6:1e-6;
num=Ca.*IP.*d2;
deno= (Ca.*IP+ IP*d2+d1*d2+Ca.*d3).*(Ca+d5);
p_open=( num./deno).^3; %this is the vector input   
dc=@(t,c) (c1.*((v1.*p_open)+v2).*(Ca_ER-c))-((v3.*(c.^2))/(c.^2+(k3^2)));
[t,c]=ode45(dc,linspace(0, 100, 1000),.19e-6);
plot(t,c);

Answer 1

错误正是如此。 ode45需要返回列向量，但变量Ca实际上是行向量。因此，在匿名函数内部，所有处理都由行向量完成，返回的是行向量，因此产生错误。

罪魁祸首是这句话：

Ca=.001e-6:.01e-6:1e-6;

以上语法声明了行向量。因此，要将其转换为列向量，只需转置：

Ca=.001e-6:.01e-6:1e-6;
Ca = Ca.'; %// Transpose

ode45中的MATLAB错误：必须返回列向量

1 个答案: