离散函数卷积的用RRRRR编程

相关词条：卷积, RRRRR, 函数, 离散

我也想创建词条赚积分
哪位有做过吗？？可以的话程序发我一段。。
比如一个简单的例子，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的分布。
本文来自: 人大经济论坛 详细出处参考：http://www.pinggu.org/bbs/viewthread.php?tid=529746&page=1&from^^uid=386674

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[i]+x2[j]
f[i,j]=f1[i]*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[i]==y)==0)
{y=c(y,vx[i])
g=c(g,vf[i])
}else
{
g=g+vf[i]*(vx[i]==y)
}
}     #卷积第二步
#y
#g
plot(y,g)
sum(g)

[biggrin]

高手啊 领教了！