生成地质体三维建模理论模型MATLAB

%%三维地质建模-曲面图
clear all
d=50%空间采样点数
T=1400%空间震荡频率
m=1000%空间区域长度 
w=2*pi/T
T2=2000
w2=2*pi/T2
f1=fopen（‘E:\gocad_install\example_data\fault_1.txt‘‘w‘）;
f2=fopen（‘E:\gocad_install\example_data\fault_2.txt‘‘w‘）;
f3=fopen（‘E:\gocad_install\example_data\h_1.txt‘‘w‘）;
f4=fopen（‘E:\gocad_install\example_data\h_2.txt‘‘w‘）;
for i=0:d
      for j=0:d
          x=i*m/d;
          y=j*m/d;
%z=200*sin（w*x）+500;
%if x<350
z1=-0.5*x+700+200*sin（w*x）;%断层1
z2=-1*x+700+100*sin（w*x）;%断层2
z3=200*sin（w2*x）+100;%地层1
z4=200*sin（w2*x）+10*rand+300;%地层2
if z1>z3 &z