MATLAB错误：必须返回一个列向量

使用错误odearguments（91行） @（T，C）（C1 *（（V1 （CA_ER-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 
[email protected](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);

来源

2015-06-05 nashynash

的错误是这一点。 ode45需要返回列向量，但您的变量Ca其实是行向量。因此，在匿名函数中，所有处理都由一个行向量完成，返回的是行向量，从而产生错误。

罪魁祸首是这样一句话：

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

以上语法声明一个行向量。因此，为了将其转化成一个列向量，简单地采取转置：

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

来源

2015-06-05 18:02:09 rayryeng

我没有转置和现在我得到这个错误：**使用错误odearguments（线93） @（T，C）（C1 *（（V1。* P_OPEN）+ V2）。*（CA_ER-C）） - （（V3。*（C.^2））/（C.^2 +（K3^2）））返回a （T，C）（C1。*（（V1。* P_OPEN）+ V2）*（CA_ER-C））的矢量。 - （（V3。*（C.^2））/（C.^2 +（K3^2））） and the initial conditions vector must have the same number of elements。** – nashynash

@nashynash - 你的第三个参数ode45。这需要成为一个向量，但是你已经把它作为一个单一的值。我不太了解你的问题发表评论......但是MATLAB的错误非常明显......很难误解他们所说的内容。 – rayryeng

第三个参数，你的意思是**（CA_ER-C）**。我的目标是在本文[文章]（http://www.pnas.org/content/89/20/9895.full.pdf）中生成图4所示的结果。 – nashynash

MATLAB错误：必须返回一个列向量

回答

相关问题