楼主: serenityz
11847 6

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

  • 0关注
  • 0粉丝

大专生

0%

还不是VIP/贵宾

-

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

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
拉您进交流群

京ICP备16021002-2号 京B2-20170662号 京公网安备 11010802022788号 论坛法律顾问:王进律师 知识产权保护声明   免责及隐私声明

GMT+8, 2024-4-30 15:43