源程序.txt

本软件研究了地基的计算方法

    %输入论文57,aiN2-1%
    %改停%
    A(i+N2,i+N2)=k2*D(i-bc0*N2,N2)+2*DL(i-bc0*N2,N2)*(1-ul(i-bc0*N2,N2))+2*DL(i-bc0*N2-1,N2)*(1-ul(i-bc0*N2-1,N2))+1/2*k1*(D(i-bc0*N2+1,column)*(1-u(i-bc0*N2+1,column)*u(i-bc0*N2+1,column))+2*EI4/h2)+1/2*k1*(D(i-bc0*N2-1,column)*(1-u(i-bc0*N2-1,column)*u(i-bc0*N2-1,column))+2*EI4/h2)+2*k1*(D(i-bc0*N2,column)*(1-u(i-bc0*N2,column)*u(i-bc0*N2,column))+2*EI4/h2)+1/2*C(i-bc0*N2,column)*h12;
    %    a(i,N2)=k2*D(i,N2-1)         +2*DL(i+1,N2)*(1-ul(i+1,N2))                  +2*DL(i,N2)*(1-ul(i,N2))                          +1/2*k1*(D(i+1,N2)*(1-u(i+1,N2)*u(i+1,N2))+2*EI4/h2)                                       +1/2*k1*(D(i-1,N2)*(1-u(i-1,N2)*u(i-1,N2))+2*EI4/h2)                                       +2*k1*(D(i,N2)*(1-u(i,N2)*u(i,N2))+2*EI4/h2)                                       +1/2*C(i,N2)*h1*h1*h2*h2%
    %输入纸17，aiN2，WiN2的系数%
    bc0=bc0+1;
end

%                            下面为地基中的第N1-2列                %
A((N1-3)*N2+N1-2,(N1-3)*N2+N1-2)=k2*D(N1-2,2)+2*DL(N1-3,1)*(1-ul(N1-3,1))+2*DL(N1-2,1)*(1-ul(N1-2,1))+1/2*k1*(D(N1-3,1)*(1-u(N1-3,1)*u(N1-3,1))+2*EI2/h2)+1/2*k1*(D(N1-1,1)*(1-u(N1-1,1)*u(N1-1,1))+2*EI2/h2)+2*k1*(D(N1-2,1)*(1-u(N1-2,1)*u(N1-2,1))+2*EI2/h2)+1/2*EI2*k1/h2+1/2*C(N1-2,1)*h12;
%                      a(N1-3,0)=k2*D(N1-3,1)+2*DL(N1-3,1)*(1-ul(N1-3,1))+2*DL(N1-2,1)*(1-ul(N1-2,1))+1/2*k1*(D(N1-4,0)*(1-u(N1-4,0)*u(N1-4,0))+2*EI2/h2)+1/2*k1*(D(N1-2,0)*(1-u(N1-2,0)*u(N1-2,0))+2*EI2/h2)+2*k1*(D(N1-3,0)*(1-u(N1-3,0)*u(N1-3,0))+2*EI2/h2)+1/2*EI2*k1/h2+1/2*C(N1-3,0)*h1*h1*h2*h2%
%输入纸18，a(N1-3)0，W(N1-3)0的系数%
A((N1-3)*N2+N1-1,(N1-3)*N2+N1-1)=k1*(D(N1-3,2)+D(N1-1,2))+k2*D(N1-2,3)+4*D(N1-2,2)*(k1+k2+2*u(N1-2,2))+2*DL(N1-3,1)*(1-ul(N1-3,1))+2*DL(N1-3,2)*(1-ul(N1-3,2))+2*DL(N1-2,1)*(1-ul(N1-2,1))+2*DL(N1-2,2)*(1-ul(N1-2,2))+C(N1-2,2)*h12;
%                      a(N1-3,1)=k1*(D(N1-4,1)+D(N1-2,1))+k2*D(N1-3,2)+4*D(N1-3,1)*(k1+k2+2*u(N1-3,1))+2*DL(N1-3,1)*(1-ul(N1-3,1))+2*DL(N1-3,2)*(1-ul(N1-3,2))+2*DL(N1-2,1)*(1-ul(N1-2,1))+2*DL(N1-2,2)*(1-ul(N1-2,2))+C(N1-3,1)*h1*h1*h2*h2%
%输入论文57，ai1,i=N1-3%
%N2为列节点末点的坐标，通式%
for j=(N1-3)*N2+N1:(N1-2)*N2+N1-4   
    A(j,j)=k1*(D(N1-3,j-(N1-3)*N2-N1+3)+D(N1-1,j-(N1-3)*N2-N1+3))+k2*(D(N1-2,j-(N1-3)*N2-N1+2)+D(N1-2,j-(N1-3)*N2-N1+4))+4*D(N1-2,j-(N1-3)*N2-N1+3)*(k1+k2+2*u(N1-2,j-(N1-3)*N2-N1+3))+2*DL(N1-3,j-(N1-3)*N2-N1+2)*(1-ul(N1-3,j-(N1-3)*N2-N1+2))+2*DL(N1-3,j-(N1-3)*N2-N1+3)*(1-ul(N1-3,j-(N1-3)*N2-N1+3))+2*DL(N1-2,j-(N1-3)*N2-N1+2)*(1-ul(N1-2,j-(N1-3)*N2-N1+2))+2*DL(N1-2,j-(N1-3)*N2-N1+3)*(1-ul(N1-2,j-(N1-3)*N2-N1+3))+C(N1-2,j-(N1-3)*N2-N1+3)*h12;
