1 2 3 4 5 6 7
//Regular convolution xn=[1 2 1 0]; yn=[1 2 1 0 0 0 0]; YDFT=dft(yn,-1) SDFT=YDFT.*YDFT sn=dft(SDFT,1) sn1=convol(xn,xn)