帮忙看看代码得到结果是什么？

锦江

1. function [M,m,df]=fftseq(m,ts,df)
   
2. %      [M,m,df]=Fftseq(m,ts,df)
   
3. %      [M,m,df]=Fftseq(m,ts)
   
4. %FFTSEQ      
   
5. fs=1/ts;
   
6. if nargin==2
   
7.    n1=0;
   
8. else
   
9.    n1=fs/df;
   
10. end
    
11. n2=length(m);
    
12. n=2^(max(nextpow2(n1),nextpow2(n2)));
    
13. M=fft(m,n);
    
14. m=[m,zeros(1,n-n2)];
    
15. df=fs/n;

复制代码

1. t0=0.15;
   
2. ts=0.001;
   
3. fc=50;
   
4. fs=1/ts;
   
5. t=[0:ts:t0];
   
6. m=[ones(1,t0/(3*ts)),zeros(1,t0/(3*ts)),ones(1,t0/(3*ts)+1)];
   
7. ft1=m.*cos(2*pi*fc*t);
   
8. c=cos(2*pi*fc*t);
   
9. ft=ft1.*c;
   
10. df=0.25;
    
11. [FT,ft,df1]=fftseq(ft,ts,df);
    
12. f=[0:df1:df1*(length(ft)-1)]-fs/2;
    
13. plot(f,fftshift(angle(FT/fs)));

复制代码

请问上面的代码最后得到的图形是什么意思？

simon21

相位角随频率的变化，相频特性

锦江

simon21 发表于 2015-10-22 09:11
相位角随频率的变化，相频特性

谢谢