1

我已经为一些运行良好的嵌套 for 循环编写了一些代码。它使用列向量、行向量和 3x2 矩阵的组合，但所有维度都是一致的。现在，当我将这些循环作为微分方程放入并使用ode15s时，突然间即使对于 t 的单次迭代它也不起作用。它说

“试图访问 MP(2,1)；索引超出范围，因为 size(MP)=[1,6]”

当实际上 MP 应该是一个 3x2 矩阵时，当然是当我在没有 ode 求解器的情况下测试代码时。有谁知道导致改变矩阵维度的 ode 求解器有什么不同？

任何帮助将不胜感激，在此先感谢！

4

1 回答 1

2

这些ode例程假定一个向量输入y0( 1xn) 并给出一个向量输出y。在调用 之后ode15s，您可以立即转换回矩阵，使用 egy = reshape(y,[3 2]);并且您的其余代码应该可以工作。

于 2012-11-01T21:28:54.237 回答