求救：有關 Copula contour plot

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

火燒眉毛了，有人可以用 matlab 畫出 Copula contour plot 嗎??

因為式子有積分一直畫不出向右圖的結果，有神人可以幫個忙嗎？？

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：contour Copula opula plot Tour 眉毛 matlab

這是我寫的程式碼和圖形，但是就無法和原圖長的一樣，有人能告訴我少了什麼步驟嗎？？

r=0.0421;
q=0;
theta1=0.1; theta2=0.15;
sigma1=0.2; sigma2=0.25;
nu=1;
omega1=(1/nu)*( log(1-theta1*nu-sigma1*sigma1*nu/2));
m1=(r-q+omega1);　　%上面程式的裡面的參數m
omega2=(1/nu)*( log(1-theta2*nu-sigma2*sigma2*nu/2));
m2=(r-q+omega2);
%%%
c=zeros(99,99);
u=0.01:0.01:0.99;
for i=1:99
    F1=@(z1)quad(@(x)normcdf((z1-m1*1-theta1.*x)./(sigma1*sqrt(x))).*gampdf(x,1/nu,nu),0.001,100)-u(i);
    x1(i)=fzero(F1,u(i));
end
for j=1:99
    F2=@(z2)quad(@(x)normcdf((z2-m2*1-theta2.*x)./(sigma2*sqrt(x))).*gampdf(x,1/nu,nu),0.001,100)-u(j);
    x2(j)=fzero(F2,u(j));
end   
for i = 1:99
    for j = 1:99
        c(i,j)=quad(@(x)normcdf((x1(i)-m1*1-theta1.*x)./(sigma1*sqrt(x))).*normcdf((x2(j)-m2*1-theta2.*x)./(sigma2*sqrt(x))).*gampdf(x,1/nu,nu),0.001,100);
    end
end
contour(u,u,c)
axis([0, 1, 0, 1]);axis normal; %axis normal; axis square;
title('\theta_1=0.1,\theta_2=0.15,\sigma=0.2,\sigma=0.25,\nu=1');