基于博弈论的认知无线电无线资源分配MATLAB程序

基于博弈论的认知无线电无线资源分配MATLAB程序

clc 
clear all 
 
B=15; 
BER=10^（-4）; 
K=1.5/（logm（0.2/BER））; 
SNR=[810]; 
k=log2（1+K.*SNR）; 
r1=10; 
r2=12; 
x=0; 
y=1; 
z=1; 
% b=[00]; 
% c=x+y*（sum（b））; 
% p=r.*k.*b-b.*c; 
 
a（1）=0.001; 
for i=1:300 
    a（2）=0.001; 
    for j=1:300 
        eigenvalue（1）=1-10/3*（2*a（1）*k（1）-a（1）*k（2）+2*a（2）*k（2）-a（2）*k（1））+10/3*（（a（1）*（k（2）-2*k（1））*a（2）*（k（1）-2*k（2））+（2*a（2）*k（2）-a（2）*k（1）-2*a（1）*k（1）+a（1）*k（2））^2）^（0.5））; 
        eigenvalue（2）=1-10/3*（2*a（1）*k（1）-a（1）*k（2）+2*a（2）*k（2）-a（2）*k（1））-10/3*（（a（1）*（k（2）-2*k（1））*a（2）*（k（1）-2*k（2））+（2*a（2）*k（2）-a（2）*k（1）-2*a（1）*k（1）+a（1）*k（2））^2）^（0.5））; 
        if （eigenvalue（1）<1）&&（eigenvalue（1）>-1）&&（eigenvalue（2）<1）&&（eigenvalue（2）>-1） 
            a（2）=a（2）+0.002; 
        end 
    end 
    if a（2）>0 
        a（2）=a（2）-0.002; 
    else 
        a（2）=0; 
    end 
    f2（i）=a（2）; 
    f1（i）=a（1）; 
    a（1）=a（1）+0.002; 
end 
 
figure（1） 
plot（f1f2‘--r‘） 
axis（[00.2500.25]） 
hold on 
 
B=15; 
BER=10^（-4）; 
K=1.5/（logm（0.2/BER））; 
SNR=[1010]; 
k=log2（1+K.*SNR）; 
r1=10; 
r2=12; 
x=0; 
y=1; 
z=1; 
% b=[00]; 
% c=x+y*（sum（b））; 
% p=r.*k.*b-b.*c; 
 
a（1）=0.001; 
for i=1:300 
    a（2）=0.001; 
    for j=1:300 
        eigenvalue（1）=1-10/3*（2*a（1）*k（1）-a（1）*k（2）+2*a（2）*k（2）-a（2）*k（1））+10/3*（（a（1）*（k（2）-2*k（1））*a（2）*（k（1）-2*k（2））+（2*a（2）*k（2）-a（2）*k（1）-2*a（1）*k（1）+a（1）*k（2））^2）^（0.5））; 
        eigenvalue（2）=1-10/3*（2*a（1）*k（1）-a（1）*k（2）+2*a（2）*k（2）-a（2）*k（1））-10/3*（（a（1）*（k（2）-2*k（1））*a（2）*（k（1）-2*k（2））+（2*a（2）*k（2）-a（2）*k（1）-2*a（1）*k（1）+a（1）*k（2））^2）^（0.5））; 
        if （eigenvalue（1）<1）&&（eigenvalue（1）>-1）&&（eigenvalue（2）<1）&&（eigenvalue（2）>-1） 
            a（2）=a（2）+0.002; 
        end 
    end 
    if a（2）>0 
        a（2）=a（2）-0.002; 
    else 
        a（2）=0; 
    end 
    f2（i）=a（2）; 
    f1（i）=a（1）; 
    a（1）=a（1）+0.002; 
end 
 
figure（1） 
plot（f1f2‘-.*b‘） 
axis（[00.2500.25]） 
 
B=15; 
BER=10^（-4）; 
K=1.5/（logm（0.2/BER））; 
SNR=[108]; 
k=log2（1+K.*SNR）; 
r1=10; 
r2=12; 
x=0; 
y=1; 
z=1; 
% b=[00]; 
% c=x+y*（sum（b））; 
% p=r.*k.*b-b.*c; 
 
a（1）=0.001; 
for i=1:300 
    a（2）=0.001; 
    for j=1:300 
        eigenvalue（1）=1-10/3*（2*a（1）*k（1）-a（1）*k（2）+2*a（2）*k（2）-a（2）*k（1））+10/3*（（a（1）*（k（2）-2*k（1））*a（2）*（k（1）-2*k（2））+（2*a（2）*k（2）-a（2）*k（1）-2*a（1）*k（1）+a（1）*k（2））^2）^（0.5））; 
        eigenvalue（2）=1-10/3*（2*a（1）*k（1）-a（1）*k（2）+2*a（2）*k（2）-a（2）*k（1））-10/3*（（a（1）*（k（2）-2*k（1））*a（2）*（k（1）-2*k（2））+（2*a（2）*k（2）-a（2）*k（1）-2*a（1）*k（1）+a（1）*k（2））^2）^（0.5））; 
        if （eigenvalue（1）<1）&&（eigenvalue（1）>-1）&&（eigenvalue（2）<1）&&（eigenvalue（2）>-1） 
            a（2）=a（2）+0.002; 
        end 
    end 
    if a（2）>0 
        a（2）=a（2）-0.002; 
    else 
        a（2）=0; 
    end 
    f2（i）=a（2）; 
    f1（i）=a（1）; 
    a（1）=a（1）+0.002; 
end 
 
figure（1） 
plot（f1f2‘-g‘） 
axis（[00.2500.25]） 
 
 
B=15; 
BER=10^（-4）; 
K=1.5/（logm（0.2/BER））; 
SNR=[77]; 
k=log2（1+K.*SNR）; 
r1=10; 
r2=12; 
x=0; 

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        4921  2020-09-23 07:56  cognitive_gaming.m