clear all
w1=3*pi/4; w2=pi/4;
w3=pi/4; w4=pi/12;
n=0:10;
n=n+0.0000001;
x=(sin(w1*n)./(w1*n)).*(sin(w2*n)./(w2*n));
%y=(sin(w3*n)./(pi*n)).*(sin(w4*n)./(pi*n));
%x=(sin(w1*n)./(w1*n));
%y=(sin(w3*n)./(pi*n)).*(sin(w4*n)./(pi*n));
z0=upfirdn(x,1,3,1);
n1=n+1/3;
v1=(sin(w1*n1)./(w1*n1)).*(sin(w2*n1)./(w2*n1));
%v1=(sin(w1*n1)./(w1*n1));
z1=upfirdn(v1,1,3,1);
n2=n+2/3;
v2=(sin(w1*n2)./(w1*n2)).*(sin(w2*n2)./(w2*n2));
%v2=(sin(w1*n2)./(w1*n2));
z2=upfirdn(v2,1,3,1);
z=[z0 0 0]+[0 z1 0]+[0 0 z2];
nu=0:32;
nu=nu+0.0000001;
xu=(sin(w3*nu)./(w3*nu)).*(sin(w4*nu)./(w4*nu));
%xu=(sin(w2*nu)./(w2*nu));
%
svd([xu' z'])