matlab生成W-M分形形面

一个二维的分形表面...

1 二维平面的 W-M 公式

$$z(x)=G^{D-1}\sum_{n=n1}^{\infty}\frac{cos(2 \pi \gamma^n x)}{\gamma ^{(2-D)n}}$$

$\gamma=1.5$；$n$ 为自然序列，通常取 $n=1,2,3…100$。

matlab 中实现，取 $G=1$，$D$ 为 $1$ 和 $1.5$:

clear
clc
G=1;
n= 0:1:100;
x = 1:1:1000;
d = [1 1.5]
f = @(x,n,d)1.5.^((d-2).*n).*cos(x.*1.5.^n.*2.*pi).*G.^(d-1)
g = @(x,d)arrayfun(@(x)trapz(n,f(x,n,d)),x)
plot(x,g(x,d(1)))
figure
plot(x,g(x,d(2)))