%a(N1-3,j)=k1*(D(N1-4,j)               +D(N1-2,j))               +k2*(D(N1-3,j-1)             +D(N1-3,j+1))+4*D(N1-3,j)*(k1+k2+2*u(N1-3,j))+2*DL(N1-3,j)*(1-ul(N1-3,j))                                                                         +2*DL(N1-3,j+1)*(1-ul(N1-3,j+1))+2*DL(N1-2,j)*(1-ul(N1-2,j))+2*DL(N1-2,j+1)*(1-ul(N1-2,j+1))+C(N1-3,j)*h1*h1*h2*h2%
    %输入论文57，aij，i=N1-3,j=2,...,N2-2的系数%
end
A((N1-2)*N2+N1-3,(N1-2)*N2+N1-3)=k1*(D(N1-3,N2)+D(N1-1,N2))    +k2*D(N1-2,N2-1)+4*D(N1-2,N2)*(k1+k2+2*u(N1-2,N2))    +2*DL(N1-3,N2-1)*(1-ul(N1-3,N2-1))+2*DL(N1-3,N2)*(1-ul(N1-3,N2))+2*DL(N1-2,N2-1)*(1-ul(N1-2,N2-1))+2*DL(N1-2,N2)*(1-ul(N1-2,N2))+C(N1-2,N2)*h12;
%                   a(N1-3,N2-1)=k1*(D(N1-4,N2-1)+D(N1-2,N2-1))+k2*D(N1-3,N2-2)+4*D(N1-3,N2-1)*(k1+k2+2*u(N1-3,N2-1))+2*DL(N1-3,N2-1)*(1-ul(N1-3,N2-1))+2*DL(N1-3,N2)*(1-ul(N1-3,N2))+2*DL(N1-2,N2-1)*(1-ul(N1-2,N2-1))+2*DL(N1-2,N2)*(1-ul(N1-2,N2))+C(N1-3,N2-1)*h1*h1*h2*h2%
%输入论文57，aiN2-1,i=N1-3%
A((N1-2)*N2+N1-2,(N1-2)*N2+N1-2)=k2*D(N1-2,N2)+2*DL(N1-3,N2)*(1-ul(N1-3,N2))+2*DL(N1-2,N2)*(1-ul(N1-2,N2))+1/2*k1*(D(N1-3,column)*(1-u(N1-3,column)*u(N1-3,column))+2*EI4/h2)+1/2*k1*(D(N1-1,column)*(1-u(N1-1,column)*u(N1-1,column))+2*EI4/h2)+2*k1*(D(N1-2,column)*(1-u(N1-2,column)*u(N1-2,column))+2*EI4/h2)+1/2*EI4*k1/h2+1/2*C(N1-2,column)*h12;
%                     a(N1-3,N2)=k2*D(N1-3,N2-1)+2*DL(N1-3,N2)*(1-ul(N1-3,N2))+2*DL(N1-2,N2)*(1-ul(N1-2,N2))+1/2*k1*(D(N1-4,N2)*(1-u(N1-4,N2)*u(N1-4,N2))+2*EI4/h2)+1/2*k1*(D(N1-2,N2)*(1-u(N1-2,N2)*u(N1-2,N2))+2*EI4/h2)+2*k1*(D(N1-3,N2)*(1-u(N1-3,N2)*u(N1-3,N2))+2*EI4/h2)+1/2*EI4*k1/h2+1/2*C(N1-3,N2)*h1*h1*h2*h2%
%输入纸19，a(N1-3)N2，W(N1-3)N2的系数%

%                            下面为地基中的第N1-1列                %

