首页 文章
  • 1 votes
     answers
     views

    scipy插值griddata极坐标的问题

    我使用scipy griddata来插入极坐标数据 . 我有一个方位角,高度和值的文件 . 首先,我将方位角和高度转换为笛卡尔坐标以与griddata一起使用 . 然后我创建从0到90的高度和0到360的方位角的步骤 . alt = np.array(alt) x,y = alt * np.cos(az), alt*np.sin(az) xi = np.linspace(0, 90, 100) y...
  • 3 votes
     answers
     views

    仿射变换与插值

    我想对一个非常低分辨率的位图进行仿射变换,我想在保留最大量信息的同时进行仿射变换 . 我的输入数据是手写字符的1位64×64像素图像,我的输出将是灰度和更高的分辨率 . 在分析图像时,我构建了一系列仿射变换(旋转,缩放,剪切,平移),我可以将其乘以单个仿射变换矩阵 . 我的问题是,给定 input image 和我的计算 affine transformation matrix ,如何以最高质量计...
  • 1 votes
     answers
     views

    使用牛顿法的非线性插值

    给定一组数据点,我试图使用MATLAB中的牛顿方法逼近函数U(x)= 8-ax ^ b中的系数a,b . x = [150 200 300 500 1000 2000]'; y = [2 3 4 5 6 7]'; a=170; b=-0.7; iter = 0; for iter=1:5 f=8-a*x.^(b) -y; J = [-x.^b -a*b*x.^(b-1)]; %J...
  • 1 votes
     answers
     views

    三次样条插值函数

    我试图使用MATLAB函数 spline 插入以下函数, 在等距点 xi = i./n , i = 0,1,...,n 和 n = 2^j , j = 4,5,...,14 . 对于每次计算,我在 x = 0:0.001:1 点记录最大误差,并使用 loglog 图对 n 绘制这些误差 . 以下是代码, index=1 for j = 4:1:14; n = 2^j; i =...
  • 1 votes
     answers
     views

    从插值数据绘制图形

    我有一个6x7矩阵和代码,如下所示 . AAA = [6.58334819836268e-08,6.58360422229042e-08,1.00000000000000e-07,6.58360422229042e-08,6.58334819836268e-08,6.52045800754157e-08,5.85884775603342e-08;6.58357144702898e-08,7.078...
  • 1 votes
     answers
     views

    拉格朗日插值图

    这是拉格朗日插值的代码 . 我正在调用一个名为 polyinterp 的函数来执行插值 . 我试图近似[-π,π]中的sin(x) . 如果我使用12个点进行插值,我得到的结果很好 . 然而,如果我尝试使用200点进行插值,我无法理解图形,因为绿线应该绘制sin(x)而蓝色圆圈表示使用的点 . x=linspace(-pi,pi,200) target=sin(x); figure('name',...
  • 0 votes
     answers
     views

    Matlab:绘制最优多项式插值的循环

    我试图编写一个循环,通过可能的多项式插值拟合进行排序,然后绘制具有最小误差值s.norm值的插值 . 到目前为止,我有: for i=1:10 [p,S] = polyfit(x,y,i); f = polyval(p,x); subplot(5,2,i) plot(x,y,'*',x,f,'-') axis([0 3 -2.5 12]) ti...
  • 0 votes
     answers
     views

    八度函数的多项式插值

    我需要使用多项式插值来插值函数,但我不确定如何编写部件来获取多项式及其系数: 我的功能: fx = @(x) sin(3 .*x) ./ (0.4 .+ (x .-2) .^2); 我的9个样本值: xs = 0:0.5:4; ys = fx(xs); My attempt: n = 2; x(1) = xs; # not sure if this is correct for k=1:...
  • 1 votes
     answers
     views

    拉格朗日多项式与Matlab

    这是我必须解决的问题: 编写一个程序来评估和绘制拉格朗日插值Iu(x)的u(x)= 1 /(1 x ^ 2),x在-5和5之间 . 为5,7,9,11,13执行此操作, 15点插值(5,7,9等数据点之间,包括-5和5) . 您的结果应显示函数和插值 . 这是我到目前为止提出的代码: int_pts = [5,7,9,11,13,15]; %various values for no. of ...
  • 0 votes
     answers
     views

    Simulink中的逐行插值

    我有一个41x8矩阵,我想在Simulink中用它作为源码 . 源的输出应该是大小为8的行 . 数据需要在行之间进行插值 . 例如: 有一个41x1单调递增的列向量,它与第一个矩阵的41行按行对应,例如: X=[0; Z=[5 20; 2.5; 10 26; 5; 18 33; 7.5; 27 41; ...
  • 1 votes
     answers
     views

    在Matlab中使用非均匀网格进行二维插值?

    在Matlab中从非均匀网格进行二维插值有什么好方法? 我的问题类似于interpolation documentation中"Grid Approximation Techniques"下的问题,但我想知道是否有更好的解决方案,而不是试图使网格有效统一,即是否有任何标准的功能/技术可以让我插入我的非直接均匀采样数据,无需插入我的采样点,以便我可以插入我的数据 . 虽然我的 X...
  • 0 votes
     answers
     views

    多时间序列内插数据的问题

    我正在处理患者数据,其中参数以不同的采样频率记录,因此具有不同的时间戳 . 我想创建一个矩阵,其中数据按“最后已知值”进行插值,直到新的原始值随时间变化 . 所以最后我有统一的矩阵,每个参数在每个时间戳都有值 . 这是我上一个问题的扩展版本, Interpolate the time stamped patient data 数据采用以下格式: Time Hear Rat...
  • 1 votes
     answers
     views

    Matlab插值一维散乱数据

    我有一些数据简单地组织为2D(x,y)坐标 . 我有大量的这些数据,而X轴的数据对于所有点都有相同的范围,它不会在数据集之间使用相同的精确X点 . 我想插入每组数据,然后为每个数据集获取相同的X点 . 每当我使用Matlab尝试插入数据时,我都会遇到问题 . 所以我的数据看起来像: x = [0 1 2 3 4 5 6 7] y = [2.2 3.7 3.9 4.1 4.2 8.9 9.1 9.3...
  • 2 votes
     answers
     views

    用非整数因子重新组合numpy数组

    我有一个模型光谱网格,它具有恒定的,非常高的光谱分辨率,我需要将它们下采样到较低的分辨率,同时保留计数的总数 . 本质上,如果前5个箱子有(标称箱中心)波长 [7.8, 7.81, 7.82, 7.83, 7.84] ,并且值 [1.01, 1.02, 1.015, 1.014, 1.02] ,我想要的箱子是某个(非整数)因子(比如2.5倍)宽,我想要我的新频谱具有标称波长 [7.81, 7.83...
  • 3 votes
     answers
     views

    Matlab重复x轴进行插值

    我必须插入风向 . 每1000 [ft]给出数据 . 例如: %winddata input in feet en degrees x=0:1000:10000; grad=[340 350 360 1 10 20 30 35 34 36 38]; 插值效果很好,我使用interp1函数 . (请参阅下面的代码 . )但是,从360度到1度的步骤是个问题 . 我希望Matlab直接(顺时针)从...
  • 3 votes
     answers
     views

    Julia vs MATLAB:为什么我的Julia代码这么慢?

    我刚开始使用Julia并将我的MATLAB代码翻译成Julia(基本上是逐行) . 我注意到Julia代码要慢得多(比如50x) . 最初的问题是一个动态编程问题,其中我插入了值函数 - 插值是代码大部分时间都在使用的位置 . 所以我尝试制作一个显示性能差异的最小示例代码 . 需要注意的重要事项是,它是插值的样条近似,并且网格最好是不规则的,即不是等间距的 . MATLAB代码: tic spa...
  • 2 votes
     answers
     views

    用于多维线性插值的预计算权重

    我有一个沿D维度的非均匀矩形网格,一个网格上的逻辑值V矩阵,以及一个查询数据点X的矩阵 . 网格点的数量在不同维度上有所不同 . 我对同一网格G和查询X多次运行插值,但对于不同的值V. 目标是预先计算插值的索引和权重并重用它们,因为它们总是相同的 . 这是一个2维的例子,我必须在循环中每次计算索引和值,但我想在循环之前只计算一次 . 我保留了我的应用程序中的数据类型(主要是单个和逻辑gpuArra...
  • 0 votes
     answers
     views

    MATLAB - 导数中带条件的多项式插值

    我有一组值(它们代表放射线的角度变化),我开始创建一个插值所有值的多项式 . time = 10; n = 10*time; t = linspace(0,time,n); %Th1 is the array of angles p1 = polyfit(t,Th1,5); %fifth-grade polynomial 然后我想要一个带边界条件的多项式:在t = 0(第一个值)和t =时间(...
  • 10 votes
     answers
     views

    使用numpy.interp进行线性插值

    我有一个浮点数的一维数组A,其中大部分是好的,但缺少一些值 . 丢失的数据被替换为nan(不是数字) . 我必须通过线性插值从附近的好值中替换数组中的缺失值 . 所以,例如: F7(np.array([10.,20.,nan,40.,50.,nan,30.])) 应该回来 np.array([10.,20.,30.,40.,50.,40.,30.]). 使用Python做到这一点最好的方法是什...
  • 0 votes
     answers
     views

    修改了缺失数据的线性插值

    想象一组具有给定x值(作为列向量)的数据和在矩阵中组合的几个y值(列向量的行向量) . 矩阵中的某些值不可用: %% Create the test data N = 1e2; % Number of x-values x = 2*sort(rand(N, 1))-1; Y = [x.^2, x.^3, x.^4, x.^5, x.^6]; % Example values Y(50:80, 4...
  • -2 votes
     answers
     views

    python:平滑插值

    我知道如何进行线性插值,这是我现在为它构建的:http://repl.it/BBL9 我没看到的是我将如何进行平滑插值 . 我正在编码一个网格变形器,由于入点和终点过渡是突然的,我有一个问题,效果会像图像顶部的线一样被切断 . 我希望它像底部的线一样平滑: 所以我想将浮点数从0.0传递到1.0并返回相应的平滑插值浮点值 . 如果我不能编码它,我正在考虑做我在节点编辑器中所做的事情:拥有2个100...
  • 1 votes
     answers
     views

    线性插值 - 制作网格

    我想在不同的模型之间进行插值 . 为方便起见,我的数据如下所示: 我有10种不同的模拟(我称之为 z ) . 对于每个 z ,我有一个 array x 和一个 array y (对于给定的 z , len(x)=len(y) ) . 例如: for z=1 : x.shape=(1200,) 和 y.shape=(1200,) for z=2 : x.shape=(1250,) 和 y.sh...
  • 4 votes
     answers
     views

    加快插值练习

    我在大约120万次观测中运行了大约45,000个局部线性回归(基本上),所以我很感激一些帮助试图加快速度,因为我很不耐烦 . 我基本上是在为一堆公司构建逐年工资 Contract - 职能工资(给予公司,年份,职位的经验) . 这是我正在使用的数据集(基本结构): > wages firm year position exp salary 1: 0007 19...
  • 3 votes
     answers
     views

    给定三维分散数据,需要在特定点插入数据并查找函数

    这可能是一个很高的顺序,但这是我需要做的......我会得到一些三维的分散数据(x,y,z) . 最终目标是能够为表面上的每个点提供f(x,y)函数 . 例如,给定一个坐标(x,y)包含在数据的凸包内,我希望程序吐出f(x,y)= ax ^ 3 bx ^ 2 cx dy ^ 3 ey ^ 2 fy g ,双三次函数,适合该点的插值数据 . 这导致我总体上探索双三次B样条和样条 . 我一直在spic...
  • 1 votes
     answers
     views

    变换后的图像应始终可见

    我正在尝试使用双线性插值来变换图像,我的输入图像是 I ,我有我的仿射矩阵 [A] ,这将给我变换图像 I' ,根据双线性插值我采用仿射矩阵 inv([A]) 的逆并且应用对于输出图像的每个点(在初始级别都是零),因为我们不能保证输出图像大小可以是任何大小,所以首先我找到了边界,所以我可以得到输出图像的大小, 现在我有输入图像,仿射矩阵和输出图像,它们具有至少可以轻松保存变换图像的大小, But ...
  • 0 votes
     answers
     views

    插值与绘制曲线不对齐(drc包)

    我试图确定蛋白质的量(TIMP)对od(光密度)与TIMP的四参数曲线 . 这里的问题在于,当我将插值绘制在它们从(图像)得到的相同标准曲线的顶部时,它们不对齐 . Misalligned interpolations-curve有任何建议吗? 谢谢 . ##STANDARDS## standards <- structure(list(con = c(0, 0, 0.156, 0.156...
  • 2 votes
     answers
     views

    符号计算插值函数的输入?

    我有一个相当复杂的函数 H(x) ,我正在尝试求解 x 的值 H(x) = constant . 我想用离散间隔生成的插值对象和H(间隔)的相应输出来做到这一点,其中其他输入保持不变 . 我表示插值对象 f . 我的问题是插值对象的调用函数接受array_like,因此将符号传递给 f(x) 以使用sage的求解器方法是不可能的 . 有关如何解决这个问题的任何想法? 我有插值函数 f . 我...
  • 6 votes
     answers
     views

    组内插值

    目标 我想在数据帧中的组内插值 . 这将为数据帧中的每个组提供任意数量的中间点 . 最小的工作示例 我有一个数据框,如: OldDataFrame <- data.frame(ID = c(1,1,1,2,2,2), time = c(1,2,3,1,2,3), Var1 = c(...
  • 0 votes
     answers
     views

    有没有办法在mixin中使用参数在mixin中使用SASS生成变量名? [重复]

    这个问题在这里已有答案: Creating or referencing variables dynamically in Sass 5个答案 关于这个问题,这里有一个我正在尝试做的例子,但我有一种感觉,使用SASS插值在使用参数在mixin中输出变量时不起作用 . $red: #f00; @mixin color-accent-class($color) { .#{...
  • 2 votes
     answers
     views

    为什么scipy.interpolate.griddata在与多处理和OpenBLAS一起使用时会挂起?

    我有3D数据(x,y,z)的随机样本,其中x和y是空间维度,z是它在平面上的位置的函数:z = f(x,y) . 我想使用样本中的插值来在常规网格上评估此函数 . 每组样本都在不同的文件中,需要单独插入(每个文件是不同的时间点) . 为了加快这些文件的处理速度,我想使用多处理(mp)模块使用mp池并行执行它们 . 但是,当我尝试在mp池调用的函数内使用mp.griddata时,进程在进入gridd...

热门问题