dectobin.m

此为基于二代提升小波的阈值压缩算法源码

%**************************************************************************
% 函数功能：十进制转化成二进制
% 输入参数：dec为输入的十进制，l指定转化成的二进制数的长度
% 返回值：转化成的二进制数s和负系数对应的位置编号ind
%**************************************************************************
function [s,ind]=DectoBin(dec,l)

ind=find(dec<0);%存储负系数的位置
dec = dec(:); % 确认dec存储为列向量
dec=abs(dec);
dec = double(dec);

if (nargin<2)
  l=1; % Need at least one digit even for 0.
else
l = double(l);	
if ~isscalar(l) || l<0, error('N must be a positive scalar.'); end
l = round(l); % 将l化为整数
end;

%
% Actual algorithm
%
[f,e]=log2(max(dec)); % How many digits do we need to represent the numbers?
s=rem(floor(dec*pow2(1-max(l,e):0)),2);