A((N1-2)*N2+N1-1,(N1-2)*N2+N1-1)=k2*D(N1-1,2)+2*DL(N1-1,1)*(1-ul(N1-1,1))+2*DL(N1-2,1)*(1-ul(N1-2,1))+1/2*k1*(D(N1,1)*(1-u(N1,1)*u(N1,1))+2*EI2/h2)+1/2*k1*(D(N1-2,1)*(1-u(N1-2,1)*u(N1-2,1))+2*EI2/h2)+2*k1*(D(N1-1,1)*(1-u(N1-1,1)*u(N1-1,1))+2*EI2/h2)+8*EI2*k1/h2+1/2*C(N1-1,1)*h12;
%                      a(N1-2,0)=k2*D(N1-2,1)+2*DL(N1-1,1)*(1-ul(N1-1,1))+2*DL(N1-2,1)*(1-ul(N1-2,1))+1/2*k1*(D(N1-1,0)*(1-u(N1-1,0)*u(N1-1,0))+2*EI2/h2)+1/2*k1*(D(N1-3,0)*(1-u(N1-3,0)*u(N1-3,0))+2*EI2/h2)+2*k1*(D(N1-2,0)*(1-u(N1-2,0)*u(N1-2,0))+2*EI2/h2)+8*EI2*k1/h2+1/2*C(N1-2,0)*h1*h1*h2*h2%
%输入纸20，a(N1-2)0，W(N1-2)0的系数%
A((N1-2)*N2+N1,(N1-2)*N2+N1)=k1*(D(N1-2,2)+D(N1,2))+k2*D(N1-1,3)+4*D(N1-1,2)*(k1+k2+2*u(N1-1,2))+2*DL(N1-2,1)*(1-ul(N1-2,1))+2*DL(N1-2,2)*(1-ul(N1-2,2))+2*DL(N1-1,1)*(1-ul(N1-1,1))+2*DL(N1-1,2)*(1-ul(N1-1,2))+C(N1-1,2)*h12;
%                  a(N1-2,1)=k1*(D(N1-3,1)+D(N1-1,1))+k2*D(N1-2,2)+4*D(N1-2,1)*(k1+k2+2*u(N1-2,1))+2*DL(N1-2,1)*(1-ul(N1-2,1))+2*DL(N1-2,2)*(1-ul(N1-2,2))+2*DL(N1-1,1)*(1-ul(N1-1,1))+2*DL(N1-1,2)*(1-ul(N1-1,2))+C(N1-2,1)*h1*h1*h2*h2%
%输入论文57，ai1,i=N1-2%
%N2为列节点末点的坐标，通式%
for j=(N1-2)*N2+row:(N1-1)*N2+N1-3   
    A(j,j)=k1*(D(N1-2,j-(N1-2)*N2-N1+2)+D(N1,j-(N1-2)*N2-N1+2))+k2*(D(N1-1,j-(N1-2)*N2-N1+1)+D(N1-1,j-(N1-2)*N2-N1+3))+4*D(N1-1,j-(N1-2)*N2-N1+2)*(k1+k2+2*u(N1-1,j-(N1-2)*N2-N1+2))+2*DL(N1-2,j-(N1-2)*N2-N1+1)*(1-ul(N1-2,j-(N1-2)*N2-N1+1))+2*DL(N1-2,j-(N1-2)*N2-N1+2)*(1-ul(N1-2,j-(N1-2)*N2-N1+2))+2*DL(N1-1,j-(N1-2)*N2-N1+1)*(1-ul(N1-1,j-(N1-2)*N2-N1+1))+2*DL(N1-1,j-(N1-2)*N2-N1+2)*(1-ul(N1-1,j-(N1-2)*N2-N1+2))+C(N1-1,j-(N1-2)*N2-N1+2)*h12;
