这些矩阵有元素由方程

1次查看(最近30天)

巴沙尔halleem 2020年5月14日

评论道: 巴沙尔halleem2020年5月21日

这些矩阵有元素通过方程:

矩阵A = 0(1728,432)，包含矩阵(12,3)作为对角线每一个元素的方程。

A = [Ax1 Bx1 Cx1 0 0…];Ay1 By1 Cy1 0 0 ....;Ax2 Bx2 Cx2 0 0 0…;Ay2 By2 Cy2 0 0…;........;Ax12 Bx12 Cx12 0 0…;Ay12 By12 Cy12 0 0…;0 0 0 Ax1 Bx1 Cx1 0 0 0…0 0 0 Ay1 By1 Cy1 0 0 0…］.

换句话说，矩阵A =[矩阵(12*3)零零零

矩阵(12*3)

矩阵(12*3)零…

... ... ... ...］

ax₁,Bx1系列,…每个元素都等于等式，它们的值在288个输入值时不同。

如果有人能安排矩阵和主矩阵A的循环。

谢谢你的帮助

谢谢你的想法

巴沙尔

5个评论
显示隐藏 4年长的评论

巴沙尔halleem 2020年5月21日

                          IWInt = EIP (17);%input('图像宽度在EIPxel');
                         
                          FWInt = EIP (1,8);%input('Format Width in mm');
                         
                          IHInt = EIP (18);%input('图像高度在EIPxel');
                         
                          FHInt = EIP (9);%input('Format Height in mm');
                         
                          PPx = EIP (10);%input('Principal Point X (mm)');
                         
                          ppy = EIP (11);%input('主点y (mm)');
                         
                          EIP (12);%的输入(k1);
                         
                          EIP (13);%的输入(k2);
                         
                          EIP (14);%的输入(“k3”);
                         
                          EIP (15);%的输入(p1);
                         
                          EIP (16);%的输入(p2);
                         
                          EIP (1,1);%input('Center X(项目单位)Or (Camera Station Xo');
                         
                          EIP(1、2);%input('Center Y(项目单位)Or (Camera Station Yo');
                         
                          EIP(1、3);%input('Center Z(项目单位)Or (Camera Station Zo');
                         
                          P3 =南(1440 1);
                         
                          为i = 1:1440
                         
                          PxP = Pxl(我,1);%input('x of point in EIPxel');
                         
                          pxmm = ((FWInt * PxP) / IWInt) -PPx;
                         
                          P3(我,1)= pxmm;
                         
                          结束
                         
                          P4 =南(1440 1);
                         
                          为i = 1:1440
                         
                          PyP = Pxl(我,2);%input('y of point in EIPxel');for j=1:2:24;k = 2:2:24;因为我= 1:1440
                         
                          pymm = ppy - ((FHInt * PyP) / IHInt);
                         
                          P4(我,1)= pymm;
                         
                          结束
                         
                          r1 =南(1440 1),这位朋友=南(1440 1);sy1 =南(1440 1);dx1 =南(1440 1);₁=南(1440 1);xmm1 =南(1440 1);ymm1 =南(1440 1);
                         
                          为i = 1:1440
                         
                          r1(我,1)= (P3(我,1)^ 2 + P4(我,1)^ 2)^ 0.5;
                         
                          朋友说(我,1)= P3(我,1)* (EIP(1、12)+ EIP (13) * r1(我,1)^ 2 + EIP (14) * r1(我,1)^ 4);
                         
                          sy1(我,1)= P4(我,1)* (EIP(1、12)+ EIP (13) * r1 (1, 1) ^ 2 + EIP (13) * r1(我,1)^ 4);
                         
                          dx1 (1) = EIP (15) * (r1(我,1)^ 2 + 2 * P3(我,1)^ 2)+ 2 * EIP (16) * P3(我,1)* P4(我,1);
                         
                          ₁(1)= EIP (16) * (r1(我,1)^ 2 + 2 * P3(我,1)^ 2)+ 2 * EIP (15) * P3(我,1)* P4(我,1);
                         
                          xmm1(我,1)= P3(我,1)+朋友(我,1)+ dx1(我,1);
                         
                          ymm1(我,1)= P4(我,1)+ sy1(我,1)+₁(我,1);
                         
                          %%%%%%%%%%%%%%%%%%%%%%%%% 计算旋转矩阵 %%%%%%%%%%%%%%%%%%%%%%%%%%
                         
                          w = EIP(我,4);φ= EIP(我,5);ka = EIP(我,6);
                         
                          (我,1)= cosd(φ)* cosd (ka);
                         
                          m (i, 2) =信德(w) *信德(φ)* cosd (ka) + cosd (w) *信德(ka);
                         
                          m (i, 3) = -cosd (w) *信德(φ)* cosd (ka) +信德(w) *信德(ka);
                         
                          m (i + 1) = -cosd(φ)*信德(ka);
                         
                          m (i + 1、2)=信德(w) *信德(φ)*信德(ka) + cosd (w) * cosd (ka);
                         
                          m (i + 1, 3) = cosd (w) *信德(φ)*信德(ka) +信德(w) * cosd (ka);
                         
                          m (i + 2, 1) =信德(φ);
                         
                          m (i + 2 2) =信德(w) * cosd(φ);
                         
                          m (i + 2、3)= cosd (w) * cosd(φ);
                         
                          % M =南(4320 3);
                         
                          M = [M (i, 1) M (i, 2) M (i, 3); M (i + 1) M (i + 1, 2) M (i + 1,3); M (i + 2, 1) M (i + 2, 2) M (i + 2,3)];
                         
                          阿喜= EIP(我,7)* m (i, 1) + xmm1(我,1)* m (i + 2, 1);Bxi = EIP(我,7)* m (i, 2) + xmm1(我,1)* m (i + 2 2);Cxi = EIP(我,7)* m (i, 3) + xmm1(我,1)* m (i + 2、3);%计算参数
                         
                          我阿姨= EIP (7) * m (i + 1, - 1) + ymm1(我,1)* m (i + 2, 1);Byi = EIP(我,7)* m (i + 1、2)+ ymm1(我,1)* m (i + 2 2);Cyi = EIP(我,7)* m (i + 1, 3) + ymm1(我,1)* m (i + 2、3);
                         
                          Lxi = Axi * EIP(我,1)+ Bxi * EIP(我,2)+ Cxi * EIP(我,3);
                         
                          Lyi =阿姨* EIP(我,1)+ Byi * EIP(我,2)+ Cyi * EIP(我,3);
                         
                          结束

先生，现在我要把结果排列成特定的矩阵

巴沙尔halleem 2020年5月21日

                          一个矩阵[Ax1 Bx1 Cx1 0 0 0 .
                         
                          Ay1By1 Cy1
                         
                          . . . . . .
                         
                          Ax12Bx12 Cx12
                         
                          Ay12By12 Cy12000
                         
                          0 0 0 Ax13Bx13 Cx13
                         
                          Ay13By13 Cy13
                         
                          ． . .
                         
                          Ax24Bx24 Cx24
                         
                          0 0 0 Ay24By24 Cy24
                         
                          ．
                         
                          ．
                         
                          ．
                         
                          ．
                         
                          ．
                         
                          Ax108Bx108 Cx108
                         
                          Ay108By108 Cy108
                         
                          ． . .
                         
                          Ax120Bx120 Cx120
                         
                          Ay120By120 Cy120]

登录评论。

登录回答这个问题。

答案(1)

Stephen Cobeldick 2020年5月20日

0
链接

直接链接到这个答案

//www.tatmou.com/matlabcentral/answers/525456-create-matrix-by-several-matrices-these-matrices-have-elements-by-equations#answer_434054

编辑:Stephen Cobeldick 2020年5月20日

拥有大量的编号变量并不是一个好的方法，你应该在创建这些子矩阵时将它们放到一个单元格数组中(使用索引是很简单的)，例如:

C = {M1, M2,．..}；

在哪里

M1 = [ax₁,Bx1系列,Cx1; Ay1 By1, Cy1,．..]；

等。当然，当你创建它们时，你应该使用索引，而不是使用编号的变量名，

那你只需要一个简单的 blkdiag 电话:

B = blkdiag (C {:});

如果实际上只有一个子矩阵重复(如你的问题所示)，那么你可以这样做:

                        C = {M};
                       
                        B = blkdiag (C {1 (1144)});

0评论
显示隐藏1以上的评论

登录评论。

登录回答这个问题。

你现在正关注这个问题

这些矩阵有元素由方程

5个评论
显示隐藏 4年长的评论

答案(1)

0评论
显示隐藏1以上的评论

另请参阅

类别

标签

社区寻宝

标签

另请参阅

社区寻宝

你现在正关注这个问题

这些矩阵有元素由方程

5个评论 显示隐藏 4年长的评论

答案(1)

0评论 显示隐藏1以上的评论

另请参阅

类别

标签

社区寻宝

标签

另请参阅

社区寻宝

一个错误发生

5个评论
显示隐藏 4年长的评论

0评论
显示隐藏1以上的评论