高手看看我这个算例，为什么 uw1，us1，uw2，us2的表达式不同，但计算结果却一样？请高手指导。
clc,clear;
%%%%%%%%%%%%%%%%%%%%%%%%% 在此处输入参数
%%%function [uw1,uw2,us1,us2]=ppcal(H)
TV=0.1;
H=168;
DWA=8.4; %%% 桩体直径
NW=0.0004; %%% kh/kw
NS=5; %%% kh/ks
SW=0.4; %%% 桩间距
Cv=0.001;
Ch=0.002;
ES=10;
EP=500;
s=3;
%%%%%%%%%%%%%%%%%%%%%%%%
q0=100000;
x=[0:0.05:1];
%%%%%%%%%%%%%%%%%%%%%%%%
G=NW*H^2/DWA^2;
Y=EP/ES;
RE=sqrt(2*sqrt(3)/pi)*SW*100/2;
n=2*RE/DWA;
N=n^2;
VH=Cv/Ch;
a=1/NS;
%%%%%%%%%%%%%%%%%%%%%%%% 第二种模式
FA=(s-1)*log(a*s)/(a*s-1)-(s-1)^2/n^2/(1-a);
FA=FA-2*(s-1)*(a*s-1)*log(a)/n^2/(1-a)^2;
FA=FA-(s-1)*(2*s^3-3*s^2+1)/3/n^4/(1-a);
FA=FA-(s-1)^2*(s-3*a*s-a+3)*(a*s-1)/2/(1-a)^3/n^4;
FA=FA+(s-1)*(a*s-1)^3*log(a)/n^4/(1-a)^4;
FA=FA-(n^2-s^2)*(1-s)^2/n^4/(1-a)+log(n/s)-3/4;
FA=FA+(4*n^2*s^2-s^4)/(4*n^4);
FA=FA*n^2/(n^2-1);
CQ=2*pi*Ch/ES/FA*(1-1/N);
%%%%%%%%%%%%%%%%%%%%%%%% 计算区
uw1=0;
us1=0;
uw2=0;
us2=0;
for m=1:1000
M=(2*m-1)*pi/2;
BM=VH*4*RE^2*M^2/H^2+8*(N-1)*VH*NW/FA+8/FA;
BM=BM*(N-1+Y)/(N-1)/(8*G/FA/M^2+1);
DCQ=pi*DWA^2/4*Ch*NW/ES*M^2/H^2;
ECQ=(N-1)*DCQ/CQ+N;
e1=-BM*TV;
uw1=uw1+2*q0*CQ/M/(CQ+DCQ)*sin(M*x)*exp(e1);
us1=us1+2*q0/M*sin(M*x)*exp(e1);
uw2=uw2+2*q0*N/M/ECQ*sin(M*x)*exp(e1);
us2=us2+2*q0*N*(1+DCQ/CQ)/M/ECQ*sin(M*x)*exp(e1);
end
disp(uw1/q0);
disp(us1/q0);
disp(uw2/q0);
disp(us2/q0);
