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)
扫码加好友，拉您进群
全部版块
我的主页
收藏
