Name: Punzalan, Harly S.
Rating:
Section: EC41FB2/ECE 401
Submission: 09/27/16
Date of
DIGITAL SIGNAL PROCESS PROBLEMS ON ZTRANSFORM
SEATWORK (FINALS)
Determine the z-transform of the
following sequences using the
definition (4.1). Indicate the region of
convergence for each sequence and
verify the z-transform expression using
MATLAB.
>> b2 = [1 -2 -3]; a2 = [1];
>>xb2 = filter(b2,a2,delta); n2 =
n;
>> [xa1,na1] =
sigadd(xb1,n1,xb2,n2);
>> xa2 = [0 0 3 2 1 -2 -3 0 0];
>> error = max(abs(xa1-xa2))
>> error =
0
>> n= [-2:1:2]; x= [3 2 1 -2 -3]
Ans. x(z)= [3 2 1 -2 -3]
2
0<z<
x ( z )=3 z + 2 z +12 z 3 z2
Consider the sequence x(n) = (0.9)n
cos(n/4)u(n). Let.
2.) >> b = [0 0 0.64];
>> a = [1 -0.8];
>> [delta,n] = impz(0,0,10);
>> xb1 = filter(b,a,delta);
>> [u,n] = stepseq(2,0,10);
>> xb2 = ((0.8).^n).*u;
>> error = max(abs(xb1-xb2))
>> error =
1.1102e-016
>> n=[0:3];x=[0.8.^n];
Ans. x(z)=[1 0.8 0.64 0.521] z>0.8
ANSWER:
4.1
1.) >> b1 = [0 2 3]; a1 = [1];
>> [delta,n] = impz(0,0,4);
>>xb1 = filter(b1,a1,delta);
>> xb1 = fliplr(xb1); n1 = -fliplr(n);
x ( z )=1+0.8 z1 0.64 z2 +0.512 z3 + + z
3.)>> b = [ 2 0.3]; a = [1 0.3 -0.4];
>> [delta,n] = impseq(0,0,7);
>> xb1 = filter(b,a,delta);
>> [u,n] = stepseq(0,0,7);
>> xb2 = (((0.5).^n).*u)+(((0.8).^n).*u);
>> error = max(abs(xb1-xb2))
>> error =
1.1102e-016
>> n=[0:3]; x=(0.5.^n)+(0.8.^n);
1.)
Y ( z ) = y (n) zn
n=
n=
Ans. x(z)= [-0.3 -0.39 -0.387 -0.347]
z>0.8
( n2 ) z
; n=0, 1, 2,
2 m
x ( z )=0.3 z10.39 z2 0.387 z3 0.3471 z 4 + +
z
= x ( m) z
n=
4.) >> b = 0.5*[2 -cos(0.4*pi)];
>> a = [1 -cos(0.4*pi) 0.25];
>> [delta,n1] = impseq(0,0,7);
>> xb1 = filter(b,a,delta);
>> xb1 = fliplr(xb1);
>> [u,n2] = stepseq(-7,-7,0);
>> xb2 =
((2.^n2).*cos(0.4*pi*n2)).*u;
>> error = max(abs(xb1-xb2))
>> error =
2.7756e-017
>> n=[4:0];x=(2.^n).*(cos(0.4.*pi.*n));
2.)
x ( z )=1+0.1545 z 0.2023 z 0.1011 z + 0.0193 z
2
x ( m) ( z )
n=
= X ( z)
Ans. x(z)= [1 0.1545 -0.2023
-0.1011 0.0193] z<0.2
1
X ( z )=Z [( 0.9 )n cos
1 ( 0.9 ) cos
5.) >> b = [1 -3]; a = [1 -9 27 -27];
>> [delta,n1] = impseq(0,0,7);
>> xb1 = filter(b,a,delta);
>> [u,n2] = stepseq(0,0,7);
>> xb2 = ((n2+1).*(3.^n2)).*u;
>> error = max(abs(xb1-xb2))
>> error =
0
>> n=[0:4];x=(n+1).*(3.^n);
( n4 ) u(n)]
12 ( 0.9 ) cos
( 4 )] z
( 4 )] z
+ ( 0.9 ) z2
10.6364 z
; / z / 0.9
1
2
11.2728 z +0.81 z
Y(z)=
Ans. x(z)= [1 6 27 108 405]
z>3
x ( z )=1+6 z1 +27 z2 +108 z3+ 405 z4 ++ z
4.2
10.6364 z2
; /z / 0.9=0.9487
11.2728 z 2 +0.81 z4
3.) >> b=[1 0 -0.9*cos(pi/4)];
>> a=[1 0 2*-0.9*cos(pi/4) 0
0.81];
>> [delta, n1]=impseq(0, 0, 13);
xb1=filter(b,a,
delta);
>> [u, n2]=stepseq(0,0,6);
x1=(((0.9).^n2).*cos(pi*n2/4)).*u;
>> xb2=zeros(1,2*length(x1));
xb2(1:2:end)=x1;
>> error=max(abs(xb1-xb2))
>> error =
1.2442e-0.16