在MatLab中绘制标准化价格的3D图

debugcn 发表于 Dev

安东尼奥·乔尔

我正在做高斯过程，我从给定的矩阵（每一行代表一年）计算每年的回归，因此代码是：

M1 = MainMatrix; %This is the given Matrix 
ker =@(x,y) exp(-1013*(x-y)'*(x-y)); 
[ns, ms] = size(M1);

for N = 1:ns
    x = M1(N,:);
    C = zeros(ms,ms);
    for i = 1:ms 
        for j = 1:ms 
            C(i,j)= ker(x(i),x(j));
        end
    end

    u = randn(ms,1); 
    [A,S, B] = svd(C);
    z = A*sqrt(S)*u; % z = A S^.5 u

我想在Graph 3D中绘制每个回归，如下所示：

在此处输入图片说明

我知道情节是一个ribbon，但我不知道该怎么做

马特

无需使用即可生成所需的图ribbon。只需surf对所有价格使用fill3-plot，对z = 0处的飞机使用-plot 。平面的边界是根据图的实际极限来计算的。因此，我们需要在绘制平面之前设置极限。然后只需进行一些调整即可生成几乎相同的外观。

这是代码：

% generate some data
days  = (1:100)';
price = days*[0.18,-0.08,0.07,-0.10,0.12,-0.08,0.05];
price = price + 0.5*randn(size(price));
years = 2002+(1:size(price,2));

% prepare plot
width = 0.6;
X = ones(size(price,1),1)*0.5;
X = [-X,X]*width;
figure; hold on;

% plot all 'ribbons'
for i = 1:size(price,2)
    h = surf([days,days],X+years(i),[price(:,i),price(:,i)]);
    set(h,'MeshStyle','column');
end

% set axis limits
set(gca,'ZLim',[-20,20]);

% plot plane at z=0
limx = get(gca,'XLim');
limy = get(gca,'YLim');
fill3(reshape([limx;limx],1,[]),[flip(limy),limy],zeros(1,4),'g','FaceAlpha',0.2)

% set labels
xlabel('Day of trading')
ylabel('Year')
zlabel('Normalized Price')

% tweak appearance
set(gca,'YTick',years);
set(gca,'YDir','reverse');
view([-38,50])
colormap jet;
grid on;
%box on;

结果如下： 结果1