用matlab 画几种copula 族函数的密度函数等高线

我用matlab 中的 contour画clayton的密度函数等高线，期望得到类似于年轮的等高线图，目前只画出了分布函数等高线，并测试根据clayton的参数alpha的不断变化，分布函数等高线由椭圆形趋于矩形，棱角逐渐鲜明。



错误的密度函数等高线

我并不是懒惰，但是几经测试依然无法得到圆环状的等高线，也再没有找出自己的问题所在，请明白的高手指点一二。

%年轮等高线测试
clear all
clc
n=700;
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密度函数表达式
contour(uu,vv,z)

应该出现的等高线

求助！如果画出，也请讲解一下，个别代码编写作用。致谢，并追加论坛币。

http://www.mathworks.com/access/ ... ats/brn2ivz-93.html
上面有多元t的，自己照猫画虎吧
good luck!

2# qwaszxgecl


感谢你的回复，但 copula 函数和多元固定函数分布是完全不同的概念。

我已经找到办法，需要在带入边际函数的时候使用循环，

for ii=1:T;
   zz(:,ii) = zz(:,ii).*normpdf(xx(ii)).*clayton_pdf(uu,uu(ii),kappa);

end

边际函数为正态分布
clayton copula

kappa = 1;
zz = normpdf(xx)*ones(1,T);  

for ii=1:T;
   zz(:,ii) = zz(:,ii).*normpdf(xx(ii)).*clayton_pdf(uu,uu(ii),kappa);
end

contour(xx,xx,zz,v,'Color',[0 0.7 0])；
    title('Clayton copula, \kappa = 1');

4# serenityz

hi, 你好！
我看到了你之前写的帖子，后来你说你解决了那个问题，可是后面给出的代码很凌乱，可以问你要一下完整的可以画出圆环状等高线的代码吗？
或者说，可以问几个问题吗？
kappa = 1;
zz = normpdf(xx)*ones(1,T);  
for ii=1:T;
   zz(:,ii) = zz(:,ii).*normpdf(xx(ii)).*clayton_pdf(uu,uu(ii),kappa);
end
contour(xx,xx,zz,v,'Color',[0 0.7 0])；
title('Clayton copula, \kappa = 1');
这段代码中，T的取值是多少？xx是什么变量？clayton_pdf()这个函数是哪里来的，可以使用copulapdf('Clayton')这个函数替代么？打扰你了，谢谢。 ：）

你好，我在人大经济论坛看到了你的关于画copula等高线的帖子，具体怎么画能指导一下吗。

10年的帖子，23年搜出来了，而今我也碰到了一样的问题，纪念一下