Matlabの当量は

Question

私はこの関数の密度プロットを作成したいdensityplot ：

はしかし、私はMATLABで同様の図をプロットするために使用するかわからないです。ここで

は私の現在のMATLABコードです：

x = [0:10:100]; 
y = [-50:10:50]; 
s = [10, 0]; 
i = [50,25]; 
for 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 

誰かがMATLABでの同等のものを提案してもらえますか？メープルにおける密度プロットのために

---

は、私はあなたがこれを達成するためにsurf（3D表面プロット）を使用することができますが、あなたはそれのための10のステップよりも細かいグリッドが必要になります

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

Answer 1

を使用しました様になる！

また、xとy座標のすべての組み合わせを取得するには、meshgridが必要です。

詳細については、コメントを参照してください。

% Set up grid points 
x = 0:0.1:100; 
y = -50:0.1:50; 
[x,y] = meshgrid(x,y); 
% Set up parameters i, s and g 
i = [50 25]; s = [10 0]; g = 9.8; 
% Work out density 
% - no need for loop if we use element-wise operations ./ and .^ 
% - power(z,2) replaced by z.^2 (same function, more concise) 
% - You forgot the sqare roots in your question's code, included using .^(1/2) 
% - line continuation with "...", could remove and have on one 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, and set to a view from above 
surf(x,y,sir,'edgecolor','none','facecolor','interp'); 
view(2); 
% Change the colour scheme 
colormap('bone') 

結果：

が

あなたの例をマッチングあなたはメープルコマンドscaletorange=-5..50を使用。これにより、スケールが-5と50（docs）に制限されるため、sirがスケール変数であるため、同じにする必要があります。 MATLABでは：今

% Restrict sir to the range [-5,50] 
sir = min(max(sir,-5),50); 
% Of course we now have to replot 
surf(x,y,sir,'edgecolor','none','facecolor','interp'); 
view(2); 

あなたが緑/黒の色を望んでいた場合は、カスタムcolormapを使用することができ、これはまた'bone'colormapのみ64色を持つことによって発生するバンディングを滑らかにします。 g=3.5（あなたが使用わからない）で

% Define the three colours to interpolate between, and n interpolation points 
black = [0 0 0]; green = [0 1 0]; white = [1 1 1]; 
n = 1000; 
% Do colour interpolation, equivalent to Maple's 'colorscheme = [black, "green", "white"]' 
% We need an nx3 matrix of colours (columns R,G,B), which we get using interp1 
colormap(interp1(1:3, [black; green; white], linspace(1,3,n))); 

、我々はほとんど同じプロットを取得