雷电感应电压matlab

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