重叠保留法.docVIP

1、已知x(n)=n，0≤n≤99，h(n)=n+1, 0≤n≤4,使用C语言编译出重叠相加法，再利用Matlab编写重叠相加法和重叠保留法程序，计算x(n)和h(n)的卷积。 (1) C语言程序及结果截图 #includestdio.h void main () { int a[105], b[5],c[104]; int i,j,k,t,r=0,s=0,p; a[0]=0; int x[10][20],xnew[5][10][20],xnewl[10][18]; for(i=0;i104;i++)//给数组赋值，生成序列x(n) { if(i4) a[i]=0; else a[i]=i-4; } for(j=0;j5;j++)//给数组赋值，生成序列h(n) { b[j]=j+1; } int n,m; for(i=0;i10;i++)//把序列x(n)分成10组 { if(i0) { n=10*i; m=10*(i+1)+3; } else { n=0; m=13; } for(j=n,k=0;j=m;j++,k++) { x[i][k]=a[j]; } } for(k=0;k5;k++)//个小段分别与序列b[n]相乘 { for(i=0;i10;i++) { for(j=0;j14;j++) { xnew[k][i][j]=x[i][j]*b[k]; } } } for(i=0;i5;i++)//把xnew乘b[n]所有段数组元素向后移动n位 { for(j=0;j10;j++) { for(k=13;k=0;k--) { xnew[i][j][k+i]=xnew[i][j][k]; } } } for(i=0;i5;i++)//移动n位后，在数组前面补n个零 { for(j=0;j10;j++) { for(k=0;ki;k++) { xnew[i][j][k]=0; } } } for(i=0;i5;i++)//各个数组元素最后一位补4-n个0 { for(j=0;j10;j++) { for(k=17,p=0;p(4-i);k--, p++) { xnew[i][j][k]=0; } } } for(i=0;i10;i++) { for(j=0;j18;j++) { xnewl[i][j]=0; } } for(i=0;i10;i++)//将各小段数组的对应各元素相加 { for(j=0;j18;j++) { for(k=0;k5;k++) { xnewl[i][j]+=xnew[k][i][j]; } } } for(i=0;i104;i++) c[i]=0; for(i=0;i10;i++)//去除重叠部分，合成最终结果 { n=i*10; m=10*(i+1)-1; for(j=n,k=4;jm,k14;j++,k++) c[j]=xnewl[i][k]; } for(i=100,j=14;i105,j18;i++,j++) { c[i]=xnewl[9][j]; } for(k=0;k104;k++)//输出最终结果 { r++; if((r-1)%10==0) { printf(\n); s++; printf(第 %d 行 ,s); } printf(%d ,c[k]); } scanf(输入数据); } (2) Matlab运行结果 在命令框输入以下命令，调取已编好的三个程序conv(hn, xn)、 overlap_add(xn, hn, N) overlap_save(xn, hn, N)，得出一下结果： n = 0:9; xn = n+1; hn = [1,0,-1]; N = 6; y1 = conv(hn, xn) y2 = overlap_add(xn, hn, N) y3 = overlap_save(xn, hn, N) ① 重叠相加法计算结果 y2 = 1.0e+003 * Colum