【频谱分配】基于频谱空洞预留算法的频谱分配算法的matlab仿真

已于 2022-03-29 17:29:36 修改
阅读量1.6k 收藏 12
点赞数 1
分类专栏： MATLAB 板块1:通信与信号处理 文章标签： 频谱空洞预留算法 频谱分配算法
于 2022-03-28 21:50:42 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/ccsss22/article/details/123807209
版权

1.软件版本

MATLAB2013a
2.本算法理论知识

        基于概率切换和最小化概率切换的频谱切换算法，由于过多的考虑了资源的合理利用，使得次用户的不切换概率并大于等于的情况，从而导致对应的次用户的服务质量显著下降，考虑到这个因素，文献[40]提出了一种基于频谱空洞预留机制的频谱分配算法，从而最大程度上保证这些次用户的服务质量。

        根据文献所述，当出现授权用户的时候，次用户必须暂停当前所占的行道，并切换到其他的频谱空洞上继续完成未完成的服务。然后，过多的频谱切换会带来很多切换时延，这样就很难保证次用户的服务质量，特别对于一些实时性较强的服务，影响更大。文献[40] 提出的频谱空洞的预留方案使不能保证服务质量的次用户直接使用预留的频谱空洞，不再进行频谱感知、频谱请求等过程，从而最大程度上降低频谱切换时延。

        根据文献所述，当空洞上进行通信时实际的不切换概率，如果 大于 ，说明用于通信的频谱空洞可以满足次用户的服务质量要求，不需要再给该次用户预留任何频谱空洞；否则，如果 ，说明次用户在实际获得的频谱空洞上进行通信时达不到自己的服务质量要求，很可能会出现中断服务等情况，需要为之预留频谱空洞，使它切换到预留频谱空洞上继续完成通信。由此可见，为次用户预留频谱空洞的条件为：

我们假设预留的频谱空洞时间至少可以刚好满足服务质量的要求，从而有：

 进一步可以得到：

 

 更进一步可以得到：

 最后可以得到：

       式子为预留因子，当预留因子为0的时候，说明不需要再为次用户的通信预留其他的频谱空洞，当预留因子为1的时候，表示在服务时间内都要在预留频谱空洞上进行通信。通过预留频谱空洞算法，就可以保证次用户通信过程中的服务质量要求，弥补频谱空洞匹配算法的不足，使次用户顺利完成自己的服务。

 3.核心代码

clc;
clear;
close all;
warning off;
RandStream.setDefaultStream(RandStream('mt19937ar','seed',1));
addpath 'My_Function\' %自定义函数路径
addpath 'My_Figure\'   %仿真结果图保存路径
addpath 'My_Result\'   %仿真结果及仿真中间结果数值保存路径
%参数化
lemda   = 0.2;
M       = 50;%用户数量
%信道数目
N       = 10;
lemda   = 1/2000 + (1/200-1/2000)*rand(N,1);%1/2000 ~ 1/200的均匀分布

Treqi   = zeros(M,1);
PQi     = zeros(M,1);
PQi2    = zeros(M,1);
PHij    = zeros(M,N);
PHij2   = zeros(M,1);
betai   = zeros(M,1);


