ind.m

来自「薛定宇高等应用数学问题的MATLAB求解例子源码」· M 代码 · 共 25 行

function y=ind(a,x)
[p,q]=size(x); [ap,aq]=size(a); z=[];
for u=1:q, z=cat(2,z,u); end
tt=setdiff(z,a); [ttp,ttq]=size(tt);
for r=ttq:-1:1, x(:,tt(r))=[-1]; end
for r=q:-1:1,
    if x(1,r)==-1, x(:,r)=[]; end
end
%i=find(x(1,q:-1:1)==-1); x(:,i)=[];
y=x;
    
for i=1:p, v(i)=codea(aq,y(i,:),10); end 
y=v'; [yy,I]=sort(y); y=[yy I];
[b,k,l]=unique(yy); y=[l I]; m=max(l); aa=zeros(m,p);
for ii=1:m
    for j=1:p
       if l(j)==ii, aa(ii,j)=[I(j)]; end
   end
end
y=aa;

function yy=codea(a,x,b)
yy=0;
for i=1:a, yy=yy+x(i)*b^(a-i); end