|
set(g3,'LineWidth',2); axis([-5,22,-1,6]);
title('3) system responce y[n]');
xlabel('n');
Результат работы данного фрагмента показан на рис. 1.6, где представлены графики решетчатых функций и
.
Рис. 1.6. Графики функций ,
и
для примеров 5_03 и 5_04.
Пример 5_04. Для исходных данных, указанных в примере 5_03, нужно построить выходной сигнал , используя DLTI представление модели системы и универсальную функцию lsim из пакета прикладных программ CST, входящего в среду MATLAB. Это же построение следует выполнить с помощью функции filter, которая является базовым инструментом пакета DSPT.
В листинге 5.4 представлен программный код на языке MATLAB, реализующий решение данного примера, с необходимыми комментариями. Естественно, что результат работы этого кода полностью совпадает с предшествующим примером (рис. 1.6).
Листинг 5.4
% Пример 5_04: Реакция на заданный входной сигнал
Clear;
% Моменты дискретного времени
n=-5:22; nu=0:8;
% Входной сигнал
u1=exp(-0.1*abs(nu)); u0=[u1 zeros(1,14)];
u=[zeros(1,5) u0];
% Импульсная характеристика
h1=[1 -1 2 3 -2 -1 -0.5 2 3 -1];
h=[zeros(1,5) h1 zeros(1,13)];
% DLTI-объект
num=h1; den=[1 zeros(1,length(num)-1)];
dz=tf(num,den,1);
% Импульсная характеристика DLTI-объекта
[yi1,ni]=impulse(dz,22);
yi=[zeros(1,5) yi1'];
% Выходной сигнал DLTI-системы
y1=lsim(dz,u0,0:22); y=[zeros(1,5) y1'];
% Формирование выходного сигнала функцией filter
yf1=filter(yi1,1,u0);
yf=[zeros(1,5) yf1];
% Графики последовательностей
Дата добавления: 2015-09-11; просмотров: 54 | Поможем написать вашу работу | Нарушение авторских прав |
<== предыдущая лекция | | | следующая лекция ==> |
Figure(1); clf reset; | | | Figure(2); clf reset; |