drawiirf.m

AFD - Advanced Filter Design using MATLABMiroslav D. Lutovac, Dejan V. Tosicversion 1.00 released 15

Clear[x0,y0,dx,ds,F]
Clear[x1,x2,x3,x4,x5,x6,x7,x8,x9,x10,x11,x12,x13,x14,x15,x16,x17,x18,x19,x20]
Clear[x21,x22,x23,x24,x25,x26,x27,x28,x29,x30,x31,x32,x33,x34,x35,x36]
Clear[y1,y2,y3,y4,y5,y6,y7,y8,y9,y10,y11,y12,y13,y14,y15,y16,y17,y18,y19,y20]
Clear[y21,y22,y23,y24,y25,y26,y27,y28,y29,y30,y31,y32,y33,y34,y35,y36]
Clear[DrawDF1,DrawDF2,DrawTDF1,DrawTDF2]

DrawDF2[x0_,y0_,dx_,ds_,F_:8] := Module[{
x1 = x0+dx/4,   x2 = x0+2*dx/4, x3 = x0+3*dx/4, x4 = x0+4*dx/4,
x5 = x0+5*dx/4, x6 = x0+6*dx/4, x7 = x0+7*dx/4, x8 = x0+8*dx/4, 
x9 = x0+9*dx/4, x10= x0+10*dx/4,x11= x0+11*dx/4,x12 =x0+12*dx/4, 
x13= x0+13*dx/4,x14= x0+14*dx/4,x15= x0+15*dx/4,x16 =x0+16*dx/4, 
x17= x0+17*dx/4,x18= x0+18*dx/4,x19= x0+19*dx/4,x20 =x0+20*dx/4, 
x21= x0+21*dx/4,x22= x0+22*dx/4,x23= x0+23*dx/4,x24 =x0+24*dx/4, 
y1 = y0-dx/4,   y2 = y0-2*dx/4, y3 = y0-3*dx/4, y4 = y0-4*dx/4,
y5 = y0-5*dx/4, y6 = y0-6*dx/4, y7 = y0-7*dx/4, y8 = y0-8*dx/4, 
y9 = y0-9*dx/4, y10= y0-10*dx/4,y11= y0-11*dx/4,y12 =y0-12*dx/4, 
y13 = y0-13*dx/4, y14= y0-14*dx/4,y15= y0-15*dx/4,y16 =y0-16*dx/4, 
y17 = y0-17*dx/4, y18= y0-18*dx/4,y19= y0-19*dx/4,y20 =y0-20*dx/4, 
y21 = y0-21*dx/4, y22= y0-22*dx/4,y23= y0-23*dx/4,y24 =y0-24*dx/4, 
y25 = y0-25*dx/4, y26= y0-26*dx/4,y27= y0-27*dx/4,y28 =y0-28*dx/4, },
Show[{DrawInput[x0,y0,"X",2,ds,F+2] 
  , DrawOutput[x24,y0,"Y",4,ds,F+2] 
  , DrawSum[x0, y0, 2, 0, 1, 1,"S1",2,ds,F]
  , DrawSum[x18, y0, 2, 0, 1, 1,"S2",2,ds,F]
  , DrawJ[x6,y0,x12,y0,0]
  , DrawAmp[x12,y0 ,x18,"b0","", 0,0.9*ds,F]
  , DrawDelay[x12,y6,y0,"","-1",13,ds,F]
  , DrawSum[x0, y6, 1, 2, 0, 1,"S3",4,ds,F]
  , DrawSum[x18, y6, 0, 2, 1, 1,"S4",0,ds,F]
  , DrawAmp[x6 ,y6 ,x12,"a1","", 6,0.9*ds,F]
  , DrawAmp[x12,y6 ,x18,"b1","", 0,0.9*ds,F]
  , DrawDelay[x12,y12,y6,"","-1",13,ds,F]
  , DrawAmp[x6 ,y12 ,x12,"a2","", 6,0.9*ds,F]
  , DrawAmp[x12,y12 ,x18,"b2","", 0,0.9*ds,F]
  , DrawJ[x6,y12,x3,y9,0]
  , DrawJ[x18,y12,x21,y9,0]
  , DrawJunction[x0 ,y0 ,"1" ,2,1,F]
  , DrawJunction[x12,y0 ,"2" ,2,1,F]
  , DrawJunction[x18,y0 ,"3" ,2,1,F]
  , DrawJunction[x24,y0 ,"4" ,2,1,F]
  , DrawJunction[x3 ,y3 ,"5" ,0,1,F]
  , DrawJunction[x21,y3 ,"6" ,0,1,F]
  , DrawJunction[x6 ,y6 ,"7" ,2,1,F]
  , DrawJunction[x12,y6 ,"8" ,1,1,F]
  , DrawJunction[x18,y6 ,"9" ,2,1,F]
  , DrawJunction[x3 ,y9 ,"10",0,1,F]
  , DrawJunction[x21,y9 ,"11",0,1,F]
  , DrawJunction[x12,y12,"12",6,1,F]
},AspectRatio -> Automatic,Axes -> False, PlotRange -> All]];

