如何在R语言中为抛物线（函数）添加随机噪声并绘图？

问题原因 #

你当前代码的问题在于rnorm(1, 100,100)只生成了单个随机数，把这个数加到整个y向量上后，所有点都偏移了相同的量，本质上还是原抛物线的平移，没有实现每个数据点独立的噪声添加。

修正后的代码 #

要给每个点都添加独立的随机噪声，需要生成和x长度一致的随机数序列：

x <- 1:100
y <- x^2
# 生成100个均值为0、标准差为500的正态分布随机噪声（参数可按需调整）
z <- y + rnorm(length(x), mean = 0, sd = 500)

plot(x, z, pch = 16, col = "#2E86AB", main = "带正态噪声的抛物线")
# 可选：叠加原抛物线作为对比
lines(x, y, lwd = 2, col = "#F24C4C")

关键说明 #

• rnorm(n, mean, sd)中第一个参数n必须等于数据点的数量（这里是100），这样每个y值都会加上一个独立的随机数。
• 噪声的均值和标准差可以根据需求调整：
  • 均值设为0时，噪声会在原函数值上下波动；若设为非零值，会同时引入偏移和波动。
  • 标准差越大，噪声波动越剧烈；越小则越接近原函数。

扩展到任意函数 #

对于任意函数（比如线性函数、正弦函数等），方法完全一致：先生成基础函数值，再添加对应长度的随机噪声。示例：

# 带噪声的正弦函数
x <- seq(0, 2*pi, length.out = 100)
y <- sin(x)
z <- y + rnorm(length(x), mean = 0, sd = 0.1)

plot(x, z, pch = 16, col = "#2E86AB", main = "带噪声的正弦函数")
lines(x, y, lwd = 2, col = "#F24C4C")

内容的提问来源于stack exchange，提问作者stats_noob