在归一化频率轴上绘出归一化的功率谱曲线

VibInfo

1.       subroutine psplot(psdr,psdi,mfre,ts)
   
2. c-----------------------------------------------------------------------
   
3. c   Routine PSPLOT: To plot the normalized power spectum curve on the
   
4. c                      normalized frequency axis from -.5 to  +.5 .
   
5. c        mfre : Points in frequency axis and must be the power of 2.
   
6. c        ts    : Sampling interval in seconds (real).
   
7. c        psdr : Real array of power spectral density values.
   
8. c        psdi : Real work array.
   
9. c                                            in Chapter 11
   
10. c-----------------------------------------------------------------------
    
11.         dimension psdr(0:mfre-1),psdi(0:mfre-1)
    
12.         m2=mfre/2
    
13.         do 40 k=0,m2-1
    
14.            psdi(k)=psdr(k)
    
15.            psdr(k)=psdr(k+m2)
    
16.            psdr(k+m2)=psdi(k)
    
17. 40      continue
    
18.         pmax=psdr(0)
    
19.         do 50 k=1,mfre-1
    
20.            if(psdr(k).gt.pmax) pmax=psdr(k)
    
21. 50      continue
    
22.         do 60 k=0,mfre-1
    
23.            psdr(k)=psdr(k)/pmax
    
24.            if(psdr(k).le.0.0)psdr(k)=.000001
    
25. 60      continue
    
26.         fs=1./ts
    
27.         fs=fs/float(mfre)
    
28.         open(3,file='psd.dat',status='new')
    
29.         do 100 k=0,mfre-1
    
30.            faxis=fs*(k-m2)
    
31.            write(3,*)faxis,10.*alog10(psdr(k))
    
32. 100      continue
    
33.         close(3)
    
34.         return
    
35.         end

复制代码

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