Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- close all
- N = 20;
- x = zeros(1,N);
- x(1,1)=1;
- a=0.5;
- w=1;
- ha = zeros(1,N);
- hb = zeros(1,N);
- hc = zeros(1,N);
- hd1 = zeros(1,N);
- hd2 = zeros(1,N);
- he1 = zeros(1,N);
- he2 = zeros(1,N);
- h = zeros(1,N);
- for n=1:N %obrobka wektorow systemu
- %ha(1,n)=(1/2)^n; %ciag dla systemu stabilnego (1/2 moze byc zastapiona inna wartoscia
- %ktora bedzie spelniac warunek |z|<1)Ten ciag bierze sie z Z-Transformaty
- %hc(1,n)=j^n;%ciag dla systemu "generator" (j moze byc zastapione inna wartoscia
- %ktora spelnia warunek |z|=1, najlepiej nie 1 bo to szczegolny
- %przypadek
- ha(1,n) = j^n; %kolo jedn
- hb(1,n) = (j*0.5)^n + (-j*0.5)^n; % bieguny sprzezone
- hc(1,n) = 25^n + 4^n; %bieguny os rzeczywista
- hd1(1,n) = 5^n+(-0.5)^n;%poza kolem jedn, blizej srodka ukl wsp
- hd2(1,n) = 10^n+(-0.5)^n;%poza kolem jedn, dalej srodka ukl wsp
- he1(1,n) = 0.1^n+0.2^n;%w kole jedn, blizej sr
- he2(1,n) = 0.4^n+0.55^n;%w kole jedn, dalej sr
- h1(1,n) = exp(a*n) * cos(w*n);
- end
- figure(1);%to dajemy po to zeby nam wykresy w oddzielnych oknach rysowalo
- stem(0:1:N-1, x(1,1:1:N), 'r*');
- title('nasz sygnal');
- figure(2);
- ya = conv(ha, x);%tu robimy splot
- stem(0:1:N-1, ya(1,1:1:N),'r*');
- title('na kole jedn');
- figure(3);
- yb = conv(hb, x);
- stem(0:1:N-1, yb(1,1:1:N),'r*');
- title('bieguny sprzezone');
- figure(4);
- yc = conv(hc, x);
- stem(0:1:N-1, yc(1,1:1:N),'r*');
- title('bieguny os rzeczywista');
- figure(5);
- yd1 = conv(hd1, x);
- yd2 = conv(hd2, x);
- semilogy(0:1:N-1, yd1(1,1:1:N),'r.');
- hold on;
- semilogy(0:1:N-1, yd2(1,1:1:N),'b.');
- title('poza kolem jedn, czer-blizej, nieb-dalej srodka ukl wsp');
- figure(6);
- ye1 = conv(he1, x);
- ye2 = conv(he2, x);
- stem(0:1:N-1, ye1(1,1:1:N),'r.');
- hold on;
- stem(0:1:N-1, ye2(1,1:1:N),'b.');
- title('w kole jedn, czer-blizej, nieb-dalej srodka ukl wsp');
- hv = fft(h1);
- % ampl = abs(hvv);
- % faza = angle(hvv);
- figure(7);
- %y = conv(h, x);
- hold on;
- stem(0:1:N-1, abs(hv),'r.');
- title('amplituda');
- figure(8);
- hold on;
- stem(0:1:N-1, faza(1,1:1:N),'r.');
- title('faza');
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement