通过功率谱密度函数还原出时域函数（公路谱）

Raspberry

频域法的核心是快速傅立叶变换，即功率谱密度函数在离散的采样点上与信号的频谱有着一个确定的关系。如果能够在功率谱密度函数上离散采样，构造出频谱 ，然后再对其进行傅立叶逆变换，即可得到时域的函数曲线。
下面以公路谱为例的matlab处理程序：

1. Gx0=256;                  %参考空间频率n0下的路面功率谱密度
   
2. n0=0.1;                   %参考空间频率n0
   
3. L=409.6;
   
4. l=0.1;
   
5. N=L/l;                     %采样点数
   
6. n1=0.01;                   %空间频率范围n1--nu
   
7. nu=3;
   
8. w=2;                       %频率指数
   
9. no=1/L;                    %空间频率间隔
   
10. Xk=[];
    
11. Xm=[];
    
12. n=linspace(0.01,3,N/2+1);
    
13. GxC=Gx0*(n/n0).^(-w);
    
14. k=0:N/2;
    
15. fik=randn(1,N)*2*pi;             %产生0到2pi的均匀分布的随机序列
    
16. pg=GxC(1:N/2+1);
    
17. Xk=sqrt((N/2+1)/(2*l)*pg).*exp(j*fik(1:N/2+1));  %调用函数GxC（n）
    
18. Xk(1)=sqrt((N/2+1)/(2*l)*pg(1)); Xk(N/2+1)=sqrt((N/2+1)/(2*l)*pg(N/2+1));
    
19. Xk=[Xk(1:2049) conj(Xk(2048:-1:2))];
    
20. Xm=ifft(Xk);            %逆傅立叶变换后得到复数形式随机序列
    
21. x=linspace(0,409.6,length(Xm));
    
22. subplot(211);
    
23. plot(x,real(Xm));                                %取实部
    
24. xlabel('行驶距离/m');
    
25. ylabel('路面不平度/mm');
    
26. Pxr=abs(fft(real(Xm))).^2/N;                     %恢复序列的功率谱
    
27. Pxr=Pxr(1:N/2+1);
    
28. subplot(212);
    
29. n=linspace(0.01,3,N/2+1);
    
30. loglog(n,GxC(1:N/2+1),'r'); hold on;
    
31. loglog(n,Pxr);                  %恢复序列的功率谱与原功率谱值比较
    
32. xlabel('空间频率n');
    
33. ylabel('功率谱密度Gx(n)'); hold off;
    
34. legend('GxC','Pxr')

复制代码

转自：http://blog.sciencenet.cn/blog-357150-279391.html