DrawTDF2[x0_,y0_,dx_,ds_,F_:8] := Module[{
x1 = x0+dx/4,   x2 = x0+2*dx/4, x3 = x0+3*dx/4, x4 = x0+4*dx/4,
x5 = x0+5*dx/4, x6 = x0+6*dx/4, x7 = x0+7*dx/4, x8 = x0+8*dx/4, 
x9 = x0+9*dx/4, x10= x0+10*dx/4,x11= x0+11*dx/4,x12 =x0+12*dx/4, 
x13= x0+13*dx/4,x14= x0+14*dx/4,x15= x0+15*dx/4,x16 =x0+16*dx/4, 
x17= x0+17*dx/4,x18= x0+18*dx/4,x19= x0+19*dx/4,x20 =x0+20*dx/4, 
x21= x0+21*dx/4,x22= x0+22*dx/4,x23= x0+23*dx/4,x24 =x0+24*dx/4, 
y1 = y0-dx/4,   y2 = y0-2*dx/4, y3 = y0-3*dx/4, y4 = y0-4*dx/4,
y5 = y0-5*dx/4, y6 = y0-6*dx/4, y7 = y0-7*dx/4, y8 = y0-8*dx/4, 
y9 = y0-9*dx/4, y10= y0-10*dx/4,y11= y0-11*dx/4,y12 =y0-12*dx/4, 
y13 = y0-13*dx/4, y14= y0-14*dx/4,y15= y0-15*dx/4,y16 =y0-16*dx/4, 
y17 = y0-17*dx/4, y18= y0-18*dx/4,y19= y0-19*dx/4,y20 =y0-20*dx/4, 
y21 = y0-21*dx/4, y22= y0-22*dx/4,y23= y0-23*dx/4,y24 =y0-24*dx/4, },
Show[{
    DrawInput[x0,y0,"X",2,ds,F+2] 
  , DrawOutput[x23,y0,"Y",0,ds,F+2] 
  , DrawAmp[x0 ,y0 ,x11,"b0","", 0,0.9*ds,F]
  , DrawAmp[x0 ,y12 ,x5,"b1","", 0,0.9*ds,F]
  , DrawAmp[x17,y12 ,x23,"a1","", 6,0.9*ds,F]
  , DrawAmp[x0 ,y24 ,x5,"b2","", 0,0.9*ds,F]
  , DrawAmp[x11,y24 ,x23,"a2","", 6,0.9*ds,F]
  , DrawDelay[x14,y9,y3,"","-1",15,ds,F]
  , DrawDelay[x8,y21,y15,"","-1",15,ds,F]
  , DrawJ[x24,y0,x17,y0,0]
  , DrawJ[x0,y0,x0,y12,0]
  , DrawJ[x0,y12,x0,y24,0]
  , DrawJ[x23,y0,x23,y12,0]
  , DrawJ[x23,y12,x23,y24,0]
  , DrawSum[x11, y0, 2, 0, 1, 1,"S1",2,ds,F]
  , DrawSum[x5, y12, 2, 0, 1, 1,"S2",2,ds,F]
  , DrawSum[x11,y12, 1, 2, 1, 0,"S3",6,ds,F]
  , DrawSum[x5 ,y24, 1, 2, 1, 0,"S4",6,ds,F]
  , DrawJunction[x0 ,y0 ,"1" ,2,1,F]
  , DrawJunction[x11,y0 , "2" ,2,1,F]
  , DrawJunction[x23,y0 , "3" ,2,1,F]
  , DrawJunction[x14,y3 , "4" ,0,1,F]
  , DrawJunction[x14,y9 , "5" ,0,1,F]
  , DrawJunction[x0 ,y12, ""  ,4,1,F]
  , DrawJunction[x5 ,y12, "6" ,2,1,F]
  , DrawJunction[x11,y12, "7" ,2,1,F]
  , DrawJunction[x17,y12, "8" ,2,1,F]
  , DrawJunction[x23,y12 ,""  ,0,1,F]
  , DrawJunction[x8 ,y15 ,"9" ,0,1,F]
  , DrawJunction[x8 ,y21 ,"10",0,1,F]
  , DrawJunction[x5 ,y24, "11",2,1,F]
  , DrawJunction[x11,y24, "12",2,1,F]
},AspectRatio -> Automatic,Axes -> False, PlotRange -> All]];