%a(N1-2,j)=k1*(D(N1-3,j)               +D(N1-1,j))             +k2*(D(N1-2,j-1)             +D(N1-2,j+1))             +4*D(N1-2,j)*(k1+k2+2*u(N1-2,j))                              +2*DL(N1-2,j)*(1-ul(N1-2,j))                              +2*DL(N1-2,j+1)*(1-ul(N1-2,j+1))                          +2*DL(N1-1,j)*(1-ul(N1-1,j))                              +2*DL(N1-1,j+1)*(1-ul(N1-1,j+1))                          +C(N1-2,j)*h1*h1*h2*h2%
%输入论文57，aij，i=N1-2,j=2,...,N2-2的系数%
end
A((N1-1)*N2+N1-2,(N1-1)*N2+N1-2)=k1*(D(N1-2,N2)+D(N1,N2))      +k2*D(N1-1,N2-1)+4*D(N1-1,N2)*(k1+k2+2*u(N1-1,N2))+2*DL(N1-2,N2-1)*(1-ul(N1-2,N2-1))+2*DL(N1-2,N2)*(1-ul(N1-2,N2))+2*DL(N1-1,N2-1)*(1-ul(N1-1,N2-1))+2*DL(N1-1,N2)*(1-ul(N1-1,N2))+C(N1-1,N2)*h12;
%                   a(N1-2,N2-1)=k1*(D(N1-3,N2-1)+D(N1-1,N2-1))+k2*D(N1-2,N2-2)+4*D(N1-2,N2-1)*(k1+k2+2*u(N1-2,N2-1))+2*DL(N1-2,N2-1)*(1-ul(N1-2,N2-1))+2*DL(N1-2,N2)*(1-ul(N1-2,N2))+2*DL(N1-1,N2-1)*(1-ul(N1-1,N2-1))+2*DL(N1-1,N2)*(1-ul(N1-1,N2))+C(N1-2,N2-1)*h1*h1*h2*h2%
%输入论文57，aiN2-1,i=N1-2%
A((N1-1)*N2+N1-1,(N1-1)*N2+N1-1)=k2*D(N1-1,N2)  +2*DL(N1-2,N2)*(1-ul(N1-2,N2))+2*DL(N1-1,N2)*(1-ul(N1-1,N2))+1/2*k1*(D(N1-2,column)*(1-u(N1-2,column)*u(N1-2,column))+2*EI4/h2)+1/2*k1*(D(N1,column)*(1-u(N1,column)*u(N1,column))+2*EI4/h2)+2*k1*(D(N1-1,column)*(1-u(N1-1,column)*u(N1-1,column))+2*EI4/h2)+8*EI4*k1/h2+1/2*C(N1-1,column)*h12;
%                     a(N1-2,N2)=k2*D(N1-2,N2-1)+2*DL(N1-2,N2)*(1-ul(N1-2,N2))+2*DL(N1-1,N2)*(1-ul(N1-1,N2))+1/2*k1*(D(N1-3,N2)*(1-u(N1-3,N2)*u(N1-3,N2))+2*EI4/h2)            +1/2*k1*(D(N1-1,N2)*(1-u(N1-1,N2)*u(N1-1,N2))+2*EI4/h2)      +2*k1*(D(N1-2,N2)*(1-u(N1-2,N2)*u(N1-2,N2))+2*EI4/h2)+8*EI4*k1/h2            +1/2*C(N1-2,N2)*h1*h1*h2*h2%
%输入纸21，a(N1-2)N2，W(N1-2)N2的系数%

%                            下面为地基中的第N1列                %
A((N1-1)*N2+N1,(N1-1)*N2+N1)=k2*D(N1,2)+2*DL(N1,1)*(1-ul(N1,1))+2*DL(N1-1,1)*(1-ul(N1-1,1))+1/2*k1*(D(N1-1,1)*(1-u(N1-1,1)*u(N1-1,1))+2*EI2/h2)+2*k1*(D(N1,1)*(1-u(N1,1)*u(N1,1))+2*EI2/h2)+25/2*EI2*k1/h2+1/2*C(N1,1)*h12;
%                  a(N1-1,0)=k2*D(N1-1,1)+2*DL(N1,1)*(1-ul(N1,1))+2*DL(N1-1,1)*(1-ul(N1-1,1))+1/2*k1*(D(N1-2,0)*(1-u(N1-2,0)*u(N1-2,0))+2*EI2/h2)+2*k1*(D(N1-1,0)*(1-u(N1-1,0)*u(N1-1,0))+2*EI2/h2)+25/2*EI2*k1/h2+1/2*C(N1-1,0)*h1*h1*h2*h2%
%输入纸22，a(N1-1)0，W(N1-1)0的系数%
A((N1-1)*N2+row,(N1-1)*N2+row)=k1*D(N1-1,2)+k2*D(N1,3)+4*D(N1,2)*(k1+k2+2*u(N1,2))+2*DL(N1-1,1)*(1-ul(N1-1,1))+2*DL(N1-1,2)*(1-ul(N1-1,2))+2*DL(N1,1)*(1-ul(N1,1))+2*DL(N1,2)*(1-ul(N1,2))+C(N1,2)*h12;
%                      a(N1-1,1)=k1*D(N1-2,1)+k2*D(N1-1,2)+4*D(N1-1,1)*(k1+k2+2*u(N1-1,1))+2*DL(N1-1,1)*(1-ul(N1-1,1))+2*DL(N1-1,2)*(1-ul(N1-1,2))+2*DL(N1,1)*(1-ul(N1,1))+2*DL(N1,2)*(1-ul(N1,2))+C(N1-1,1)*h1*h1*h2*h2%
%输入论文57，a(N1-1)1,i=N1-2%
%N2为列节点末点的坐标，通式%
for j=(N1-1)*N2+N1+2:N1*N2+N1-2   
    A(j,j)=k1*D(N1-1,j-(N1-1)*N2-N1+1)+k2*(D(N1,j-(N1-1)*N2-N1)+D(N1,j-(N1-1)*N2-N1+2))+4*D(N1,j-(N1-1)*N2-N1+1)*(k1+k2+2*u(N1,j-(N1-1)*N2-N1+1))+2*DL(N1-1,j-(N1-1)*N2-N1)*(1-ul(N1-1,j-(N1-1)*N2-N1))+2*DL(N1-1,j-(N1-1)*N2-N1+1)*(1-ul(N1-1,j-(N1-1)*N2-N1+1))+2*DL(N1,j-(N1-1)*N2-N1)*(1-ul(N1,j-(N1-1)*N2-N1))+2*DL(N1,j-(N1-1)*N2-N1+1)*(1-ul(N1,j-(N1-1)*N2-N1+1))+C(N1,j-(N1-1)*N2-N1+1)*h12;
