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

li205212021

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

0 r/ _2 E2 f  m4 f' \( d4 _理论上不管num 为几 第一个值不会变的
$ a3 b3 v( U0 q" b但是
num = 1 时得到第一个值
( q9 Z$ Y" F! [; b: t
; P" C5 w2 c7 p0 @num = 2时得到两个值 但是，第一个值和前面num=1的值不一样是怎么回事
0 ~, ~* S+ ]$ c            
" |: M8 `- V( o8 o/ x
各位大大帮忙看下，不然睡不成觉了
3 @" T. S7 ]$ b/ J- w4 b' ^1 c部分程序如下num = 2; %设置层数
. n( J3 k/ g0 C1 C5 L( }
3 ?6 h- n0 {. d- x6 ?7 z2 N1 |        g=linspace(1,2.4,200);
% G: B7 e: c+ @7 Q        for x=g
/ ]5 o4 U' b+ ^, \/ T4 H9 Z                %循环镶嵌处
            for len =1:1:num  %len为厚度
) @8 D2 g" b0 t, q' x                F1=2*pi*na*d1*cos(c2)/x;
2 M" ^: y, Z+ A  E; M: a  E                M1=[cos(F1),i*sin(F1)/p1;i*p1*sin(F1),cos(F1)];

                d2 = d2 - 0.005*(len-1);% 厚度

                F2{len,1}=2*pi*nb*d2*cos(c3)/x; %第len层 角度
) [( K9 b; [. o* {  A9 B$ }  `                M2{len,1}=[cos(F2{len,1}),i*sin(F2{len,1})/p2;i*p2*sin(F2{len,1}),cos(F2{len,1})]; %len层 矩阵表示

                temp{len,1} = M1 * M2{len,1}; % 矩阵A 乘 矩阵 B(len)  得到num个 A B 相乘的矩阵
                result = [1,0;0,1]; %初值为单位矩阵
, d  a  I  ?1 u) i. P+ U                result = result * temp{len,1};%得到中间所有矩阵的乘积
9 f# K& j+ W4 e          end
# j9 u5 g6 r; w0 p' b* `+ W" L
1 _, G7 L+ K& \5 `8 f

- R: F8 X6 X/ x+ ~+ v6 [( N- @
" T& ~5 C9 m- p2 J5 d& g6 J

Tdongxue

你的M2{len,1}是由F2{len,1}计算出来的，而
F2{len,1}=2*pi*nb*d2*cos(c3)/x;
7 u9 D* N! D8 h9 w& w你应该分析F2在num发生变化时状况，例如你这里的d2，在num不同时
7 N  a0 s9 C3 |: ]2 ^d2 = d2 - 0.005*(len-1);
恐怕d2得到的值也会不一样，另外还有你没有给出的nb和c3
/ A. Y# s% s9 I; S6 M' s0 M这些都要你一步步分析的，而不是盯着M2奇怪就能解决问题的

li205212021

Tdongxue 发表于 2023-2-8 14:57
你的M2{len,1}是由F2{len,1}计算出来的，而
F2{len,1}=2*pi*nb*d2*cos(c3)/x;
你应该分析F2在num发生变化 ...

5 `' S" B5 E4 H找到个错误 就是 result = [1,0;0,1]; 初值应该在第一个循环里面第二个循环外面
具体不知道哪儿错了 变的就只有d2 从而试F2 M2 变化其余都是已知的
5 b- u" y+ _0 t: u3 }% T! `  Q) |