作者:懒羊羊
时间:2020.9.16 19:24 耗时2h30min原创
Lanzhou University of Technology
‘’个人代码备忘录‘’

  • ✪实验作业1

  • ✪实验目的

  • 1、熟悉MATLAB程序编辑与开发环境,掌握使用该软件解决移动通信问题的具体方法;

  • 2、掌握无线电波传播损耗计算的基本知识;

  • 3、掌握运用Hata-Okumura模型预测的方法实现大城市、中小城市、郊区、开阔农村地带的无线传播损耗预测系统;

  • 4、通过无线传播损耗预测仿真加深对Hata-Okumura无线传播损耗预测模型的理解。

  • 详细代码如下:

clc;clear;
%“author:<十三月凉>”
%“time:2020.9.16”
%“Lanzhou University of Technology”
%Hata-Okumura模型
%% 公式说明
%f:中心频率; hb:基站天线有效高度; hm:移动台天线有效高度,单位:m; d:通信距离,单位:km; a_hm:移动台天线高度校正;
%市区传播公式损耗公式:L=69.55+26.16*log10(f)-13.82*log10(hb)+(44.9-6.55*log10(hb))*log10(d)-a_hm;
%    其中高度校正a_hm: 1.大城市:a_hm=8.29*(log10(1.54*hm))^2-1.1; f<=200MHz
%                              a_hm=3.2*(log10(11.75*hm))^2-4.97; f>=400MHz
%                     2.中小城市:a_hm=(1.1*log10(f)-0.7)*hm-(1.56*log10(f)-0.8);
%郊区:Lj=L-2*(log10(f/28))^2-5.4;
%开阔农村地带:Ln=L-4.78*(log10(f))^2-18.33*log10(f)-40.98;
hm=3;
hb=200;d1=1:0.1:10;
figure()
subplot(221);
for f1=300:200:1500a_hm=3.2*(log10(11.75*hm))^2-4.97; Ld1=69.55+26.16*log10(f1)-13.82*log10(hb)+(44.9-6.55*log10(hb))*log10(d1)-a_hm;plot(d1,Ld1,'rx-');hold on;
end
xlabel('传播距离(单位km)');ylabel('传输损耗L(单位dB)');
title('(大城市)传输损耗随传播距离的变化曲线');
hold off;subplot(222);
for f1=300:200:1500a_hm=(1.1*log10(f1)-0.7)*hm-(1.56*log10(f1)-0.8);Lx1=69.55+26.16*log10(f1)-13.82*log10(hb)+(44.9-6.55*log10(hb))*log10(d1)-a_hm;plot(d1,Lx1,'yx-');hold on;end
xlabel('传播距离(单位km)');
ylabel('传输损耗L(单位dB)');title('(中小城市)传输损耗随传播距离的变化曲线');hold off;
subplot(223);
for f1=300:200:1500a_hm=3.2*(log10(11.75*hm))^2-4.97;Ld=69.55+26.16*log10(f1)-13.82*log10(hb)+(44.9-6.55*log10(hb))*log10(d1)-a_hm;Lj1=Ld-2*(log10(f1/28))^2-5.4;plot(d1,Lj1,'bx-');hold on;
end
xlabel('传播距离(单位km)');ylabel('传输损耗L(单位dB)');
title('(郊区)传输损耗随传播距离的变化曲线');
hold off;
subplot(224);
for f1=300:200:1500a_hm=3.2*(log10(11.75*hm))^2-4.97;Ld=69.55+26.16*log10(f1)-13.82*log10(hb)+(44.9-6.55*log10(hb))*log10(d1)-a_hm;Ln1=Ld-4.78*(log10(f1))^2-18.33*log10(f1)-40.98;plot(d1,Ln1,'cx-');hold on;
end
xlabel('传播距离(单位km)');ylabel('传输损耗L(单位dB)');title('(开阔农村地带)传输损耗随传播距离的变化曲线');
hold off;
hm=3;hb=200;
f2=300:10:1500;
figure()
subplot(221);
for d2=1:0.5:10a_hm=3.2*(log10(11.75*hm))^2-4.97;Ld2=69.55+26.16*log10(f2)-13.82*log10(hb)+(44.9-6.55*log10(hb))*log10(d2)-a_hm;    plot(f2,Ld2,'rx-');hold on;
end
xlabel('工作频率(单位Hz)');ylabel('传输损耗L(单位dB)');title('(大城市)传输损耗随工作频率的变化曲线');
hold off;subplot(222);for d2=1:0.5:10a_hm=(1.1*log10(f2)-0.7)*hm-(1.56*log10(f2)-0.8); Lx2=69.55+26.16*log10(f2)-13.82*log10(hb)+(44.9-6.55*log10(hb))*log10(d2)-a_hm;plot(f2,Lx2,'yx-');hold on;end
xlabel('工作频率(单位Hz)');
ylabel('传输损耗L(单位dB)');title('(中小城市)传输损耗随工作频率的变化曲线');
hold off;subplot(223);
for d2=1:0.5:10a_hm=3.2*(log10(11.75*hm))^2-4.97;Ld=69.55+26.16*log10(f2)-13.82*log10(hb)+(44.9-6.55*log10(hb))*log10(d2)-a_hm;Lj2=Ld-2*(log10(f2/28)).^2-5.4;plot(f2,Lj2,'bx-');hold on;
end
xlabel('工作频率(单位Hz)');ylabel('传输损耗L(单位dB)');title('(郊区)传输损耗随工作频率的变化曲线');
hold off;subplot(224);
for d2=1:0.5:10a_hm=3.2*(log10(11.75*hm))^2-4.97;Ld=69.55+26.16*log10(f2)-13.82*log10(hb)+(44.9-6.55*log10(hb))*log10(d2)-a_hm;Ln2=Ld-4.78*(log10(f2)).^2-18.33*log10(f2)-40.98;    plot(f2,Ln2,'cx-');hold on;
end;
xlabel('工作频率(单位Hz)');ylabel('传输损耗L(单位dB)');
title('(开阔农村地带)传输损耗随工作频率的变化曲线');
hold off;

