MATLAB案例 | 基于Copula的可靠度分析

半撅儿

已于 2024-10-17 00:25:56 修改

阅读量1.4k

点赞数 22

于 2024-09-26 00:36:08 首次发布

本文链接：https://blog.csdn.net/ikhui7/article/details/142533121

版权

本文详细介绍了Copula函数的绘制及在可靠度分析中的应用

各种类型Copula函数绘图

在这里插入图片描述

完整代码

clear
clc

y_gaussian = copularnd('gaussian', 0.9, 1000);
y_t = copularnd('t', 0.91, 17.53,1000);
y_Gumbel = copularnd('Gumbel', 3.22,1000);
y_Clayton = copularnd('Clayton', 3.06,1000);
y_Frank = copularnd('Frank', 12.04,1000);

subplot(2,3,1.4)
plot(y_gaussian(:,1), y_gaussian(:,2),'.')
xlabel('u')
ylabel('v')
title('Gaussian(\theta = 0.90)')
set(gca,'FontName','Times New Roman','FontSize',14);
subplot(2,3,2.6)
plot(y_t(:,1), y_t(:,2),'.')
xlabel('u')
ylabel('v')
title('t(\theta = 0.90,  \nu =17.53)')
set(gca,'FontName','Times New Roman','FontSize',14);
subplot(2,3,4)
plot(y_Gumbel(:,1), y_Gumbel(:,2),'.')
xlabel('u')
ylabel('v')
title('Gumbel(\theta = 3.22)')
set(gca,'FontName','Times New Roman','FontSize',14);
subplot(2,3,5)
plot(y_Clayton(:,1), y_Clayton(:,2),'.')
xlabel('u')
ylabel('v')
title('Clayton(\theta = 3.06)')
set(gca,'FontName','Times New Roman','FontSize',14);
subplot(2,3,6)
plot(y_Frank(:,1), y_Frank(:,2),'.')
xlabel('u')
ylabel('v')
title('Frank(\theta = 12.04)')
set(gca,'FontName','Times New Roman','FontSize',14);