过滤数据误差
4视图(30天)
显示旧的评论
我想过滤sla数据2191数据点(每日数据从2013年到2018年)。我用下面的代码。当我运行这段代码,我得到一个错误。
在filtfilt > getCoeffsAndInitialConditions(第200行)
在filtfilt(第97行)
在filter_matlab(第8行)
警告:矩阵接近奇异或严重了。结果可能是不准确的。RCOND = 1.419047 e-16。
它的意义是什么?
我的代码是:
sla = ncread (sla_ave。数控”、“SLA_AVE”);
highcut = (2191/100) / (2191/2);%定义highcut频率100天弧度
lowcut = (2191/250) / (2191/2);%定义lowcut频率250天弧度
[b] =黄油(6日[lowcut highcut]);%创建6阶巴特沃斯带通滤波器
freqz (b);%看看你的过滤特性
FilteredData = filtfilt (b, a, sla);%不相移过滤数据
t = datenum (datestr (ncread (sla_ave。数控”、“时间”)+ datenum(1950年,1,1)));
情节(t, sla,“线宽”,1.5);
抓住;
情节(t, FilteredData“线宽”,1.5);
轴紧;
xtck = t (1:100:2191);
集(gca, XTick, xtck);
集(gca, XTickLabel, datestr (t (1:100:2191),“m-yy”),“FontWeight”,“大胆”,“字形大小”,8);
此外,的意义是什么定义100天的弧度。为什么它是必要的吗?
