%計算初始設計標高
%適合規則場地的計算
clear?allclc;
a=input（‘請輸入場地分塊正方形邊長:‘）;
disp（‘請輸入A矩陣‘）;??%可以直接在command窗口輸入A=[]來建立A矩陣從而輸入初始參數
%?A=xlsread（‘E:\shumo\A.xlsx‘）;
hang=size（A1）;
lie=size（A2）;
N=（hang-1）*（lie-1）;
Sum1=A（11）+A（1lie）+A（hang1）+A（hanglie）;
S=sum（A2）;
sum1=S（1）-A（11）-A（1lie）;
sum2=S（hang）-A（hang1）-A（hanglie）;
s=sum（A1）;
sum3=s（1）-A（11）-A（hang1）;
sum4=s（lie）-A（1lie）-A（hanglie）;
Sum2=2*（sum1+sum2+sum3+sum4）;
sum5=zeros（hang1）;
for?ii=2:hang-1
????sum5（ii）=S（ii）-A（ii1）-A（iilie）;
end
Sum3=4*sum（sum5）;
sigma=Sum1+Sum2+Sum3;
H0=sigma/（4*N）
[xy]=meshgrid（0:a:（（lie-1）*a）0:a:（（hang-1）*a））;
surfc（xyA）;
%?[xyZ]=griddata（xyAlinspace（min（x）max（x））‘linspace（min（y）max（y））‘cubic‘）;%三維樣條插值
%?meshc（xyA）;hold?on;
%?surf（xyA）;
%?shading?interp;



%不考慮土體的可松性影響，余土加寬邊坡，計算場地的施工高度
ix=input（‘請輸入場地的X方向坡度：‘）;
iy=input（‘請輸入場地的Y方向坡度：‘）;
x0=（lie+1）/2;
y0=（hang+1）/2;
Hn=zeros（hanglie）;
for?ii=1:fix（y0）
????for?jj=1:fix（x0）
????????Hn（iijj）=H0-