问题标签 [piecewise]

我有两条分段线性曲线 c1 和 c2，我想要一条新的分段线性曲线 c3 作为 c1 和 c2 的较小部分。

有没有一种简洁的算法来获得 c3？

他们的观点是：

我正在上一门关于模糊系统的课程，我在电脑上做笔记。这意味着我必须不时在我的电脑上绘制图表。由于这些图的定义非常明确，我觉得用它来绘制它们numpy是一个好主意（我用 LaTeX 做笔记，而且我在 python shell 上很快，所以我想我可以摆脱这个）。

模糊隶属函数的图是高度分段的，例如：

为了绘制这个，我尝试了以下代码numpy.piecewise（这给了我一个神秘的错误）：

在这一点上，我对如何绘制这个函数感到相当困惑。我不太了解错误消息，这进一步阻碍了我调试它的努力。

最终，我希望将此函数绘制并导出到 EPS 文件中，因此我也将不胜感激这些方面的任何帮助。

有谁知道为什么 MATLAB 的“分段三次插值”会在二维曲面图拟合中给出对角化图数据？附件是使用此三次拟合的代码的二维曲面图（在代码和数据下方）。可以看出，图像/图中的峰值在五个孤立的数据点处对角化（从左上角指向右下角）。而不是预期的对称高斯分布。为什么三次插值会这样做？

三次拟合显示在： http: //www.mathworks.co.uk/help/curvefit/fit.html

下面是带有“分段三次插值”的 matlab 代码，用于拟合来自以下的数据：

“helm2Coils126.txt”和“positionsData2_20x20_2.txt”。

对此的任何帮助将不胜感激。

提前谢谢你，

布伦丹·达雷尔

================================CODE================= =====================

====================== 代码结束========================= ====

%====== helm2Coils126.txt ==== 复制并粘贴到 .txt 文件中 ========
%=================== （10 列，40 行）======================

=========================结束第一个文本文件===================== =======

======= positionData2_20x20_2.txt ======== 复制并粘贴到 .txt 文件中 =========

===================（2 列，400 行）======================= ======

===================结束第二个文本文件=========================== ===

当我进行信号处理时，我有一个特殊的分段函数：

如果有任何连续函数可以近似地表达这个分段函数，那么它对于任何 x 是连续可微的？

帮助，

我对 Matlab 中的 Simulink 环境有点陌生。所以，我试图在其中实现一个分段函数，但我不知道该怎么做。我要实现的功能如下：

在这些方程中，唯一随时间变化的是x。那么有人可以帮助我吗？

提前致谢。

编辑：提供有关我的方法的更多信息，我想补充一点，我的“第一个解决方案”是尝试使用Matlab Functionam304 建议的 bloc 来实现它。但是，当我尝试编译它时遇到了问题。所以我移动了一个并开发了以下解决方案，该解决方案有效，但我发现它有点粗糙。

Switch 1：Vref用作阈值，Switch 2：-Vref用作阈值，Switch 3：-Vref用作阈值，Switch 4：Vref用作阈值。

我有一个有界函数。

When x > 1, y = 1;

when 0 <= x <= 1, y = x;

when x < 0, y = 0.

有谁知道任何简单的数学函数来表示 y（不要使用分段函数）。

例如，Heaviside 函数就是一种类似的方法。

谢谢

我需要构建一个具有任意数量的间隔和函数的分段函数，能够对一个 numpy 输入数组进行操作。

我可以使用 for 循环和指标数组来做到这一点，如下面的代码片段所示，但有没有更 Pythonic 的方式来做到这一点？

我尝试使用 numpy.piecewise 但据我所知，段和函数的数量需要在源代码中静态定义。

我正在尝试了解分段混合效果模型的摘要输出，并且可以使用一些见解。具体来说，我想知道如何获得断点左右线的回归截距和斜率。据我了解，下面输出中给出的截距是断点左侧的回归线，而 I(Days * (Days < 6.07)) 给出的值是该线的斜率。但是，我不认为 I(Days * (Days >= 6.07)) 是断点右侧线的斜率，也不是两个斜率的差异。

我从上一个线程中提取了断点：https ://stats.stackexchange.com/questions/19772/estimating-the-break-point-in-a-broken-stick-piecewise-linear-model-with-rando

我试图通过消除随机效应来简化：当 I() 包含在 lm 模型中时，斜率/截距与上面的混合模型非常相似，我仍然感到困惑。

mod_lm<-lm(反应 ~ I(Days*(Days < 6.07)) + I(Days*(Days>= 6.07)), data = sleepstudy) 摘要(mod_lm)

但是，当 I() 从 lm 公式中删除时，我理解输出，并且结果是有意义的。

mod_lm<-lm（反应〜天*（天< 6.07）+天*（天> = 6.07），数据= sleepstudy）摘要（mod_lm）

当 I() 项从 lmer 公式中删除时， lmer 将不会运行。

有人可以告诉我在模型预测器上使用 I() 时如何解释 lmer() 输出，或者告诉我如何在没有模型预测器上的 I() 的情况下运行 lmer() 模型？

我很感激任何可用的指导，因为我无法在 R 帮助页面上找到任何关于这个的指导！

谢谢你。

我在 R 中定义的分段函数有一些问题。

我将函数定义为：

其中 gs1,gs2,gs3 是我之前定义的函数。

为了绘制函数 g1，我尝试了：

但 R 显示以下内容：

我认为问题可能是 R 需要函数的参数作为数字而不是向量？我不确定我是否正确。

我通过使用以下方法找出了绘图问题：

但是我仍然想问有没有办法处理这种问题，因为我发现我定义的分段函数也不适用于许多其他命令。例如，假设我想将 g1 积分超过 0 到 1，如果我只是在 R 中尝试“积分”功能，则会出现相同的警告消息。

由于广播错误，我无法np.piecewise进行多维绘图。

有没有人有办法解决这个问题？

这是我在简化的可执行脚本中的内容：

这给出了错误：