 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 3 | 1/1 | 返回列表 |
| 查看: 744 | 回復(fù): 2 | |||
keekeelock銀蟲 (小有名氣)
|
[求助]
matlab編程循環(huán)語句 跑兩個(gè)理論相等矩陣 求助
|
各位大神求助，用matlab，循環(huán)求解最佳值，初始矩陣通過一定旋轉(zhuǎn)和平移后變成另外一個(gè)結(jié)束矩陣，現(xiàn)在已經(jīng)知道初始的矩陣和結(jié)束矩陣，求解過程的一個(gè)問題。我現(xiàn)在自己的想法是，設(shè)定幾個(gè)變量，構(gòu)造關(guān)系，讓初始矩陣進(jìn)行關(guān)于變量的運(yùn)算，最后形成一個(gè)關(guān)于變量矩陣，通過這兩個(gè)矩陣?yán)碚撓嗟鹊年P(guān)系，選取一個(gè)精度值，小于這個(gè)精度時(shí)跳出這樣子，循環(huán)語句一直有問題，新手求助��！ %a=5,b=20,轉(zhuǎn)動(dòng)角=30,dx=10,dy=10,dz=10 %M1=[75 0 75;0 0 0;200 0 -200;1 1 1;]; clc; M(1,1)=431.38;M(2,1)=449.42;M(3,1)=10.711; M(1,2)=318.18;M(2,2)=531.11;M(3,2)=-150.96; M(1,3)=329.28;M(2,3)=546.19;M(3,3)=-363.74; M(4,1)=1;M(4,2)=1;M(4,3)=1; jingdu=20;gama=97.27; d=gama*pi/180; G=sin(d); H=cos(d); Td=[H -G 0 0;G H 0 0;0 0 1 0;0 0 0 1]; M1=Td*M; N(1,1)=75;N(2,1)=0;N(3,1)=200; N(1,2)=0;N(2,2)=0;N(3,2)=0; N(1,3)=75;N(2,3)=0;N(3,3)=-200; N(4,1)=1;N(4,2)=1;N(4,3)=1; tr1=[1 0 0 0;0 1 0 -11;0 0 1 0;0 0 0 1;]; tr2=[1 0 0 25;0 1 0 -100;0 0 1 200;0 0 0 1;]; tr3=[1 0 0 -210;0 1 0 880;0 0 1 40;0 0 0 1;]; tr4=[1 0 0 185;0 1 0 -160;0 0 1 -400;0 0 0 1;]; delta1=10;delta2=10;delta3=10; dx=0;dy=0;dz=0; for alfa=-10:0.1:10 for bata=-20:0.1:20 a=alfa*pi/180; A=sin(a); B=cos(a); Ta=[1 0 0 0;0 B -A 0;0 A B 0;0 0 0 1]; b=bata*pi/180; C=sin(b); D=cos(b); Tb=[D 0 C 0;0 1 0 0;-C 0 D 0;0 0 0 1]; N1=(tr4*(tr3*(tr2*(Ta*(tr1*(Tb*N)))))); delta11=abs(M1(1,1)-N1(1,1)-(M1(1,2)-N1(1,2)));delta12=abs(M1(1,1)-N1(1,1)-(M1(1,3)-N1(1,3))); delta21=abs(M1(2,1)-N1(2,1)-(M1(2,2)-N1(2,2)));delta22=abs(M1(2,1)-N1(2,1)-(M1(2,3)-N1(2,3))); delta31=abs(M1(3,1)-N1(3,1)-(M1(3,2)-N1(3,2)));delta32=abs(M1(3,1)-N1(3,1)-(M1(3,3)-N1(3,3))); dx=((M1(1,1)-N1(1,1))+(M1(1,2)-N1(1,2))+(M1(1,3)-N1(1,3)))/3; dy=((M1(2,1)-N1(2,1))+(M1(2,2)-N1(2,2))+(M1(1,3)-N1(1,3)))/3; dz=((M1(3,1)-N1(3,1))+(M1(2,2)-N1(2,2))+(M1(2,3)-N1(2,3)))/3; sprintf('%8.4f %8.4f %8.4f %8.4f %8.4f',dx,dy,dz,alfa,bata) i=i+1; if delta11<jingdu&&delta12<jingdu&&delta21<jingdu&&delta22<jingdu&&delta31<jingdu&&delta32<jingdu sprintf('無法輸出解') sprintf('%8.4f %8.4f %8.4f %8.4f %8.4f',dx,dy,dz,alfa,bata) break; end end end |
回復(fù)此樓» 猜你喜歡
287求調(diào)劑
已經(jīng)有5人回復(fù)
-
22408 344分 求調(diào)劑 一志愿 華電計(jì)算機(jī)技術(shù)
已經(jīng)有3人回復(fù)
-
289求調(diào)劑
已經(jīng)有3人回復(fù)
-
317求調(diào)劑
已經(jīng)有10人回復(fù)
-
廣西大學(xué)家禽遺傳育種課題組2026年碩士招生（接收計(jì)算機(jī)專業(yè)調(diào)劑）
已經(jīng)有3人回復(fù)
-
學(xué)校已經(jīng)提交到NSFC，還能修改嗎？
已經(jīng)有8人回復(fù)
-
08工學(xué)調(diào)劑
已經(jīng)有5人回復(fù)
-
281求調(diào)劑（0805）
已經(jīng)有25人回復(fù)
-
085600材料與化工
已經(jīng)有6人回復(fù)
-
265求調(diào)劑
已經(jīng)有8人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助:
- matlab循環(huán)語句怎么讀取excel第二列并存入新excel
已經(jīng)有6人回復(fù)
- matlab如何以循環(huán)語句形式將三維數(shù)組存儲成多個(gè)asc格式的文件
已經(jīng)有4人回復(fù)
- matlab循環(huán)矢量化求助
已經(jīng)有4人回復(fù)
- matlab循環(huán)計(jì)算轉(zhuǎn)為矩陣 30金幣求助
已經(jīng)有7人回復(fù)
- matlab如何實(shí)現(xiàn)矩陣多行數(shù)據(jù)的循環(huán)調(diào)用
已經(jīng)有5人回復(fù)
- matlab中如何通過for語句得到已有矩陣的部分?jǐn)?shù)據(jù)
已經(jīng)有3人回復(fù)
- 求助，matlab公式繪圖，公式里面有兩個(gè)數(shù)組，編出小程序，但不出圖像
已經(jīng)有4人回復(fù)
- 跪求matlab的for循環(huán)語句的編寫
已經(jīng)有20人回復(fù)
- matlab編程大家看這兩個(gè)式子該如何編寫呢？
已經(jīng)有3人回復(fù)
- 【求助】matlab矩陣中尋找距離最近的行向量
已經(jīng)有4人回復(fù)
- 請幫忙用matlab求矩陣特征值
已經(jīng)有8人回復(fù)
- 求助用matlab表示一個(gè)矩陣
已經(jīng)有4人回復(fù)
- 【求助】請教在在matlab里面怎么輸入含參數(shù)矩陣指數(shù)形式-----------【已完結(jié)】
已經(jīng)有5人回復(fù)
csgt0
榮譽(yù)版主 (著名寫手)
彩色掛圖
-
專家經(jīng)驗(yàn): +2 - 應(yīng)助: 367 (碩士)
- 貴賓: 0.457
- 金幣: 7807
- 散金: 8897
- 紅花: 20
- 沙發(fā): 7
- 帖子: 2809
- 在線: 515.8小時(shí)
- 蟲號: 449046
- 注冊: 2007-11-02
- 專業(yè): 金屬非晶態(tài)、準(zhǔn)晶和納米晶
- 管轄: 計(jì)算模擬
2樓2013-12-02 16:56:05
keekeelock
銀蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 193.3
- 帖子: 109
- 在線: 63.9小時(shí)
- 蟲號: 2734574
- 注冊: 2013-10-18
- 性別: GG
- 專業(yè): 零件加工制造
3樓2013-12-02 21:11:59
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 keekeelock 的主題更新
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 289求調(diào)劑 +3 | 懷瑾握瑜l 2026-03-20 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 廣西大學(xué)家禽遺傳育種課題組2026年碩士招生（接收計(jì)算機(jī)專業(yè)調(diào)劑） +3 | 123阿標(biāo) 2026-03-17 | 3/150 |
|
|
|
[考研] 281求調(diào)劑（0805） +14 | 煙汐憶海 2026-03-16 | 25/1250 |
|
|
|
[考研] 298-一志愿中國農(nóng)業(yè)大學(xué)-求調(diào)劑 +9 | 手機(jī)用戶 2026-03-17 | 9/450 |
|
|
|
[考研] 一志愿西南交通 專碩 材料355 本科雙非 求調(diào)劑 +4 | 西南交通專材355 2026-03-19 | 4/200 |
|
|
|
[考研] 085410人工智能專碩317求調(diào)劑（0854都可以） +4 | xbxudjdn 2026-03-18 | 4/200 |
|
|
|
[考研] 329求調(diào)劑 +6 | 想上學(xué)吖吖 2026-03-19 | 6/300 |
|
|
|
[考博] 申博26年 +3 | 八6八68 2026-03-19 | 3/150 |
|
|
|
[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑 +11 | yangfz 2026-03-17 | 11/550 |
|
|
|
[考研] 一志愿 西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。 +3 | 晨昏線與星海 2026-03-19 | 3/150 |
|
|
|
[考研] 0703化學(xué) 305求調(diào)劑 +4 | FY_yy 2026-03-14 | 4/200 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +3 | 妮妮ninicgb 2026-03-17 | 3/150 |
|
|
|
[考博] 26申博 +4 | 八6八68 2026-03-16 | 4/200 |
|
|
|
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6 | 薛云鵬 2026-03-15 | 6/300 |
|
|
|
[考研] 283求調(diào)劑 +3 | 聽風(fēng)就是雨； 2026-03-16 | 3/150 |
|
|
|
[基金申請]
今年的國基金是打分制嗎？
 50+3
|
zhanghaozhu 2026-03-14 | 3/150 |
|
|
|
[考研] 318求調(diào)劑 +3 | Yanyali 2026-03-15 | 3/150 |
|
|
|
[考研] 求老師收留調(diào)劑 +4 | jiang姜66 2026-03-14 | 5/250 |
|
|
|
[考研] 26考研一志愿中國石油大學(xué)(華東)305分求調(diào)劑 +3 | 嘉年新程 2026-03-15 | 3/150 |
|
|
|
[考研] 22408總分284求調(diào)劑 +3 | InAspic 2026-03-13 | 3/150 |
|