Я пытаюсь построить контуры для переменной "T", и должно получиться что-то вроде этого рисунка:
ИЗОБРАЖЕНИЕ ОБРАЗЦОВ СОВМЕСТНОГО СОЕДИНЕНИЯ Сам T является функцией WL и Hs. Ниже мой код:
clc
clear all
WL = gevrnd(-0.4176,0.6703,1.4837,100000,1);
Hs = gprnd(-0.0960,0.2152,1.7001,100000,1);
U = cdf('generalized extreme value',WL,-0.4176,0.6703,1.4837);
V = cdf('Generalized Pareto',Hs,-0.0960,0.2152,1.7001);
mu = 0.08;
for i=1:length(WL)
for j=1:length(Hs)
Cc(i,j) = exp(1-((((1-log(U(i))).^0.9838-1).^1.0075+((1-log(V(j))).^0.9838-1)...
.^1.0075).^0.9926+1).^1.0165);
test(i,j)=1-U(i)-V(j)+Cc(i,j);
T(i,j)=mu/test(i,j);%year
end
end
v=[1 25];
[C,h]=contour(WL,Hs,T,v);
clabel(C,h);
Кто-нибудь имел опыт такого анализа?