我想创建变量概率密度函数的三维图。我有一个矩阵,大小为189x10000,其中行对应于时间,列是模拟结果。有人可以帮助我创建一个随时间变化的密度图吗?我希望我的图像如下所示:
使用以下代码可以得到下面的效果:
![enter image description here](https://istack.dev59.com/LH9zA.webp)
A = [1:185]'; % substitute for date vector
K = linspace( -20, 20, 100);
f = zeros(185,100);
xi = zeros(185,100);
r = normrnd(0,1,[185,10000]);
for i=1:185
[f(i,:),xi(i,:)] = ksdensity(r(I,:));
end
a = figure;
meshc(A, K', f')
datetick('x', 'yyyy')
view(85, 50)
set(gca, 'YLim', [-15, 10])
set(gca, 'XLim', [A(1), A(end)])
xlabel('Time')
使用以下代码可以得到下面的效果: