加速度功率谱

加速度功率

跪求输出信号的功率谱matlab计算程序，在simulink中将信号x输出到工作空间后，如何求出该信号的功率普？

NVHht

给你两个例子参考参考(具体可参考《Matlab 7辅助信号处理技术与应用》)：
x(n)的功率谱估计：

1. clear all;
   
2. Fs=1000;%采样频率
   
3. %产生含有噪声的序列
   
4. var=sqrt(1/exp(1.0));
   
5. n=0:1/Fs:1;
   
6. N=length(n);
   
7. e=var*randn(1,N);
   
8. w0=100*pi;
   
9. w1=50*pi;
   
10. xn=exp(j*w0*n-j*pi)+exp(j*w1*n-j*0.7*pi)+e;
    
11. 
    
12. %绘制信号波形
    
13. subplot(311)
    
14. plot(n,abs(xn))
    
15. xlabel('n')
    
16. title('x(n)=exp(j*w0*n-j*pi)+exp(j*w1*n-j*0.7*pi)+e(n)')
    
17. 
    
18. %计算序列的DFT
    
19. nfft=1024;
    
20. Xk=fft(xn,nfft);
    
21. %计算序列的PSD
    
22. Pxx1=abs(Xk).^2/N;
    
23. %绘制功率谱图形
    
24. index=0:round(nfft/2-1);
    
25. k=index*N/nfft;
    
26. plot_Pxx1=10*log10(Pxx1(index+1));
    
27. subplot(312)
    
28. plot(k,plot_Pxx1)
    
29. ylabel('公式直接计算的功率谱')
    
30. 
    
31. %periodogram函数计算的功率谱
    
32. window=boxcar(length(xn));
    
33. [Pxx2,f]=periodogram(xn,window,nfft,Fs);
    
34. plot_Pxx2=10*log10(Pxx2(index+1));
    
35. subplot(313)
    
36. plot(k,plot_Pxx2)
    
37. xlabel('periodogram函数计算的功率谱')
    
38. 
    
39. x(n)的间接法功率谱估计：
    
40. clear all;
    
41. Fs=1000;% 采样频率
    
42. % 产生含有噪声的序列
    
43. var=sqrt(1/exp(1.0));
    
44. n=0:1/Fs:1;
    
45. N=length(n);
    
46. e=var*randn(1,N);
    
47. w0=100*pi;
    
48. w1=50*pi;
    
49. xn=exp(j*w0*n-j*pi)+exp(j*w1*n-j*0.7*pi)+e;
    
50. 
    
51. % 绘制信号波形
    
52. subplot(311)
    
53. plot(n,abs(xn))
    
54. xlabel('n')
    
55. title('x(n)=exp(j*w0*n-j*pi)+exp(j*w1*n-j*0.7*pi)+e(n)')
    
56. 
    
57. % 计算序列的自相关函数
    
58. m=-500:500
    
59. [r,lag]=xcorr(xn,500,'biased')
    
60. subplot(312)
    
61. hndl=stem(m,r);
    
62. set(hndl,'Marker','.')
    
63. set(hndl,'MarkerSize',2);
    
64. ylabel('自相关函数R(m)')
    
65. 
    
66. % 利用间接法计算功率谱
    
67. k=0:1000;
    
68. w=(pi/500)*k;
    
69. M=k/500;
    
70. X=r*(exp(-j*pi/500).^(m'*k));
    
71. magX=abs(X);
    
72. subplot(313)
    
73. plot(M,10*log10(magX));
    
74. xlabel('功率谱的BT法估计')

复制代码

finland

请教各位高手，为什么对于同一个信号，采用不同的方法，其功率谱的幅植相差很大？？？那么功率谱的幅值的物理意义究竟是什么呢

simon21

请教各位高手，为什么对于同一个信号，采用不同的方法，其功率谱的幅植相差很大？？？那么功率谱的幅值的物理意义究竟是什么呢

http://forum.vibunion.com/thread-17307-1-1.html
已作回复