版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

Kevin1987s

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 7
帖子: 15
在線: 7.3小時(shí)
蟲號(hào): 1405380
注冊(cè): 2011-09-17
專業(yè): 機(jī)械動(dòng)力學(xué)

[求助] 求大神幫我改一下下面的程序！matlab迭代法的不知道錯(cuò)在哪里！

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;       %%設(shè)定初值
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的表達(dá)式

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的表達(dá)式

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

end

x y z p初值已知，迭代100次上面這個(gè)程序運(yùn)算出來(lái)貌似不是我要的結(jié)果。。。迭代循環(huán)是不是有問(wèn)題�。恐x謝大神！

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

回復(fù)此樓

» 猜你喜歡

085601材料工程找調(diào)劑已經(jīng)有11人回復(fù)
總分293求調(diào)劑已經(jīng)有9人回復(fù)
334分一志愿武理材料求調(diào)劑已經(jīng)有5人回復(fù)
327求調(diào)劑已經(jīng)有6人回復(fù)
329求調(diào)劑已經(jīng)有10人回復(fù)
275求調(diào)劑已經(jīng)有15人回復(fù)
291求調(diào)劑已經(jīng)有6人回復(fù)
309求調(diào)劑已經(jīng)有7人回復(fù)
297求調(diào)劑已經(jīng)有11人回復(fù)
一志愿南航 335分 | 0856 | GPA 4.07 | 有科研經(jīng)歷已經(jīng)有6人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

求助��！求大神幫忙分析一下鎂氧水泥的xrd圖，萬(wàn)分感謝！已經(jīng)有0人回復(fù)
投期刊居然讓我格式按照手冊(cè)，怎么辦，求助已經(jīng)有5人回復(fù)
求大神們幫我翻譯一下下面的內(nèi)容，謝謝了！已經(jīng)有2人回復(fù)
求大神幫我改一下下面的程序！matlab迭代法的已經(jīng)有2人回復(fù)

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

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

csl_2012

捐助貴賓 (初入文壇)

應(yīng)助: 6 (幼兒園)
金幣: 838.3
散金: 300
紅花: 1
帖子: 50
在線: 113小時(shí)
蟲號(hào): 1734989
注冊(cè): 2012-04-03
性別: GG
專業(yè): 動(dòng)力學(xué)與控制

【答案】應(yīng)助回帖

★
感謝參與，應(yīng)助指數(shù) +1
xzhdty: 金幣+1, 謝謝參與 2012-11-15 21:03:40
jjdg: 編輯內(nèi)容 2012-11-16 10:05

你好，我剛看了你修改過(guò)的程序，你的目的應(yīng)該是一直x,y,z,p的初值，然后經(jīng)過(guò)迭代(經(jīng)過(guò)A-W的中間處理)，最后得到迭代100次后的x,y,z,p值，如果正常的話，計(jì)算結(jié)果應(yīng)該逐漸趨于某一組值，即收斂。
你的程序中a(i-1)=-0.1096+(i-1)*0.002648是不是應(yīng)該修改為：
a(i-1)=-0.1096+(i-2)*0.002648，因?yàn)閕=2時(shí)，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;          %%設(shè)定初值

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的表達(dá)式

    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的表達(dá)式

    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)

計(jì)算結(jié)果：

x2=

  196.4266

y2=

  -92.4946

z2=

 -619.1816

p2=

   -1.2086

希望對(duì)你有用。

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

贊一下(1人)

回復(fù)此樓

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

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

Kevin1987s

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 7
帖子: 15
在線: 7.3小時(shí)
蟲號(hào): 1405380
注冊(cè): 2011-09-17
專業(yè): 機(jī)械動(dòng)力學(xué)

引用回帖:

