使用matlab，对以上程序中y函数进行加窗校正，数字滤波，时频分析

Richard0824

1. fs=180;
   
2. N=1024;
   
3. n=0:N-1;
   
4. t=n/fs;
   
5. y=0.68*sin(2*pi*10*t+pi/2)+0.45*cos(2*pi*20*t+1/4*pi)+0.28*randn(size(t));
   
6. figure(1)
   
7. subplot(3,1,1)
   
8. plot(t,y)
   
9. xlabel('time/s')
   
10. ylabel('amplitude/1')
    
11. axis([0,1,-1.5,2])
    
12. title('0.68*sin(2*pi*10*t+pi/2)+0.45*cos(2*pi*20*t+1/4*pi)+randn')
    
13. a=fft(y,N);
    
14. b=abs(a/N);
    
15. b1=b(1:N/2+1);
    
16. b1(2:end-1)=2*b1(2:end-1);
    
17. f=(0:(N/2))*fs/N;
    
18. subplot(3,1,2)
    
19. plot(f,b1)
    
20. xlabel('frequency/Hz')
    
21. ylabel('amplitude/1')
    
22. title('A  0.68   0.45    fs  10   20')
    
23. axis([0,50,0,1])
    
24. xlabel('frequency/Hz')
    
25. ylabel('amplitude/1')
    
26. title('corrected')
    
27. axis([0,50,0,1])
    

复制代码

使用matlab，对以上程序中y函数进行加窗校正，数字滤波，时频分析