sd_round.m

我认为很不错的语音处理的ｍａｔｌａｂ源代码

function [A2, A_str, real_digitsL, real_digitsR, imag_digitsL, imag_digitsR]=sd_round(A, N, flag, mult)
% % sd_round: Rounds an array to a specified number of Significant Digits, significant figures, digits of precision
% % 
% % *********************************************************************
% % 
% % Syntax;
% % 
% % [A2, A_str, real_digitsL, real_digitsR, imag_digitsL,...
% % imag_digitsR]=sd_round(A, N, flag);
% % 
% % *********************************************************************
% % 
% % Description
% % 
% % sd_round stands for "Significant Digits Round".
% %
% % This program rounds a 2-d matrix of numbers to a specified number
% % of significant digits.  
% % 
% % This program support five different styles of rounding the last digit:
% % to the nearest integer, up, down, toward zero, and away from zero.  
% % 
% % This program supports real and complex numbers.  
% % 
% % The program outputs the rounded array, a cell string of the 
% % rounded matrix the number of digits, to the left and right of the
% % decimal place.  
% % 
% % This program is useful for presenting scientific data that 
% % requires rounding to a specifid number of significant digits 
% % for publication.  
% % 
% % Significant digits are counted starting from the first non-zero 
% % digit from the left.  
% % 
% %
% % *********************************************************************
% % 
% % Input variables
% %
% % A is the input matrix of number to be rounded.  
% %      default is empty array A=[];.
% % 
% % N is the number of significant digits. 
% %      default is N=3;
% % 
% % flag specifiies the style of rounding.  
% %      This program supports four different styles of rounding.
% %      flag == 1 rounds to the nearest integer 
% %      flag == 2 rounds up   
% %      flag == 3 rounds down 
% %      flag == 4 rounds toward zero
% %      flag == 5 rounds away from zero
% %      otherwise round to the nearest integer
% %      default is round to the nearest integer
% % 
% % mult is a whole number.  The program rounds the last digit to mult.  
% %      It is preferred that mult be between 1 and 9; however, all whole 
% %      numberS >= 1 are valid input.  The prorgam rounds mult to the 
% %      nearest integer and makes sure the value is at least 1.
% %      default is mult=1;
% % 
% % *********************************************************************
% %
% % Output variables
% % 
% % A2 is the rounded array.
% % 
% % A_str           % The rounded array is converted to a cell
% %                 % string format with the specified rounding and showing
% %                 %  the trainling zeros.  
% %                 % This is convenient for publishing tables in a tab 
% %                 % delimited string format
% % 
% % real_digitsL    % The number of real digits to the left of the decimal 
% %                 % point
% % 
% % real_digitsR    % The number of real digits to the right of the decimal
% %                 % point
% % 
% % imag_digitsL    % The number of imaginary digits to the left of the 
% %                 % decimal point
% % 
% % imag_digitsR 	% The number of imaginary digits to the right of the 
% %                 % decimal point
% %
% % *********************************************************************
% %
% 
% Example1=''; 
% D1=pi;        % Double or Complex two dimensional array of numbers
% N=3;          % Number of significant digits.  3 is the default
% 
% [P1, P1_str]=sd_round(pi, N);
% 
% % P1 should be 3.14 which has 3 significant digits.
%    
% Example2=''; 
% D1=pi/1000000;    % Double much smaller than 1
% N=3;              % Number of significant digits.  3 is the default
% flag=1;           % round to the nearest digit
% mult=5;           % round to a multiple 5
% 
% [P1, P1_str]=sd_round(D1, N, 1, 5);
% 
% % P1_str should be 0.00000315 which has 3 significant digits.
% % and the last digit is rounded to the nearest multiple of 5.  
% 
% Example3=''; 
% N=4;                          % N is the number of significant digits
% D2=10.^randn(10, 100);        % D2 is the unrounded array 
% [P2, P2_str]=sd_round(D2, N); % P2 is the rounded array 
%                               % of real numbers 
%                               % P2_str is the cell array of strings of  
%                               % rounded real numbers
% Example4=''; 
% D3=log10(randn(10, 100));     % D3 is an unrounded array of complex 
%                               % numbers
% [P3, P3_str]=sd_round(D3, 4); % P3 is the rounded array of 
%                               % complex numbers
%                               % P3_str is the cell array of strings of  
%                               % rounded complex numbers
%                         
% % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% % 
% % Program Written by Edward L. Zechmann
% %    
% %      date 23 November 2007
% %  
% %  modified 26 November 2007   updated comments
% %  
% %  modified 17 December 2007   added outputs
% %                              real_digitsL
% %                              real_digitsR
% %                              imag_digitsL
% %                              imag_digitsR
% %  
% %  modified 28 December 2007   added string output
% %                              fixed bug in real_digitsR
% %                              fixed bug in imag_digitsR
% %  
% %  modified  7 January  2008   fixed bug in real_digitsR
% %                              fixed bug in imag_digitsR
% %                              sped up program by only
% %                              converting the array to strings
% %                              if it is an output argument
% %     
% %  modified  1 March    2008   Added support for rounding 
% %                              to nearest integer, up, down, 
% %                              and toward zero
% % 
% %  modified  3 March    2008   updated comments
% % 
% %  modified 16 March    2008   Changed Program name from
% %                              p_round to sd_round.  
% %                              Added another rounding style 
% %                              flag =5; (away from 0).  
% %                              Updated comments.
% %                                         
% %  modified 18 August   2008   Added option to round last digit to a 
% %                              multiple of a given number. 
% %                              Fixed a bug in rounding powers of 10.
% %                              Improved examples.
% %  
% %  modified 21 August   2008   Fixed a bug in rounding numbers less 
% %                              than 1.  Added an example.  
% %  
% %  modified 25 August   2008   Modified program to recalculate the 
% %                              number of digits after rounding,
% %                              because rounding can change the number 
% %                              of digits to the left and right of the 
% %                              decimal place. Updated Comments
% %  
% % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% % 
% % Please Feel Free to Modify This Program
% %    
% % See Also: pow10_round, round, ceil, floor, fix, fix2, round2, roundd
% %