%a(N1-1,j)=k1*D(N1-2,j)               +k2*(D(N1-1,j-1)+D(N1-1,j+1))                    +4*D(N1-1,j)*(k1+k2+2*u(N1-1,j))                          +2*DL(N1-1,j)*(1-ul(N1-1,j))                          +2*DL(N1-1,j+1)*(1-ul(N1-1,j+1))                          +2*DL(N1,j)*(1-ul(N1,j))+2*DL(N1,j+1)*(1-ul(N1,j+1))+C(N1-1,j)*h1*h1*h2*h2%   
    %输入论文57，a(N1-1)j,j=2,...,N2-2的系数%
end
A(N1*N2+N1-1,N1*N2+N1-1)=k1*D(N1-1,N2)+k2*D(N1,N2-1)+4*D(N1,N2)*(k1+k2+2*u(N1,N2))+2*DL(N1-1,N2-1)*(1-ul(N1-1,N2-1))+2*DL(N1-1,N2)*(1-ul(N1-1,N2))+2*DL(N1,N2-1)*(1-ul(N1,N2-1))+2*DL(N1,N2)*(1-ul(N1,N2))+C(N1,N2)*h12;
%           a(N1-1,N2-1)=k1*D(N1-2,N2-1)+k2*D(N1-1,N2-2)+4*D(N1-1,N2-1)*(k1+k2+2*u(N1-1,N2-1))+2*DL(N1-1,N2-1)*(1-ul(N1-1,N2-1))+2*DL(N1-1,N2)*(1-ul(N1-1,N2))+2*DL(N1,N2-1)*(1-ul(N1,N2-1))+2*DL(N1,N2)*(1-ul(N1,N2))+C(N1-1,N2-1)*h1*h1*h2*h2%
%输入论文57，a(N1-1)N2-1%
A(N1*N2+N1,N1*N2+N1)=k2*D(N1,N2)+2*DL(N1,N2)*(1-ul(N1,N2))+2*DL(N1-1,N2)*(1-ul(N1-1,N2))+1/2*k1*(D(N1-1,column)*(1-u(N1-1,column)*u(N1-1,column))+2*EI4/h2)+2*k1*(D(N1,column)*(1-u(N1,column)*u(N1,column))+2*EI4/h2)+25/2*EI4*k1/h2+1/2*C(N1,column)*h12;
%         a(N1-1,N2)=k2*D(N1-1,N2-1)+2*DL(N1,N2)*(1-ul(N1,N2))+2*DL(N1-1,N2)*(1-ul(N1-1,N2))+1/2*k1*(D(N1-2,N2)*(1-u(N1-2,N2)*u(N1-2,N2))+2*EI4/h2)+2*k1*(D(N1-1,N2)*(1-u(N1-1,N2)*u(N1-1,N2))+2*EI4/h2)+25/2*EI4*k1/h2+1/2*C(N1-1,N2)*h1*h1*h2*h2%
%输入纸23，a(N1-1)N2，W(N1-1)N2的系数%

