乘以两个矩阵

2视图(30天)

显示旧的评论

shiv白肢野牛 2021年9月26日

0
链接

这个问题直接联系

https://de.mathworks.com/matlabcentral/answers/1461104-multiply-by-two-matrix

回答: 沃尔特·罗伯森 2021年9月27日

                          乔丹= [cos (kj hj)我* sin (kj hj) / kj;我* kj * sin (kj hj)因为(kj hj)];
                         
                          kj =√nj ^ 2 * k02-x ^ 2);

取

n1 = 1.521; n2 = 2.66; k0 = 1; h1 = 1.5 e-6; h2 = 1 e-6

M1乘*平方米

4评论
显示隐藏 3年长的评论

shiv白肢野牛 2021年9月26日

编辑:沃尔特·罗伯森 2021年9月27日

                                %这是程序
                               
                                函数kp1
                               
                                关闭所有;
                               
                                清晰的所有;
                               
                                函数y = F (x)
                               
                                为t2 = 1
                               
                                k0 = 1;
                               
                                n1 = 1.521;
                               
                                % n2 = 2.66;
                               
                                n2 = 4.1 - -0.211 * 1我;
                               
                                ns = 1.512;
                               
                                % 1数控= 0.15 * 3.2;
                               
                                数控= 1;
                               
                                k1 =√n1。^ 2 * k0。^ 2 x ^ 2);
                               
                                k2 =√n2。^ 2 * k0。^ 2 x ^ 2);
                               
                                t1 = 1.5;
                               
                                m11公路= cos (t1 * k1) * cos (t2 * k2) - (k2 / k1) * sin (t1 * k1) *罪(t2 * k2);
                               
                                m12 = (1 / k2) * (cos (t1 * k1) * sin (t2 * k2) * 1 i) + (1 / k1) * (cos (t2 * k2) * sin (t1 * k1) * 1我);
                               
                                m21 = (k1) * cos (t2 * k2) * sin (t1 * k1) * 1 + (k2) * cos (t1 * k1) * sin (t2 * k2) * 1我;
                               
                                锰= cos (t1 * k1) * cos (t2 * k2) - (k1和k2) * sin (t1 * k1) * sin (t2 * k2);
                               
                                一个= [m11公路m12;m21锰];
                               
                                disp (eig (A))
                               
                                gs = x。^ ^ 2 * k0 2-ns。^ 2;
                               
                                gc = x。^ ^ 2 * k0 2-nc。^ 2;
                               
                                y = 1 * (gs * m11公路+ gc *锰)-m21 + gc * g * m12;
                               
                                % y = x ^ 2 + x + 1;
                               
                                结束
                               
                                结束
                               
                                p0 = 1;
                               
                                p1 = 1.5;
                               
                                p2 = 2;
                               
                                托尔= 10 ^ 5;
                               
                                n = 50;
                               
                                h1 = p1 - p0;
                               
                                h2 = p2 - p1;
                               
                                del1 = (F (p1) - F (p0)。/ h1;
                               
                                del2 = (F (p2) - F (p1)。/ h2;
                               
                                d = (del2-del1)。/ (h2 + h1);
                               
                                我= 3;
                               
                                %步骤2
                               
                                而我< = n
                               
                                %步骤3
                               
                                b = del2 + h2。* d;
                               
                                D =√6 ^ 2 - 4 (b。* F (p2)。* D);%可能是复杂的
                               
                                %步骤4
                               
                                如果abs (b - D) < abs (b + D)
                               
                                E = b + D;
                               
                                其他的
                               
                                E = b - D;
                               
                                结束
                               
                                %步骤5
                               
                                h = 2。* F (p2)。/ E;
                               
                                p = p2 + h;
                               
                                如果我= = 3
                               
                                表{1}=穆勒”年代的迭代方法;
                               
                                表{2}=“我f P (P)”;
                               
                                表{3}=“- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -”;
                               
                                结束
                               
                                str = sprintf (% 3 u: % 6.6 + % 6.6 fi % 6.6 f + % 6.6 fi ',我,真正的(p),图像放大(p),真正的(F (p)),图像放大(F (p)));
                               
                                表{I + 1} = str;% #好< * AGROW >
                               
                                %步骤6
                               
                                如果abs (h) <托尔
                               
                                val = p;
                               
                                表= char(表);
                               
                                打破
                               
                                结束
                               
                                p0 = p1;
                               
                                p1, p2;
                               
                                p2 = p;
                               
                                h1 = p1 - p0;
                               
                                h2 = p2 - p1;
                               
                                del1 = (F (p1) - F (p0)。/ h1;
                               
                                del2 = (F (p2) - F (p1)。/ h2;
                               
                                d = (del2-del1)。/ (h2 + h1);
                               
                                我=我+ 1;
                               
                                结束
                               
                                disp (p)
                               
                                结束