高斯型随机粗糙面MATLAB仿真

一、概要
一般可认为随机粗糙表面高度可由许多不同频率的谐波叠加而成，因此可采用线性滤波法来生成随机粗糙面：先在频域谱进行滤波，然后进行逆傅里叶变换得到粗糙面高度，来生成模拟随机粗糙面模型。本文将基于上述原理，给出仿真步骤以及仿真结果。
二、仿真过程

clear
clc
%% 初始参数设置
omega=0.5*10^-6;%高度起伏均方根为0.1064um
kxi=4.75*10^-6;%横向相干长度为4.75um

M=180;%x方向采样点数
N=180;%y方向采样点数

Lx=M*kxi;%二维随机粗糙表面x方向长度
Ly=N*kxi;%二维随机粗糙表面y方向长度

ux=-M/2:M/2;%X方向
uy=-N/2:N/2;%Y方向
Ux=2*pi*ux/Lx;%X方向离散波数
Uy=2*pi*uy/Ly;%Y方向离散波数

vx=-M/2:M/2;%X方向
vy=-N/2:N/2;%Y方向
[vx,vy]=meshgrid(vx,vy);
deltax