ssa

1
function newpw = TrendSubtraction(pw)
2
    smoot_pw = smoothdata(pw,"movmean",32); 
3
%     smooth_pw = smooth(pw,128); 
4
    newpw=pw-smoot_pw;
5
    
6
    figure();
7
    subplot(2,1,1); 
8
    plot(pw); hold on; plot(smoot_pw); grid on; hold off; 
9
    legend('pw','smoothdata');
10
%     plot(pw); hold on; plot(smoot_pw); grid on; plot(smooth_pw); hold off; 
11
%     legend('pw','smoothdata','smooth');
12
    
13
    subplot(2,1,2);
14
    plot(newpw); grid on;
15
    
16
end