function 中23行,括号不匹配,修改为:MS = ( ( BS(:,1)/toa(1) + BS(:,2)/toa(2) + BS(:,3)/toa(3) + BS(:,4)/toa(4) ) / ( 1/toa(1) +1/toa(2)+1/toa(3)+1/toa(4) ) )';44行,矩阵维数不匹配,修改为:EP = iEP + delt; 运行结果:iEP = 1.0e-007 * 0.1179 0.2357