如何同时获取一个矩阵中包含的多个时间序列的斜率？

Question

A是180x156矩阵，B是180x156矩阵。在B内是180个向量（对应于每个纬度的平均温度），跨越156个时间值（或几个月）。

A是一系列时间值（从时间= 1到时间= 156）。我通过repmat创建了A（1：时间，180,1）。

但是情节（polyfit（A，B），1）* [1; 0]）则不输出任何实质内容。

Answer 1

您不清楚您想要线条的系数，还是想要绘制线条。如果您只想绘制我们需要创建一些虚拟数据的行：

x = 1:156;
m = repmat((0:179)/179, 156, 1)';
A = repmat(x, 180, 1);
B = (m.*A)+randn(180, 156)+179*m;

我们需要用符号绘制它

plot(A', B', '*')

然后lsline函数可以发挥其魔力

lsline

您可以看到lsline的来源，这样您就可以创建自己的函数，返回系数矩阵而不是绘制线的句柄。