求问各位大佬这个程序应该怎么改

等式坐标的数组维度与右边不匹配，建议将Feats改为cell数组，即17行改为：

Feats{k} = sort( diag(D),'desend' );