%                            下面为地基中的第N1+1列(最后一列)                %
A(N1*N2+row,N1*N2+row)=2*DL(N1,1)*(1-ul(N1,1))+1/2*k2*(D(row,2)*(1-u(row,2)*u(row,2))+2*EI3/h1)+1/2*k1*(D(N1,1)*(1-u(N1,1)*u(N1,1))+2*EI2/h2)+2*EI3*k2/h1+2*EI2*k1/h2+1/4*C(row,1)*h12;
%                a(N1,0)=2*DL(N1,1)*(1-ul(N1,1))+1/2*k2*(D(N1,1)*(1-u(N1,1)*u(N1,1))+2*EI3/h1)+1/2*k1*(D(N1-1,0)*(1-u(N1-1,0)*u(N1-1,0))+2*EI2/h2)+2*EI3*k2/h1+2*EI2*k1/h2+1/4*C(N1,0)*h1*h1*h2*h2%
%输入纸24，aN1.0，WN1.0的系数%
A(N1*N2+N1+2,N1*N2+N1+2)=k1*D(N1,2)+2*DL(N1,1)*(1-ul(N1,1))+2*DL(N1,2)*(1-ul(N1,2))+1/2*k2*(D(row,3)*(1-u(row,3)*u(row,3))+2*EI3/h1)+2*k2*(D(row,2)*(1-u(row,2)*u(row,2))+2*EI3/h1)+25/2*EI3*k2/h1+1/2*C(row,2)*h12;
%                a(N1,1)=k1*D(N1-1,1)+2*DL(N1,1)*(1-ul(N1,1))+2*DL(N1,2)*(1-ul(N1,2))+1/2*k2*(D(N1,2)*(1-u(N1,2)*u(N1,2))+2*EI3/h1)+2*k2*(D(N1,1)*(1-u(N1,1)*u(N1,1))+2*EI3/h1)+25/2*EI3*k2/h1+1/2*C(N1,1)*h1*h1*h2*h2%
%输入纸25，aN1.1，WN1.1的系数%
A(N1*N2+N1+3,N1*N2+N1+3)=k1*D(N1,3)+2*DL(N1,2)*(1-ul(N1,2))+2*DL(N1,3)*(1-ul(N1,3))+1/2*k2*(D(row,4)*(1-u(row,4)*u(row,4))+2*EI3/h1)+2*k2*(D(row,3)*(1-u(row,3)*u(row,3))+2*EI3/h1)+1/2*k2*(D(row,2)*(1-u(row,2)*u(row,2))+2*EI3/h1)+8*EI3*k2/h1+1/2*C(row,3)*h12;
%                a(N1,2)=k1*D(N1-1,2)+2*DL(N1,2)*(1-ul(N1,2))+2*DL(N1,3)*(1-ul(N1,3))+1/2*k2*(D(N1,3)*(1-u(N1,3)*u(N1,3))+2*EI3/h1)+2*k2*(D(N1,2)*(1-u(N1,2)*u(N1,2))+2*EI3/h1)+1/2*k2*(D(N1,1)*(1-u(N1,1)*u(N1,1))+2*EI3/h1)+8*EI3*k2/h1+1/2*C(N1,2)*h1*h1*h2*h2%
%输入纸26，aN1.2，WN1.2的系数%
A(N1*N2+N1+4,N1*N2+N1+4)=k1*D(N1,4)+2*DL(N1,4)*(1-ul(N1,4))+2*DL(N1,3)*(1-ul(N1,3))+1/2*k2*(D(row,5)*(1-u(row,5)*u(row,5))+2*EI3/h1)+1/2*k2*(D(row,3)*(1-u(row,3)*u(row,3))+2*EI3/h1)+2*k2*(D(row,4)*(1-u(row,4)*u(row,4))+2*EI3/h1)+1/2*EI3*k2/h1+1/2*C(row,4)*h12;
%                a(N1,3)=k1*D(N1-1,3)+2*DL(N1,4)*(1-ul(N1,4))+2*DL(N1,3)*(1-ul(N1,3))+1/2*k2*(D(N1,4)*(1-u(N1,4)*u(N1,4))+2*EI3/h1)+1/2*k2*(D(N1,2)*(1-u(N1,2)*u(N1,2))+2*EI3/h1)+2*k2*(D(N1,3)*(1-u(N1,3)*u(N1,3))+2*EI3/h1)+1/2*EI3*k2/h1+1/2*C(N1,3)*h1*h1*h2*h2%
%输入纸27，aN1.3，WN1.3的系数%
%N2为列节点末点的坐标，通式%
for j=N1*N2+N1+5:row*N2+N1-3   
    A(j,j)=k1*D(N1,j-N1*N2-N1)+2*DL(N1,j-N1*N2-N1)*(1-ul(N1,j-N1*N2-N1))+2*DL(N1,j-N1*N2-N1-1)*(1-ul(N1,j-N1*N2-N1-1))+1/2*k2*(D(row,j-N1*N2-N1+1)*(1-u(row,j-N1*N2-N1+1)*u(row,j-N1*N2-N1+1))+2*EI3/h1)+1/2*k2*(D(row,j-N1*N2-N1-1)*(1-u(row,j-N1*N2-N1-1)*u(row,j-N1*N2-N1-1))+2*EI3/h1)+2*k2*(D(row,j-N1*N2-N1)*(1-u(row,j-N1*N2-N1)*u(row,j-N1*N2-N1))+2*EI3/h1)+1/2*C(row,j-N1*N2-N1)*h12;
  %a(N1,j)=k1*D(N1-1,j)       +2*DL(N1,j+1)*(1-ul(N1,j+1))+2*DL(N1,j)*(1-ul(N1,j))+1/2*k2*(D(N1,j+1)*(1-u(N1,j+1)*u(N1,j+1))+2*EI3/h1)                                                                     +1/2*k2*(D(N1,j-1)*(1-u(N1,j-1)*u(N1,j-1))+2*EI3/h1)+2*k2*(D(N1,j)*(1-u(N1,j)*u(N1,j))+2*EI3/h1)+1/2*C(N1,j)*h1*h1*h2*h2%
    %输入纸28，aN1j,j=4,...,N2-4的系数%
