数学建模实录（3）

1。函数名：CrtPower(a)作用：计算矩阵中的向量权优点：直接由原始数据可以换算出矩阵中的向量权源代码：function rstPower=CrtPower(a)%**********************************************************% This is Help Information About Power() Functi

孤剑

1832人浏览 · 2004-08-30 16:55:00

孤剑 · 2004-08-30 16:55:00 发布

1。函数名：CrtPower(a)

作用：计算矩阵中的向量权

优点：直接由原始数据可以换算出矩阵中的向量权

源代码：

function rstPower=CrtPower(a)
%**********************************************************
% This is Help Information About Power() Function.
%           Find the min number in the matrix.
%         Verison:1.2.0 Finish Date:30/08/2004
%   Usage:
%       Power(a)      %a is matrix .
%       return a rstPower.
%***********************************************************

if nargout>1
    error('Too many output arguments!');
else
if (nargin<1 | nargin>=2)
        error('Too many input arguments!');
    else
        SumCols=length(a(1,:));%Get the Cols of the matrix.
  if SumCols>=3
            SumRows=length(a(:,1));%Get the Rows of the matrix.
   rstPower=zeros(SumRows);
      for Col=2:SumCols-1
                ParCol=a(:,Col);
                BefCol=a(:,Col-1);
                AftCol=a(:,Col+1);
    for RowCount=1:SumRows;
                    for tmpRow=1:SumRows;
                        if RowCount==tmpRow
                            rstPower(RowCount,tmpRow)=0;
                        else
                            if ParCol(RowCount)-BefCol(tmpRow)>=0
                                tmpM=ParCol(RowCount)-BefCol(tmpRow);
                                if AftCol(RowCount)-ParCol(tmpRow)>=0;
                                    tmpN=AftCol(RowCount)-ParCol(tmpRow);
                                else
                                    tmpN=AftCol(RowCount)-ParCol(tmpRow);
                                end;
                                rstPower(RowCount,tmpRow)=rstPower(RowCount,tmpRow)+abs(tmpM)+abs(tmpN)-(Col-2)*abs(tmpM);
                            else
                                tmpM=(SumCols-Col+1)*(ParCol(RowCount)-BefCol(tmpRow));
                                if AftCol(RowCount)-ParCol(tmpRow)>=0;
                                    tmpN=AftCol(RowCount)-ParCol(tmpRow);
                                    rstPower(RowCount,tmpRow)=rstPower(RowCount,tmpRow)+abs(tmpM+tmpN)-(Col-2)*abs(tmpM);
                                else
                                    tmpN=AftCol(RowCount)-ParCol(tmpRow);
                                    rstPower(RowCount,tmpRow)=rstPower(RowCount,tmpRow)+abs(tmpM)+abs(tmpN)-(Col-2)*abs(tmpM);
                                end;
                            end;
                        end;%End if
                     end;%End tmpRow for
                end;%End RowCount for
   end;%End Col For
        else
            error('Length of the martrix is too short!The minium number is 3!');
        end; %End if SumCols>=3
    end;%End If (nargin<1 | nargin>=2)
end;

2。函数名：FindMinParam(a,iCount,Diff)

作用：查找矩阵中指定的iCount个最小数。查找同行中（Diff=0，默认值）、不同行（Diff=1）中最小值

优点：直接给出数据就可以查找出矩阵中较小元素所在的行和列

源代码：

function rst=FindMinParam(a,iCount,Diff)
%**********************************************************
% This is Help Information About FindMinParam() Function.
% Find the min number in defferent Rows and Cols the matrix.
%           Verison:1.1.2 Finish Date:28/08/2004
%   Usage:
%       FindMinParam(a,iCount,Diff)
%   a is matrix .
%   iCount is Counter.
%   Diff is the parame to Find the mininum in Different Row.
%***********************************************************

if nargout>1
    error('Too many output arguments!');
else
if (nargin==0 | nargin>3)
        error('Too many input arguments!');
else
        Cols=length(a(1,:));
  Rows=length(a(:,1));
        if nargin==1
            iCount=Rows * Cols;
            Diff=0;
        end;
        if nargin==2 | nargin==3
            if iCount>Rows * Cols
                error('The search number is too big!');
                iCount=Rows * Cols;
            elseif iCount<1
                error('The mininum is 1');
                iCount=1;
            else
                iCount=iCount;
            end;
            if nargin==3
                if Diff==1
                    Diff=1;
                else
                    Diff=0;
                end;
            else
                Diff=0;
            end;
        end;

        rst=zeros(iCount,3);
  for Count=1:iCount
            Succ=0;
            for RowCount=1:Rows;
                for ColCount=1:Cols;
           if (min(min(a))==a(RowCount,ColCount));
                        tmpMin=min(min(a));
            tmpRow=RowCount;
            tmpCol=ColCount;
                        if Diff==1
                            a(RowCount,:)=inf;
                a(:,ColCount)=inf;
                        else
                            a(RowCount,ColCount)=inf;
                        end;
                        Succ=1;
                        break;
                    end;
                end     %End For
                if Succ==1
                    break;
                end;
            end;
            if Succ==1
                rst(Count,1)=tmpMin;
                rst(Count,2)=tmpRow;
                rst(Count,3)=tmpCol;
            end;
  end;
  disp(rst);
end;
end;