2樓: Originally posted by csl_2012 at 2012-11-15 17:54:28
你好，我剛看了你修改過(guò)的程序，你的目的應(yīng)該是一直x,y,z,p的初值，然后經(jīng)過(guò)迭代(經(jīng)過(guò)A-W的中間處理)，最后得到迭代100次后的x,y,z,p值，如果正常的話，計(jì)算結(jié)果應(yīng)該逐漸趨于某一組值，即收斂。
你的程序中a(i-1)= ...

太感謝您了！可以求disp(x2)，怎么求disp(x100)呢？把里面的2改了就算不出來(lái)了，小弟新手，剛用ml，謝謝指點(diǎn)！

[ 發(fā)自手機(jī)版 http://www.gaoyang168.com/3g ]

贊一下

回復(fù)此樓

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

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

csl_2012

捐助貴賓 (初入文壇)

應(yīng)助: 6 (幼兒園)
金幣: 838.3
散金: 300
紅花: 1
帖子: 50
在線: 113小時(shí)
蟲號(hào): 1734989
注冊(cè): 2012-04-03
性別: GG
專業(yè): 動(dòng)力學(xué)與控制

★
jjdg: 金幣+1, 感謝參與 2012-11-16 10:05:13

我計(jì)算的就是x100,在迭代過(guò)程中，每次x2在x1的基礎(chǔ)上計(jì)算，你要的迭代100次，我最后輸出的就是迭代100次后的結(jié)果。x1,x2只是變量而已，你好好看看程序就明白了。

贊一下

回復(fù)此樓

4樓2012-11-15 22:02:46

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

Kevin1987s

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 7
帖子: 15
在線: 7.3小時(shí)
蟲號(hào): 1405380
注冊(cè): 2011-09-17
專業(yè): 機(jī)械動(dòng)力學(xué)

引用回帖:

4樓: Originally posted by csl_2012 at 2012-11-15 22:02:46
我計(jì)算的就是x100,在迭代過(guò)程中，每次x2在x1的基礎(chǔ)上計(jì)算，你要的迭代100次，我最后輸出的就是迭代100次后的結(jié)果。x1,x2只是變量而已，你好好看看程序就明白了。

謝謝您！那就是這個(gè)方法不能得到我預(yù)計(jì)的結(jié)果。。。還得換個(gè)思路。。。

[ 發(fā)自手機(jī)版 http://www.gaoyang168.com/3g ]

贊一下

回復(fù)此樓

5樓2012-11-15 22:38:14

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

Kevin1987s

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 7
帖子: 15
在線: 7.3小時(shí)
蟲號(hào): 1405380
注冊(cè): 2011-09-17
專業(yè): 機(jī)械動(dòng)力學(xué)

引用回帖:

我想問(wèn)一下，您幫我改的這個(gè)程序是不是類似于 x2=x1+t1(根隨i變化的方程)，x3=x2+t2...
x100=x99+t99, 是這個(gè)意思么？

[ 發(fā)自手機(jī)版 http://www.gaoyang168.com/3g ]

贊一下

回復(fù)此樓

6樓2012-11-15 22:58:53

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

csl_2012

捐助貴賓 (初入文壇)

應(yīng)助: 6 (幼兒園)
金幣: 838.3
散金: 300
紅花: 1
帖子: 50
在線: 113小時(shí)
蟲號(hào): 1734989
注冊(cè): 2012-04-03
性別: GG
專業(yè): 動(dòng)力學(xué)與控制

★
jjdg: 金幣+1, 感謝參與 2012-11-16 10:05:00

是的，你可以找本數(shù)值分析的書看看，有關(guān)于迭代法的詳細(xì)介紹。

贊一下

回復(fù)此樓

7樓2012-11-15 23:18:41

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

Kevin1987s

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 7
帖子: 15
在線: 7.3小時(shí)
蟲號(hào): 1405380
注冊(cè): 2011-09-17
專業(yè): 機(jī)械動(dòng)力學(xué)

引用回帖:

7樓: Originally posted by csl_2012 at 2012-11-15 23:18:41
是的，你可以找本數(shù)值分析的書看看，有關(guān)于迭代法的詳細(xì)介紹。

好的，太感謝了！嘻嘻~

[ 發(fā)自手機(jī)版 http://www.gaoyang168.com/3g ]

贊一下

回復(fù)此樓

8樓2012-11-16 00:10:17

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 Kevin1987s 的主題更新

返回列表

普通表情龍兔虎貓

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 一志愿：西北大學(xué)，英一數(shù)一408-284分求調(diào)劑 +3	12.27 2026-03-27	3/150	2026-03-29 12:53 by 熱情沙漠
[考研] 求收留 +5	1943443204 2026-03-28	5/250	2026-03-29 12:11 by 無(wú)際的草原
[考研] 各位老師好，我的一志愿為北京科技大學(xué)085601材料專碩 +8	Koxui 2026-03-28	8/400	2026-03-29 09:50 by laoshidan
[碩博家園] 求調(diào)劑有機(jī)化學(xué)考研356分 +3	Nadiums 2026-03-25	4/200	2026-03-29 00:50 by 544594351
[考研] 生物學(xué)學(xué)碩，一志愿湖南大學(xué)，初試成績(jī)338 +6	YYYYYNNNNN 2026-03-26	7/350	2026-03-28 20:52 by 唐沐兒
[考研] 304求調(diào)劑 +6	曼殊2266 2026-03-27	6/300	2026-03-28 14:10 by 唐沐兒
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +8	邱gl 2026-03-27	8/400	2026-03-28 12:42 by 唐沐兒
[考研] 0856，材料與化工321分求調(diào)劑 +12	大饞小子 2026-03-27	13/650	2026-03-28 10:56 by self2008
[考研] 274求調(diào)劑 +17	顧九笙要謙虛 2026-03-24	23/1150	2026-03-27 15:16 by caszguilin
[考研] 一志愿陜師大生物學(xué)071000，298分，求調(diào)劑 +5	SYA！ 2026-03-23	5/250	2026-03-27 09:29 by 不吃魚的貓
[考研] 284求調(diào)劑 +11	junqihahaha 2026-03-26	12/600	2026-03-27 04:37 by wxiongid
[考研] 341求調(diào)劑 +7	青檸檬1 2026-03-26	7/350	2026-03-27 00:19 by wxiongid
[考研] 327求調(diào)劑 +7	prayer13 2026-03-23	7/350	2026-03-26 20:48 by 不吃魚的貓
[考研] 材料調(diào)劑 5+4	想要一壺桃花水 2026-03-25	10/500	2026-03-26 19:56 by 不吃魚的貓
[考研] 263求調(diào)劑 +6	yqdszhdap－ 2026-03-22	10/500	2026-03-26 13:11 by 公瑾逍遙
[考研] 303求調(diào)劑 +6	藍(lán)山月 2026-03-25	6/300	2026-03-25 22:47 by 418490947
[考研] 26考研-291分-廈門大學(xué)（085601）-柔性電子學(xué)院材料工程專業(yè)求調(diào)劑 +3	min3 2026-03-24	4/200	2026-03-25 18:22 by xcjcqu
[考研] 0854電子信息求調(diào)劑 +7	α____ 2026-03-22	9/450	2026-03-25 13:37 by α____
[考研] 0854電子信息求調(diào)劑 324 +4	Promise-jyl 2026-03-23	4/200	2026-03-25 11:36 by Sugarlight
[考研] 086003食品工程求調(diào)劑 +6	淼淼111 2026-03-24	6/300	2026-03-25 10:29 by 3Strings

24小時(shí)熱門版塊排行榜

[求助] 求大神幫我改一下下面的程序！matlab迭代法的 不知道錯(cuò)在哪里！

» 猜你喜歡

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

【答案】應(yīng)助回帖

[求助] 求大神幫我改一下下面的程序！matlab迭代法的不知道錯(cuò)在哪里！