clear
a=?linspace（0?2*?pi200）?;
b=?linspace（-pi?pi200）?;
L=1.25*pi;??%L=kl=2l/入

%計算歸一化方向圖函數
f?=?abs（（cos（L.*（cos（b）））-cos（L））./sin（b））;
F=f/max（f）;

%畫出e面方向圖
subplot（121）;
polar（b+0.5*pi?F?）hold?on;
polar（b+0.5*piones（1200）*0.707‘k--‘）t