✪需要注意的一点是subplot语句必须要放在循环之前,否则绘图会出现错误!

仿真结果如下:

and

✪一点其他的想法:
✪以上代码还可以进行优化,代码优化思路如下:
1.利用case选择语句进行模型选择,即将循环嵌套进判断语句里;
2. 利用position语句进行图注以对代码行数进行压缩

matlab Hata Okumura模型相关推荐

  1. 基于MATLAB的COST231 Hata信道模型的仿真与理论分析

    COST231 Hata模型应用频率在1 500MHz到2 000MHz之间,适用于小区半径大于1km的宏蜂窝系统,发射有效天线高度在30m到200m之间,接收有效天线高度在1m到10m之间.路径损耗 ...

  2. matlab求迟滞,MATLAB PI迟滞模型问题

    之前我用matlab进行迟滞模型建模,误差非常大,我之前怀疑是自己的参数辨识出问题了,所以我用论文中的阈值向量和权重向量进行建模,输入数据也是用论文数据,但是出来的拟合的线误差非常大,和文献中相差太大 ...

  3. matlab 传感器的迟滞,MATLAB PI迟滞模型问题

    之前我用matlab进行迟滞模型建模,误差非常大,我之前怀疑是自己的参数辨识出问题了,所以我用论文中的阈值向量和权重向量进行建模,输入数据也是用论文数据,但是出来的拟合的线误差非常大,和文献中相差太大 ...

  4. 用matlab简单电路模型,基于MATLAB的电路模型仿真应用

    基于MATLAB的电路模型仿真应用实验指导书 一.实验目的 1.掌握采用M文件及SIMULINK对电路进行仿真的方法. 2.熟悉POWERSYSTEM BLOCKSET 模块集的调用.设置方法. 3. ...

  5. matlab传函零极点形式,2013实验一 MATLAB 中控制系统模型的建立与仿真

    实验一 MATLAB 中控制系统模型的建立与仿真 一. 实验目的 (1)熟悉MATLAB 控制系统工具箱中线性控制系统传递函数模型的相关函数. (2)熟悉SIMULINK 模块库,能够使用SIMULI ...

  6. matlab加热模型,基于MATLAB钢丝感应加热模型及其温度场特性的研究

    <工业加热>第 卷 年第 期 11 热能工程 前 言 自 20 世纪 30 年代,感应加热技术大量应用于工业生产以来,由于感应加热自身的优点和感应加热理论和技术的不断进步,迄今已成为一项节 ...

  7. Matlab - Solidworks 机器人建模(4)—— 如何把SolidWorks模型导入到Matlab (Simscape模型)

    Matlab - Solidworks 机器人建模(4)-- 如何把SolidWorks模型导入到Matlab (Simscape模型) 1.说在前面 本文会介绍怎么直接把solidworks的模型转 ...

  8. 汽车matlab小论文,基于matlab汽车道路模型研究论文

    基于matlab汽车道路模型研究论文 武汉工业学院 毕业设计(论文) 2011届 毕业设计(论文)题目:基于MATLAB的汽车道路模型研究 院 (系) 机械工程学院 专业名称 机械设计制造及其自动化 ...

  9. matlab获得模型切面,基于MATLAB的STL模型切片分层算法

    吴建 吴婷 陈廷豪 包涵 摘  要:为提高STL模型的切片效率,提出一种基于MATLAB的STL模型切片分层新算法.首先对STL模型进行数据预处理,筛选出只与切平面相交的三角片集合,然后利用相邻三角面 ...

