R：如何产生嘈杂的正弦函数

Question

我对整个R-thing仍然很陌生。

我有以下目标;我有一个正弦函数描述了一段时间内的钙粒子数： 类似于y = a * sin（b * t）+ c

实际上，在随机事件中描述了钙的产生和去除，我想在我的函数中添加一个随机噪声项（最好是平均噪声幅度可扩展）。

像z = y + random * Amplitude

之类的东西 你可以帮帮我吗？

最佳

Answer 1

这是我将要使用的一种方法 - 我提供了两种关于如何生成错误的选项（均匀分布与高斯分布）：

### Equation: y=a*sin(b*t)+c.unif*amp
# variables
n <- 100 # number of data points
t <- seq(0,4*pi,,100)
a <- 3
b <- 2
c.unif <- runif(n)
c.norm <- rnorm(n)
amp <- 2

# generate data and calculate "y"
set.seed(1)
y1 <- a*sin(b*t)+c.unif*amp # uniform error
y2 <- a*sin(b*t)+c.norm*amp # Gaussian/normal error

# plot results
plot(t, y1, t="l", ylim=range(y1,y2)*c(1,1.2))
lines(t, y2, col=2)
legend("top", legend=c("y1", "y2"), col=1:2, lty=1, ncol=2, bty="n")

Answer 2

使用y <- jitter(a*sin(b*t) + c)函数的

jitter()会为您的函数添加随机噪音。你可以指定＆＃34;金额＆＃34; jitter（）内部的参数来控制振幅。