end
A(row*N2+N1-2,row*N2+N1-2)=k1*D(N1,N2-2)+2*DL(N1,N2-3)*(1-ul(N1,N2-3))+2*DL(N1,N2-2)*(1-ul(N1,N2-2))+1/2*k2*(D(row,N2-3)*(1-u(row,N2-3)*u(row,N2-3))+2*EI3/h1)+1/2*k2*(D(row,N2-1)*(1-u(row,N2-1)*u(row,N2-1))+2*EI3/h1)+2*k2*(D(row,N2-2)*(1-u(row,N2-2)*u(row,N2-2))+2*EI3/h1)+1/2*EI3*k2/h1+1/2*C(row,N2-2)*h12;
%                     a(N1,N2-3)=k1*D(N1-1,N2-3)+2*DL(N1,N2-3)*(1-ul(N1,N2-3))+2*DL(N1,N2-2)*(1-ul(N1,N2-2))+1/2*k2*(D(N1,N2-4)*(1-u(N1,N2-4)*u(N1,N2-4))+2*EI3/h1)+1/2*k2*(D(N1,N2-2)*(1-u(N1,N2-2)*u(N1,N2-2))+2*EI3/h1)+2*k2*(D(N1,N2-3)*(1-u(N1,N2-3)*u(N1,N2-3))+2*EI3/h1)+1/2*EI3*k2/h1+1/2*C(N1,N2-3)*h1*h1*h2*h2%
%输入纸29，aN1(N2-3)，WN1(N2-3)的系数%
A(row*N2+N1-1,row*N2+N1-1)=k1*D(N1,N2-1)+2*DL(N1,N2-2)*(1-ul(N1,N2-2))+2*DL(N1,N2-1)*(1-ul(N1,N2-1))+1/2*k2*(D(row,N2-2)*(1-u(row,N2-2)*u(row,N2-2))+2*EI3/h1)+1/2*k2*(D(row,N2)*(1-u(row,N2)*u(row,N2))+2*EI3/h1)+2*k2*(D(row,N2-1)*(1-u(row,N2-1)*u(row,N2-1))+2*EI3/h1)+8*EI3*k2/h1+1/2*C(row,N2-1)*h12;
%                     a(N1,N2-2)=k1*D(N1-1,N2-2)+2*DL(N1,N2-2)*(1-ul(N1,N2-2))+2*DL(N1,N2-1)*(1-ul(N1,N2-1))+1/2*k2*(D(N1,N2-3)*(1-u(N1,N2-3)*u(N1,N2-3))+2*EI3/h1)+1/2*k2*(D(N1,N2-1)*(1-u(N1,N2-1)*u(N1,N2-1))+2*EI3/h1)+2*k2*(D(N1,N2-2)*(1-u(N1,N2-2)*u(N1,N2-2))+2*EI3/h1)+8*EI3*k2/h1+1/2*C(N1,N2-2)*h1*h1*h2*h2%
%输入纸30，aN1(N2-2)，WN1(N2-2)的系数%
A(row*N2+N1,row*N2+N1)=k1*D(N1,N2)+2*DL(N1,N2)*(1-ul(N1,N2))+2*DL(N1,N2-1)*(1-ul(N1,N2-1))+1/2*k2*(D(row,N2-1)*(1-u(row,N2-1)*u(row,N2-1))+2*EI3/h1)+2*k2*(D(row,N2)*(1-u(row,N2)*u(row,N2))+2*EI3/h1)+25/2*EI3*k2/h1+1/2*C(row,N2)*h12 ;
%                 a(N1,N2-1)=k1*D(N1-1,N2-1)+2*DL(N1,N2)*(1-ul(N1,N2))+2*DL(N1,N2-1)*(1-ul(N1,N2-1))+1/2*k2*(D(N1,N2-2)*(1-u(N1,N2-2)*u(N1,N2-2))+2*EI3/h1)+2*k2*(D(N1,N2-1)*(1-u(N1,N2-1)*u(N1,N2-1))+2*EI3/h1)+25/2*EI3*k2/h1+1/2*C(N1,N2-1)*h1*h1*h2*h2%
%输入纸31，aN1(N2-1)，WN1(N2-1)的系数%
A(row*column,row*column)=2*DL(N1,N2)*(1-ul(N1,N2))+1/2*k2*(D(row,N2)*(1-u(row,N2)*u(row,N2))+2*EI3/h1)+1/2*k1*(D(N1,column)*(1-u(N1,column)*u(N1,column))+2*EI4/h2)+2*EI4*k1/h2+2*EI3*k2/h1+1/4*C(row,column)*h12;
%               a(N1,N2)=2*DL(N1,N2)*(1-ul(N1,N2))+1/2*k2*(D(N1,N2-1)*(1-u(N1,N2-1)*u(N1,N2-1))+2*EI3/h1)+1/2*k1*(D(N1-1,N2)*(1-u(N1-1,N2)*u(N1-1,N2))+2*EI4/h2)+2*EI4*k1/h2+2*EI3*k2/h1+1/4*C(N1,N2)*h1*h1*h2*h2%
%输入纸32，aN1N2，WN1N2的系数%
%....................对A稀疏化，节省内存...................
% A=sparse(A);
 %.......................................................
 %                          验证A 对称不？？？？？？？？？？%
