uvod - Laboratorij za obdelavo signalov in daljinska vodenja

8.9 Uporaba programa MATLAB pri računanju konvolucije 217
ny = [0:10]
ker so okrogli oklepaji v zgornjem ukazu preprosto za izračunati. Zaporedje x[n],
ker je dovolj kratko, lahko zapišemo direktno z:
x=[1 1 1 1 1 1]
lahko pa uporabimo funkcijo ones(n,m), kjer sta n in m dimenziji matrike enica.
Ker rabimo vrstični vektor s 6 enicami, je n = 1 in m = 6, oziroma:
x=ones(1,6)
Zaporedje y[n], ki ga izračunamo s konvolucijo, narišemo s funkcijo stem(ny,y).
Torej nalogo rešimo z naslednjim preprostim programom:
MATLAB 8.1: Izračun konvolucije y[n] = x[n] ∗ x[n]
ny=[0:10]; % interval y[n]: n je element [0,10]
x=ones(1,6);
% funkcija x[n]=1, n=[1,6]
y=conv(x,x);
% konvolucija y[n]=x[n]*x[n]
figure; stem(ny,y);
% izris zaporedja y[n]
xlabel(’n’);ylabel(’x[n]*x[n]’)
nx=[0:5];
figure; stem(nx,x);
% izris zaporedja x[n]
axis ([0 5 0 max(y)])
xlabel(’n’);ylabel(’x[n]’)
Da sta sliki 8.27b in 8.27a v enakem merilu, smo ukazu za naris slike 8.27a dodali še
ukaz:
axis ([0 6 0 max(y)])
s katerim smo drugi sliki določili enako skalo amplitudnega razmaha kot ga ima
slika 8.27b.
♦
ZGLED 8.9.2 (izračun konvolucije dveh različnih pulzov)
Imejmo sistem s prehodnim impulznim odzivom:
x[n] =
{
n , 0 n 5
0 , sicer
in s pomočjo programa MATLAB izračunajmo odziv tega sistema na (i) zaporedje iz
predhodnega zgleda in (ii) na zaporedje:
datoteka: signal_A
x[n] =
{
1 , 0 n 25
0 , sicer