[ZX, ZY, ZXX, ZYY, ZXY] = trigradient2(X, Y, Z, T, M) 使用最小二乘线性回归计算函数 Z(X,Y) 的导数。 方程组是用泰勒级数从每个点到相邻顶点建立的。 如果一个顶点连接到少于五个顶点，则也使用两条边距离内的顶点。 这种推导方法比一阶方法提供了更好的结果。 特别是计算出的二阶场导数的误差明显小于用一阶函数推导两倍场的误差。 输入： X= 带有 x 坐标的向量。 Y= 带有 y 坐标的向量。 Z = 矩阵，每个点都有函数值。 如果 Z 有多个列 计算每一列的导数。 可选参数： T = 三角剖分（具有多边形顶点的 Nx3 矩阵）。 如果没有给出 X,Y 的 delaunay 三角剖分，则使用。 M= 用于计算的方法。 默认值为 0。 0：一个大方程组。 快速地。 1：多个小方程系统。 较慢但取决于输入值更准确。
输出： ZX=dz/dx
2021-10-21 15:37:44
6KB
1
三阶样条插值的matlab程序代码
function m=naspline(x,y,dy0,dyn,xx)
%三阶样条插值
%格式：m=naspline(x,y,dy0,dyn,xx) x为结点向量，y为数据……
大学计算方法课作业
2021-10-20 10:54:39
1KB
1
计算贝塞尔函数的一阶导数的零点。 使用以下更改更新 BessDerivZerosBisect.m： 1. 允许 m = 0 2. 允许用户指定所需的特定 m 和 k 值。 3. 允许容差输入参数4. 使用查表为 m 和 k 的小值获得更接近的初始化。 5. 增加错误检查6. 计算导数 wrt 到 x 而不是 wrt m（虽然旧方法也适用。）
2011/05/11 使用了 Vincent 的改进。
2021-10-08 14:55:34
4KB
1
计算相关勒让德多项式的精确导数可能很棘手。 即使在高级文本中，它们通常也被写成递归关系和/或包含阶乘的（归一化）因子。 因此，一个简单的软件实现将很快遇到 IEEE754 双精度的限制，导致 NaN/inf 或显着的精度损失，已经在非常低的 N 度。 LEGENDRE_DERIVATIVE 是导数计算的完全矢量化、数值稳定且经过稳健验证的实现。 它允许快速准确地计算任何 N 次的导数。
它的工作原理与 MATLAB 自己的 LEGENDRE 相同，不同之处在于它不计算多项式值，而是计算导数的值。
2021-09-08 15:43:14
6KB
1