以下为第一问理想反射面的求解代码,完整题目的视频讲解和代码请联系微信:18655167262
clear
close
clc
YMax=zeros(1,1200); %建立零矩阵,用来存放每个l所对应的最大l_max
YSum=zeros(1,1200); %建立零矩阵,用来存放每个l所对应的最大l_sum
Theta_i=zeros(1,1200); %%建立零矩阵,用来存放每个l所对应的最大theta范围
for i=1:1200 %l从-0.6~0.6,步长为0.001,所以是1200
syms theta
R=300.4; %基准状态的半径
F_l=0.466*R-0.6+0.001*i; %焦距,相当于F-l
CP=R-0.466*R; %馈源仓所在圆的半径
rho=2.*F_l./(1-cos(theta)); %抛物线的极坐标方程,主要通过这个来计算焦距
r=rho.*sin(theta); %理想抛物面上任意一点到对称轴上的距离
deltaTheta=abs(double(solve(r-150,theta))); %计算理想抛物面口径为300时theta为多