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

li205212021

通过改变num的值计算不同层数的值，需要记忆前一个值
' X! ^& @3 ]. g$ V
理论上不管num 为几 第一个值不会变的
但是
/ K* k$ O% N  h1 ~. unum = 1 时得到第一个值
+ I) Y* M, L8 A# v# H
num = 2时得到两个值 但是，第一个值和前面num=1的值不一样是怎么回事
) P; @4 a: o  o7 k            
5 F) I+ h- f' H; @- ^! z) K/ b
  Q3 l* F, d8 \2 `各位大大帮忙看下，不然睡不成觉了
% j, z+ h/ C7 i! h3 b部分程序如下num = 2; %设置层数
: g/ m  f; M5 M+ I6 c% Z* k  v
6 p) \) G& |2 S  T        g=linspace(1,2.4,200);
- _/ K1 V! V5 u- r% X2 ^        for x=g
5 B/ Z) X6 a8 s8 l, L$ t; N' x2 M                %循环镶嵌处
( m( k8 ~  }; F+ a% ^            for len =1:1:num  %len为厚度
( U# v2 p3 h5 k9 i- a                F1=2*pi*na*d1*cos(c2)/x;
                M1=[cos(F1),i*sin(F1)/p1;i*p1*sin(F1),cos(F1)];
- X6 V$ k  w# H0 k% n
- Q4 z. P" |9 C5 c; G% J1 H: S                d2 = d2 - 0.005*(len-1);% 厚度

                F2{len,1}=2*pi*nb*d2*cos(c3)/x; %第len层 角度
1 I* p; l* ^8 n- R, I* U7 [. l                M2{len,1}=[cos(F2{len,1}),i*sin(F2{len,1})/p2;i*p2*sin(F2{len,1}),cos(F2{len,1})]; %len层 矩阵表示

( G- u! i, O! y; r8 W& l: ?, K                temp{len,1} = M1 * M2{len,1}; % 矩阵A 乘 矩阵 B(len)  得到num个 A B 相乘的矩阵
+ L! @1 H8 x4 t! l8 e                result = [1,0;0,1]; %初值为单位矩阵
& V4 K/ X$ Y2 A" N( L                result = result * temp{len,1};%得到中间所有矩阵的乘积
          end
$ i# l# c& V/ |  C0 {
& ]4 ]1 n9 y0 b- |2 p0 T* G
' t2 T3 \  H& ~' @# N/ h7 |  e

* E* R1 _0 G3 b' |6 E+ j

Tdongxue

你的M2{len,1}是由F2{len,1}计算出来的，而
: B" W- E3 C/ p+ j& l2 _F2{len,1}=2*pi*nb*d2*cos(c3)/x;
你应该分析F2在num发生变化时状况，例如你这里的d2，在num不同时
* Z+ K0 D' ]; m) @6 M  O8 s# a3 I2 `* T, Wd2 = d2 - 0.005*(len-1);
# e- U, P% z6 v+ h( R  f恐怕d2得到的值也会不一样，另外还有你没有给出的nb和c3
这些都要你一步步分析的，而不是盯着M2奇怪就能解决问题的

li205212021

Tdongxue 发表于 2023-2-8 14:57
' ~$ [" {, c# c# J0 O  Z你的M2{len,1}是由F2{len,1}计算出来的，而
$ e0 J/ C3 ?$ _; u7 R5 [, x: S5 RF2{len,1}=2*pi*nb*d2*cos(c3)/x;
你应该分析F2在num发生变化 ...

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