Diferencie e integre datos interpolados con cantidades 

Lleve a cabo una interpolación tridimensional en una lista de datos de cantidad.

In[1]:=

X data = N@Flatten[ Outer[Quantity[Sin[#1 + #2^2], "Pounds"] &, Range[10], Range[10]], 1]; grid = Flatten[ Outer[List, Quantity[Range[10], "Seconds"], Quantity[Range[10], "Feet"]], 1]; data = Join[grid, Transpose[{data}], 2] /. {x_, y_, z_} :> {{x, y}, z};

In[2]:=

X if = Interpolation[data]

Out[2]=

Grafique la solución.

muestre la entrada completa de Wolfram Languageoculte la entrada

In[3]:=

X fun[t_?NumericQ, p_?NumericQ] := QuantityMagnitude[ Quiet[if[Quantity[t, "Seconds"], Quantity[p, "Feet"]] ] ]; Plot3D[fun[tt, p], {tt, 1, 10}, {p, 1, 10}, AxesLabel -> {"Seconds", "Feet", "Pounds"}, PlotRange -> All, ImageSize -> Medium]

Out[3]=

Las cantidades de entrada y salida asociadas con la función de interpolación pueden ser extraídas utilizando la propiedad "QuantityUnits".

In[4]:=

X if["QuantityUnits"]

Out[4]=

Diferencie la función de interpolación.

In[5]:=

X Derivative[#, 0][if][Quantity[3, "Seconds"], Quantity[1, "Meters"]] & /@ Range[3]

Out[5]=

In[6]:=

X Derivative[#, #][if][Quantity[3, "Seconds"], Quantity[1, "Meters"]] & /@ Range[3]

Out[6]=

Integre la función de interpolación.

In[7]:=

X Integrate[ if[x, y], {x, Quantity[1, "Seconds"], Quantity[10, "Seconds"]}, {y, Quantity[1, "Feet"], Quantity[10, "Feet"]}]

Out[7]=