plot - Matlab equivalent of Maple densityplot -


i create density plot of function:

enter image description here

in maple, 1 use densityplot function achieve (code @ end), gives:

density plot

however, not sure use plotting similar figure in matlab.

here current matlab code:

x = [0:10:100]; y = [-50:10:50]; s = [10, 0]; = [50,25]; ii = 1 : length(x)     sir(ii) = -10 * 9.8 * log10((power((x(ii) - s(1)),2) + power((y(ii) - s(2)),2)) / (power((x(ii) - i(1)),2) + power((y(ii) - i(2)),2)));   end

could suggest equivalent in matlab?

for density plot in maple, used 

densityplot(sir(x,y), x=0..100, y=-50..50, axes=boxed, style=patchnogrid, scaletorange=-5..50, colorscheme = [black, "green", "white"])

you can use surf (a 3d surface plot) achieve this, need finer grid steps of 10 good!

also need meshgrid combinations of x , y coordinates.

please see comments further details.

% set grid points x = 0:0.1:100; y = -50:0.1:50; [x,y] = meshgrid(x,y); % set parameters i, s , g = [50 25]; s = [10 0]; g = 9.8; % work out density   % - no need loop if use element-wise operations ./ , .^ % - power(z,2) replaced z.^2 (same function, more concise) % - forgot sqare roots in question's code, included using .^(1/2) % - line continuation "...", remove , have on 1 line sir = -10*g*log10( ((x-s(1)).^2 + (y-s(2)).^2).^(1/2) ./ ...                    ((x-i(1)).^2 + (y-i(2)).^2).^(1/2)        );  % plot, , set view above surf(x,y,sir,'edgecolor','none','facecolor','interp'); view(2); % change colour scheme colormap('bone')

result:

plot

matching example

you used maple command scaletorange=-5..50. limits scale between -5 , 50 (docs), since sir our scale variable, should limit same. in matlab:

% restrict sir range [-5,50] sir = min(max(sir,-5),50); % of course have replot surf(x,y,sir,'edgecolor','none','facecolor','interp'); view(2);

now, if wanted black/green colours, can use custom colormap, smooth out banding caused 'bone' colormap having 64 colours.

% define 3 colours interpolate between, , n interpolation points black = [0 0 0]; green = [0 1 0]; white = [1 1 1];  n = 1000;  % colour interpolation, equivalent maple's 'colorscheme = [black, "green", "white"]' % need nx3 matrix of colours (columns r,g,b), using interp1 colormap(interp1(1:3, [black; green; white], linspace(1,3,n)));

with g=3.5 (not sure used), identical plot

compare


Comments

Post a Comment

Popular posts from this blog

Is there a better way to structure post methods in Class Based Views -

i'm using class based views in django , never specify url in ajax calls , omit action parameter in form because know these requests go through post method of class based view. notice view - profile page make many requests depending on user updating on page - post method become convoluted , if/else mess example: def post(self, request, *args, **kwargs): # handle user changing available date if request.post.get('availabledate') != none: self.updateavailabledate(request) return jsonresponse({'result':'success'}) if request.post.get('newprojectname') != none: creative_user = creativeuserprofile.objects.get(id = request.user.id) project = project.create_project(creative_user, request.post['newprojectname']) project.save() return jsonresponse({'projectid': project.id }) if request.post.get('projectimage') != none , request.files['file']: proj...
Read more

performance - Why is XCHG reg, reg a 3 micro-op instruction on modern Intel architectures? -

i'm doing micro-optimization on performance critical part of code , came across sequence of instructions (in at&t syntax): add %rax, %rbx mov %rdx, %rax mov %rbx, %rdx i thought had use case xchg allow me shave instruction , write: add %rbx, %rax xchg %rax, %rdx however, dimay found agner fog's instruction tables , xchg 3 micro-op instruction 2 cycle latency on sandy bridge, ivy bridge, broadwell, haswell , skylake. 3 whole micro-ops , 2 cycles of latency! 3 micro-ops throws off 4-1-1-1 cadence , 2 cycle latency makes worse original in best case since last 2 instructions in original might execute in parallel. now... cpu might breaking instruction micro-ops equivalent to: mov %rax, %tmp mov %rdx, %rax mov %tmp, %rdx where tmp anonymous internal register , suppose last 2 micro-ops run in parallel latency 2 cycles. given register renaming occurs on these micro-architectures, though, doesn't make sense me done way. why wouldn't register renamer ...
Read more

c# - Asp.net web api : redirect unauthorized requst to forbidden page -

im trying redirect unauthorized request forbidden page instead i'm getting forbidden page in response body , how can fix ? here's startup class : app.createperowincontext(storecontext.create); app.createperowincontext<applicationusermanager>(applicationusermanager.create); app.createperowincontext<applicationsigninmanager>(applicationsigninmanager.create); app.usecookieauthentication(new cookieauthenticationoptions { authenticationtype = defaultauthenticationtypes.applicationcookie, expiretimespan = timespan.fromdays(30), }); app.useexternalsignincookie(defaultauthenticationtypes.externalcookie); app.usetwofactorsignincookie(defaultauthenticationtypes.twofactorcookie, timespan.fromminutes(5)); app.usetwofactorrememberbrowsercookie(defaultauthenticationtypes.twofactorrememberbrowsercookie); the method im trying reach : [httpget] [authorize(roles = "admin")] public string getcurrentusername() { return userman...
Read more