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

li205212021

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

理论上不管num 为几 第一个值不会变的
但是
+ F/ ]/ v! H! j+ `9 R5 c2 N; ]9 ?6 d, Inum = 1 时得到第一个值

# x& m; ]+ P) jnum = 2时得到两个值 但是，第一个值和前面num=1的值不一样是怎么回事
           

各位大大帮忙看下，不然睡不成觉了
部分程序如下num = 2; %设置层数
* i* C4 V) I2 l; j4 F+ Z% [$ [
$ u$ M1 B+ b! _' F% v) }% `0 l% v        g=linspace(1,2.4,200);
        for x=g
( `2 P% S, F2 t/ e                %循环镶嵌处
            for len =1:1:num  %len为厚度
# O5 U; i, c2 a9 z) K. y& b) S" S                F1=2*pi*na*d1*cos(c2)/x;
                M1=[cos(F1),i*sin(F1)/p1;i*p1*sin(F1),cos(F1)];
3 G- }/ |- [$ t" }; `/ w$ l: X
1 H" N0 o! J0 y. c; n                d2 = d2 - 0.005*(len-1);% 厚度

" n1 X6 E* e' x1 @9 G( M# x1 M                F2{len,1}=2*pi*nb*d2*cos(c3)/x; %第len层 角度
                M2{len,1}=[cos(F2{len,1}),i*sin(F2{len,1})/p2;i*p2*sin(F2{len,1}),cos(F2{len,1})]; %len层 矩阵表示
) f5 g' ?# B8 y- U: s
                temp{len,1} = M1 * M2{len,1}; % 矩阵A 乘 矩阵 B(len)  得到num个 A B 相乘的矩阵
. x) [+ e1 M' A; Q                result = [1,0;0,1]; %初值为单位矩阵
                result = result * temp{len,1};%得到中间所有矩阵的乘积
5 j& P% e, O. T. A" c' V          end
  s, Y. d. ?- T  j& y( k5 l& E! ?+ H

5 g, Y) f: `4 K( i1 `

/ m( ?( |+ v% _: B+ _6 t6 q0 d

Tdongxue

你的M2{len,1}是由F2{len,1}计算出来的，而
! T4 q8 Z0 n+ v) g9 ]( q! mF2{len,1}=2*pi*nb*d2*cos(c3)/x;
( }2 u; q9 a$ c# e) ^4 T+ {( S/ B你应该分析F2在num发生变化时状况，例如你这里的d2，在num不同时
+ f1 l" _8 j; [1 x1 K+ Dd2 = d2 - 0.005*(len-1);
恐怕d2得到的值也会不一样，另外还有你没有给出的nb和c3
这些都要你一步步分析的，而不是盯着M2奇怪就能解决问题的

li205212021

Tdongxue 发表于 2023-2-8 14:57
3 a7 A7 m& _) T. n你的M2{len,1}是由F2{len,1}计算出来的，而
, u$ w: Q% j4 }2 K" }9 L' vF2{len,1}=2*pi*nb*d2*cos(c3)/x;
你应该分析F2在num发生变化 ...

, Q* F, s( `: n; {. c5 a找到个错误 就是 result = [1,0;0,1]; 初值应该在第一个循环里面第二个循环外面
具体不知道哪儿错了 变的就只有d2 从而试F2 M2 变化其余都是已知的
* G/ D/ k& M) y) R/ ]