%计算方程： W=10*10^4; %单位N a=1250; %单位mm b=1500; h=1330; L=4*b; C=h/L; A=asin(C); L1=2*b+a; L2=2*b-a; D=L1*sin(A)/L2/cos(A); B=atan (D); P=L*cos(A)/[ L1*sin(B+A)]*W; P0=1.5*P/2; %单位N