-
1 votesanswersviews
在3D表面上使用分析Perlin噪声导数
我已经在试验Catlike Coding's tutorials中的一个代码了，虽然我想到假设'up'轴是1并且在2D平面上采样时进行标准化，我找不到任何有关在3D表面上采样时应该做什么的信息 . 这就是我所追求的（这使用中央差异法线）： 这就是那些法线在世界空间中的样子： 但我得到的衍生物看起来像这样： 看起来他们可能是相对于未移位球体的表面，所以这是否意味着我需要用未移位的球体法线重新定向它们... -
16 votesanswersviews
在MATLAB中计算数值导数的最佳方法是什么？
（注意：这是一个社区Wiki . ） Suppose I have a set of points xi = {x0,x1,x2,...xn} and corresponding function values fi = f(xi) = {f0,f1,f2,...,fn}, where f(x) is, in general, an unknown function. （在某些情况下，我们可能提前... -
10 votesanswersviews
3D Perlin噪声分析导数
我目前正在使用Shader Model 4（DirectX 10 HLSL）实现3D Perlin噪声凹凸贴图 . 生成噪声本身并不是一个大问题（有大量的教程和代码），但我没有找到3D Perlin噪声的分析衍生物 . 考虑衍生品的唯一网站是Ińigo Quilez's site和相关的GameDev.net discussion . 问题是在第一个链接中噪声是基于值的，而不是基于梯度的（这是我的... -
2 votesanswersviews
Perlin Noise：为什么二阶插值导数会产生正常/阴影伪影？
我正在研究/实现Perlin Noise和Improved Perlin Noise的版本 . Perlin在他的论文中说他取代了smoothstep函数 3t^2 - 2t^3 他曾经使用以下功能在网格单元的角上插入8个线性函数： 6t^5 - 15t^4 + 10t^3 因为smoothstep函数的二阶导数是不连续的 . 他说（并且在他所显示的图像中清晰可见），这会导致一些视觉假象，因... -
1 votesanswersviews
3D perlin噪声的分析导数产生线伪影
当我试图在噪声计算中利用分析导数时，我会得到线条伪影 . 我是_789976的方法：3D Perlin noise analytical derivative 例如，在尝试使用IQ噪声时： float IQturbulence(float3 p, int octaves, float freq, float amp, float gain, float lacunarity) { floa... -
0 votesanswersviews
多项式曲面导数的置信区间
我的问题如下： 我已经将曲面拟合到某些xyz坐标数据以获得多项式曲面 . 这是变量x和y中的多项式，给出了z值的表面 . 我知道如何在平均响应和新的预测间隔（标准的东西）上计算95％的置信区间 . 获得表面函数后，也可以很容易地确定表面上某个点的导数，但我正在寻找的是： How do I compute confidence levels on the derivative value dS/dx... -
1 votesanswersviews
如何在Simulink中使用高阶限速器平滑矩形信号？
想象一下，我有一个位置/位移 x 的矩形参考值，我需要平滑它 . 翻译运动的数学很简单： speed: v = x' acceleration: a = v' = x'' jerk. j = a' = v'' = x''' 我需要限制所有这些值 . 所以我考虑在Simulink中使用速率限制器：这种方法适用于斜坡信号，如下面的输出所示： 但是， x 的参考信号没有斜坡，它们是矩形/台阶 . 因此... -
0 votesanswersviews
求解Maple上求和的导数
有没有办法解决枫木上的总和衍生物 . 例如， y := qr->(1-qn-(sum(qn[j]-delta*qr, j))-c[r])*qr+(wn-c)*(qn+sum(qn[j], j)) 我拿了第一个偏导数w.r.t qr，我认为结果是正确的 diff(y(qr), qr) 但是，当我解决衍生物的结果时，枫树没有给我任何结果 . solve(diff(y(qr), qr) = 0... -
1 votesanswersviews
在pyomo中如何从目标函数中提取二阶导数
我正在使用pyomo，并且已经定义了一个模型，具有与之相关的目标函数 . 模型解决后，目标函数附加了某些参数 . 所以如果我有一个多索引变量 [x1, x2, x3] ，我的二次目标函数会看起来像这样： (x1^2 + 13*x2^2 + 10*x3^2) + (2*x1 +......) . 我的问题是：鉴于我实际上可以从目标中以字符串格式访问此表达式，有没有办法获得关于所有变量的此函数的二阶... -
1 votesanswersviews
在方案中绘制衍生物
我写了一个衍生函数如下： (define (der f h) (lambda (x) (/ (- (f(+ x h)) (f x)) h))) 以及寻找第二和第三衍生物的另一个功能： (define (der-multiple f h n) (if (= 0 n) 0 (der (f h))) (- n 1)) 我试图将sin（x）的导数绘制如下，但它不起作用 . ... -
1 votesanswersviews
如何应用Newton-Raphson方法来查找五次函数的根
说明 我开发了一种实现Newton-Raphson方法的算法来找到一个五次函数的根 . 我必须反映的结果是 303.6 . 但是，我的实施无法衡量 . 数据 Parameters g = 9.81;Ds = 0.198;uj = 805.9;W = 0.0557; Equation 0.024 *（（gDs / uj ^ 2）^（1/3））（Y ^（5/3））0.2 *（Y ^（2/3））... -
0 votesanswersviews
给定点处的非线性pyomo约束的梯度
I（重复地）需要在给定点处的非线性pyomo约束 con 的数值梯度信息（即，相应的pyomo模型的变量都被设置为特定值） . 我已阅读this post并决定（稍加修改）行 from pyomo.core.base.symbolic import differentiate var_list = list(model.component_objects(Var, active=True)) g... -
3 votesanswersviews
模拟终止是因为导数不是有限的
我尝试运行我的Simulink应用程序时遇到此错误 . 时间9.6046876340724416E-7中块'proiect / Filtru / Integrator'中状态'1'的导数不是有限的 . 模拟将停止 . 解决方案中可能存在单一性 . 如果没有，请尝试减小步长（通过减小固定步长或通过拧紧误差容差） 我的Simulink应用程序是这样的： 我得到错误的块看起来像这样： 在"... -
3 votesanswersviews
如何使用Sympy导出Matrix元素
给定矩阵和向量的乘积 A.v 对于形状（m，n）和v的dim n，其中m和n是符号，我需要相对于矩阵元素计算导数 . 我还没有找到使用正确向量的方法，所以我从2 MatrixSymbol 开始： n, m = symbols('n m') j = tensor.Idx('j') i = tensor.Idx('i') l = tensor.Idx('l') h = tensor.Idx('h')... -
2 votesanswersviews
在matlab中绘制差异误差
我写了一个matlab代码，它使用中心差异来近似 cos(x) 的导数 . 意思 cos(x)'=(cos(x+h)-cos(x-h))/2h 大约 . 因为h从0.1,0.01,0.001开始，依此类推，直到10 ^ -8 这是我的代码 function [pos,neg,D]=shifted_cos(x) for i=1:8 var=rand(1)*0.5*10^(-5... -
0 votesanswersviews
采用Tensorflow中的渐变，tf.gradient
我正在使用tensorflow的这个函数来获取我的函数jacobian . 遇到两个问题： 如果我没有弄错的话，张量流文档在以下两段中与自身相矛盾： gradient（）将ops添加到图形中以输出ys相对于xs的偏导数 . 它返回长度为len（xs）的Tensor列表，其中每个张量是y中y的和（dy / dx） . Blockquote Blockquote返回：xs中每个x的sum（dy...