if A==A'
    fprintf(1,'A 是对称矩阵\n');
end 
%系数矩阵A赋值结束%
%    计算方程的常数项b，标号从(1)算起  %
 %                          下面对b进行赋值            %
b(1)=(1/4)*h12*q(1);%1对应地基中的点(0,0)的常数项%
b(column)=(1/4)*h12*q(column);%N2+1对应地基中的点(0,N2)的常数项%
b(N1*column+1)=(1/4)*h12*q(N1*column+1);%N1*(N2+1)+1对应地基中的点(N1，0)的常数项%
b(row*column)=(1/4)*h12*q(row*column);%对应地基中的点(N1,N2)的常数项%
for j=2:N2
    b(j)=(1/2)*h12*q(j);%对应地基中第一条梁EI1的点(0,1)到(0,N2-1)的常数项%
end
for j=N1*column+2:row*column-1
    b(j)=(1/2)*h12*q(j);%对应地基中第三条梁EI3的点(N1,1)到(N1,N2-1)的常数项%
end
for i=N2+2:column:(N1-1)*N2+N1
    b(i)=(1/2)*h12*q(i);%对应地基中第二条梁EI2的点(0,N2)的常数项%
end
for i=2*N2+2:column:N1*N2+N1
    b(i)=(1/2)*h12*q(i);%N2+1对应地基中第四条梁EI4的点(0,N2)的常数项%
end
for i=N2+3:column:(N1-1)*N2+column
    for j=i:i+N2-2
        b(j)=h12*q(j);%对应地基中的内部的点，i=1,2,...,N1-1,j=1,2,...,N2-1的常数项%
    end
end
%常数项b赋值结束%
%解方程得扰度，单位为毫米
%...........无迭代求解....................................
X=A\b;%或者inv(A)*b
X=X*10^3;
sj1=zeros(N1+1,N2+1);
sj1=reshape(X,N1+1,N2+1);
%输出结果%
fprintf(' \n.....................程序的运行结果；地基上(%d*%d)各点扰度值,单位为毫米(从左下角开始，按纵方向，列排)...............\n ',row,column);
for i=1:column
    fprintf(1,' 第 %d 列  ',i);%可以改进disp();不需循环
end
fprintf(1,'\n');
for i=1:1+N2
    for j=i:column:column*N1+i
        fprintf(1,'  %5.4f  ',X(j,1));
    end
    fprintf(1,'\n');
end
fprintf(' \n...........程序运行结束！！...............\n ');
%.....................