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

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

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