gaussian beam height plot

hi there, is there any matlab code that return the value of a 3d plot? for example: mu = [1 2;-3 -5]; sigma = cat(3,[2 0;0 .5],[1 0;0 1]); p = ones(1,2)/2; obj = gmdistribution(mu,sigma,p);

ezsurf(@(x,y)pdf(obj,[x y]),[-10 10],[-10 10])

so i got this graph so simulate my gaussian beam for nanomachining, the problem is by using this m file, i cant specify the height of the beam. i think the only one way is to return the z value of the graph, and multiply by certain coefficient that reach the height that i want then regraph it. any help would be much appreciated.

thank you

[Oleg Komarov]

Can you format the code?

Answer 1

[Oleg Komarov]

Can you be more specific? What do you want to do with the graph?

Also, pdf map the domain into the [0 1] interval. You can't have more than 100% probability.

If you are just trying to set values for the z axis then:

set(gca,'zlim',[0,1])

or

inspect(gca)

[Win Thomas Halim]

yes, i know i cant get more than 100% probability, thats why i need the values of z from my graph. from there i will be able to edit it at certain height.

here my task is to plot a 3d graph projecting a focused ion beam milling a material. the milling should be looks like gaussian distribution for each point. therefore i find the gmdistribution command come in handy to project the milling, but the problem is i cant modify the height. the result should be like:

http://img8.imageshack.us/i/fibn.jpg

i dont know how to format the code, therefore i just copy and paste it here:

mu = [1 2;-3 -5];

sigma = cat(3,[2 0;0 .5],[1 0;0 1])

p = ones(1,2)/2;

obj = gmdistribution(mu,sigma,p);

ezsurf(@(x,y)pdf(obj,[x y]),[-10 10],[-10 10])

just need to know how to set the height.

thank you very much oleg.

[Oleg Komarov]

To get the values of z call:

h = ezsurf(@(x,y)pdf(obj,[x y]),[-10 10],[-10 10])

get(h, 'zdata')

Answer 2

[Win Thomas Halim]

thank you Mr. Komarov, thats very helpfull!