Learn how easy it is to sync an existing GitHub or Google Code repo to a SourceForge project! See Demo

Close

Diff of /wavelets/wfilt_remez.m [a9068d] .. [4530c5] Maximize Restore

  Switch to unified view

a/wavelets/wfilt_remez.m b/wavelets/wfilt_remez.m
...
...
85
   AK=AK/sum(AK);
85
   AK=AK/sum(AK);
86
end
86
end
87
%%%%%%%%%%%%%%%%%%%%%%%%%%% STEP 2' %%%%%%%%%%%%%%%%%%%%%%%%%%%
87
%%%%%%%%%%%%%%%%%%%%%%%%%%% STEP 2' %%%%%%%%%%%%%%%%%%%%%%%%%%%
88
% Daubechies factor
88
% Daubechies factor
89
% PK(z)=((1+z^(-1))/2)^2*K QK(z)
89
% PK(z)=((1+z^(-1))/2)^2*K QK(z)
90
if computeroots & K>0
90
if computeroots && K>0
91
   QK=binom(2*K,1:K);
91
   QK=binom(2*K,1:K);
92
   QK=QK.*abs(QK);
92
   QK=QK.*abs(QK);
93
   QK=cumsum(QK);
93
   QK=cumsum(QK);
94
   QK=QK./abs(binom(2*K-1,1:K));
94
   QK=QK./abs(binom(2*K-1,1:K));
95
   QK=[QK QK(K-1:-1:1)];
95
   QK=[QK QK(K-1:-1:1)];