piecewise

0熱度

2回答

誠然，v，a是Eigen::VectorXd矢量具有N維的，我想作以下分段操作：的v的分段乘以a，即，矢量(a[1]*v[1], ..., a[n]*v[n])，和的v的分段方形，即，矢量(v[1]*v[1], ..., v[n]*v[n])。 EigenEigen是否提供上述操作的方法，還是我需要手動實現它們？這當然很簡單，但我希望他們儘可能快地跑。

-1熱度

1回答

分段線性插值R

我知道你可以使用近似函數進行線性插值。但是，我無法找到任何地方如何分段線性內插這對點。例如，假設x < -c（1,4,3,2）y < -c（4,3,2,1）有人可以給我一個答案嗎？提前致謝

1熱度

2回答

分段函數不工作

我定義了一個分段函數是這樣的： def C1(d,x): if ((x<d/2) and (x>-d/2)): return (1/d) else: return 0 但它總是返回我0 In [23]: C1(1000,1) Out[23]: 0

2熱度

1回答

在MATLAB中添加分段多項式

我需要添加從多個數據集派生的分段多項式。有沒有一種簡單的方法來將分段多項式添加到一起而不進行插值？換句話說，給定PP1和PP2，是否有辦法生成PP3（其中PP3保持分段多項式形式）？ e.g ... t1 = linspace(0,1,5); t2 = linspace(0,1,7); pp1 = spline(t1,sin(pi*t1)); pp2 = spline

4熱度

2回答

用一條直線和一條水平線在斷點處進行逐段迴歸

我想要做一個具有一個斷點的分段線性迴歸，其中迴歸線的後半部分有slope = 0。有一些如何進行分段線性迴歸的例子，如here。我遇到的問題是我不清楚如何解決模型的一半的斜率爲0 我試圖 lhs <- function(x) ifelse(x < k, k-x, 0) rhs <- function(x) ifelse(x < k, 0, x-k) fit <- lm(y ~ lhs(x) +

2熱度

2回答

Python/Scipy中的分段恆定或0度樣條插值

儘管scipy的文檔指出scipy.interpolate.UnivariateSpline將在訂單k < = 5上運行，但在引擎蓋下還有一個額外的約束條件，用於order> = 1。有沒有人知道用Numpy/Scipy實現0度樣條或分段常量插值的方法？

1熱度

2回答

循環不會爲分段函數返回正確的值

我試圖編寫代碼，它是較大程序的一部分，它將在s的每個點返回值z。但是，當我運行代碼時，我只能得到z=0，或者如果最後else被忽略，代碼將返回零向量。有人有線索，我犯了錯誤嗎？我用這個source的方法1。任何幫助將不勝感激，我已經嘗試了幾個月現在這個工作。 % clc;close all; %// not generally appreciated %initial values b=1.

1熱度

1回答

的Python：np.vectorize返回「浮動」

運行下面的代碼： import matplotlib.pyplot as plt import numpy as np def xon (ton, t): if ton <= t: return (t-ton)/5 else: return 0 vxon = np.vectorize(xon) t = np.linspace(0, 49, 5

0熱度

1回答

在分段包中使用davies.test需要幫助

對不起，提出這樣一個愚蠢的問題，但我真的很困惑！我有一個斷點的分段線性迴歸，現在我想在包segmented要測試與davies.test斜坡顯著差異：現在哪一個 lin.mod <- lm(y~x) segmented.mod <- segmented(lin.mod, seg.Z=~x) 是戴維斯測試是否正確？ davies.test(lin.mod, seg.Z=~x) 或 dav

2熱度

1回答

Sympy分段表達式偶數和奇數

的目的是實現一個分段表達式，給出0，當n爲偶數時，以及1當n爲奇數。一個辦法做到這一點是使用地板之類的函數如下： from sympy import * from sympy.abc import n f = Lambda((n,), Piecewise((0, Eq(n, floor(n/S(2)))), (1, Eq(n, floor(n/S(2))+1))))