最新文章

  1. Ubuntu编译安装VASP
  2. mapreduce编程实例(2)-求最大值和最小值
  3. 采用Flume实时采集和处理数据
  4. oracle查询哪个sql占用资源,查询oracle最耗资源的sql语句
  5. 安装head插件依赖包grunt-cli
  6. 随记(电脑硬件对比)
  7. 用c++编写的植物大战僵尸
  8. 光影精灵usb安装linux,惠普光影精灵5笔记本怎么装win10系统(uefi+gpt)
  9. python教学反思_Python语言教学反思
  10. Pyhton 兔子繁衍问题(后有递归求法)(设有一对新生的兔子从第4个月开始他们每个月月初都生一对新兔子, 新生的兔子从第4个月开始又每个月月初生一对兔子,求n个月后兔子的总对数)
  11. 第一啪,第一啪电影网,第一啪电影网用的哪里的模板diyipa.cc
  12. 区块链之发行以太坊 Token
  13. C++ QT有道翻译 爬虫 分析 破解
  14. OS-实现一个RR调度算法
  15. vim列删除、列修改
  16. 合泰单片机做电压表_合泰单片机基础教程
  17. tf.keras.layers.MaxPooling2D
  18. window 命令行bat脚本软连接执行exe
  19. access的未来发展 ieee_IEEE Access 期刊怎么样,在国外认可度高吗?
  20. 群论在物理学中的应用

热门文章

  1. 硬件工程师如何选择一款合适的霍尔电流传感
  2. Java内置队列和高性能队列Disruptor
  3. android 高德地图 3d,在React Native中使用 高德地图组件react-native-amap3d
  4. 图形设计六字诀影叠置断悖质优漫动
  5. source和resource的区别
  6. 月薪10K!汽修专业小伙转行5G网络优化工程师,零基础跨行真的很难嘛?
  7. C++学习 1: 两个无符号数相减
  8. 电话面试常见问题及应答技巧
  9. python基础教程期末考试题_智慧树Python程序设计基础期末教程考试完整答案
  10. 视频直播源码,插入图片、删除图片、设置图片大小、提取图片