循环后画曲面

zw123zw

x1=[-0.486 -0.377 -0.298 -0.192 -0.057 -0.052 0.162 0.271 0.364 0.486];
y1=[0.065 0.061 0.025 0.087 0.087 0.078 0.018 0.092 0.065 0.048];
t=1;
for deltax=-100:100;
    for deltay=-100:100;
        x=x1+deltax;
        y=y1+deltay;
        n=10;
        for i=1:n;
           XX(i)=x(i)^2;
           YY(i)=y(i)^2;
           XY(i)=x(i)*y(i);
        end
        xx=sum(XX);
        yy=sum(YY);
        xy=sum(XY);
        xt=sum(x);
        yt=sum(y);
        h=[xx xy xt;xy yy yt;xt yt n];
        H=2*h;            
        CN=cond(H,2);
        LG=log10(CN);
        A(t)=LG;
          t=t+1;
           end
          end
surf(deltax,deltay,A)
无法生成曲面 ，请指教啊

winner245

你应该把A存成一个二维矩阵，而不是一个向量
surf 要求 A 是一个矩阵

zw123zw

为什么是二维不是一维？

winner245

zw123zw 发表于 2014-1-5 12:42
为什么是二维不是一维？

因为 surf 默认是对一个矩形网格画曲面图（对于非矩形网格曲面，需要变成映射实现），所以，默认 A 是一个矩阵，矩阵的行对应横轴方向，列对应纵轴方向。这个在help里也介绍了，建议读读help。

直接修改你的代码如下

x1=[-0.486 -0.377 -0.298 -0.192 -0.057 -0.052 0.162 0.271 0.364 0.486];
y1=[0.065 0.061 0.025 0.087 0.087 0.078 0.018 0.092 0.065 0.048];
deltax=-100:100;
deltay=-100:100;
A = zeros(length(deltay),length(deltax));
for m = 1:length(deltay)
    for n = 1:length(deltax)
        x=x1+deltax(n);
        y=y1+deltay(m);
        n=10;
        for i=1:n;
            XX(i)=x(i)^2;
            YY(i)=y(i)^2;
            XY(i)=x(i)*y(i);
        end
        xx=sum(XX);
        yy=sum(YY);
        xy=sum(XY);
        xt=sum(x);
        yt=sum(y);
        h=[xx xy xt;xy yy yt;xt yt n];
        H=2*h;
        CN=cond(H,2);
        LG=log10(CN);
        A(m,n)=LG;
    end
end
surf(deltax,deltay,A)

zw123zw

哦 我就是看help A要是201×201矩阵就可以了  还是谢谢你