重叠保留法可编辑.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

voidmain()

{

inta[105],b[5],c[104];

inti,j,k,t,r=0,s=0,p;

a[0]=0;

intx[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;

}

intn,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*

Columns1through10