楼主: statchao
6582 14

怎么画出这个函数的图像? [推广有奖]

  • 1关注
  • 2粉丝

副教授

80%

还不是VIP/贵宾

-

威望
0
论坛币
1133 个
通用积分
0.0084
学术水平
3 点
热心指数
2 点
信用等级
0 点
经验
13500 点
帖子
337
精华
0
在线时间
1424 小时
注册时间
2005-4-26
最后登录
2023-2-12

楼主
statchao 发表于 2010-12-23 11:14:03 |AI写论文

+2 论坛币
k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群
赵安豆老师微信:zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

感谢您参与论坛问题回答

经管之家送您两个论坛币!

+2 论坛币
  1. x=rnorm(1000)
  2. sigma=3
  3. f <- function(x, sigma) {
  4. if (any(x < 0)) return (0)
  5. stopifnot(sigma > 0)
  6. return((x / sigma^2) * exp(-x^2 / (2*sigma^2)))
  7. }  
复制代码



请教怎么画出这个函数的图像呢?谢谢啊!
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

关键词:RETURN Sigma turn Norm GMA 函数 图像

沙发
qoiqpwqr 发表于 2010-12-23 11:50:19
你的f(x, sigma) = 0,只有一个点,怎么画图啊?

藤椅
statchao 发表于 2010-12-23 12:09:15
这是定义的一个函数,如何画出这个定义的函数的图像。

板凳
qoiqpwqr 发表于 2010-12-23 12:34:01
函数里面的x是一个数还是一组数?

如果x是一个数的话,curve(f(x, sigma = 3), -5, 5)

报纸
statchao 发表于 2010-12-23 14:27:09
不好意思,可能我的表述不清楚。我的意思是,现在定义了一个函数,
f <- function(x, sigma) {
if (any(x < 0)) return (0)
stopifnot(sigma > 0)
return((x / sigma^2) * exp(-x^2 / (2*sigma^2)))
}
希望画出这个函数的图像。谢谢啊!

地板
shenbaiseshatan 在职认证  发表于 2010-12-23 15:36:42
  1. plot(1:20,f(1:20,2),type='l',col=1)
  2. lines(1:20,f(1:20,5),type='l',col=2)
复制代码
胜人者有力,自胜者强!

7
statchao 发表于 2010-12-23 23:09:18
谢谢楼上热心的朋友!
假设x=rnorm(100)该如何做呢?我用下面的程序得到的图像很怪。
plot(order(x),f(order(x),2),type='l',col=1)

8
Epicureanism 发表于 2010-12-24 01:18:26
order返回的是序,去掉order...

9
shenbaiseshatan 在职认证  发表于 2010-12-24 10:26:49
如楼上所述,使用
  1. plot(x,f(x,2),type='l',col=1)
复制代码
此外,为了结果的可重复性,建议设置随机种子set.seed()
胜人者有力,自胜者强!

10
statchao 发表于 2010-12-24 12:55:22
不好意思啊,这个还是不行,如下提示
> set.seed(12345)
> x=rnorm(100)
> plot(x,f(x,2),type='l',col=1)
错误于xy.coords(x, y, xlabel, ylabel, log) :
  'x' and 'y' lengths differ

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注cda
拉您进交流群
GMT+8, 2025-12-30 08:32