路面不平度与功率谱密度

nocturne

这个是我生成路面不平度的程序：
clear all
% ]  Y! q7 p& Y* [nn=0.0048
Gxn0=262144
N=5200
ll=0.04
2 `3 w  F3 o' o3 @for ii=1:N/2
3 y: Y! D' Q) Z' E: ~+ Z! k  nk=ii*nn
4 W$ z2 D  e" X6 u8 _0 |Gx(ii)=Gxn0*(nk/0.1)^(-2)
Xk(ii)=sqrt(N*Gx(ii)/2/ll)
8 A0 K# |; k4 C% \4 Nend
%R=2*pi*normrnd(0,1,1,N/2)
2 X; ~* y1 a9 g( p. _5 R$ BR=2*pi*rand(1,N/2)
* P6 K) {: I1 R7 c8 S%x2=sqrt(2)*exp(i*(3*pi/4))
. a( U3 |. K) l0 \# _9 b( a0 Bfor ii=1:N/2
  Xkf(ii)=Xk(ii)*exp(i*R(ii))
   Xkf(N-ii)=Xk(ii)*exp(-i*R(ii))
" }. E) ^# h( C0 R) z3 s# kend
3 g, j8 i! n  Z# V9 J. Q2 `Xkf(N/2)=0
1 b# `  f; m1 [0 U
! t" v6 ]( l% p$ B%tt=Xk(2)*sin(R(2))
for jj=1:N
    ii=(1:1:N-1)
& N5 M- Z1 |( [# W4 {1 x$ w- y    bb=exp(i*2*pi*ii*(jj-1)/N)
    cc=Xkf.*bb
* }8 n/ _9 C2 I+ i    aa=sum(cc)
: u0 C- X. A/ J" Z7 q* ]%aa=0
%for ii=1:N-1
8 ?) w8 i, P% S9 k% aa=aa+Xkf(ii)*exp(i*2*pi*ii*(jj-1)/N)
%end
; C, F+ W& o4 J0 g2 b0 Uxmm(jj)=aa/N
end
%bb
tt=[ll:ll:N*ll]
# \3 R: l9 r% I8 r( ^" Iplot(tt,xmm)
8 M' b* [4 g  Y! C" t+ ?/ w怎么用welch法生成功率谱密度，坐标为双对数坐标。
+ V+ m6 A- l* ?

hfiwioq

感觉你的原始信号的生成有点问题，注意复数计算和for循环中i，j的使用，i，j当作复数单位是可以写成1i，1j，
clear all
clc
n=0.0048;
Gxn0=262144;
* M8 B. H+ E8 m( [- `3 [N=5200;
ll=0.04;
for k=1:N/2
0 m) Z( O; u% B6 I6 Z) [    nk=k*n;
3 ^' U. _6 M4 }+ j: V2 Z" b    Gx(k)=Gxn0*(nk/0.1)^(-2);
    Xk(k)=sqrt(N*Gx(k)/2/ll);
end
. k/ M# E6 D$ d5 MR=2*pi*rand(1,N/2);
( A0 p2 J8 v* i8 p" V3 _for k=1:N/2
( L. X  W7 ~! W4 s8 D. K, p    Xkf(k)=Xk(k)*exp(1i*R(k));
    Xkf(N-k)=Xk(k)*exp(-1i*R(k));
- P/ B6 d& L+ I0 d! a  S% U. [end
Xkf(N/2)=0;
5 c6 Z8 j) I6 a2 afor j=1:N
    k=(1:1:N-1);
    bb=exp(1i*2*pi*k*(j-1)/N);
& p. Z7 I6 c; q' v. K    cc=Xkf.*bb;
0 m8 s3 F* T& R, D4 E7 l  u    aa=sum(cc);
    Xm(j)=abs(aa/N);                    %%此处修改了！%%%
: S# Y- }2 P3 O0 C4 ~end
t=(1:N)*ll;
1 e. w$ h  V+ k$ D4 v4 _. M; Efigure(1)
plot(t,Xm)
L=length(Xm);
nfft=2^nextpow2(round(length(Xm)/4));
/ [1 R" ]: S# I: D8 G9 Gwindow=chebwin(nfft);
overlap=round(length(Xm)/8);
ns=ceil(N/L);
7 u2 g0 v8 K# U- w( U[pxx,f] = pwelch(Xm,window,overlap,nfft,ns,'oneside');
figure(2)
4 W  E3 r) j' t' Z. q; Zplot(log10(f),log10(pxx))

S2_hard

你画的图就默认把虚部去掉了，如果单纯想把虚部去掉，real（xmm）即可，但是不一定合理，是不是取模更合适，abs（xmm）

nocturne

怎么生成这个图形。