%迭代法
clear
clc

x1=210;
y1=-90;
z1=-620;
p1=0.2988;
% a1=-0.1096;
% b1=0.5996;
% c1=0.3507;
% d1=-0.3882;          %%设定初值

for  i=1:100      
    a=-0.1096+(i-1)*0.002648;
    b=0.5996+(i-1)*0.000269;
    c=0.3507+(i-1)*0.004468;
    d=-0.3882+(i-1)*0.007478;        %%第（i-1）的角度值

    A1=x1+122.9075-48*sin(p1+0.87)+62*sin(p1)-(2^(-1/2))*(265+375*sin(a));
    B1=y1+48.55-48*cos(p1+0.87)+62*cos(p1)-(2^(-1/2))*(285+375*sin(a));
    C1=z1-375*cos(a);

    D1=x1-7.9075-48*sin(p1+0.87)+62*sin(p1)+(2^(-1/2))*(265+375*sin(b));
    E1=y1+48.55-48*cos(p1+0.87)+62*cos(p1)-(2^(-1/2))*(285+375*sin(b));
    F1=z1-375*cos(b);

    G1=x1-7.9075-48*sin(p1+0.87)+(2^(-1/2))*(265+375*sin(c));
    H1=y1-48.55-48*cos(p1+0.87)+(2^(-1/2))*(285+375*sin(c));
    K1=z1-375*cos(c);

    L1=x1+122.9075-48*sin(p1+0.87)-(2^(-1/2))*(265+375*sin(d));
    M1=y1-48.55-48*cos(p1+0.87)+(2^(-1/2))*(285+375*sin(d));
    N1=z1-375*cos(d);

    P1=375*(2^(1/2)*A1*cos(a)+2^(1/2)*B1*cos(a)-2*C1*sin(a));
    Q1=375*(-2^(1/2)*D1*cos(b)+2^(1/2)*E1*cos(b)+2*F1*sin(b));
    R1=-375*(2^(1/2)*G1*cos(c)+2^(1/2)*H1*cos(c)+2*K1*sin(c));
    S1=375*(2^(1/2)*L1*cos(d)-2^(1/2)*M1*cos(d)-2*N1*sin(d));

    T1=-(A1*(48*cos(p1+0.87)-62*cos(p1))+B1*(48*sin(p1+0.87))-62*sin(p1));
    U1=-(D1*(48*cos(p1+0.87)-62*cos(p1))+E1*(48*sin(p1+0.87))-62*sin(p1));
    V1=-G1*48*cos(p1+0.87)+H1*48*sin(p1+0.87);
    W1=-L1*48*cos(p1+0.87)+M1*48*sin(p1+0.87);  % A-W的表达式

    x2=x1+2*A1*0.151796/P1+2*B1*0.01542/P1+2*C1*0.256127/P1+2*T1*0.428675/P1;
    y2=y1+2*D1*0.151796/Q1+2*E1*0.01542/Q1+2*F1*0.256127/Q1+2*U1*0.428675/Q1;
    z2=z1+2*G1*0.151796/R1+2*H1*0.01542/R1+2*K1*0.256127/R1+2*V1*0.428675/R1;
    p2=p1+2*L1*0.151796/S1+2*M1*0.01542/S1+2*N1*0.256127/S1+2*W1*0.428675/S1;      %%所求的x y z p的表达式
  
    x1=x2;  %迭代
    y1=y2;
    z1=z2;
    p1=p2;
%     [x(i),y(i),z(i),p(i)]      %%输出x y z p
end
disp('x2=');
disp(x2)
disp('y2=')
disp(y2)
disp('z2=')
disp(z2)
disp('p2=')
disp(p2)

计算结果：
x2=
  196.4266

y2=
  -92.4946

z2=
-619.1816

p2=
   -1.2086

2楼: Originally posted by csl_2012 at 2012-11-15 17:54:28
你好，我刚看了你修改过的程序，你的目的应该是一直x,y,z,p的初值，然后经过迭代(经过A-W的中间处理)，最后得到迭代100次后的x,y,z,p值，如果正常的话，计算结果应该逐渐趋于某一组值，即收敛。
你的程序中a(i-1)= ...

4楼: Originally posted by csl_2012 at 2012-11-15 22:02:46
我计算的就是x100,在迭代过程中，每次x2在x1的基础上计算，你要的迭代100次，我最后输出的就是迭代100次后的结果。x1,x2只是变量而已，你好好看看程序就明白了。

4楼: Originally posted by csl_2012 at 2012-11-15 22:02:46
我计算的就是x100,在迭代过程中，每次x2在x1的基础上计算，你要的迭代100次，我最后输出的就是迭代100次后的结果。x1,x2只是变量而已，你好好看看程序就明白了。

7楼: Originally posted by csl_2012 at 2012-11-15 23:18:41
是的，你可以找本数值分析的书看看，有关于迭代法的详细介绍。