��Byczek Marta Chrupczalski Karol Grabowska Patrycja 26.03.2015r. Identyfikacja ukB
ad�w mechatroniki Projekt 1  Splot sygnaB
u 1. Wykresy a) Rys.1 Wykres otrzymany dla obliczeo podpunktu a) b) Rys.2 Wykres otrzymany dla obliczeo podpunktu b) c) Rys.2 Wykres otrzymany dla obliczeo podpunktu c) 2.Wnioski Powy|
sze wykresy pokazuj
, |
e wyniki obliczeo wykonanie przy u|
yciu funkcji conv, jak i te otrzymane na podstawie definicji splotu s
jednakowe. Dla n-pr�bek sygnaB
�w x[n] i h[n] ilo[
d pr�bek sygnaB
u y[n] wynosi 2n-1. 3. Kod Matlaba clc; clear all; close all; %podpunkt a x1 = [ ones(1,1) zeros(1,9) ]; u1 = [ ones(1,4) zeros(1,6) ]; for i=1:10 h1(i) = sin((2*pi*i)/8)*u1(i); end y1=zeros(size(h1)); for n=1:10 for k=1:10 if (n-k+1)<�=0 y1(n)=y1(n); else y1(n)=y1(n)+x1(k)*h1(n-k+1); end end end % wykresy figure(1) subplot(4,1,1) stem(x1) xlabel('Pr�bka'); ylabel('Warto[

'); title('SygnaB
wej[
ciowy x[n] '); subplot(4,1,2) stem(h1) xlabel('Pr�bka'); ylabel('Warto[

'); title('Odpowiedz
impulsowa h[n]'); subplot(4,1,3) konw1=conv(h1,x1); stem(konw1(1:10)) xlabel('Pr�bka'); ylabel('Warto[

'); title('Odpowiedz
y[n]- funkcja conv'); subplot(4,1,4) stem(y1) xlabel('Pr�bka'); ylabel('Warto[

'); title('Odpowiedz
y[n] obliczona z definicji'); %podpunkt b x2=[ ones(1,4) zeros(1,6)]; u2=[ ones(1,4) zeros(1,6)]; for i=1:10 h2(i) = sin((2*pi*i)/8)*u2(i); end y2=zeros(size(h2)); % funkcja splotu wg definicji for n=1:10 for k=1:10 if (n-k+1)<�=0 y2(n)=y2(n); else y2(n)=y2(n)+x2(k)*h2(n-k+1); end end end % wykresy figure(2) subplot(4,1,1) stem(x2) xlabel('Pr�bka'); ylabel('Warto[

'); title('SygnaB
wej[
ciowy x[n]'); subplot(4,1,2) stem(h2) xlabel('Pr�bka'); ylabel('Warto[

'); title('Odpowiedz
impulsowa h[n]'); subplot(4,1,3) konw2=conv(h2,x2); stem(konw2(1:10)) xlabel('Pr�bka'); ylabel('Warto[

'); title('Odpowiedz
y[n]- funkcja conv'); subplot(4,1,4) stem(y2) xlabel('Pr�bka'); ylabel('Warto[

'); title('Odpowiedz
y[n] obliczona z definicji'); %podpunkt c u3h=[ ones(1,5) zeros(1,5)]; u3x=[ ones(1,3) zeros(1,7)]; for i=1:10 h3(i) = sin((2*pi*i)/8)*u3h(i); x3(i) = sin((2*pi*i)/8)*u3x(i); end y3=zeros(size(h3)); % funkcja splotu wg definicji for n=1:10 for k=1:10 if (n-k+1)<�=0 y3(n)=y3(n); else y3(n)=y3(n)+x3(k)*h3(n-k+1);  end end end % wykresy figure(3) subplot(4,1,1) stem(x3) xlabel('Pr�bka'); ylabel('Warto[

'); title('SygnaB
wej[
ciowy x[n] '); subplot(4,1,2) stem(h3) xlabel('Pr�bka'); ylabel('Warto[

'); title('Odpowiedz
impulsowa h[n]'); subplot(4,1,3) konw3=conv(h3,x3); stem(konw3(1:10)) xlabel('Pr�bka'); ylabel('Warto[

'); title('Odpowiedz
y[n]- funkcja conv '); subplot(4,1,4) stem(y3) xlabel('Pr�bka'); ylabel('Warto[

'); title('Odpowiedz
y[n] obliczona z definicji');