contour plot for different values

6 次查看（过去 30 天）
Beena
Beena 2014-7-31
回答： Mike Garrity 2014-8-1
Hi I am using this code to plot contours of normal marginals
A=[1 2;3 4];
B=[3 8;5 9];
MA=mean(A,2);
MB=mean(B,2);
rho=corrcoef(A,B);
MU1 = MA;
SIGMA1 = rho;
MU2 = MB;
SIGMA2 = [1 0; 0 1];
u=linspace(-8,8,1000);
[X,Y]= meshgrid(u); % mashgrid of u
F1 = mvnpdf([X(:),Y(:)],MU1',SIGMA1);%pdf of marginals
F2 = mvnpdf([X(:),Y(:)],MU2',SIGMA2);%pdf of marginals
F1 = reshape(F1,length(Y),length(X));
F2 = reshape(F2,length(Y),length(X));
Fa=mvncdf([X(:),Y(:)],MU1',SIGMA1); %cdf of marginals
Fb=mvncdf([X(:),Y(:)],MU1',SIGMA1); %cdf of marginal
F = copulapdf('Gaussian',[Fa(:), Fb(:)],0);% Gaussian copula of marginals
F = reshape(F,length(Y),length(X));
mvd = reshape(F,1000,1000).*F1.*F2; % computing multivariate density of copula and marginals
figure;
contour(X,Y,mvd);
I want to plot a contour plot of different values of F1 and F2 but when I am using contour (X,Y,mvd,[val1,val2]); it just gives one line of the contour plot. But I want the complete contour plot of using just two values of coordinates.Please help with this Thanks in advance for your help.
kind regards Ben

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

回答（1 个）

Mike Garrity
Mike Garrity 2014-8-1
What you're doing looks fine, but I don't know what val1 and val2 were. If I do this:
contour(X,Y,mvd,[.04e-3, .12e-3])
I get two levels.
If you do this:
[~,h] = contour(X,Y,mvd)
then h will be a handle to a contour object. You can then ask it what its doing and try different levels.
For example,
get(h,'LevelList')
ans =
1.0e-03 *
0.02 0.04 0.06 0.08 0.10 0.12 0.14 0.16 0.18
Do those values look like the ones you were passing in?
You can then change the levels like this:
set(h,'LevelList',[.03 .09]*1e-3)

类别

在 Help Center 和 File Exchange 中查找有关 Contour Plots 的更多信息

标签

产品

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by