matlab方格网法土方量计算程序

%计算初始设计标高
%适合规则场地的计算
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-