for nn = 1:M
    nn
    TIME    = 100;
    NUM     = 0;
    %中继信道lemda
    lemda2  =  1/1000 + (1/100-1/1000)*rand(N,1);%1/1000 ~ 1/100的均匀分布
    for t = 1:TIME
        for i = 1:nn
            %步骤1：对次用户业务服务质量所要求的PQi等进行初始化，即最小的不发生切换的概率
            PQi(i)  = rand(1,1)/2;

            %步骤2：次用户SUi的目标切换概率1-PQi
            PQi2(i) = 1 - PQi(i);

            %步骤3：计算次用户在不同频谱空洞上通信的实际切换概率Phij
            Treqi(i) = 10 + (300-10)*rand(1,1);
            for j = 1:N
                PHij(i,j) = 1-exp(-lemda(j)*Treqi(i));
            end

            %步骤4：根据Hi，找到切换概率和目标切换概率的差值最小的频谱空洞j
            for j = 1:N
                Hs1(j) = abs(PQi(i) - exp(-lemda(j)*Treqi(i)));
            end
            %考虑中继协作，CR用户和中继进行竞争，获得最佳的II值。
            for j = 1:N
                PHij(i,j) = 1-exp(-lemda2(j)*Treqi(i));
            end
            for j = 1:N
                Hs2(j) = abs(PQi(i) - exp(-lemda2(j)*Treqi(i)));
            end
            Hss = [Hs1,Hs2];
            [VV,II] = min(Hss);
            
            %步骤5：把频谱控制j分配给对应的次用户i
            lemdas   = [lemda;lemda2];
            H(i)     = Hss(II);%将最佳的空洞j分配给次用户i
            %步骤6：计算次用户SUi在频谱空洞j上通信的不切换概率
            PHij2(i) =1-exp(-lemdas(II)*Treqi(i));


            %步骤7：比较1-Phij和PQi的大小，1-Phij > PQi，beta=0，则不需要预留频谱空洞，1-Phij <=
            %PQi，beta不等于0，则需要频谱预留
            if PHij2(i) > PQi
               betai(i) = 0;
            else
               betai(i) = 1+log(PQi(i))/lemdas(II)/Treqi(i); 
            end
            %步骤8：计算对应的切换次数
            if betai(i) ~=0;
               NUM = NUM + 1;
            end    
        end
    end
    PRO(nn) = NUM/TIME;

end
hold on
plot(1:M,PRO,'k','Linewidth',2);
PRO1 = PRO;



Treqi   = zeros(M,1);
PQi     = zeros(M,1);
PQi2    = zeros(M,1);
PHij    = zeros(M,N);
PHij2   = zeros(M,1);
betai   = zeros(M,1);


for nn = 1:M
    nn
    TIME    = 100;
    NUM     = 0;
    for t = 1:TIME
        for i = 1:nn
            %步骤1：对次用户业务服务质量所要求的PQi等进行初始化，即最小的不发生切换的概率
            PQi(i)  = rand(1,1)/2;

            %步骤2：次用户SUi的目标切换概率1-PQi
            PQi2(i) = 1 - PQi(i);

            %步骤3：计算次用户在不同频谱空洞上通信的实际切换概率Phij
            Treqi(i) = 10 + (300-10)*rand(1,1);
            for j = 1:N
                PHij(i,j) = 1-exp(-lemda(j)*Treqi(i));
            end

            %步骤4：根据Hi，找到切换概率和目标切换概率的差值最小的频谱空洞j
            for j = 1:N
                Hs(j) = abs(PQi(i) - exp(-lemda(j)*Treqi(i)));
            end
            %步骤5：把频谱控制j分配给对应的次用户i
            [VV,II] = min(Hs);
            H(i) = Hs(II);%将最佳的空洞j分配给次用户i


            %下面的步骤引入频谱预留的思想，进行切换的改进
            %步骤6：计算次用户SUi在频谱空洞j上通信的不切换概率1-Phij
            PHij2(i) = 1-(1-exp(-lemda(II)*Treqi(i)));


            %步骤7：比较1-Phij和PQi的大小，1-Phij > PQi，beta=0，则不需要预留频谱空洞，1-Phij <=
            %PQi，beta不等于0，则需要频谱预留
            if PHij2(i) > PQi
               betai(i) = 0;
            else
               betai(i) = 1+log(PQi(i))/lemda(II)/Treqi(i); 
            end
            %步骤8：计算对应的切换次数
            if betai(i) ~=0;
               NUM = NUM + 1;
            end    
        end
    end
    PRO(nn) = NUM/TIME;

end
hold on
plot(1:M,PRO,'b','Linewidth',2);
PRO2 = PRO;

