求助：关于粗糙表面的分形理论

wuchengyun9226

有大神对粗糙表面的分形理论有研究吗？最近看关于摩擦学的分形，里面有一个W-M函数及其曲线图的表示，看不懂，想请教大神帮忙看看下面的图怎么画！相当于是这个公式

以下是我用MATLAB写的程序，但画出来的图跟这个相差太大，请高手帮忙看看：
clear all
clc
lammada=1.5;
D=1.5;
g1=[];g2=[];
for x=linspace(0,0.0001,200)
    for n=-200:200;
    g=lammada^((D-2)*n)*(1-cos(lammada^n*x));
    g1=[g1;g];
    end
    g2=[g2;sum(g1)];
end
x=linspace(0,1,200);
plot(x,g2)
hold on

liuyalong008

1. n=-1000:1000;
   
2. x = 0:0.001:1;
   
3. d = [1.2 1.5 1.8]
   
4. f = @(x,n,d)1.5.^((d-2).*n).*(1-cos(x.*1.5.^n))
   
5. g = @(x,d)arrayfun(@(x)trapz(n,f(x,n,d)),x)
   
6. plot(x,[g(x,d(1));g(x,d(2));g(x,d(3))])
   

复制代码

caoyaocui

楼主这个问题很感兴趣！
刚刚用Mathcad试了下，没画出来。
坐等Matlab板块的大侠出手啊！{:3_63:}

TBE_Legend

111111

TBE_Legend

TBE_Legend 发表于 2014-12-18 20:11
111111

fasdfasdfa 222

wuchengyun9226

TBE_Legend 发表于 2014-12-18 20:47
fasdfasdfa 222

大神，我想问下这个是MATLAB仿出来的吗，我看不太懂您的这个程序啊

wuchengyun9226

TBE_Legend 发表于 2014-12-18 20:47
fasdfasdfa 222

请问一下在MATLAB里这个命令应该怎么表达

wuchengyun9226

liuyalong008 发表于 2014-12-19 10:54

太感谢了 ！

wuchengyun9226

TBE_Legend 发表于 2014-12-18 20:11
111111

很感谢大神的帮助

wuchengyun9226

wuchengyun9226 发表于 2014-12-19 13:25
很感谢大神的帮助

大神有时间能帮我看看这个帖子吗，还是那个分形模拟的问题，只是函数稍微有点改变：http://forum.simwe.com/thread-1118296-1-1.html