Wolfram
Mathematica
8의 신기능: 새롭게 향상된 과학 및 정보 시각화 기능
◄
이전
|
다음
►
그래픽 및 시각화
Riemann 합의 시각화
Riemann 합은 곡선 아래의 면적이나 곡면 아래의 부피를 근사하기 위해서 별도의 사각형을 사용합니다.
In[1]:=
X
RiemannBar[box : {{x0_, x1_}, {y0_, y1_}}, {x_, y_}, _] := Block[{area = (y1 + y0) (x1 - x0)}, Sow[area]; {ChartElementData["Rectangle"][box, {}, {}], {Opacity[1], Black, Point[{x, y}]}, If[Abs[area] < 10^-2, {}, {Opacity[1], Black, Text[Rotate[NumberForm[area, {Infinity, 2}], Pi/2], Mean /@ box]}]}]
In[2]:=
X
RiemannPlot[f_, {x_, x0_, x1_, dx_: 1.}, opts___] := Block[{plot, areas, extent, points, curve}, extent = OptionValue[Flatten[{opts}], ExtentSize]; points = If[IntegerQ[dx], RandomReal[{x0, x1}, dx], Switch[extent, Full, Range[x0, x1, dx], Left, Range[x0 + dx, x1, dx], Right, Range[x0, x1 - dx, dx]]]; {plot, areas} = Reap[DiscretePlot[f, Evaluate@{x, points}, opts, ImageSize -> 275, ExtentElementFunction -> RiemannBar, FillingStyle -> Opacity[0.5], PlotStyle -> PointSize[Medium]]]; Show[plot, Plot[f, {x, x0, x1}], PlotLabel -> Row[{"Estimated Area: ", Total[Flatten@areas], Spacer[10], "Actual Area: ", NIntegrate[f, {x, x0, x1}]}], Frame -> True, PlotRange -> All, Axes -> {True, False}] ]
In[3]:=
X
Grid[{{RiemannPlot[BesselJ[5, x], {x, 0, 10, 0.5}, ExtentSize -> Full, PlotStyle -> ColorData["Crayola"]["Asparagus"]], RiemannPlot[BesselJ[5, x], {x, 0, 10, 0.5}, ExtentSize -> Left, PlotStyle -> ColorData["Crayola"]["TealBlue"]]}, {RiemannPlot[ BesselJ[5, x], {x, 0, 10, 0.5}, ExtentSize -> Right, PlotStyle -> ColorData["Crayola"]["RawSienna"]], RiemannPlot[BesselJ[5, x], {x, 0, 10, 20}, ExtentSize -> Full, PlotStyle -> ColorData["Crayola"]["Sunglow"]]}}]
Out[3]=
