一维高斯随机粗糙面建模.zip

clear all close all clc %N： 表面点数 %h： 均方根高度 %l： 相关长度 %L： 表面方向长度 N = 128; l = 1; h = 0.05; L = 102.4; [x,f] = cucaomian1D(N,l,h,L); f1 = f; h = 0.2; [x,f] = cucaomian1D(N,l,h,L); f2=f; figure(1); plot(x,f1,x,f2,'k:');legend('h=0.05w','h=0.2w'); axis([-51 51 -1.0 1.0]); text(-15,-0.8,'l=1.0w , L=102.4w','fontname','Times New Roman','FontSize',12); xlabel('x/w','fontname','Times New Roman','FontSize',14); ylabel('f/w','fontname','Times New Roman','FontSize',14); l = 0.5; h = 0.1; [x,f] = cucaomian1D(N,l,h,L); f3 = f; l = 1.5; [x,f] = cucaomian1D(N,l,h,L); f4 = f; figure(2); plot(x,f3,x,f4,'k:');legend('l=0.5w','l=1.5w'); axis([-51 51 -1.0 1.0]); text(-15,-0.8,'h=0.1w , L=102.4w','fontname','Times New Roman','FontSize',12); xlabel('x/w','fontname','Times New Roman','FontSize',14); ylabel('f/w','fontname','Times New Roman','FontSize',14);