DrawDF1[x0_,y0_,dx_,ds_,F_:8] := Module[{
x1 = x0+dx/4,   x2 = x0+2*dx/4, x3 = x0+3*dx/4, x4 = x0+4*dx/4,
x5 = x0+5*dx/4, x6 = x0+6*dx/4, x7 = x0+7*dx/4, x8 = x0+8*dx/4, 
x9 = x0+9*dx/4, x10= x0+10*dx/4,x11= x0+11*dx/4,x12 =x0+12*dx/4, 
x13= x0+13*dx/4,x14= x0+14*dx/4,x15= x0+15*dx/4,x16 =x0+16*dx/4, 
x17= x0+17*dx/4,x18= x0+18*dx/4,x19= x0+19*dx/4,x20 =x0+20*dx/4, 
x21= x0+21*dx/4,x22= x0+22*dx/4,x23= x0+23*dx/4,x24 =x0+24*dx/4, 
y1 = y0-dx/4,   y2 = y0-2*dx/4, y3 = y0-3*dx/4, y4 = y0-4*dx/4,
y5 = y0-5*dx/4, y6 = y0-6*dx/4, y7 = y0-7*dx/4, y8 = y0-8*dx/4, 
y9 = y0-9*dx/4, y10= y0-10*dx/4,y11= y0-11*dx/4,y12 =y0-12*dx/4, 
y13 = y0-13*dx/4, y14= y0-14*dx/4,y15= y0-15*dx/4,y16 =y0-16*dx/4, 
y17 = y0-17*dx/4, y18= y0-18*dx/4,y19= y0-19*dx/4,y20 =y0-20*dx/4, 
},
Show[{
    DrawInput[x0,y0,"X",2,ds,F+2] 
  , DrawOutput[x21,y0,"Y",4,ds,F+2] 
  , DrawSum[x0, y0, 2, 0, 1, 1,"S1",2,ds,F]
  , DrawSum[x15, y0, 2, 0, 1, 1,"S2",2,ds,F]
  , DrawJ[x6,y0,x9,y0,0]
  , DrawAmp[x9,y0 ,x15,"b0","", 0,0.9*ds,F]
  , DrawDelay[x9,y6,y0,"","-1",13,ds,F]
  , DrawAmp[x3,y6 ,x9 ,"a1","", 6,0.9*ds,F]
  , DrawAmp[x9,y6 ,x15,"b1","", 0,0.9*ds,F]
  , DrawJ[x3,y3,x3,y6,1]
  , DrawJ[x15,y6,x18,y3,0]
  , DrawJunction[x0 ,y0 ,"1" ,2,1,F]
  , DrawJunction[x9 ,y0 ,"2" ,2,1,F]
  , DrawJunction[x15,y0 ,"3" ,2,1,F]
  , DrawJunction[x21,y0 ,"4" ,2,1,F]
  , DrawJunction[x3 ,y3 ,"5" ,0,1,F]
  , DrawJunction[x18,y3 ,"6" ,0,1,F]
  , DrawJunction[x9 ,y6 ,"7" ,6,1,F]
},AspectRatio -> Automatic,Axes -> False, PlotRange -> All]];

DrawTDF1[x0_,y0_,dx_,ds_,F_:8] := Module[{
x1 = x0+dx/4,   x2 = x0+2*dx/4, x3 = x0+3*dx/4, x4 = x0+4*dx/4,
x5 = x0+5*dx/4, x6 = x0+6*dx/4, x7 = x0+7*dx/4, x8 = x0+8*dx/4, 
x9 = x0+9*dx/4, x10= x0+10*dx/4,x11= x0+11*dx/4,x12 =x0+12*dx/4, 
x13= x0+13*dx/4,x14= x0+14*dx/4,x15= x0+15*dx/4,x16 =x0+16*dx/4, 
x17= x0+17*dx/4,x18= x0+18*dx/4,x19= x0+19*dx/4,x20 =x0+20*dx/4, 
y1 = y0-dx/4,   y2 = y0-2*dx/4, y3 = y0-3*dx/4, y4 = y0-4*dx/4,
y5 = y0-5*dx/4, y6 = y0-6*dx/4, y7 = y0-7*dx/4, y8 = y0-8*dx/4, 
y9 = y0-9*dx/4, y10= y0-10*dx/4,y11= y0-11*dx/4,y12 =y0-12*dx/4, 
y13 = y0-13*dx/4, y14= y0-14*dx/4,y15= y0-15*dx/4,y16 =y0-16*dx/4, 
},
Show[{
    DrawInput[x0,y0,"X",2,ds,F+2] 
  , DrawOutput[x18,y0,"Y",0,ds,F+2] 
 , DrawAmp[x0 ,y0 ,x6,"b0","", 0,0.9*ds,F]
 , DrawAmp[x0 ,y12,x6,"b1","", 0,0.9*ds,F]
  , DrawAmp[x12,y12 ,x18,"a1","", 6,0.9*ds,F]
  , DrawDelay[x9,y9,y3,"","-1",15,ds,F]
  , DrawJ[x18,y0,x12,y0,0]
  , DrawJ[x0,y0,x0,y12,0]
  , DrawJ[x18,y0,x18,y12,0]
  , DrawSum[x6, y0, 2, 0, 1, 1,"S1",2,ds,F]
  , DrawSum[x6,y12, 1, 2, 1, 0,"S2",6,ds,F]
 , DrawJunction[x0 ,y0 ,"1" ,2,1,F]
 , DrawJunction[x6,y0 , "2" ,2,1,F]
 , DrawJunction[x18,y0 , "3" ,2,1,F]
 , DrawJunction[x9,y3 , "4" ,0,1,F]
 , DrawJunction[x9,y9 , "5" ,0,1,F]
 , DrawJunction[x6,y12, "6" ,2,1,F]
 , DrawJunction[x11,y12, "7" ,2,1,F]
},AspectRatio -> Automatic,Axes -> False, PlotRange -> All]];