用matlab计算离轴三反初始结构出错

佞臣

参考(建议收藏)几种离轴三反初始结构计算汇总方法二和附件pdf编程计算求解a1时，解出来是复数计算不正确，只取实数计算也有问题，，，有木有佬知道的可指点一二

1. %%
   
2. %%方法二求初始结构
   
3. clear;
   
4. clc;
   
5. D1=-500
   
6. D2=500
   
7. L3=-500
   
8. f=-1500
   
9. %输入初始参数
   
10. A=D2*L3/(f^2)-D1/f;
    
11. B=2*L3*(D1-D2)/(f^2)+D1*D2/(f^2);
    
12. C=L3/f*(D2/f-D1*L3/(f^2));
    
13. a1=((-B+sqrt(B*B-4*A*C))/(2*A));
    
14. 
    
15. % p1=[A B C];
    
16. % y1=poly2sym(p1)
    
17. % a12=real(solve(y1,0))
    
18. % a1=double(a12(1))
    
19. a2=L3/(a1*f);
    
20. b2=a1*(1-a2)/D2*f;
    
21. b1=(1-a1)*f/(D1*b2);
    
22. r1=2*f/(b1*b2);
    
23. r2=2*a1/(b2*(1+b1))*f;
    
24. r3=2*a1*a2*f/(1+b2);
    

复制代码

小黄爱学习

您好 您找到错误了吗

Nayuta

我也是这个问题，发现直接带入数据会出现复数。

Nayuta

找到问题了，原公式B应该是 - d1d2/f^2不是+，这样算出来的R1R2R3就和原文一样了，应该是论文有误？