用多项式拟合作图矩阵

Question

我有一个由我用plotmatrix创建的子图组成的矩阵，并且希望为每个子图添加最适合的行，但是我想不出怎么做。是否有一种方法分别适用于每个子图？

下面是我到目前为止所得到的一个简化的例子

x = randn(50,3);
y = x*[-1 2 1;2 0 1;1 -2 3;]';
[H,AX,BigAx,P,PAx] = plotmatrix(x,y);

我想我需要学习如何以某种方式对子图矩阵进行索引。

Answer 1

我不知道如何覆盖plotmatrix生成的现有矩阵，但是您可以自己创建矩阵，并用polyfit的结果覆盖每个子图

figure;
x = randn(50,3);
y = x*[-1 2 1;2 0 1;1 -2 3;]';

degree=4;

rows = size(x,2);
cols = size(y,2);

for k=1:rows
    for m=1:cols
        subplot(rows, cols, (k-1)*rows+m);
        hold all;
        scatter(x(:,k),y(:,m),'.');
        p = polyfit(x(:,k),y(:,m),degree);
        x_p = linspace(min(x(:,k)), max(x(:,k)));
        y_p = polyval(p,x_p);
        plot(x_p, y_p,'LineWidth',2);
        hold off;
    end
end

​

​