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

) C7 M! `9 h9 p; j* d% A' n$ T找到个错误 就是 result = [1,0;0,1]; 初值应该在第一个循环里面第二个循环外面
- [+ l9 M1 a6 h/ J具体不知道哪儿错了 变的就只有d2 从而试F2 M2 变化其余都是已知的

你的M2{len,1}是由F2{len,1}计算出来的，而
& m% U5 C/ W/ Y2 G' b. JF2{len,1}=2*pi*nb*d2*cos(c3)/x;
2 w* k* U4 l3 M你应该分析F2在num发生变化时状况，例如你这里的d2，在num不同时
d2 = d2 - 0.005*(len-1);
恐怕d2得到的值也会不一样，另外还有你没有给出的nb和c3
7 D! q- R: e4 V! V+ t这些都要你一步步分析的，而不是盯着M2奇怪就能解决问题的