Fixes example by defining IIR Butterworth filter for comparison with …

…non-causal method.
sync2brain · Jan 16, 2025 · f3ef50a · f3ef50a
1 parent 4645e5a
commit f3ef50a
Showing 1 changed file with 5 additions and 3 deletions.
diff --git a/toolbox/examples/demo_phase_prediction_error_simple.m b/toolbox/examples/demo_phase_prediction_error_simple.m
@@ -57,11 +57,13 @@
 %% Phase error using standard non-causal methods
 disp('Determining phase using standard non-causal methods...');
 
-% Obtain applied filter
-oscBPFcoeffs = bd.getparam('OSC/alpha', 'bpf_fir_coeffs');
+% Prepare IIR Butterworth filter
+peakFrequency = 10;
+oscBPFfilter = designfilt('bandpassiir','FilterOrder',12,'HalfPowerFrequency1',peakFrequency-2,...
+    'HalfPowerFrequency2',peakFrequency+2,'SampleRate',osc_alpha_ip.Properties.SampleRate,'DesignMethod','butter');
 
 % Compute phase prediction error
-[phaseError, meanError, meanDev] = bossapi.boss.computePhasePredictionError(oscBPFcoeffs,...
+[phaseError, meanError, meanDev] = bossapi.boss.computePhasePredictionError(oscBPFfilter,...
                         syncedData.spf_sig_idx1(1+numSamples:end-1), syncedData.osc_idx1(2:end-numSamples));
 
 disp('Done.');