mmdigit.m

基本的源程序

function y=mmdigit(x,n,b,t)
% 演示输入输出参数个数判断函数
% 舍入函数，以b为底，将x舍入到n
% 若b未给定，那么b=10，以十进制舍入
% 如果x为复数，那么实部和虚部分别进行舍入
% 其他的舍入选项，fix, ceil, floor

% 判断输入参数的个数
if nargin<2
    error('输入参数个数太少')
elseif nargin==2
    b=10;
    t='round';
elseif nargin==3
    t='round';
end

n=round(abs(n(1)));
if isempty(b), b=10;
else b=round(abs(b(1)));
end

if isreal(x)
    y=abs(x)+(x==0);
    e=floor(log(y)./log(b)+1);
    p=repmat(b,size(x)).^(n-e);
    if strncmpi(t,'round',1)
        y=round(p.*x)./p;
    elseif strncmpi(t,'fix',2)
        y=fix(p.*x)./p;
    elseif strncmpi(t,'ceil',2)
        y=ceil(p.*x)./p;        
    elseif strncmpi(t,'floor',2)
        y=floor(p.*x)./p; 
    else
        error('未知舍入错误');
    end
else
    y=complex(mmdigit(real(x),n,b,t),mmdigit(imag(x),n,b,t));
end