MATLAB案例 | Copula的密度函数和分布函数图

本文介绍各种类型(Gaussian、t、Gumbel、Clayton、Frank)Copula的密度函数和分布函数图的绘制

完整代码

复制代码
clc
close all
clear

%% ********************计算Copula的密度函数和分布函数图************************
[Udata,Vdata] = meshgrid(linspace(0,1,31));  % 为绘图需要,产生新的网格数据
% 调用copulapdf函数计算网格点上的二元正态Copula密度函数值
Cpdf_norm = copulapdf('Gaussian',[Udata(:), Vdata(:)],0.7);
% 调用copulacdf函数计算网格点上的二元正态Copula分布函数值
Ccdf_norm = copulacdf('Gaussian',[Udata(:), Vdata(:)],0.7);

% 调用copulapdf函数计算网格点上的二元t-Copula密度函数值
Cpdf_t = copulapdf('t',[Udata(:), Vdata(:)],0.7,5);
% 调用copulacdf函数计算网格点上的二元t-Copula分布函数值
Ccdf_t = copulacdf('t',[Udata(:), Vdata(:)],0.7,5);

% 调用copulapdf函数计算网格点上的二元Gumbel-Copula密度函数值
Cpdf_Gumbel = copulapdf('Gumbel',[Udata(:), Vdata(:)],1.5);
% 调用copulacdf函数计算网格点上的二元Gumbel-Copula分布函数值
Ccdf_Gumbel = copulacdf('Gumbel',[Udata(:), Vdata(:)],1.5);

% 调用copulapdf函数计算网格点上的二元Clayton-Copula密度函数值
Cpdf_Clayton = copulapdf('Clayton',[Udata(:), Vdata(:)],1);
% 调用copulacdf函数计算网格点上的二元Clayton-Copula分布函数值
Ccdf_Clayton = copulacdf('Clayton',[Udata(:), Vdata(:)],1);

% 调用copulapdf函数计算网格点上的二元Frank-Copula密度函数值
Cpdf_Frank = copulapdf('Frank',[Udata(:), Vdata(:)],2);
% 调用copulacdf函数计算网格点上的二元Gumbel-Copula分布函数值
Ccdf_Frank = copulacdf('Frank',[Udata(:), Vdata(:)],2);

%% 绘制二元正态Copula的密度函数和分布函数图
% 绘制二元Gaussian-Copula的密度函数和分布函数图
figure(1);  % 新建图形窗口
subplot(1,2,1)
surf(Udata,Vdata,reshape(Cpdf_norm,size(Udata)));  % 绘制二元正态Copula密度函数图
xlabel('U');  % 为X轴加标签
ylabel('V');  % 为Y轴加标签
zlabel('c(u,v)');  % 为z轴加标签
title('二元正态Copula的密度函数')
subplot(1,2,2)
surf(Udata,Vdata,reshape(Ccdf_norm,size(Udata)));  % 绘制二元正态Copula分布函数图
xlabel('U');  % 为X轴加标签
ylabel('V');  % 为Y轴加标签
zlabel('C(u,v)');  % 为z轴加标签
title('二元正态Copula的分布函数')

% 绘制二元t-Copula的密度函数和分布函数图
figure(2);  % 新建图形窗口
subplot(1,2,1)
surf(Udata,Vdata,reshape(Cpdf_t,size(Udata)));  % 绘制二元t-Copula密度函数图
xlabel('U');  % 为X轴加标签
ylabel('V');  % 为Y轴加标签
zlabel('c(u,v)');  % 为z轴加标签
title('二元t-Copula的密度函数')
subplot(1,2,2)
surf(Udata,Vdata,reshape(Ccdf_t,size(Udata)));  % 绘制二元t-Copula分布函数图
xlabel('U');  % 为X轴加标签
ylabel('V');  % 为Y轴加标签
zlabel('C(u,v)');  % 为z轴加标签
title('二元t-Copula的分布函数')

% 绘制二元Gumbel-Copula的密度函数和分布函数图
figure(3);  % 新建图形窗口
subplot(1,2,1)
surf(Udata,Vdata,reshape(Cpdf_Gumbel,size(Udata)));  % 绘制二元t-Copula密度函数图
xlabel('U');  % 为X轴加标签
ylabel('V');  % 为Y轴加标签
zlabel('c(u,v)');  % 为z轴加标签
title('二元Gumbel-Copula的密度函数')
subplot(1,2,2)
surf(Udata,Vdata,reshape(Ccdf_Gumbel,size(Udata)));  % 绘制二元t-Copula分布函数图
xlabel('U');  % 为X轴加标签
ylabel('V');  % 为Y轴加标签
zlabel('C(u,v)');  % 为z轴加标签
title('二元Gumbel-Copula的分布函数')

% 绘制二元Gumbel-Copula的密度函数和分布函数图
figure(4);  % 新建图形窗口
subplot(1,2,1)
surf(Udata,Vdata,reshape(Cpdf_Clayton,size(Udata)));  % 绘制二元t-Copula密度函数图
xlabel('U');  % 为X轴加标签
ylabel('V');  % 为Y轴加标签
zlabel('c(u,v)');  % 为z轴加标签
title('二元Clayton-Copula的密度函数')
subplot(1,2,2)
surf(Udata,Vdata,reshape(Ccdf_Clayton,size(Udata)));  % 绘制二元t-Copula分布函数图
xlabel('U');  % 为X轴加标签
ylabel('V');  % 为Y轴加标签
zlabel('C(u,v)');  % 为z轴加标签
title('二元Clayton-Copula的分布函数')

% 绘制二元Gumbel-Copula的密度函数和分布函数图
figure(5);  % 新建图形窗口
subplot(1,2,1)
surf(Udata,Vdata,reshape(Cpdf_Frank,size(Udata)));  % 绘制二元t-Copula密度函数图
xlabel('U');  % 为X轴加标签
ylabel('V');  % 为Y轴加标签
zlabel('c(u,v)');  % 为z轴加标签
title('二元Frank-Copula的密度函数')
subplot(1,2,2)
surf(Udata,Vdata,reshape(Ccdf_Frank,size(Udata)));  % 绘制二元t-Copula分布函数图
xlabel('U');  % 为X轴加标签
ylabel('V');  % 为Y轴加标签
zlabel('C(u,v)');  % 为z轴加标签
title('二元Frank-Copula的分布函数')

绘制结果





相关推荐
liulilittle1 小时前
HTTP简易客户端实现
开发语言·网络·c++·网络协议·http·编程语言
程序员皮皮林1 小时前
Java jar 如何防止被反编译?代码写的太烂,害怕被人发现
java·开发语言·jar
qczg_wxg1 小时前
高阶组件介绍
开发语言·javascript·react native·ecmascript
CHANG_THE_WORLD2 小时前
C++ 并发编程指南 实现无锁队列
开发语言·c++·缓存·无锁队列·无锁编程
这里没有酒2 小时前
[C语言] 结构体 内存对齐规则 内存大小计算
c语言·开发语言
CHANG_THE_WORLD2 小时前
C++ 内存模型:用生活中的例子理解并发编程
开发语言·c++·生活
闯闯桑2 小时前
toDF(columns: _*) 语法
开发语言·前端·spark·scala·apache
滴滴滴嘟嘟嘟.2 小时前
Qt图表功能学习
开发语言·qt·学习
charlie1145141912 小时前
Android开发——初步了解AndroidManifest.xml
android·xml·开发语言·学习·安卓·安全架构
Edward.W2 小时前
用 Go + HTML 实现 OpenHarmony 投屏(hdckit-go + WebSocket + Canvas 实战)
开发语言·后端·golang