用matlab画Clayton Copula的密度函数图

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

我用matlab 中的 mesh画clayton的密度函数图，期望得到类似于帽子的密度函数图。
错误的密度函数图：

代码如下：
n=70;
x=linspace(-2,2,n);
y=linspace(-2,2,n);
u=normpdf(x');
v=normpdf(y');
[uu,vv]=meshgrid(u,v);
onen=ones(n);
r=1;
z=(1+r).*((uu.*vv).^(-r-1)).*((uu.^(-r))+(vv.^(-r))-onen).^((-2-r)/r); % Clayton copula密度函数表达式
mesh(uu,vv,z)


因为之前在论坛看到一位前辈画过Clayton Copula的 密度函数等高线图如下：

这是代码：
n=70;T = 100;
tic;
xx = (-2:4/(T-1):2)';
uu = normcdf(xx);
v = (0.02:0.03:0.2);
kappa = 1;
zz = normpdf(xx)*ones(1,T);% this the the part of the pdf from the X variable
tic;
for ii=1:T;
zz(:,ii) = zz(:,ii).*normpdf(xx(ii)).*clayton_pdf(uu,uu(ii),kappa);
x=linspace(-2,2,n);
y=linspace(-2,2,n);
end
u=normpdf(x');
v=normpdf(y');
[uu,vv]=meshgrid(u,v);
onen=ones(n);
r=1;
z=(1+r).*((uu.*vv).^(-r-1)).*((uu.^(-r))+(vv.^(-r))-onen).^((-2-r)/r); % Clayton copula密度函数表达式
contour(xx,xx,zz,v,'Color',[0 0.7 0]);
title('Clayton copula, \kappa = 1');


可以知道Clayton Copula是可以画出帽子那样形状的密度函数图的，大概形状应该是这样（用的二元正态分布作比较。。）：


求助！！！怎样画出Clayton Copula 帽子形状的密度函数图。致谢，并追加论坛币。

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏2 回帖

关键词：二元正态分布 函数表达式 密度函数 正态分布 等高线图

请问楼主解决了吗，我也遇到了同样的问题？