| 查看: 808 | 回复: 2 | ||
keekeelock银虫 (小有名气)
|
[求助]
matlab编程循环语句 跑两个理论相等矩阵 求助
|
各位大神求助，用matlab，循环求解最佳值，初始矩阵通过一定旋转和平移后变成另外一个结束矩阵，现在已经知道初始的矩阵和结束矩阵，求解过程的一个问题。我现在自己的想法是，设定几个变量，构造关系，让初始矩阵进行关于变量的运算，最后形成一个关于变量矩阵，通过这两个矩阵理论相等的关系，选取一个精度值，小于这个精度时跳出这样子，循环语句一直有问题，新手求助！！ %a=5,b=20,转动角=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 |
回复此楼» 猜你喜欢
论文撤稿了
已经有8人回复
-
博士申请
已经有7人回复
-
青B发送上会通知了吗
已经有8人回复
-
化学专业申博
已经有4人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有3人回复
-
评审有感
已经有27人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab循环语句怎么读取excel第二列并存入新excel
已经有6人回复
- matlab如何以循环语句形式将三维数组存储成多个asc格式的文件
已经有4人回复
- matlab循环矢量化求助
已经有4人回复
- matlab循环计算转为矩阵 30金币求助
已经有7人回复
- matlab如何实现矩阵多行数据的循环调用
已经有5人回复
- matlab中如何通过for语句得到已有矩阵的部分数据
已经有3人回复
- 求助，matlab公式绘图，公式里面有两个数组，编出小程序，但不出图像
已经有4人回复
- 跪求matlab的for循环语句的编写
已经有20人回复
- matlab编程大家看这两个式子该如何编写呢？
已经有3人回复
- 【求助】matlab矩阵中寻找距离最近的行向量
已经有4人回复
- 请帮忙用matlab求矩阵特征值
已经有8人回复
- 求助用matlab表示一个矩阵
已经有4人回复
- 【求助】请教在在matlab里面怎么输入含参数矩阵指数形式-----------【已完结】
已经有5人回复
2楼2013-12-02 16:56:05
keekeelock
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 193.3
- 帖子: 109
- 在线: 63.9小时
- 虫号: 2734574
- 注册: 2013-10-18
- 性别: GG
- 专业: 零件加工制造
3楼2013-12-02 21:11:59
10