循环问题，两次循环不一致

li205212021

通过改变num的值计算不同层数的值，需要记忆前一个值

理论上不管num 为几 第一个值不会变的
5 O5 ~: i" O, f但是
+ G" _7 f0 j: o3 Wnum = 1 时得到第一个值
* X9 D! Y8 I! K/ g* Y+ i2 h& ?
+ L& F" \8 ]  knum = 2时得到两个值 但是，第一个值和前面num=1的值不一样是怎么回事
           
0 A% D8 Q9 I& s0 P8 l/ a( P
$ a9 D8 W, o% y" }( T各位大大帮忙看下，不然睡不成觉了
部分程序如下num = 2; %设置层数

+ \' D2 b& n; Y! Q1 N0 u3 G/ k        g=linspace(1,2.4,200);
/ J  I; A; E- i        for x=g
                %循环镶嵌处
            for len =1:1:num  %len为厚度
                F1=2*pi*na*d1*cos(c2)/x;
                M1=[cos(F1),i*sin(F1)/p1;i*p1*sin(F1),cos(F1)];
1 B- D* z2 I4 W1 J  B4 z& A
- }) d% m3 q4 \0 J7 c8 _7 G$ z& b! V; \                d2 = d2 - 0.005*(len-1);% 厚度
  m  p! a9 x, {7 @% s
0 {; i/ N: }2 B5 o                F2{len,1}=2*pi*nb*d2*cos(c3)/x; %第len层 角度
2 ]. Z( o% R# h- U- ?6 [. j                M2{len,1}=[cos(F2{len,1}),i*sin(F2{len,1})/p2;i*p2*sin(F2{len,1}),cos(F2{len,1})]; %len层 矩阵表示
* W$ ?. X# w# O' d) r# v. Y2 K, ?, e
: ], ]  n  O7 O$ \- G6 X& k) Y                temp{len,1} = M1 * M2{len,1}; % 矩阵A 乘 矩阵 B(len)  得到num个 A B 相乘的矩阵
                result = [1,0;0,1]; %初值为单位矩阵
                result = result * temp{len,1};%得到中间所有矩阵的乘积
$ h+ [) J- n# J% s          end

. ?# g  D, |6 ?
; O2 i& N5 n- E5 c5 u

Tdongxue

你的M2{len,1}是由F2{len,1}计算出来的，而
F2{len,1}=2*pi*nb*d2*cos(c3)/x;
7 m  k7 H" [! r0 t) O! E  M你应该分析F2在num发生变化时状况，例如你这里的d2，在num不同时
d2 = d2 - 0.005*(len-1);
恐怕d2得到的值也会不一样，另外还有你没有给出的nb和c3
这些都要你一步步分析的，而不是盯着M2奇怪就能解决问题的

li205212021

Tdongxue 发表于 2023-2-8 14:57
& {5 G4 P4 I* h+ j+ E& g: P你的M2{len,1}是由F2{len,1}计算出来的，而
F2{len,1}=2*pi*nb*d2*cos(c3)/x;
/ Q" n7 m# ~: L0 b1 {4 x2 o9 U你应该分析F2在num发生变化 ...

找到个错误 就是 result = [1,0;0,1]; 初值应该在第一个循环里面第二个循环外面
具体不知道哪儿错了 变的就只有d2 从而试F2 M2 变化其余都是已知的