直接按定义求序列x(n)，h(n) 的线性卷积y(n)

VibInfo

直接按定义求序列x(n)，h(n) 的线性卷积y(n)

主程序：

1. C----------------------------------------------------------------------
   
2. C main program HCONVO1:To test subroutine CONVO1
   
3. c Please link subroutine CONVO1
   
4. C----------------------------------------------------------------------
   
5.         dimension h(0:7),x(0:7),y(0:7)
   
6.         data m/4/,n/4/
   
7.         l=n+m-1
   
8. C
   
9.         do 10 i=0,3
   
10.            h(i)=1.
    
11.            x(i)=i+1.
    
12. 10      continue
    
13.         call convo1(x,h,y,n,m,l,Ierror)
    
14.         write(*,*)'   Ierror=',Ierror
    
15.         do 2 k=0,l-1
    
16. 2          write(*,*)k,y(k)
    
17.         stop
    
18.         end

复制代码

子程序：

1.         subroutine convo1(x,h,y,n,m,L,Ierror)
   
2. c--------------------------------------------------------------------
   
3. c Routine CONVO1: To implement Linear Convolution y(n)=x(n)*h(n)
   
4. c input parameters:
   
5. c x(n):L dimensioned real array,signal data is stored in x(0) to x(n-1).
   
6. c h(n):L dimensioned real array,impulse response is stored in h(0) to h(m-1).
   
7. c n   : the data length of x.
   
8. c m   : the data length of h.
   
9. c L   : the data length of y, L must be >=n+m-1
   
10. c output parameters:
    
11. c y(n):L dimensioned real array, y(n)=x(n)*h(n),n=0,...L-1.
    
12. c                                      in Chapter 1
    
13. c--------------------------------------------------------------------
    
14.         dimension x(0:L-1),h(0:L-1),y(0:L-1)
    
15.         Ierror=0
    
16.         if(l.ge.(m+n-1))goto 1
    
17.         Ierror=1
    
18.         return
    
19. 1       do 2 i=n,L-1
    
20. 2          x(i)=0.0
    
21.         do 3 i=m,L-1
    
22. 3          h(i)=0.0
    
23. c
    
24.         do 11 k=0,L-1
    
25.            sum=0.0
    
26.            do 12 i=0,k
    
27.               sum=sum+x(i)*h(k-i)
    
28. 12          continue
    
29.            y(k)=sum
    
30. 11       continue
    
31.         return
    
32.         end

复制代码