数学分析实验new.pptVIP

偏导数与全微分 多元函数的极值 平面上的几何变换 几个练习 机翼轮廓线的加工 3.4 数控机床的刀具补偿 三次样条插值 三次样条函数 满足以下： (a) 在每一区间段[xi-1,xi]上都是三次多项式，记为Si(x). (b) S(x)在区间[a,b]上有二阶连续导数. (c) 在插值节点处满足插值条件S(xi)=yi, i=0,1,2,…,n. 使用方法 NumericalMath`SplineFit` SplineFit[Data,Cubic] sindata=Table[{t,Sin[t]},{t,0,6Pi,0.3Pi}]; splinesin=SplineFit[sindata,Cubic] ParametricPlot[splinesin[t],{t,0,20},Compiled-False] fff=SplineFit[t,Cubic]; ParametricPlot[fff[t],{t,0,12}, AspectRatio-1/3] 实际问题：刀具补偿的估计 在数控机床加工零件时，由于刀具磨损会影响加工精度，要对刀具的磨损进行补偿. 为测定刀具的磨损速度，实验室每隔一小时测量一次刀具的厚度（单位: mm），得到实验数据如下： 25.3 6 25.7 5 24.8 26.1 26.3 26.5 26.8 27.0 厚度 7 4 3 2 1 0 时间 t={0,1,2,3,4,5,6,7}; y={27.0,26.8,26.5,26.3,26.1,25.7,25.3,24.8}; data=Transpose[{t,y}] ph0=ListPlot[data,PlotStyle- {RGBColor[1,0,0],PointSize[0.02]] * 计算f 关于x1,x2,…的全导数 Dt[f,x1,x2,…] 计算f 关于x的全导数，其中c1,c2,…为常数 Dt[f,x,Constants- {c1,c2,…}] 计算f 关于x的全导数 Dt[f,x] 求f 的全微分 Dt[f[x,y]] 计算f 关于xi的n阶偏导数 D[f,{xi,n}] 计算f 关于x的偏导数，y是x的函数 D[f,x,NonConstants-{y}] 计算f 关于x1,x2,…的高阶偏导数 D[f,x1,x2,…] 计算f 关于x的偏导数 D[f,x] 意义 格式 D[Exp[x+y+z^2],z] D[x^2+y^2,x,y] D[Exp[x+2y],{y,5}] D[x[1]^2+x[2]^2,x[1]] D[g[x^2,y^2],x] D[g[x,y^2],x,x,y] D[x^2+y^2,x,NonConstants-{y}] D[x^2+y[x]^2,x] D[f[Sin[x],y[x^2]],x] Dt[x^2+y^2,x] Dt[x^2+y^2+z^2,x] Dt[x^2+y^2+z^2,x,Constants-{y}] Dt[ArcSin[x/y]] t1={{D[Exp[x+y+z^2],z], D[x^2+y^2,x,y], D[Exp[x+2y],{y,5}]}, D[x[1]^2+x[2]^2,x[1]], D[g[x^2,y^2],x], D[g[x,y^2],x,x,y], D[x^2+y^2,x,NonConstants-{y}], D[x^2+y[x]^2,x], D[f[Sin[x],y[x^2]],x], Dt[x^2+y^2,x], Dt[x^2+y^2+z^2,x], Dt[x^2+y^2+z^2,x,Constants-{y}], Dt[ArcSin[x/y]]}}; t2={{D[Exp[x+y+z^2],z], D[x^2+y^2,x,y], D[Exp[x+2y],{y,5}], D[x[1]^2+x[2]^2,x[1]], D[g[x^2,y^2],x], D[g[x,y^2],x,x,y], D[x^2+y^2,x,NonConstants-{y}], D[x^2+y[x]^2,x], D[f[Sin[x],y[x^2]],x], Dt[x^2+y^2,x], Dt[x^2+y^2+z^2,x], Dt[x^2+y^2+z^2,x,Constants-{y}], Dt[ArcSin[x/y]]}}; TableForm[Transpose[Join[t1,t2]]] 隐函数求导 例 已知x2+y2+z2-4z=0, 求 . eq=(x^2+y^2+z^2-4z==0); y/:Dt[y,x]=0; deq1=Dt[eq,x] deq2=Dt[eq,x,