Warning: Rank deficient, rank = 0,问题

飞云

下面是我的函数文件：
function f=trytwo(A,x)
a=A(1);b=A(2);
f=1/3.*(3/5000.*x/b^4*exp(-1)-3/1250.*x/b^4/a*(-4/a-1)*exp(-1)+3/625.*x/b^4/a*exp(-1)).*((b.^(-a)).^3-((x/b).^a).^3);
然后在MATLAB中运行：
x=[...];
y=[...];
A0=[0.00001,0.00001];
[A,resnorm]=nlinfit(x,y,'trytwo',A0)
运行后就显示下面的警告：
>> [A,resnorm]=nlinfit(x,y,'trytwo',A0)
Warning: Rank deficient, rank = 0,  tol =  -1.#INDe+000.
> In nlinfit at 161
把后面的 .*((b.^(-a)).^3-((x/b).^a).^3)除去后就能进行计算，但在计算公式中它又是必要部分，请问能如何解决？谢谢！

[ 本帖最后由 sigma665 于 2008-11-6 09:45 编辑 ]

科技在线

rank = 0说明的这个矩阵是零秩的，还能算？

sigma665

x=[...];
y=[...];

这两行是啥意思？