Treqi   = zeros(M,1);
PQi     = zeros(M,1);
PQi2    = zeros(M,1);
PHij    = zeros(M,N);
PHij2   = zeros(M,1);
betai   = zeros(M,1);
for nn = 1:M
    nn
    TIME    = 100;
    NUM     = 0;
    for t = 1:TIME
        for i = 1:nn
            %步骤1：对次用户业务服务质量所要求的PQi等进行初始化，即最小的不发生切换的概率
            PQi(i)  = rand(1,1)/2;

            %步骤2：次用户SUi的目标切换概率1-PQi
            PQi2(i) = 1 - PQi(i);

            %步骤3：计算次用户在不同频谱空洞上通信的实际切换概率Phij
            Treqi(i) = 10 + (300-10)*rand(1,1);
            for j = 1:N
                PHij(i,j) = 1-exp(-lemda(j)*Treqi(i));
            end

            %步骤4：根据Hi，找到切换概率和目标切换概率的差值最小的频谱空洞j
            for j = 1:N
                Hs(j) = abs(PQi(i) - exp(-lemda(j)*Treqi(i)));
            end
            %步骤5：把频谱控制j分配给对应的次用户i
            [VV,II] = min(Hs);
            H(i) = Hs(II);%将最佳的空洞j分配给次用户i
            %步骤6：计算次用户SUi在频谱空洞j上通信的切换概率Phij
            PHij2(i) = 1-(1-exp(-lemda(II)*Treqi(i)));
            %步骤6：计算对应的切换次数
            NUM = NUM+PHij2(i);   
        end
    end
    PRO(nn) = NUM/TIME;
end
hold on
plot(1:M,PRO,'r','Linewidth',2);
PRO3 = PRO;
xlabel('次用户个数');
ylabel('切换次数');
legend('中继协作基于最小概率的频谱空洞预留算法','基于最小概率的频谱空洞预留算法','基于最小概率的频谱空洞匹配算法');

save result1.mat M PRO1 PRO2 PRO3

4.操作步骤与仿真结论

 

5.参考文献

[01]Wei Dang Lu, Student Member, IEEE, Yi Gong, Senior Member, IEEE, See Ho Ting, Member, IEEE, Xuan Li Wu, Member, IEEE, and Nai Tong Zhang. Cooperative OFDM Relaying for Opportunistic Spectrum Sharing: Protocol Design and Resource Allocation. IEEE TRANSACTIONS ON WIRELESS COMMUNICATIONS, VOL. 11, NO. 6, JUNE 2012

[02]Li-Chun Wang, Fellow, IEEE, Chung-Wei Wang, Student Member, IEEE, and Chung-Ju Chang, Fellow, IEEE. Modeling and Analysis for Spectrum Handoffs in Cognitive Radio Networks IEEE TRANSACTIONS ON MOBILE COMPUTING, VOL. 11, NO. 9, SEPTEMBER 2012

[03]Y. Shi and Y. T. Hou. Optimal power control for multi-hop software defined radio networks[R]. IEEE INFOCOM 2007.2007:1694-1702.

A01-84

6.完整源码获得方式

方式1：微信或者QQ联系博主

方式2：订阅MATLAB/FPGA教程，免费获得教程案例以及任意2份完整源码

