shiftright.m

来自「intlab 工具用于快速计算 各各层的倒数等等」· M 代码 · 共 26 行

function m = shiftright(m,r)
%SHIFTRIGHT   Shift array m right by r bits
%
%   m = shiftright(m,r)
%
%last r bits of m(:,end) are zero, otherwise lost
%

% written  12/30/98     S.M. Rump
%

global INTLAB_LONG_BETA

  % array m divided by 2^r
  q = m .* ( 2.^(-r)*ones(1,size(m,2)) );

  % array m shifted right by r bits (last r bits lost)
  floorq = floor(q);

  % first component of m
  m(:,1) = floorq(:,1);

  % shifted array plus last r bits of previous array component
  m(:,2:end) = ( q(:,1:end-1) - floorq(:,1:end-1) ) * INTLAB_LONG_BETA + ...
               floorq(:,2:end);