I have two scripts: one preparing/analysing data and the other plotting it. The analysis part filters signals; I would like to run this twice with the use of different filters and settings separately to see the change in the resulting plots. I know I can run the second script twice using a for-loop and inserting an IF/switch statement in the analsying part, but I wondered if something like this can be done from the 'calling' script, so i dont have to temper the other.
This is a snip of the analysis part (1):
switch runFilt
case 1
fcut = 400;
[B,A] = butter(1, fcut / (fsample/2),'low');
case 2
fcut = 10;
[B, A] = cheby1(5,0.9,fcut/(fsample/2),'low');
otherwise
return
end
matrix(5,:) = filtfilt(B, A, matrix(5,:));
matrix(6,:) = filtfilt(B, A, matrix(6,:));
processing{end+1} = ['Low-pass ', num2str(fcut), ' Hz'];
Snip of the calling function part (2):
for runFilt = 1:2
step1_prepdata
step2_analyze_data
line06_m1 = nan(1,4);
line10_m1 = nan(1,4);
line14_m1 = nan(1,4);
line06_m2 = nan(1,4);
line10_m2 = nan(1,4);
line14_m2 = nan(1,4);
for ii = 1:4
switch runFilt
case 1
line06_m1(1,ii) = mean(m1(1).area(:,ii,2));
line10_m1(1,ii) = mean(m1(1).area(:,ii,3));
line14_m1(1,ii) = mean(m1(1).area(:,ii,4));
ButterM1 = [line06_m1; line10_m1; line14_m1];
case 2
line06_m1(1,ii) = mean(m1(1).area(:,ii,2));
line10_m1(1,ii) = mean(m1(1).area(:,ii,3));
line14_m1(1,ii) = mean(m1(1).area(:,ii,4));
ChebyM1 = [line06_m1; line10_m1; line14_m1];
otherwise
return
end
end
If anyone knows how to process part(1) into part(2) that would be very helpful. The reason is, I need to do several of these scenario-changes and I don't want to change a lot in such a large code. Thanks in advance.
