Emd 노이즈 제거 프로그램 MATLAB 코드를 무릎을 꿇고 중국어 해석이 있어 이해하기 쉽다.

함수 IMF = emd(x, n); %% 은 (는) Grilling 의 EMD 와의 충돌을 피하기 위해 함수 이름을 EMD 1 등으로 변경하는 것이 좋습니다.

%%n 은 원하는 IMF 의 숫자입니다.

C = x ('; 신호의% 사본을 라인 벡터로 입력합니다

N = 길이 (x); --

% 루프는 입력 신호를 n 개의 연속 IMF 로 분해합니다

Imf = []; 연속 IMF 및 t= 1:n 의 나머지를 포함하는% 매트릭스

연속 IMF 의% 사이클

%-

% 내부 루프는 각 IMF 를 찾습니다

H = c；; % 필터링 프로세스가 시작될 때 h 는 신호입니다

Sd =1; 필터링 프로세스를 중지하는 데 사용되는% 표준 편차

반면 SD & gt0.3%, 표준 편차가 0.3 (일반)%% 이상 필터링 중지 기준

% 로컬 최대/최소 점 찾기

D = diff (h); % 근사 파생 상품%% 각 점의 파생 상품을 찾습니다

Maxmin = []; % 를 입력하여 최적 값을 저장합니다. 지금까지는 최소값과 최대값이 다르지 않았습니다

I= 1:N-2 의 경우

D(i)==0% 0% 의 0% 에서 도수가 0 인 점이' 주재점' 인 경우, 주둔점이 항상 극한점은 아닙니다 (예: x=0, 여기서 y = x = 3).

Sign (d (I-1)) ~ = sign (d (I+1))% 가 최대값인 경우%% 상주 점의 양쪽에 있는 파생 기호가 다를 경우 (

Maxmin = [maxmin, I]; %% 는 신호의 극점 좌표를 찾습니다 (최대값 및 최소값에 관계없이)

끝

Elseif sign (d (I)) ~ = sign (d (I+1))% 0 에 걸쳐 있으므로 %% If y=|x| x=0 에 있습니다

인터럽트 방법

Maxmin = [maxmin, I+1]; % 정의 0 은 at I+ 1 (비 I)%% 입니다. 극한점을 판단하는 또 다른 방법입니다.

끝

끝

If size(maxmin, 2) < 2% 인 경우 신호가 잔차 컴포넌트의 정의를 충족하는지 여부를 판단하는 잔차%% 입니다.

파열

끝

% maxmin 을 maxes 와 mins 로 분할%% 최대값과 최소값을 분리합니다.

If maxmin( 1)>maxmin(2)% 첫 번째는 최소값이 아닌 최대값입니다

Maxes = maxmin (1:2: length (maxmin));

Mins = maxmin (2: 2: length (maxmin));

Else% 는 반대입니다

Maxes = maxmin (2: 2: length (maxmin));

Mins = maxmin (1:2: length (maxmin));

End% 는 끝점을 최대값과 최소값으로 만듭니다

최대 es = [1최대 es n];

분 = [1 분 n];

%-%최대 및 최소 포락선을 얻기 위한 스플라인 보간; 전 국제 통화 기금

Maxenv = spline (최대 es, h (최대 es),1:n); %% 스플라인 함수 보간을 사용하여 모든 최대값 점을 맞춥니다.

Minenv = spline(mins, h(mins),1:n); %% 스플라인 함수 보간을 사용하여 모든 최소 점을 맞춥니다.

M = (maxenv+minenv)/2; % 최대 및 최소 포락선의 평균%% 위/아래 포락선의 평균을 구합니다.

Prevh = h；; % 수정 전 이전 값의 사본 %% h 은 (는) 분해 전 신호입니다.

H = h-m; %-평균-h %%-랩 평균을 뺍니다.

% 표준 편차 계산

Eps = 0.0000001; 0 값을 방지하려면%

SD = sum (((prevh-h)) 입니다. 2)./(prevh. 2+EPS)); %% 계산 중지 기준

끝

Imf = [IMF; H]; % 추출된 IMF 를 매트릭스 IMF 에 저장합니다

% if size(maxmin, 2) < 2, 그럼 h 는 잔기지다

알고리즘의% 정지 기준. 만약 우리가 N 전에 종점에 도착한다면.

If size(maxmin, 2) < 2

파열

끝

C = c-h; % 신호에서 추출된 IMF 를 뺍니다

끝

돌아오다

참조 데이터

/링크? Url = dv2 ef 87 togx 8 ZC BCT1ujsz 2 Kut wrm 4 fut 5 kbm zy 5 maan 5 yv7 API bq1y8 fsag 5 jvb F2 fkl i5 jhgptxu 95