Skip to main content

plotting - How to 2D-plot function in two varibles


I try to make a plot of a function in multi variables as in that paper arXiv:1312.1935, FIG. 2 .


Mathematica graphics


I tried some thing like:


k[s_,f_] = s + f;

ListPlot[Table[k[s, f], {s, -1, 1, 0.5}, {f, -1, 1, 0.5}]]

But it plotted k[s, f] on the y- axis. While I'd like to have s and f on the x and y axises.



There is also ContourPlot, or PlotRegion, but to my knowledge the function k[s,f] will be plotted as continuous regions, while i'd like to present it as points (with known values like in the FIG).


So any help ?



" If the data is 3D and the third entry is obtained by applying a function like k[s,f] to the first two entries (like data set dt3d below), then the function we use to style the data is slightly different:"


Actually I don't understand from here. I understand in the example of td and styleddt that in Style[{##} and k[##], that ## refers to the two variables which k is function of them.


But now I try to plot another function, like Y[s,f,d]= s+f+d;, with -1 < s < 1, f= 0.5, and -0.5 < d < 0.5, and I want to plot Y[s,f,d] only at -2 < Y < 0, or we can use Piecewise as before to know Y values.



Answer



Update: I guess you want to color a list of 2D points using a function like k[s,f]. You can Style the original data and use the resulting data with ListPlot to get something like Figure 2 in the linked paper.


dt = RandomReal[{-1, 1}, {1000, 2}];


The first 10 rows of dt:


TableForm[dt[[;; 10]], TableHeadings -> {None, {"s", "f"}}]

Mathematica graphics


