msort.sci

小波分解源代码

function y=msort(x)
//  Copyright Aldo I Maalouf
m=size(x,1);
n=size(x,2);
if m~=1&n~=1
  for j=1:n
    for i=1:m-1;
     if x(i,j)>x(i+1,j)
       b=x(i,j);
       x(i,j)=x(i+1,j);
       x(i+1,j)=b;
     end
    end
  end
  for i=1:m
    for j=1:n
       y(i,j)=x(i,j);
    end
  end
  
elseif m==1 
    for i=1:n
      for j=(i+1):n
        if x(i)>x(j)
           b=x(i);
           x(i)=x(j);
           x(j)=b
        end
      end
     end
     
     for i=1:n
       y(i)=x(i);
     end
       y=y';
elseif n==1
     for i=1:m
       for j=i:m
         if x(i)>x(j)
            b=x(i);
            x(i)=x(j);
            x(j)=b
          end
       end
     end
end     
     for i=1:m
       y(i)=x(i);
     end
endfunction