[交流] 【分享】线性卷积c程序 已有3人参与

#include /*  y(n)=x(n)*h(n)      m--length of x(n);      n--length of h(n);      l=m+n-1  length of y(n) */ void conv(float x[],int m,float h[],int n,float y[],int l) {         int i,j,k;         for(i=0; i         {                 y=0.0;                 for(j=0; j                 {   k=i-j;                     if(k>=0 && k                         y += x[j]*h[k];                 }         } } void main() {              //参数m,n,a1,a2,b1,b2,x[],h[]可任意改动         const int m=3;//length of x(n)         const int n=3;//length of h(n)         const int l=m+n-1;//length of y(n)         const int a1=0;//x(n)左下标,表示x(n)={x(0),x(1),x(2)}         const int a2=2;//x(n)右下标         const int b1=4;//h(n)左下标,h(n)={h(4),h(5),h(6)}         const int b2=6;//h(n)右下标         int c1;//卷积结果序列的左下标         int c2;//卷积结果序列的右下标         int i;         float x[m]={1,1,1};         float h[n]={1,1,1};         float y[l];         c1=a1+b1;         c2=a2+b2;         conv(x,m,h,n,y,l);                 printf("result:\n";         for(i=0; i              printf("y(%d)=%f\n",c1+i,y); }

回复此楼