%%三維地質(zhì)建模-曲面圖
clear?all
d=50%空間采樣點(diǎn)數(shù)
T=1400%空間震蕩頻率
m=1000%空間區(qū)域長度?
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