flag1=0;

if (nargin < 1 || isempty(A)) || ~isnumeric(A)
    flag1=1;
    A=[];
    A2=[];
    warningdlg('Error in p_round did not Input Matrix A');
end

if isempty(A)
    flag1=1;
    A=[];
    A2=[];
    warningdlg('Error in p_round Matrix A is Empty');
end

if (nargin < 2 || isempty(N)) || ~isnumeric(N)
    N=3;
end

if (nargin < 3 || isempty(flag)) || ~isnumeric(flag)
    flag=1;
end

if (nargin < 4 || isempty(mult)) || ~isnumeric(mult)
    mult=1; 
end

mult=round(mult);
if mult < 1
    mult=1;
end
    
% Number of digits to keep
N=round(N);

if ~isempty(A)
    real_digitsL=zeros(size(A));
    real_digitsR=zeros(size(A));
    imag_digitsL=zeros(size(A));
    imag_digitsR=zeros(size(A));
else
    real_digitsL=[];
    real_digitsR=[];
    imag_digitsL=[];
    imag_digitsR=[];
end

if isequal(flag1, 0)
    if isnumeric(A)

        if isreal(A)
            
            % Digit furthest to the left of the decimal point
            D1=ceil(log10(abs(A)));
            buf1=D1( abs(A)-10.^D1 == 0)+1;
            D1( abs(A)-10.^D1 == 0)=buf1;
            
            % Number of digits to the left of the decimal place
            real_digitsL=max(D1, 0);
            real_digitsL(A==0)=0;
            
            % rounding factor
            dec=10.^(N-D1);
            
            % Number of digits to the right of the decimal place
            real_digitsR=max(N-D1, 0);
            real_digitsR(A==0)=N;
            
            % Rounding Computation
            % This program supports five different styles of rounding.
            % flag == 1 rounds to the nearest integer