Читайте также: |
|
Clear;
% Базовые моменты времени
n1=-10:10;
% Заданные отсчеты сигнала
u1=sin(0.5*n1).*exp(-0.4*abs(n1));
u=[0 0 0 u1 0 0 0]; n=[-13 -12 -11 n1 11 12 13];
% График
figure(1); clf reset; h2=stem(n,u,'bo','fill');
set(h2,'LineWidth',2); grid on;
xlabel('n'); ylabel('u[n]');
% Отсчеты единичного импульса
n1=-10; n2=10; n0=0;
nn=[n1:n2]; xx=[(nn-n0)>=0]; %!!!
% График
figure(2); clf reset;
h2=stem(nn,xx,'ro','fill');
set(h2,'LineWidth',2); grid on;
axis([-10,10,0,1.2]);
xlabel('n'); ylabel('u[n]');
Рис. 1.5. График функции для примера 5_01.
Пример 5_02. Необходимо построить конечную импульсную характеристику DLTI-системы в виде последовательности
,
заданной в моменты дискретного времени . Для применения формулы свертки нужно сгенерировать две вспомогательные последовательности: отраженную последовательность
и сдвинутую последовательность
(
).
Все указанные последовательности должны быть представлены на отдельных графиках в едином окне функции plot среды MATLAB.
Данный пример можно выполнить с помощью программного кода на языке MATLAB, представленного в листинге 5.2.
Листинг 5.2
% Пример 5_02: Импульсная характеристика и свертка
Clear;
% Моменты дискретного времени
n=-8:1:8;
% Импульсная характеристика
h=[0 0 0 0 0 0 0.1 0.2 0.4 0.5 0.6 0.8 0.7 0.5 0.2 0 0];
% Отражение от начала координат
h1=flipud(h'); h1=h1'; nn=length (h);
% Сдвиг вправо на 3 отсчета
h2=[0 0 0 h1(1:nn-3)];
% Графики последовательностей
Дата добавления: 2015-09-11; просмотров: 76 | Поможем написать вашу работу | Нарушение авторских прав |