freq.m

the file contain many matlab signal source code and many example .

w=0:0.001:pi;
n=length(w);
for n1=1:n;
    if  w(n1)<=0.3*pi
        Hd(n1)=1;
    else
         Hd(n1)=0;
     end
 end
 k=10;
 [A,fai]=IIR_dirFreq(Hd,k)
 for n1=1:k
     m1=4*(n1-1);
     b(n1,1)=1;b(n1,2)=fai(m1+1);b(n1,3)=fai(m1+2);
     a(n1,1)=1;a(n1,2)=fai(m1+3);a(n1,3)=fai(m1+4);
 x=zeros(1,2*k);
 y=x;
end
for n1=1:k
    x=conv(x,b(n1));
    y=conv(y,a(n1));
end