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]

高手啊 领教了！