设计巴特沃斯模拟低通滤波器求转移函数

VibInfo

设计巴特沃斯模拟低通滤波器，求出转移函数G(p)

1.       subroutine butwcf(l,k,ln,d,c,ierror)
   
2. c----------------------------------------------------------------------
   
3. c    routine BUTWCF: To design low-pass Butterworth analog filter:
   
4. c                    H(s)=D(s)/C(s) ,
   
5. c       If IERROR=0: no errors detected
   
6. c                =1: invalid filter order l and k
   
7. c                                       in Chapter 7
   
8. c----------------------------------------------------------------------
   
9.         dimension d(0:4),c(0:4)
   
10.         pi=4.*atan(1.)
    
11.         ierror=1
    
12.         if(l.le.0.or.k.gt.int((l+1)/2)) return
    
13.         ierror=0
    
14.         d(0)=1.
    
15.         c(0)=1.
    
16.         do 10 i=1,ln
    
17.            d(i)=0.
    
18.            c(i)=0.
    
19. 10       continue
    
20.         orderk=k-(l+1.)/2.
    
21.         if(orderk.eq.0) goto 20
    
22.         c(1)=(-2.)*cos((2*k+l-1)*pi/(2*l))
    
23.         c(2)=1.
    
24.         return
    
25. 20      c(1)=1.
    
26.         return
    
27.         end

复制代码

[ 本帖最后由 VibInfo 于 2006-8-6 07:26 编辑 ]