Matlab:绘制简单能量的接收机工作特性曲线(Energy_detection_simulation_ok)-阿里云开发者社区

开发者社区> 一个处女座的程序猿> 正文

Matlab:绘制简单能量的接收机工作特性曲线(Energy_detection_simulation_ok)

简介: Matlab:绘制简单能量的接收机工作特性曲线(Energy_detection_simulation_ok)
+关注继续查看

输出结果

image.png


实现代码

% This code is to plot receiver operating characteristic curve for simple energy

%绘制简单能量的接收机工作特性曲线

% detection, when the primary signal is real Gaussian signal and noise is

% addive white real Gaussian. Here, the threshold is available

% analytically.

% Code written by: Sanket Kalamkar, Indian Institute of Technology Kanpur,

% India.


%% 以下代码绘制在虚警概率一定时,检测概率和信噪比之间的关系曲线称为检测器的检测性能曲线

clc

close all

clear all

L = 1000;                                             % The number of samples

snr = 0.01:0.01:10;

Pf = 10e-4;                                           % Pf = Probability of False Alarm 虚警概率确定

%% Simulation to plot SNR vs.Probability of Detection (Pd)

for m = 1:length(snr)

   i = 0;

   thresh = (qfuncinv(Pf)./sqrt(L))+ 1;              % Theoretical value of Threshold, refer, Sensing-Throughput Tradeoff for Cognitive Radio Networks, Y. C. Liang

   for kk = 1:5000                                   % Number of Monte Carlo Simulations(https://cn.mathworks.com/discovery/monte-carlo-simulation.html)

       n = randn(1,L);                               % AWGN noise with mean 0 and variance 1

       s = sqrt(snr(m)).*randn(1,L);                 % Real valued Gaussina Primary User Signal

       y = s + n;                                    % Received signal at SU(认知用户接收到的信号)

       energy = abs(y).^2;                           % Energy of received signal over N samples

       energy_fin =(1/L).*sum(energy);               % Test Statistic for the energy detection

       if(energy_fin >= thresh)                      % Check whether the received energy is greater than threshold, if so, increment Pd (Probability of detection) counter by 1

           i = i+1;

       end

   end

   Pd(m) = i/kk;

end

plot(10*log(snr), Pd,  'r')

xlabel('信噪比,单位db');

ylabel('检测概率');

title('能量感知检测性能曲线');

grid on

hold on

%% Theroretical expression of Probability of Detection; refer above reference.

thresh = (qfuncinv(Pf)./sqrt(L))+ 1;

%Pd_the = qfunc(((thresh - (snr + 1)).*sqrt(L))./(sqrt(2).*(snr + 1))); % 原来代码中的表达与论文中不一致

for k = 1:length(snr)

   Pd_the(k) = qfunc(((thresh - (snr(k) + 1)).*sqrt(L))./(sqrt(2).*snr(k) + 1)); % 与论文中的方程式保持一致

end

plot(10*log(snr), Pd_the, 'b')

legend('实际检测概率', '理论检测概率', 'Location', 'SouthEast');

hold on

grid on


版权声明:本文内容由阿里云实名注册用户自发贡献,版权归原作者所有,阿里云开发者社区不拥有其著作权,亦不承担相应法律责任。具体规则请查看《阿里云开发者社区用户服务协议》和《阿里云开发者社区知识产权保护指引》。如果您发现本社区中有涉嫌抄袭的内容,填写侵权投诉表单进行举报,一经查实,本社区将立刻删除涉嫌侵权内容。

相关文章
人工智能和自动化在工作场所中的重要性
人工智能为解决问题增加了新的维度,使一系列行业的企业受益。由于AI识别和预测大量数据模式的能力越来越可靠,因此AI正在帮助企业领导者事半功倍。
167 0
【对讲机的那点事】当对讲机不工作时如何处理?
对讲机属于电子类产品,是一种双向移动通信工具,在不需要任何网络支持的情况下,就可以通话,没有话费产生,适用于相对固定且频繁通话的场合,对讲机提供一对一,一对多的通话方式,一按就说,操作简单,令沟通更自由,尤其是紧急调度和集体协作工作的情况下,尤为彰显对...
942 0
阿里云服务器端口号设置
阿里云服务器初级使用者可能面临的问题之一. 使用tomcat或者其他服务器软件设置端口号后,比如 一些不是默认的, mysql的 3306, mssql的1433,有时候打不开网页, 原因是没有在ecs安全组去设置这个端口号. 解决: 点击ecs下网络和安全下的安全组 在弹出的安全组中,如果没有就新建安全组,然后点击配置规则 最后如上图点击添加...或快速创建.   have fun!  将编程看作是一门艺术,而不单单是个技术。
3956 0
【对讲机的那点事】450MHz模拟无线列调的工作原理(连载二)
在前文中我们谈到了450MHz模拟无线列调的工作原理,接下来我们来聊聊B制式模拟无线列调: B制式模拟车站台的功能框图: B制式工作原理: 车站对讲机应根据系统功能、设备制式的不同和线路条件,采用交替配置或单一配置方式。
1051 0
从RNN到LSTM,性能良好的神经网络到底是如何工作的?
在2016年Google率先发布上线了机器翻译系统后,神经网络表现出的优异性能让人工智能专家趋之若鹜。本文将借助多个案例,来带领大家一同探究RNN和以LSTM为首的各类变种算法背后的工作原理。
14093 0
CSDN社区分享面试经历活动作品11——我的第一份实习工作
我的第一次工作献给了她,虽然只是我的第一份实习工作,但是仍然让我至今难以忘怀~按时间顺序来:第一帖http://topic.csdn.net/u/20100915/13/1dbfadcd-e5a4-462f-bb01-29be96a10792.html貌似是我人生的第一次面试。
925 0
阿里云服务器如何登录?阿里云服务器的三种登录方法
购买阿里云ECS云服务器后如何登录?场景不同,阿里云优惠总结大概有三种登录方式: 登录到ECS云服务器控制台 在ECS云服务器控制台用户可以更改密码、更换系.
5676 0
【对讲机的那点事】450MHz模拟无线列调的工作原理(连载四)
在前文中我们谈到了450MHz模拟无线列调的B制式设备、C制式设备和车次号信息发送系统频点,接下来我们来聊聊无线列调调度通信业务的几个场景: 无线列调调度通信业务的场景: ●场景一:调度员呼叫司机: ●场景二:...
877 0
+关注
一个处女座的程序猿
国内互联网圈知名博主、人工智能领域优秀创作者,全球最大中文IT社区博客专家、CSDN开发者联盟生态成员、中国开源社区专家、华为云社区专家、51CTO社区专家、Python社区专家等,曾受邀采访和评审十多次。仅在国内的CSDN平台,博客文章浏览量超过2500万,拥有超过57万的粉丝。
1701
文章
0
问答
文章排行榜
最热
最新
相关电子书
更多
文娱运维技术
立即下载
《SaaS模式云原生数据仓库应用场景实践》
立即下载
《看见新力量:二》电子书
立即下载