不同状态角 θz 下的 K2eff 和 Kav 计算结果

 

不同状态角 θ z 下的 K 2 eff 和 K av 计算结果

clear 
clc
close all
syms theta
for theta_z=[pi/5,pi/3,pi/2,2*pi/3,4*pi/5,pi]  %
theta_zd=theta_z*180/pi
theta_1=-theta_z/2;
theta_2=theta_z/2;
ftheta=cos(theta);
ftheta2=cos(theta)^2;
Keff2=int(ftheta2,theta_1,theta_2)/theta_z;
K_eff2=vpa(Keff2,4)  %化简为小数显示
Kav=int(ftheta,theta_1,theta_2)/theta_z;
K_av=vpa(Kav,4)  %化简为小数显示
Keff2Kav=Keff2/Kav;
Keff2_Kav=vpa(Keff2Kav,4)

end