Diferencie e integre dados intepolados com quantidades 

Realize uma interpolação tridimensional em uma lista de dados de quantidade.

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]=

Crie o gráfico da solução.

mostre o input completo de Wolfram Languageoculte o input

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]=

As quantidades de entrada e de saída associadas com a função de interpolação podem ser extraídas usando a propriedade "QuantityUnits".

In[4]:=

X if["QuantityUnits"]

Out[4]=

Diferencie a função interpolada.

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 a função interpolada.

In[7]:=

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

Out[7]=