楼主: serenityz
12936 6

用matlab 画几种copula 族函数的密度函数等高线 [推广有奖]

  • 0关注
  • 0粉丝

大专生

0%

还不是VIP/贵宾

-

威望
0
论坛币
3618 个
通用积分
0.0000
学术水平
0 点
热心指数
0 点
信用等级
0 点
经验
812 点
帖子
30
精华
0
在线时间
41 小时
注册时间
2007-12-4
最后登录
2023-11-22

楼主
serenityz 发表于 2009-12-18 03:34:37 |AI写论文
6论坛币
我用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)

应该出现的等高线

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

关键词:MATLAB Copula atlab opula matla Copula Density contour
Never Leave that till tomorrow.

沙发
qwaszxgecl 发表于 2009-12-21 02:53:10
http://www.mathworks.com/access/ ... ats/brn2ivz-93.html
上面有多元t的,自己照猫画虎吧
good luck!

藤椅
serenityz 发表于 2010-1-2 04:14:50
2# qwaszxgecl


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

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

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

end

边际函数为正态分布
clayton copula
Never Leave that till tomorrow.

板凳
serenityz 发表于 2010-1-2 04:17:12
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');
Never Leave that till tomorrow.

报纸
nacsgunt 发表于 2010-12-5 10:19:27
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')这个函数替代么?打扰你了,谢谢。 :)

地板
zhangwang1943 发表于 2013-6-13 12:09:45
你好,我在人大经济论坛看到了你的关于画copula等高线的帖子,具体怎么画能指导一下吗。

7
乌昂王0 发表于 2023-3-7 15:45:43
10年的帖子,23年搜出来了,而今我也碰到了一样的问题,纪念一下

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注cda
拉您进交流群
GMT+8, 2025-12-25 22:55