二维粗糙面的模拟的matlab代码

以下是一个使用Matlab模拟二维粗糙面的代码示例：

% 设置参数
N = 256; % 图像大小
p = 3; % 噪声维度
nOctaves = 6; % Octaves数量
persistence = 0.5; % Persistence参数

% 生成Perlin噪声
noise = zeros(N);
for octave = 1:nOctaves
    freq = 2^(octave-1);
    amp = persistence^(octave-1);
    [X,Y] = meshgrid(1:N);
    noise = noise + amp*interp2(rand(p,N/freq), X/freq, Y/freq, 'spline');
end

% 生成Fractal Brownian Motion
fbm = zeros(N);
for octave = 1:nOctaves
    freq = 2^(octave-1);
    amp = persistence^(octave-1);
    [X,Y] = meshgrid(1:N);
    fbm = fbm + amp*interp2(rand(p,N/freq), X/freq, Y/freq, 'spline');
end

% 显示结果
figure;
subplot(1,2,1);
imagesc(noise);
colormap(gray);
title('Perlin Noise');
axis equal;
axis off;
subplot(1,2,2);
imagesc(fbm);
colormap(gray);
title('Fractal Brownian Motion');
axis equal;
axis off;

该代码首先使用Perlin噪声生成随机的二维图像。然后使用Fractal Brownian Motion技术增加表面的粗糙度。最后，将结果显示在两个子图中，分别显示Perlin噪声和Fractal Brownian Motion的结果。