通过多个矩阵创建矩阵这些矩阵具有通过方程式的元素

1次查看(最近30天)

Bashar Halleem. 于2020年5月14日

评论： Bashar Halleem.2020年5月21日

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

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

A = [AX1 BX1 CX1 0 0 0 ...;AY1 BY1 CY1 0 0 0 ....;AX2 BX2 CX2 0 0 0 ...;AY2 BY2 CY2 0 0 0 ...;........; AX12 BX12 CX12 0 0 0 ...;AY12 BY12 CY12 0 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）零zeros

Zeros Zeros矩阵（12 * 3）零..

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

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

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

谢谢你的帮助

谢谢你的任何想法

Bashar.

5点评论
显示隐藏＆nbsp4旧的评论

Bashar Halleem. 2020年5月21日

                          IWInt = EIP (17);%input('图像宽度在EIPxel');
                         
                          Fwint = EIP（1,8）;％输入（'mm'的格式宽度）;
                         
                          IHInt = EIP (18);%input('图像高度在EIPxel');
                         
                          fhint = EIP（1,9）;％输入（'mm'的格式高度）;
                         
                          ppx = eip（1,10）;％输入（'主点x（mm）'）;
                         
                          ppy = EIP (11);％输入（'主点Y（mm）'）;
                         
                          EIP (12);%的输入(k1);
                         
                          EIP（1,13）;％输入（'K2'）;
                         
                          EIP（1,14）;％输入（'k3'）;
                         
                          EIP（1,15）;%的输入(p1);
                         
                          EIP（1,16）;%的输入(p2);
                         
                          EIP (1,1);％输入（'Center X（项目单位）或（相机站XO'）;
                         
                          EIP（1,2）;%input('Center Y(项目单位)Or (Camera Station Yo');
                         
                          EIP（1,3）;％输入（'Center Z（项目单位）或（相机站ZO'）;
                         
                          p3 = nan（1440,1）;
                         
                          为了我= 1：1440
                         
                          PxP = Pxl(我,1);%input('x of point in EIPxel');
                         
                          pxmm =（（fwint * pxp）/ iwint）-ppx;
                         
                          P3(我,1)= pxmm;
                         
                          结尾
                         
                          p4 = nan（1440,1）;
                         
                          为了我= 1：1440
                         
                          PYP = PXL（I，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);
                         
                          为了我= 1：1440
                         
                          R1（I，1）=（p3（i，1）^ 2 + p4（i，1）^ 2）^ 0.5;
                         
                          SX1（I，1）= P3（I，1）*（EIP（1,12）+ EIP（1,13）* R1（I，1）^ 2 + EIP（1,14）* R1（I，1）^ 4）;
                         
                          sy1(我,1)= P4(我,1)* (EIP(1、12)+ EIP (13) * r1 (1, 1) ^ 2 + EIP (13) * r1(我,1)^ 4);
                         
                          DX1（I，1）= EIP（1,15）*（R1（I，1）^ 2 + 2 * P3（I，1）^ 2）+ 2 * EIP（1,16）* P3（I，1）* p4（i，1）;
                         
                          ₁(1)= EIP (16) * (r1(我,1)^ 2 + 2 * P3(我,1)^ 2)+ 2 * EIP (15) * P3(我,1)* P4(我,1);
                         
                          XMM1（I，1）= P3（I，1）+ SX1（I，1）+ DX1（I，1）;
                         
                          ymm1(我,1)= P4(我,1)+ sy1(我,1)+₁(我,1);
                         
                          %%%%%%%%%%%%%%%%%%%%%%%%%数值remix %%%%%%%%%%%%%%%%%%%%%%%%%%
                         
                          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）* sind（phi）* cosd（ka）+ sind（w）* sind（ka）;
                         
                          m (i + 1) = -cosd(φ)*信德(ka);
                         
                          m（i + 1,2）=  -  Sind（w）* sind（phi）* sind（ka）+ cosd（w）* cosd（ka）;
                         
                          m（i + 1,3）= Cosd（w）* sind（phi）* sind（ka）+ sind（w）* cosd（ka）;
                         
                          m (i + 2, 1) =信德(φ);
                         
                          m（i + 2,2）=  -  sind（w）* cosd（phi）;
                         
                          m (i + 2、3)= cosd (w) * cosd(φ);
                         
                          ％m = nan（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);％compute参数
                         
                          Ayi = EIP（I，7）* M（I + 1,1）+ YMM1（I，1）* M（i + 2,1）;BYI = EIP（i，7）* m（i + 1,2）+ ymm1（i，1）* m（i + 2,2）;Cyi = EIP（I，7）* M（I + 1,3）+ YMM1（I，1）* M（i + 2,3）;
                         
                          LXI = AXI * EIP（I，1）+ BXI * EIP（I，2）+ CXI * EIP（I，3）;
                         
                          Lyi = Ayi * EIP（I，1）+ BYI * EIP（I，2）+ CYI * EIP（I，3）;
                         
                          结尾

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

Bashar Halleem. 2020年5月21日

                          一种矩阵[AX1 BX1 CX1 0 0 0
                         
                          AY1By1 Cy1
                         
                          。。
                         
                          AX12BX12 CX12.
                         
                          AY12By12 Cy120.0.0.
                         
                          0 0 0 Ax13Bx13 Cx13
                         
                          Ay13截至13 Cy13.
                         
                          。
                         
                          AX24Bx24 Cx24
                         
                          0 0 0 Ay24BY24 CY24
                         
                          。
                         
                          。
                         
                          。
                         
                          。
                         
                          。
                         
                          Ax108Bx108 Cx108
                         
                          Ay108By108 Cy108
                         
                          。
                         
                          Ax120BX120 CX120
                         
                          Ay120By120 Cy120]

登录评论。

登录来回答这个问题。

答案（1）

斯蒂芬 2020年5月20日

0.
关联

直接链接到这个答案

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

编辑:斯蒂芬 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点评论
显示隐藏＆nbsp4旧的评论

答案（1）

0评论
显示隐藏-1旧的评论

也可以看看

类别

标签

社区寻宝

你现在关注这个问题

通过多个矩阵创建矩阵这些矩阵具有通过方程式的元素

5点评论 显示隐藏＆nbsp4旧的评论

答案（1）

0评论 显示隐藏-1旧的评论

也可以看看

类别

标签

社区寻宝

发生错误

5点评论
显示隐藏＆nbsp4旧的评论

0评论
显示隐藏-1旧的评论