调用ode45求解微分方程时，经常报错左右侧元素数目对不上，反复数输入输出向量的长度，看起来完全一致，却始终无法运行。这类问题不是简单的计数错误，几个隐性细节没处理好就会触发维度冲突。

1. 检查微分方程返回格式

ode45要求方程函数必须返回列向量，不能是行向量。

如果函数内部计算后输出的是行向量，即使元素总数正确，也会触发维度不匹配报错，用y(:)强制转为列向量即可修复。

2. 初值向量维度对齐

求解时传入的初始条件向量，必须和微分方程的状态变量数目完全一致。

比如3阶微分方程需要传入3个初值，少传或多传一个元素，都会直接触发维度不匹配，用length函数核对初值长度即可排查。

3. 避免隐式矩阵运算

方程内部不小心用了矩阵乘法，而非逐点运算。

比如用代替.，会把一维向量运算扩展成二维矩阵，输出元素数目直接翻倍，替换成逐点运算符号就能解决。

4. 分步调试定位

在微分方程函数开头插入size(y)语句，打印每一步输入的维度。

运行时就能直观看到哪一步维度发生了异常变化，不用盲目数元素数目，快速定位冲突点。

本文凡亿教育原创文章，转载请注明来源！