function?[]=distancebd（tem）
%該函數是你的畫圖程序，將此程序放到你的matlab文件庫
%調用的格式為distancebd（參數值）
%以下是加裝耦合線時候導線上面感應雷過電壓百分比的變化情況???
%先畫出耦合線距離桿塔中心線距離dm為0.4m，然后是0.5m
%參數意義如下h：耦合線的高度
%u-per:表示Ud‘/Ud即加裝耦合線前后的導線上感應雷過電壓百分比，其中uperAuperB，uperC分別表示ABC三項情況
%最終的圖形反應的就是h與u-per的關系，h的范圍從圖中來看是0-15m
%所有數據根據圖5.3，下面就是程序：
%數據初始化
if?（tem>0.69）|（tem<0）
????disp（‘You?are?a?fool!‘）;
else
dm=tem;
disab=0.69;%ab相距離
da=disab-dm;%耦合線與a相的水平距離
ha=10;%ac相的高度
Rb=0.0033;%耦合線半徑
Rg=10;
%下面先畫出0.4m時候的A相，然后BC相只需要復制就可以，這里A相是圖5.3中從左至右的第一個
htem=0:0.01:15;
upertemp1=zeros（size（htem））;
i=1;
for?h=0:0.01:15
??????Dba=sqrt（（ha+h）^2+da^2）;
??????dba=sqrt（abs（ha-h）^2+da^2）;
??????Zba=60*log（Dba/dba）;
??????Zbb=60*log（2*h/Rb）;
??????uperA=1-（Zba/（Zbb+2*Rg））*（h/ha）;
??????upertemp1（i）=uperA;
??????if?h==9
??????????temp2=100*（1-upertemp1（i））;
??????end
??????i=i+1;???
end
upertemp1=100*（1-upertemp1）;%文中是電壓降低百分比而且這里將其放大100倍
%去掉百分號
plot（htemupertemp1‘r‘）;%到此為止在hm=0.4時候的A相圖就出來了
h=9;
text（htemp2‘\rightarrow相導線1‘）
hold?on;%這是為了在同一幅圖中繼續下一相圖
%以下的程序就是復制
%畫B相
da=dm;%此時耦合線距離中間線即B相的距離為dm=4；所以只要將上面的da賦值為4即可
ha=10.95;%此時B相的高度為10.95m
%下面一切如上
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
htem=0:0.01:15;
upertemp1=zeros（size（htem））;
i=1;
for?h=0:0.01:15
??????Dba=sqrt（（ha