clear
r = 0:0.1:1;
phi = -5:0.0001:25;
for i1 = 1:11
y1 = 1./(1+4*r(i1)*sin(phi/2).^2/(1-r(i1))^2);%反射光强
plot(phi,y1)
hold on
end
text(-3,0.1,'r = 0.8')
text(5,0.2,'r = 0.5')
text(7,0.44,'r = 0.6')
text(12,0.5,'r = 0.01')
text(16,0.65,'r = 0.4')
text(20,0.89,'r = 0.23')
xlabel('角变量')
ylabel('反射光强')
figure(2)
for i1 = 1:11
y2 = 1./(1+(1-r(i1))^2./(4*r(i1)*sin(phi/2).^2));%反射光强
plot(phi,y2);
hold on
end
text(-3,0.1,'r = 0.8')
text(5,0.2,'r = 0.5')
text(7,0.44,'r = 0.6')
t
