Лаба 12 / Listbox

%Объявление функции
function Listbox()
global Haxes_2 Haxes_3 Haxes_5 Haxes_6 Hbtn Hlist t c EEG dy fname
r = get(Hlist,'Value')-1;
if r==0
axes(Haxes_2)
hold off
for i=1:c
plot(t, EEG(:,i)-i*dy);
xlim([0,5]);
title(fname)
xlabel('t,c')
ylabel('U, мкВ')
hold on
end
cla(Haxes_3);
cla(Haxes_5);
cla(Haxes_6);
return
else
%Рассчет фильтра
Fs=250;
Fc=[1 4;4 8;8 13;14 30]; % Граничные частоты ритмов
Ord=fix(4*Fs/(Fc(r,1)+3)); % Выбор порядка фильтра
df1=Fc(r,1)/5; % Первая (левая) частота среза АЧХ
Fstop1 = Fc(r,1)-0.7*df1; % Конец первой полосы задержки
Fpass1 = Fc(r,1)+0.3*df1; % Начало полосы пропускания
df2=Fc(r,2)/5; % Вторая (правая) частота среза АЧХ
Fpass2 = Fc(r,2)-0.3*df2; % Конец полосы пропускания
Fstop2 = Fc(r,2)+0.7*df2; % Начало второй полосы задержки
Wstop1 = 1; % Коэффициент АЧХ фильтра
Wpass = 1; % Коэффициент АЧХ фильтра
Wstop2 = 1; % Коэффициент АЧХ фильтра
b=firls(Ord, [0 Fstop1 Fpass1 Fpass2 Fstop2 Fs/2]/(Fs/2),...
[0 0 1 1 0 0], [Wstop1 Wpass Wstop2]); % Расчет фильтра
%Построение в поле 5 амплитудно-частотную характеристику фильтра
df=0.2; % Шаг по частоте
Nf=fix((Fs/2)/df); % Число отсчетов для расчета АЧХ
[h,w]=freqz(b,1,Nf); % Расчет АЧХ
Fmax=50; % Максимальная частота АЧХ
Nmax=fix(Fmax/df); % Число отсчетов для вывода АЧХ
axes(Haxes_5) % Активизация поля 5
plot(w(1:Nmax)/(2*pi)*Fs, abs(h(1:Nmax)))% Построение АЧХ
title('АЧХ фильтра')
xlabel('f,Гц')
ylabel('|H(f)|')
%Постройте в поле 6 импульсную характеристику фильтра
axes(Haxes_6) % Активизация поля 6
stem(b,'.') % Построение импульсной характеристики
title('Импульсная характеристика')
xlabel('n')
ylabel('x(n)')
%Выполните фильтрацию ЭЭГ
EEG_F=filter(b,1,EEG);
%Выполните компенсацию вносимой фильтром задержки:
N=length(EEG_F(:,1)); % Число отсчётов ЭЭГ
D=fix(Ord/2); % Определение числа тактов задержки
EEG_F(1:N-D,:)=EEG_F(D+1:N,:); % Начальный участок
for i=1:c
EEG_F(N-D+1:N,i)=EEG_F(N,i); % Компенсация задержки
end
%Очистите поле 2, заново выведите исходные ЭЭГ и на том же графике
cla(Haxes_2);
axes(Haxes_2)
for i=1:c
%plot(t, EEG(:,i)-i*dy,t,EEG_F(:,i)-i*dy,'Color','r');
plot(t, EEG(:,i)-i*dy,'Color','b');
hold on
plot(t,EEG_F(:,i)-i*dy,'Color','r');
xlim([0,5]);
title(fname)
xlabel('t,c')
ylabel('U, мкВ')
end

Var(1,:)=var(EEG); % Дисперсия исходной ЭЭГ
Var(2,:)=var(EEG_F); % Дисперсия фильтрованной ЭЭГ
axes(Haxes_3); % Активизация поля 3
barh(Var','group'); % Построение столбиковой диаграммы
title('Оценка дисперсии ЭЭГ до и после фильтрации')
set(Haxes_3,'YDir','Reverse') % Изменение порядка по вертикали

end