Okumura-Hata模型的计算机仿真.doc

实验基本原理及公式移动通信中电波传播的实际情况是复杂多变的。实践证明，任何试图使用一个或几个理论公式计算的结果，都将引入较大误差。甚至与实测结果相差甚远。为此，人们通过大量的实地测量和分析，总结归纳了多种经验模型。通常在一定情况下，使用这些模型对移动通信电波传播特性进行估算，都能获得比较准确的预测结果。目前应用较为广泛的是 OM 模型（Okumura 模型），为了在系统设计时，使Okumura 模型能采用计算机进行预测，Hata 对 Okumura 模型的基本中值场强曲线进行了公式化处理，所得基本传输损耗的计算公式如下：式中， d 为收发天线之间的距离，km；f 为工作频率，MHz；hb为基站天线有效高度，m； (hm) 为移动台天线高度校正因子，hm为移动台天线高度（m）。 (hm) 由下式计算：此公式适用范围为：150MHzf1920MHz 30mhb200m，1mhm10m，1kmd20km，准平滑地形。实验内容1、比较不同地形下的衰减情况：针对市区、郊区和开阔区，比较在相同条件（基站天线高度、移动台天线高度、传播距离）下，随着频率的改变，电波传播衰减的情况。2、比较相同地形不同城市类型下的衰减情况：针对市区、郊区或者开阔区任一种地形，比较在相同条件（基站天线高度、移动台天线高度、传播距离）下，不同城市类型（大城市或者中、小城市），随着频率的改变，电波传播衰减的情况。实现代码：#include#includefloat computer1(float *p,float *q,float h1,float h2);float computer2(float *p,float *q,float h1,float h2);float computer3(float *p,float *q,float h1,float h2);float computer5(float *p,float *q,float h1,float h2);float computer4(float *p,float *q,float h1,float h2);float computer6(float *p,float *q,float h1,float h2);void main()float frequency,distan,height1,height2;int eara,type;for(;)printf(choose your eara:nt1.city earant2.suburbannt3.open earan);scanf(%d,&eara);if(eara=1)printf(choose your city:nt1.bignt2.small and middn );scanf(%d,&type);printf(input the frequencyn);scanf(%f,&frequency);printf(input your distancen);scanf(%f,&distan);printf(input your height of jizhann);scanf(%f,&height1);printf(input your height of mobilen);scanf(%f,&height2);if(frequency1920|frequency200|height110|height220|distan1)printf(your distance is illegaln);continue;switch(eara)case 1:if(type=1)printf(%.3f,computer1(&frequency,&distan,height1,height2);else printf(%.3f,computer4(&frequency,&distan,height1,height2); break;case 2:printf(%.3f,computer2(&frequency,&distan,height1,height2);break;case 3:printf(%.3f,computer3(&frequency,&distan,height1,height2);break;default:printf(some input of yours is error);float computer1(float *p,float *q,float h1,float h2)float m,n;if(*p=300)m=computer5(p,q,h1,h2);else m=computer6(p,q,h1,h2);n=69.55+26.16*log10(*p)-13.82*log10(h1)-m+(44.9-6.55*log10(h2)*log10(*q);return n;float computer2(float *p,float *q,float h1,float h2)return (computer1(p,q,h1,h2)-2*pow(log10(*p/28),2)-5.4);float computer3(float *p,float *q,float h1,float h2)return (computer1(p,q,h1,h2)-4.78*pow(log10(*p),2)+18.33*log10(*p)-40.98);float computer5(float *p,float *q,float h1,float h2)float m;m=8.29*pow(2,log10(1.54*h2)-1.1;return m;float computer6(float *p,float *q,float h1,float h2)return (3.2*pow(2,log10(11.75*h2)-4.97);float computer4(float *p,float *q,float h1,f