НИС Франции. 1 лаба_7. Цель работы Синтез пф по заданной форме их. Описание
Скачать 68.29 Kb.
|
Цель работы: Синтез ПФ по заданной форме ИХ. Описание Задается порядок ЦФ 4 и форма его ИХ (вид №4) на конечном интервале. Определить коэффициенты фильтра ak и bk с помощью функции prony(...). Найти реальную ИХ синтезированного фильтра, используя функции filter (...) и impz(...). Графики Листинг N = 7; X = [-0.03 -0.4 -1 -0.8 -0.5 -0.3 -0.05 -0.02 0.03 0.03 -0.03 -0.02];%Импульсионная характеристика while 1 upr = menu ('Лабораторная работа','Первая - вариант 10','Выход'); if upr == 1 subplot(3,1,1); %создание графического окна 3 строки 1 столбец, и выводить впервый plot(X); %вывод графика по индексам от 1 до количества элементов xlabel('Номер отсчета');%подписка координат ylabel('Амплитуда'); axis([1 12 -1 1]);%масштабирование осей и вывод на экран title 'Исходная ИХ тип №7' [bk,ak] = prony (X,N,N);%метод Прони IMP = [1 0 0 0 0 0 0 0 0 0 0 0]; X1 = filter(bk,ak,IMP);%фильтрует сигнал, заданный в виде одномерногомассива x, используя дискретный фильтр, описываемый конечно-разностнымиуравнениями subplot(3,1,2);%вывод на второй экран plot(X1); xlabel('Номер отсчета'); ylabel('Амплитуда'); axis([1 12 -1 1]); title 'Синтезированная ИХ функцией filter(...)' X2 = impz(bk,ak,length(X));%Вычисляет импульсную характеристику фильтра,заданного векторами коэффициентов числителя (b) и знаменателя (a) функциипередачи. Length – размерность ИХ subplot(3,1,3); plot(X2); xlabel('Номер отсчета'); ylabel('Амплитуда'); axis([1 12 -1 1]); title 'Синтезированная ИХ функцией impz(...)' else disp ('Выход из программы'); break; end end Вывод: Синтезированные функции, созданные разными способами, совпадают с исходной импульсной характеристикой, значит методы, что используются в функциях filter и impz для синтезирования результативно правильные. |