%********************绘制 Copula 的密度函数和分布函数图************************
[Udata,Vdata] = meshgrid(linspace(0,1,31)); % 为绘图需要,产生新的网格数据
% 调用 copulapdf 函数计算网格点上的二元正态 Copula 密度函数值
Cpdf_norm = copulapdf('Gaussian',[Udata(:), Vdata(:)],rho_norm);
% 调用 copulacdf 函数计算网格点上的二元正态 Copula 分布函数值
Ccdf_norm = copulacdf('Gaussian',[Udata(:), Vdata(:)],rho_norm);
% 调用 copulapdf 函数计算网格点上的二元 t-Copula 密度函数值
Cpdf_t = copulapdf('t',[Udata(:), Vdata(:)],rho_t,nuhat);
% 调用 copulacdf 函数计算网格点上的二元 t-Copula 分布函数值
Ccdf_t = copulacdf('t',[Udata(:), Vdata(:)],rho_t,nuhat);
% 绘制二元正态 Copula 的密度函数和分布函数图
figure; % 新建图形窗口
surf(Udata,Vdata,reshape(Cpdf_norm,size(Udata))); % 绘制二元正态 Copula 密度函数图
xlabel('U'); % 为 X 轴加标签
ylabel('V'); % 为 Y 轴加标签
zlabel('c(u,v)'); % 为 z 轴加标签
figure; % 新建图形窗口
surf(Udata,Vdata,reshape(Ccdf_norm,size(Udata))); % 绘制二元正态 Copula 分布函数图
xlabel('U'); % 为 X 轴加标签
ylabel('V'); % 为 Y 轴加标签
zlabel('C(u,v)'); % 为 z 轴加标签
% 绘制二元 t-Copula 的密度函数和分布函数图
figure; % 新建图形窗口
surf(Udata,Vdata,reshape(Cpdf_t,size(Udata))); % 绘制二元 t-Copula 密度函数图
xlabel('U'); % 为 X 轴加标签
ylabel('V'); % 为 Y 轴加标签
zlabel('c(u,v)'); % 为 z 轴加标签
figure; % 新建图形窗口
surf(Udata,Vdata,reshape(Ccdf_t,size(Udata))); % 绘制二元 t-Copula 分布函数图
xlabel('U'); % 为 X 轴加标签
ylabel('V'); % 为 Y 轴加标签
zlabel('C(u,v)'); % 为 z 轴加标签
%**************求 Kendall 秩相关系数和 Spearman 秩相关系数***********************
% 调用 copulastat 函数求二元正态 Copula 对应的 Kendall 秩相关系数
Kendall_norm = copulastat('Gaussian',rho_norm)
% 调用 copulastat 函数求二元正态 Copula 对应的 Spearman 秩相关系数
Spearman_norm = copulastat('Gaussian',rho_norm,'type','Spearman')
% 调用 copulastat 函数求二元 t-Copula 对应的 Kendall 秩相关系数
Kendall_t = copulastat('t',rho_t)
% 调用 copulastat 函数求二元 t-Copula 对应的 Spearman 秩相关系数
评论1