流场高通滤波求助

Zahor

怎么做一个三维流场lon*lat*time的七天高通滤波呢，求大神

Zahor

clc;close all;clear
file='201601uvgos.nc'
ncdisp(file);

ugos=ncread(file,'ugos');
data=nan(size(ugos));
[b,a]=butter(3,2/7,'high');
for ii=1:3600;
    for jj=1:1702;
            data(ii,jj,:)=filtfilt(b,a,squeeze(ugos(ii,jj,:)));
        end
    end
参照论坛写了个butter高通滤波，但是报错：
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
错误使用 filtfilt
输入 应为 有限。

出错 filtfilt>efiltfilt (第 118 行)
    validateattributes(x,{'double'},{'finite','nonempty'},'filtfilt');

出错 filtfilt (第 93 行)
        y=efiltfilt(b,a,x);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
是因为数据中的缺测数据导致的吗

FDU.EDU

fh=7;
n=3;   
fs=1;
wn_pass=1/fh/(1/(2*fs));
%高通滤波
[B,A]=butter(n,wn_pass,'high');

Zahor

FDU.EDU 发表于 2022-10-14 12:14
fh=7;
n=3;   
fs=1;

好的谢谢你，如果数据有缺测值是不是会导致出错呢