【matlab】三维路面谱生成代码

生成A/B/C/D级三维路面谱
主函数

function main()
clc
clear
close all
%% 参数设置
x=0:1:500;%米%x方向总距离
y=0:1:50;%米%y方向总距离
[X,Y]=meshgrid(x,y);
Gq0=64*10^(-6);%Gq0为Gq(n0)%%%%%%%%%%%%%%%%%%%%在这里修改道路等级
z=q_gen(X,Y,Gq0);%生成三维路面谱
%% 画图
figure(1)
surf(x,y,z);
xlabel('x/m')
ylabel('y/m')
zlabel('q(x,y)/m')

q_gen.m

function z=q_gen(X,Y,Gq0)%生成三维路面谱
% Gq0=256*10^(-6);%
n0=0.1;%参考空间频率
n1=0.011;%频率下限
n2=2.83;%频率上限
m=200;%m为频率范围[n1,n2]划分的m个小区间
W=2;
%%
z=0;%道路振幅
for i=1:m                                               n=n1+(i-0.5)*(n2-n1)/m; % n为每个小区间内中心频率Gq=Gq0*(n/n0)^(-W);%Gq=Gq0*(n/n0)^(-W)  W=2,n0=0.1z=z+sqrt(2*Gq*(n2-n1)/m)*sin(2*pi*n*(X+Y)+2*pi*rand(1));
end

效果