我有两个时间序列数据,每个数据间隔10分钟,每个8280X1,我想简单地计算两个时间序列在连续时间间隔(每24小时)的偏移量(y-intercept)和斜率(M)。有没有这方面的matlab代码。到目前为止,我这样做了
polyfit(wl1, wl2, 1) % wl1 and wl2 is matrix column of time series
ans =
1.0184 -4.3500
但这需要完整的数据。每72行执行一次手动操作可能需要数周时间。
帮助
发布于 2018-06-05 06:35:44
使用你的解决方案,这样con就是一个矩阵。
for i=1:115
[con(i,1:2)]= [polyfit(B(:,i),BB(:,i),1)]
end
或者跳过重塑部分,执行以下操作:
for i=1:115
temp = (i-1)*72;
[con(i,1:2)] = polyfit(wl1(temp+(1:72)),wl2(temp+(1:72)),1);
end
发布于 2018-06-05 06:10:59
不管怎样,
我试着用它工作。!但如果它以矩阵格式给出答案,我会非常高兴,我想我必须努力解决它。
clc
A = (wl1)'; %wl1 is 8280x1 dim
B = reshape(A,72,115);
AA= (wl2)'; %wl2 is 8280x1 dim
BB= reshape(AA,72,115);
for i=1:115
[con]= [polyfit(B(:,i),BB(:,i),1)]
end
https://stackoverflow.com/questions/50687587
复制相似问题