MATLAB编制路面谱程序

完整代码如下：

％ ==================================================％==x==(0:

0.1:

409.6);％离散路面点 n0= 0.1;

Gq_n0=e-6; delt_n= 0.0024;

thta=2*pi*(rand(1,3/ 0.0024)); q=0; for j=1:3/ 0.0024 n= delt_n *j;

Gq_n二Gq_nO* nO 八2/n^2; Ai=sqrt(2*Gq_n*delt_n); q= q + Ai*sin( 2*pi*n*x+thta(j) ); end

1/ 4

% plot(x,q); %可画出空间路面谱图 ％ hold on;

% ====================================================%=%==求= 功率谱密度

Fs=10;

Hs=spectrum.welch; hpsd=psd(Hs,q,Fs,Fs); Pw = hpsd.Data; Fw = hpsd.Frequencies; h1=log(Fw,Pw,color,blue); hold on;

% ================================================%====

Gq_n0=32e-6;

n=( 0.01: 0.1:10);

Gq_n=Gq_n0* 门0八 2. /n.A2;

log(n,Gq_n,color,red); hold on; Gq_n0=e-6; n=(

2/ 4

0.01: 0.1:10);

Gq_n二Gq_nO* 门0八 2./n.A2;

log(n,Gq_n,color,green); hold on; Gq_n0=128e-6; n=( 0.01: 0.01:10);

Gq_n=Gq_n0*n0A 2./n.A2;

log(n,Gq_n,color,red);

％ =================================================％==从=结

= 果上 可以看出，所模拟的路面波在高频段和理论上的是一致的，而在低频段有很大

误差。不知是在求功率谱时有问题，还是在其他过程中不对，希望大家探讨一 下。谢谢。

2.游客 2010-4-26 23:10:31 Says:

我以前也做过，我的理解误差的原因主要在谐波叠加法和 matlab 的一些函 数的计算上面(尤其是PSD函数)，不过这样的路面谱用来做仿真可以了

3. 游客 2009-1-9 2:19:56 Says:

3/ 4

另外，对以下求PSD寸的函数不是很清楚:

hpsd=psd(Hs,q,Fs,Fs); Pw = hpsd.Data; Fw = hpsd.Frequencies； h1=log(Fw,Pw,color,blue);

为什么Pw和Fw的维数是129x1，有办法设置这两个参数吗？

直接用psd(Hs,q,Fs,Fs也能画出图形来，但结果与 psd(Fw,Pw不同，还不知具体原因 ;

4/ 4

道