function y=quling(a)     % 該程序把矩陣中全為0的行和列刪掉
[m,n]=size(a);
jihao1=[];
jihao2=[];
for k=n:-1:1
    if a(:,k)==0
        jihao2=[jihao2,k];
    else
        break;
    end
end
for k=m:-1:1
    if a(k,:)==0
        jihao1=[jihao1,k];
    else
        break;
    end
end
a(:,jihao2)=[];
y=a;
y(jihao1,:)=[];