fpga和matlab
关注
  • 1
    点赞
  • 踩
  • 12
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
专栏目录
动态频谱分配算法
FPGA/MATLAB学习教程/源码/项目合作开发
02-16 5992
clc clear M=4; % Size of The constellation sayan=0; % Count...
【老生谈算法】matlab实现频谱分析算法源码——频谱分析
m0_53407570的博客
10-20 1159
这类方法对线谱（正弦信号的谱）最合适，对检测噪声下的正弦信号很有效，特别是低信噪比的情况。因此，矩阵可逆且方程一定有解。这类方法使用不同的方式来估计频谱：不是试图直接从数据中估计PSD，而是将数据建模成一个由白噪声驱动的线性系统的输出，并试图估计出该系统的参数。Burg法最首要的优势在于解决含有低噪声的间隔紧密的正弦信号，并且对短数据的估计，在这种情况下AR功率谱密度估计非常逼近与真值。Spectral estimation（谱估计）的目标是基于一个有限的数据集合描述一个信号的功率（在频率上的）分布。
全球各地5G频谱分配情况如何？最全无线通信频率分配表
热门推荐
成功不必在我，而功力必不唐捐！
12-12 2万+
先看看无线电信号的频谱如何划分：   1、 5G NR    3GPP已指定5G NR 支持的频段列表，5G NR频谱范围可达100GHz，指定了两大频率范围： ①Frequency range 1 （FR1）：就是我们通常讲的6GHz以下频段 •频率范围：450MHz - 6.0GHz •最大信道带宽100MHz ②Frequency range 2 （FR2）：就...
RSA整理--频谱路由分配算法
MineralterMan的专栏
09-06 8267
0.背景介绍 随着互联网的高速时代来临，智能设备、消费类电子等对于数据流量的告诉需求，迫使数据带宽呈几何级数增长。据Cisco公司2015年1月份的财务报告指出，2016年，全球IP流量将达到每年1.1皆字节，即月均88.7艾字节（10亿吉字节[GB]）。到2020年，全球IP流量将达到每年2.3皆字节，月均194艾字节。全球IP流量将在今后5年增长近3倍。 图1.Cisco固网业务流量规
Matlab中的RandStream的应用场景
zhshow的专栏
11-26 8013
一，matlab中生成随机数主要有三个函数：rand, randn,randi 1，rand  生成均匀分布的伪随机数。分布在（0~1）之间   主要语法：rand(m,n)生成m行n列的均匀分布的伪随机数             rand(m,n,'double')生成指
频谱细化matlab程序,频谱细化-----Zoom-FFT算法介绍及MATLAB实现
weixin_39716877的博客
03-18 614
栅栏效应若长度为N=T/dt=T⋅fsN=T/dt=T\cdot {{f}_{s}}N=T/dt=T⋅fs​，其中TTT为采样时间，dtdtdt为采样间隔，fsf_sfs​为采样频率，通过傅里叶变换后得到X(fi)X\left( {{f}_{i}} \right)X(fi​)，fi=i⋅fs/N,i=0,1,2,...,N/2{{f}_{i}}=i\cdot {{f}_{s}}/N,i=0,1,2...
基于遗传算法的认知无线电频谱分配算法
04-11
基于遗传算法的认知无线电频谱分配算法，matlab代码
基于CVX工具箱凸优化的动态频谱分配算法matlab仿真,含仿真操作录像
03-29
1.版本：matlab2021a，包含...3.内容：基于CVX工具箱凸优化的动态频谱分配算法matlab仿真，包括CVX工具箱的调用和使用。 4.注意事项：注意MATLAB左侧当前文件夹路径，必须是程序所在文件夹位置，具体可以参考视频录。
基于遗传算法和粒子群算法的认知无线电频谱分配算法
08-22
针对认知无线电中空闲频谱资源的最优分配问题，分别采用了遗传算法和粒子群算法进行求解。该代码是利用遗传算法和粒子群算法求解该问题的matlab仿真代码。
CR系统中基于博弈论的频谱分配算法及仿真
01-31
在分析了无线电频谱管理现状及频谱资源紧张的前提下，引出频谱管理中的认知无线电技术，同时提出了基于博弈论的电磁频谱管理的模型化方法，深入地分析了博弈论算法的收敛性及效用函数的选择，并且通过Matlab仿真以...
动态频谱分配算法代码
08-09
动态频谱分配算法代码，写的真的非常好，希望大家使用，很好很好很好很好1DFRelayPA
认知无线电频谱分配算法
03-08
认知无线电的静态频谱分配的一种匹配博弈的算法 难得一见哦 是我师兄自己编的哦
matlab频谱分析代码-WavePart:Matlab（r）函数集，用于将定向波谱分配到其风和不同的膨胀分量。首先使用分水岭定义算法识别分区
05-24
matlab频谱分析代码WavePart Matlab（r）函数集，用于将定向波谱分配到其风和不同的膨胀分量。 分区首先使用分水岭定义算法进行标识，然后主要按照Hanson和Phillips（2001）中描述的方法进行修改。 作者： 道格拉斯·卡尔（Douglas Cahl）和乔治·沃加利斯（George Voulgaris） 地球，海洋与环境学院南卡罗来纳大学美国南卡罗来纳州哥伦比亚市29205。 代码更新： 2020年1月22日-waveparams.m-平均频率（fm）估算不正确； 它已得到纠正。 2020年1月22日-更新了partition.m，以捕获将平坦频谱作为输入时的情况。 2020年1月22日-master_partition.m已更新，因此它以正确的输出数量调用waveparams.m。 2020年1月21日-waveparams.m-Hrms估计不正确； 现在的平均方向在-180至+180度范围内。 引文： Douglas，C.和Voulgaris，G.，2019，WavePART：用于定向海浪频谱划分的MATLAB（r）软件。：Zenodo，doi：10.528
国内无线频谱资源分配图集[图]
08-09
国内无线频谱资源分配图集。
全球频谱资源（9K~30GHz)分配详情高清图
08-01
做通信相关的童鞋需要了解，关于0~30GHz在全球的分布情况详情，哪些频段用于什么功能一张图显示的很清楚。
认知无线电中基于干扰图模型的动态频谱分配的MATLAB仿真
03-29
认知无线电中基于干扰图模型的动态频谱分配。 M=5; %频谱池中信道的个数 N=Nc/2; %发送节点的个数 r=10.^(15/20); %最小SINR值，r=15dB a=4; %衰减因子 V=[rand(N,1)*10,rand(N,1)*10]; %产生10*10范围随机分布的N个...
基于GA优化的竞价博弈频谱分配算法的matlab仿真
FPGA/MATLAB学习教程/源码/项目合作开发
10-26 1274
目前，现有的关于频谱分配的方法主要有基于竞价的分配方法，基于博弈论的分配方法。再比如，有2个频率资源，5个用户，在发送业务类似的情况下，那么一个频率分配2个用户，一个频率分配3个用户的排队时间，肯定优于1个资源分配一个用户，其余一个频率分配四个用户。：用户的需求情况进行分配，如果当前用户的通信属于实时性需求较强的情况，则优先分配频谱资源（这里涉及到具体的业务模型，比如视频，电话等各种业务）而对于实时性要求较低的情况（如短信，文字信息等各种业务），则可以分配同一频率，做类似同频不同时的通信。
matlab设置随机流的种子
weixin_30432579的博客
11-09 2372
matlab2007之后改变了设置方式: RandStream.setDefaultStream(RandStream('mt19937ar','seed',sum(100*clock))); normrnd(mean,var); 如此就可取到不同初始值的正态随机数 之前的方式类似于 normrnd('seed', tseed) ; 转载于:https://www.cnblogs....
【老生谈算法】matlab信号频谱仿真源码.doc
11-05
MATLAB 信号频谱仿真源码详解 本文档提供了 MATLAB 信号频谱仿真源码，涵盖了信号时域波形和频谱图的生成。通过对代码的解读，我们可以了解到信号处理的基本概念和 MATLAB 编程技术。 信号时域波形 信号时域波形...
matlab写一个基于二机制遗传算法的频谱资源分配的仿真代码
最新发布
05-25
以下是一个基于二进制遗传算法的频谱资源分配的 MATLAB 仿真代码： ```matlab clc; clear; close all; %% 参数设置 % 频带宽度 bw = 20e6; % 子载波数量 nSubcarriers = 512; % 子载波带宽 subcarrierBw = bw / ...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

fpga和matlab

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值