Drexel dragonThe Math ForumDonate to the Math Forum



Search All of the Math Forum:

Views expressed in these public forums are not endorsed by Drexel University or The Math Forum.


Math Forum » Discussions » Software » comp.soft-sys.matlab

Topic: help with special deconvolution
Replies: 4   Last Post: Jan 17, 2013 8:06 PM

Advanced Search

Back to Topic List Back to Topic List Jump to Tree View Jump to Tree View   Messages: [ Previous | Next ]
Jonathon Cheah

Posts: 5
Registered: 8/21/12
help with special deconvolution
Posted: Jan 6, 2013 1:39 AM
  Click to see the message monospaced in plain text Plain Text   Click to reply to this topic Reply

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 de-conv 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 32-bit 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';



Point your RSS reader here for a feed of the latest messages in this topic.

[Privacy Policy] [Terms of Use]

© Drexel University 1994-2014. All Rights Reserved.
The Math Forum is a research and educational enterprise of the Drexel University School of Education.