styleddt = Style[{##}, PointSize[.02],
     Piecewise[{{Orange, -1 <= k[##] <= 0}, {Green, .5 <=  k[##] <= .9}, 
      {Red, 1.2 <= k[##] <= 1.6}}, Blue]] & @@@ dt;
labels = {"-1<=k[s,f]<=0", ".5<=k[s,f]<=.9", "1.2<=k[s,f]<=1.6", "otherwise"};
colors = {Orange, Green, Red, Blue};
legend = Row[Style[##, "Panel", 18] & @@@ Transpose[{labels, colors}], Spacer[5]];


ListPlot[styleddt, DataRange -> {{-1, 1}, {-1, 1}}, Frame -> True,
 ImageSize -> 500, AspectRatio -> 1, PlotLabel -> legend]

Mathematica graphics


If the data is 3D and the third entry is obtained by applying a function like k[s,f] to the first two entries (like data set dt3d below), then the function we use to style the data is slightly different:


dt3d = {##, k@##} & @@@ dt;

The first 10 rows of dt3d:


TableForm[dt3d[[;; 10]], TableHeadings -> {None, {"s", "f", "k[s, f]"}}]


Mathematica graphics


styleddata = Style[{#, #2}, PointSize[.02],
     Piecewise[{{Orange, -1 <= #3 <=  0}, 
      {Green, .5 <= #3 <= .9}, {Red, 1.2 <= #3 <= 1.6}},  Blue]] & @@@ dt3d;

We get the same picture as above using:


ListPlot[styleddata, DataRange -> {{-1, 1}, {-1, 1}}, Frame -> True,
         ImageSize -> 500, AspectRatio -> 1, PlotLabel -> legend]


Original post:


Here are few alternative ways to use 2D plots / charts to visualize your data.


Graphics:


Graphics[{PointSize[Large],
  {ColorData[{"TemperatureMap", {-2, 2}}][#3], Disk[{#, #2}, .1], 
     Black, Text[Style[#3, 14, Bold], {#, #2}]} & @@@ (Join @@ 
     Table[{s, f, k[s, f]}, {s, -1, 1, 0.5}, {f, -1, 1, 0.5}])}, 
 Frame -> True]

Mathematica graphics



MatrixPlot


MatrixPlot[Table[k[s, f], {s, -1, 1, 0.5}, {f, -1, 1, 0.5}], 
 DataRange -> {{-1, 1}, {-1, 1}}, ColorFunction -> "TemperatureMap", 
 DataReversed -> True, PlotLegends -> Automatic]

Mathematica graphics


BubbleChart


BubbleChart[Labeled[Style[{#, #2, Abs@#3 /. 0. -> 1}, 
     ColorData[{"TemperatureMap", {-2, 2}}][#3]], #3] & @@@ (Join @@ 
    Table[{s, f, k[s, f]}, {s, -1, 1, 0.5}, {f, -1, 1, 0.5}]), 

 BubbleScale -> (1 &)]

Mathematica graphics


Comments

Post a Comment

Popular posts from this blog

front end - keyboard shortcut to invoke Insert new matrix

I frequently need to type in some matrices, and the menu command Insert > Table/Matrix > New... allows matrices with lines drawn between columns and rows, which is very helpful. I would like to make a keyboard shortcut for it, but cannot find the relevant frontend token command (4209405) for it. Since the FullForm[] and InputForm[] of matrices with lines drawn between rows and columns is the same as those without lines, it's hard to do this via 3rd party system-wide text expanders (e.g. autohotkey or atext on mac). How does one assign a keyboard shortcut for the menu item Insert > Table/Matrix > New... , preferably using only mathematica? Thanks! Answer In the MenuSetup.tr (for linux located in the $InstallationDirectory/SystemFiles/FrontEnd/TextResources/X/ directory), I changed the line MenuItem["&New...", "CreateGridBoxDialog"] to read MenuItem["&New...", "CreateGridBoxDialog", MenuKey["m", Modifiers-...
Post a Comment
Read more

How to thread a list

I have data in format data = {{a1, a2}, {b1, b2}, {c1, c2}, {d1, d2}} Tableform: I want to thread it to : tdata = {{{a1, b1}, {a2, b2}}, {{a1, c1}, {a2, c2}}, {{a1, d1}, {a2, d2}}} Tableform: And I would like to do better then pseudofunction[n_] := Transpose[{data2[[1]], data2[[n]]}]; SetAttributes[pseudofunction, Listable]; Range[2, 4] // pseudofunction Here is my benchmark data, where data3 is normal sample of real data. data3 = Drop[ExcelWorkBook[[Column1 ;; Column4]], None, 1]; data2 = {a #, b #, c #, d #} & /@ Range[1, 10^5]; data = RandomReal[{0, 1}, {10^6, 4}]; Here is my benchmark code kptnw[list_] := Transpose[{Table[First@#, {Length@# - 1}], Rest@#}, {3, 1, 2}] &@list kptnw2[list_] := Transpose[{ConstantArray[First@#, Length@# - 1], Rest@#}, {3, 1, 2}] &@list OleksandrR[list_] := Flatten[Outer[List, List@First[list], Rest[list], 1], {{2}, {1, 4}}] paradox2[list_] := Partition[Riffle[list[[1]], #], 2] & /@ Drop[list, 1] RM[list_] := FoldList[Transpose[{First@li...
Post a Comment
Read more

functions - Get leading series expansion term?

Given a function f[x] , I would like to have a function leadingSeries that returns just the leading term in the series around x=0 . For example: leadingSeries[(1/x + 2)/(4 + 1/x^2 + x)] x and leadingSeries[(1/x + 2 + (1 - 1/x^3)/4)/(4 + x)] -(1/(16 x^3)) Is there such a function in Mathematica? Or maybe one can implement it efficiently? EDIT I finally went with the following implementation, based on Carl Woll 's answer: lds[ex_,x_]:=( (ex/.x->(x+O[x]^2))/.SeriesData[U_,Z_,L_List,Mi_,Ma_,De_]:>SeriesData[U,Z,{L[[1]]},Mi,Mi+1,De]//Quiet//Normal) The advantage is, that this one also properly works with functions whose leading term is a constant: lds[Exp[x],x] 1 Answer Update 1 Updated to eliminate SeriesData and to not return additional terms Perhaps you could use: leadingSeries[expr_, x_] := Normal[expr /. x->(x+O[x]^2) /. a_List :> Take[a, 1]] Then for your examples: leadingSeries[(1/x + 2)/(4 + 1/x^2 + x), x] leadingSeries[Exp[x], x] leadingSeries[(1/x + 2 + (1 - 1/x...
Post a Comment
Read more