Search All of the Math Forum:
Views expressed in these public forums are not endorsed by
Drexel University or The Math Forum.



help with special deconvolution
Posted:
Jan 6, 2013 1:39 AM


Below is a code for special convolution of (32 bits) x, (50 bits) h giving (82 + 6 padded zeros )Y. I cannot find a way to deconv Y and x to get back h exactly. The code below gives a slightly different value than Matlab conv(x,h) so Matlab deconv does not return h. Tried ifft(fft(Y)*1./fft(x)) and other common methods with no luck. Any suggestions?  clc clear all % filter: x=dec2bin(4073739089.0); % signal: h=dec2bin( 60036946411689664.0); %convolution: m=length(x); n=length(h); X=[x,zeros(1,m)]; H=[h,zeros(1,(n+m))]; T=zeros(1,m); Y=zeros(1,(n+m)); % shift H into 32bit T, a bit at a time. for i=1:m+n T=[T(2:end) 0]; if str2num(H(i))>0 T(m)=1; else T(m)=0; end for j=1:m Y(i)=Y(i)+sum(str2num(X(j))& T(j)); end Y(i)=mod(Y(i),2); end dec2hex(bi2de(fliplr(Y))) % Answer Y='DBFB31F057E34000000000';



