24小时热门版块排行榜

返回列表

Kevin1987s

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 7
帖子: 15
在线: 7.3小时
虫号: 1405380
注册: 2011-09-17
专业: 机械动力学

[求助] 求大神帮我改一下下面的程序！matlab迭代法的不知道错在哪里！

x(1)=210;
y(1)=-90;
z(1)=-620;
p(1)=0.2988;
a(1)=-0.1096;
b(1)=0.5996;
c(1)=0.3507;
d(1)=-0.3882;       %%设定初值
for    i=2:101;
a(i-1)=-0.1096+(i-1)*0.002648;
b(i-1)=0.5996+(i-1)*0.000269;
c(i-1)=0.3507+(i-1)*0.004468;
d(i-1)=-0.3882+(i-1)*0.007478;       %%第（i-1）的角度值

A(i-1)=x(i-1)+122.9075-48*sin(p(i-1)+0.87)+62*sin(p(i-1))-(2^(-1/2))*(265+375*sin(a(i-1)));
B(i-1)=y(i-1)+48.55-48*cos(p(i-1)+0.87)+62*cos(p(i-1))-(2^(-1/2))*(285+375*sin(a(i-1)));
C(i-1)=z(i-1)-375*cos(a(i-1));

D(i-1)=x(i-1)-7.9075-48*sin(p(i-1)+0.87)+62*sin(p(i-1))+(2^(-1/2))*(265+375*sin(b(i-1)));
E(i-1)=y(i-1)+48.55-48*cos(p(i-1)+0.87)+62*cos(p(i-1))-(2^(-1/2))*(285+375*sin(b(i-1)));
F(i-1)=z(i-1)-375*cos(b(i-1));

G(i-1)=x(i-1)-7.9075-48*sin(p(i-1)+0.87)+(2^(-1/2))*(265+375*sin(c(i-1)));
H(i-1)=y(i-1)-48.55-48*cos(p(i-1)+0.87)+(2^(-1/2))*(285+375*sin(c(i-1)));
K(i-1)=z(i-1)-375*cos(c(i-1));

L(i-1)=x(i-1)+122.9075-48*sin(p(i-1)+0.87)-(2^(-1/2))*(265+375*sin(d(i-1)));
M(i-1)=y(i-1)-48.55-48*cos(p(i-1)+0.87)+(2^(-1/2))*(285+375*sin(d(i-1)));
N(i-1)=z(i-1)-375*cos(d(i-1));

P(i-1)=375*(2^(1/2)*A(i-1)*cos(a(i-1))+2^(1/2)*B(i-1)*cos(a(i-1))-2*C(i-1)*sin(a(i-1)));
Q(i-1)=375*(-2^(1/2)*D(i-1)*cos(b(i-1))+2^(1/2)*E(i-1)*cos(b(i-1))+2*F(i-1)*sin(b(i-1)));
R(i-1)=-375*(2^(1/2)*G(i-1)*cos(c(i-1))+2^(1/2)*H(i-1)*cos(c(i-1))+2*K(i-1)*sin(c(i-1)));
S(i-1)=375*(2^(1/2)*L(i-1)*cos(d(i-1))-2^(1/2)*M(i-1)*cos(d(i-1))-2*N(i-1)*sin(d(i-1)));

T(i-1)=-(A(i-1)*(48*cos(p(i-1)+0.87)-62*cos(p(i-1)))+B(i-1)*(48*sin(p(i-1)+0.87))-62*sin(p(i-1)));
U(i-1)=-(D(i-1)*(48*cos(p(i-1)+0.87)-62*cos(p(i-1)))+E(i-1)*(48*sin(p(i-1)+0.87))-62*sin(p(i-1)));
V(i-1)=-G(i-1)*48*cos(p(i-1)+0.87)+H(i-1)*48*sin(p(i-1)+0.87);
W(i-1)=-L(i-1)*48*cos(p(i-1)+0.87)+M(i-1)*48*sin(p(i-1)+0.87);  && A-W的表达式

x(i)=x(i-1)+2*A(i-1)*0.151796/P(i-1)+2*B(i-1)*0.01542/P(i-1)+2*C(i-1)*0.256127/P(i-1)+2*T(i-1)*0.428675/P(i-1);
y(i)=y(i-1)+2*D(i-1)*0.151796/Q(i-1)+2*E(i-1)*0.01542/Q(i-1)+2*F(i-1)*0.256127/Q(i-1)+2*U(i-1)*0.428675/Q(i-1);
z(i)=z(i-1)+2*G(i-1)*0.151796/R(i-1)+2*H(i-1)*0.01542/R(i-1)+2*K(i-1)*0.256127/R(i-1)+2*V(i-1)*0.428675/R(i-1);
p(i)=p(i-1)+2*L(i-1)*0.151796/S(i-1)+2*M(i-1)*0.01542/S(i-1)+2*N(i-1)*0.256127/S(i-1)+2*W(i-1)*0.428675/S(i-1);    %%所求的x y z p的表达式

[x(i),y(i),z(i),p(i)]    %%输出x y z p

end

x y z p初值已知，迭代100次上面这个程序运算出来貌似不是我要的结果。。。迭代循环是不是有问题啊？谢谢大神！

[ Last edited by Kevin1987s on 2012-11-15 at 11:03 ]

回复此楼

» 猜你喜欢

职称评审没过，求安慰已经有34人回复
垃圾破二本职称评审标准已经有17人回复
回收溶剂求助已经有6人回复
投稿Elsevier的Neoplasia杂志，到最后选publishing options时页面空白，不能完成投稿已经有22人回复
申请26博士已经有5人回复
EST投稿状态问题已经有7人回复
毕业后当辅导员了，天天各种学生超烦已经有4人回复
聘U V热熔胶研究人员已经有10人回复
求助文献已经有3人回复
投稿返修后收到这样的回复，还有希望吗已经有8人回复

» 本主题相关价值贴推荐，对您同样有帮助:

1楼 2012-11-15 10:21:40

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

csl_2012

捐助贵宾 (初入文坛)

应助: 6 (幼儿园)
金币: 838.3
散金: 300
红花: 1
帖子: 50
在线: 113小时
虫号: 1734989
注册: 2012-04-03
性别: GG
专业: 动力学与控制

【答案】应助回帖

★
感谢参与，应助指数 +1
xzhdty: 金币+1, 谢谢参与 2012-11-15 21:03:40
jjdg: 编辑内容 2012-11-16 10:05

你好，我刚看了你修改过的程序，你的目的应该是一直x,y,z,p的初值，然后经过迭代(经过A-W的中间处理)，最后得到迭代100次后的x,y,z,p值，如果正常的话，计算结果应该逐渐趋于某一组值，即收敛。
你的程序中a(i-1)=-0.1096+(i-1)*0.002648是不是应该修改为：
a(i-1)=-0.1096+(i-2)*0.002648，因为i=2时，a(1)=-0.1096.
修改后的程序如下：

CODE:

%迭代法

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

希望对你有用。

[ Last edited by jjdg on 2012-11-16 at 10:05 ]

赞一下(1人)

回复此楼

2楼2012-11-15 17:54:28

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

Kevin1987s

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 7
帖子: 15
在线: 7.3小时
虫号: 1405380
注册: 2011-09-17
专业: 机械动力学

引用回帖:

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)= ...

太感谢您了！可以求disp(x2)，怎么求disp(x100)呢？把里面的2改了就算不出来了，小弟新手，刚用ml，谢谢指点！

[ 发自手机版 http://muchong.com/3g ]

赞一下

回复此楼

3楼2012-11-15 19:47:03

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