问题标签 [piecewise]

所以我试图在 MATLAB 中使用 heaviside 和 ezplot 函数绘制 f(t) 分段函数。现在我对 MATLAB 一点也不熟悉。如果有人知道我为什么会收到此错误，那将很有帮助。

f = '12+(-2t+8)*heaviside(t-2)+(2t-12)*heaviside(t-6)'

f = 12+(-2t+8)*heaviside(t-2)+(2t-12)*heaviside(t-6)

ezplot(f)

使用 inlineeval 时出错（第 15 行） 内联表达式中的错误 ==> 12+(-2t+8).*heaviside(t-2)+(2t-12).*heaviside(t-6) 错误：意外的 MATLAB 表达式。

inline/feval 中的错误（第 34 行）INLINE_OUT_ = inlineeval(INLINE_INPUTS_, INLINE_OBJ_.inputExpr, INLINE_OBJ_.expr);

ezplotfeval 中的错误（第 52 行）z = feval(f,x(1));

ezplot>ezplot1 中的错误（第 469 行）[y, f, loopflag] = ezplotfeval(f, x);

ezplot 错误（第 145 行）[hp,cax] = ezplot1(cax, f{1}, vars, labels, args{:});

我一直在研究R非线性模型，比如我们：

我会确定一个伯努利变量是否会影响Alpha1*time和截距。这样的伯努利变量可以是：

有

我发现它nls可以拟合这个模型，但没有这个分类变量的影响，我还发现这个nlme包可以估计变量的影响，而不是非线性模型的项。我的问题是：这些软件包中的哪一个可以帮助我？而且，我怎么能在代码中添加这个伯努利变量？问候和感谢。

如果我有一个像

有没有人在情节中添加标签？

例如，如果我想将位置 x=1 标记为“进入磁场”之类的东西，有什么办法可以做到吗？

我正在尝试对我的数据进行分段拟合，但我需要使用任意数量的线段来进行。有时有三个部分；有时有两个。我将拟合系数存储在 actable 中，并将段的边界存储在 btable 中。

以下是我的界限的示例值：

以下是我的系数的示例值：

这是我的代码的样子：

问题是 lambda 将自己设置为列表的最后一个实例，因此它使用所有段的最后一个段的系数。我不知道如何在不使用 lambda 的情况下执行分段函数。

目前，我这样做是在作弊：

但这让我内心感到恶心。我知道必须有更好的解决方案。有人可以帮忙吗？

我使用 sympy 创建了一个分段函数：

它可以计算0 到 2 之间XX的任何值x（仅感兴趣的范围）。

当我绘制它时，它返回一个错误，即不能将负数提升到分数幂，这永远不会发生，因为只有当x大于 1 时才会应用分数幂。

它是否出于某种原因在整个范围内处理这两个部分？有什么想法吗？

假设您有一个由区间定义的函数，例如

当我们将其与

我们得到

而我想得到

有没有办法获得这样的结果？

我感谢任何帮助segmented.lm（或任何其他函数）在此示例中找到明显断点的帮助：

它返回以下错误：

solve.default(crossprod(x1), crossprod(x1, y1)) 中的错误：系统在计算上是奇异的：倒数条件数 = 1.51417e-20

如果我改变 K：

我得到另一个错误：

segmented.lm(lm(y ~ x, data = data), seg.Z = ~x, psi = NA, control = seg.control(K = 1)) 中的错误：间隔中只有 1 个数据：断点（s ) 在边界或彼此太近

我正在尝试在 R 中使用偏移量的泊松回归（版本 3.1.1）使用分段函数。我的分段包 i 版本 0.4-0.0。

在较早的帖子中，给出了一个可重现的示例：

但是我无法重现该示例 - 我收到一条错误消息：

其他人有同样的问题吗？甚至更好 - 解释？

谢谢

/基拉

我想从分段指数分布中生成一个随机数，它在不同的时间间隔具有不同的危险率。

我的理解是可以应用标准指数分布的无记忆特性。

任何人都知道以下代码是否适用于此目的？特别是，向量“S”包含与相应危险率相关的区间上限（包含在向量“lambda”中。

我希望这个问题很清楚。

我一直在寻找一个解决方案，但仍然没有得到它的工作。我非常感谢每一个有用的建议！

我试图在 R 中定义一个分段线性函数，这意味着，我想定义一个在不同间隔上具有不同线性表达式的函数。这些间隔取决于 xWerte 和 yWerte 中给出的数据点（以及这些数据点的数量），因此切割点的数量和不同的函数规则并没有固定（这就是为什么我用 for 循环而不是用ifelse 构造...)

到目前为止，我一直在尝试定义切割（间隔的切割点）和这些间隔中的每一个的值 = 函数规则。我写下代码以使您了解我所做的事情，尽管在我看来这部分正在工作并且我的问题与我的代码的这一部分无关......

这一切似乎都在做我想做的事。 现在，这是我的问题：

这告诉我

这是 vals[length(vals)]=vals[6]。但由于 searchInt(0.8) = 5，它应该显示 vals[5]=0.63806。

我知道我成功地向 vals 添加了不同的函数规则，因为如果我这样做了

R告诉我

所以，不同区间的正确函数规则是存在的，但它们不能被 fx 正确调用。在我看来，R 似乎“忘记”了 vals 中曾经有不同的函数规则。它只是调用添加到 vals 的最后一个函数规则，不管 searchInt(x) 的值是多少。

我发现这可能是局部和全局变量的问题，但我不知道如何解决它:(

非常感谢您阅读并尝试帮助我！

PS：如有错误请见谅。我不是母语人士，但我已经尽力了；）