Contour plots for random variables

Question

soheil radfar 2020-1-3

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/498846-contour-plots-for-random-variables

评论： soheil radfar 2020-1-4

Hi everybody,

I am trying to plot contours for variable "T". T itself is a function of X and Y. Below is my code, which is clearly wrong!

clc
clear all
X = gevrnd(-0.4176,0.6703,1.4837,100000,1); %X follows GEV distribution
Y = gprnd(-0.0960,0.2152,1.7001,100000,1); %Y follows GPD distribution
U = cdf('generalized extreme value',X,-0.4176,0.6703,1.4837);
V = cdf('Generalized Pareto',Y,-0.0960,0.2152,1.7001);
C = exp(1-((((1-log(U)).^0.9838-1).^1.0075+((1-log(V)).^0.9838-1).^1.0075).^0.9926+1).^1.0165); % Copula
% Return period: T = mu / (1+C-U-V)
mu = 0.08;
T = @(U,V) mu./(1+C-U-V);
fc = fcontour(T,U,V)
fc.LevelList = [1 5 10 25]; %Joint contours for T = 1; 5; 10 and 25 years

How shouhd I do this?

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

per isakson 2020-1-3

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/498846-contour-plots-for-random-variables#answer_408609

编辑：per isakson 2020-1-3

在 MATLAB Online 中打开

The documentation on fcontour, Plot contours says:

Specify a function of the form z = f(x,y). The function must accept two matrix input arguments and return a matrix output argument of the same size.

The vector C is the problem

>> T = @(U,V) mu./(1-U-V);
>> fcontour( T, [2.5,3.5,1.8,3.2] )

produces a plot without throwing an error.

Try this script

%%
Z = T(U,V);
Finterp = scatteredInterpolant( X, Y, Z );
xmin = 2.5; % min(X);
xmax = 3.5; % max(X);
ymin = min(Y);
ymax = max(Y);
N = 1000;
xrange = linspace(xmin, xmax, N);
yrange = linspace(ymin, ymax, N);
[XGrid, YGrid] = meshgrid(xrange, yrange);
ZGrid = Finterp(XGrid, YGrid);
contourf(XGrid, YGrid, ZGrid);
xlabel('X');
ylabel('Y');
zlabel('T');