离散函数卷积的用RRRRR编程

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

哪位有做过吗？？可以的话程序发我一段。。

比如一个简单的例子，X的分布列为：1，3，5；0.5，0.3，0.2.
Y的分布列为1，2，3，4，5；0.4，0.2，0.2，0.1，0.1
求X+Y的分布。

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：R编程 函数 离散 卷积 RRRRR

x1=1:100
f1=rep(0.1,100)
x2=seq(1,199,by=2)
f2=(1:100)/10000
    #建立两个函数


n1=length(x1)
n2=length(x2)
x=matrix(0,n1,n2)
f=x
g=0              #赋值

for(i in 1:n1)
{for(j in 1:n2)
{
x[i,j]=x1+x2[j]
f[i,j]=f1*f2[j]
}
}                #卷积第一步

vx=as.vector(x)    #可以用c(x)
vf=as.vector(f)
y=vx[1]
j=1
g=vf[1]     #赋值
n=n1*n2


for (i in 2:n)
{
if (sum(vx==y)==0)
{y=c(y,vx)
g=c(g,vf)
}else
{
g=g+vf*(vx==y)
}
}     #卷积第二步
#y
#g
plot(y,g)
sum(g)