Hi there,
I have this figure that you can see in the picture. Is it possible to change the coloring in the legend so that you can see the difference between the two like in the picture?
Thanks.
MATLAB: Coloring in plots,legend
coloringfigureMATLABplot
Related Solutions
Hello Aishwarya,
The problem here is that the array x is 2x31 and does not have every possible combination of an element of X1 and an element of X2, which is 2x961. The method 2 code below fixes that by implementing the same idea as in method 1,
X = [X1(:) X2(:)];For method 2 to be correct you need the factor in front to be to the +1/2 power, since you have already put the factor into the denominator.
mu = [0,0]; sigma = [0.5 0.8; 0.8 2.0]; x1 = -3:0.2:3; x2 = -3:0.2:3; X = [x1; x2] [X1,X2] = meshgrid(x1,x2); X = [X1(:) X2(:)]; y = mvnpdf(X,mu,sigma); y = reshape(y,length(x1),length(x2)); figure(1) surf(X1,X2,y) xlabel('x') ylabel('y') zlabel('Bivariate Gaussian Distributions') % method 2
mu = [0;0]; covar = [0.5 0.8; 0.8 2.0]; x1 = -3:0.2:3; x2 = -3:0.2:3; [X1,X2] = meshgrid(x1,x2); X = [X1(:), X2(:)]'; arg = sum((X-mu).*(pinv(covar)*(X-mu))); y2 = 1/(2*pi*(det(covar))^(1/2)) * exp(-(1/2)*arg); y2 = reshape(y2,size(X1)); figure(2) surf(X1,X2,y2) xlabel('x') ylabel('y') zlabel('Bivariate Gaussian Distributions') max(max(abs(y-y2))) ans = 1.8041e-16% FUNCTION Y = MVNPDF(X, [MU], [SIGMA])
%
% X & MU: vectors of same size
% sigma: square matrix
%% default MU = zeros(size(X))
% SIGMA = diag(ones(size(X)));
function y = mvnpdf(x,mu,sigma)if nargin < 2 mu = zeros(size(x));endif nargin < 3 sigma = diag(ones(size(x)));endx = x(:); % Column vector
mu = mu(:); % Column vector% Check dimensions
n = length(x);if size(mu,1) ~= n | size(mu,2) ~= 1 | size(sigma,1) ~= n | ... size(sigma,2) ~= n error('Parameter dimensions must agree');enda = (2*pi)^-(n/2) * det(sigma)^-.5;b = exp(-.5 * (x-mu).' * inv(sigma) * (x-mu));y = a * b;Just check by your own which Expression (a or b) becomes complex in your application.
Best wishes
Torsten.
Best Answer