%?Program?
%
%?enfor.m
%
%?Simulation?program?to?realize?enforcement?algorithm
%
%?Programmed?by?Kang
%

%%%%%%%%%%%%%%%%?preparation?part?%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clc;
clear;
M?=?4;???????%信道數
T?=?400;???????%迭代次數
S?=?8;???????%第二用戶數
e=2.71828;
sum=0;????????%存放和
%ch?=?zeros（1S）;?????????%每步信道選擇矩陣
n0?=0.01;???????????%白噪聲干擾
r?=?5;??????%覆蓋半徑
ts=20;
p=0.5;%power?!!!!!!!!
u?=?zeros（ST）;?%?存放各信道各次迭代的結果
w?=?zeros（ST）;??%信道i在t時刻對用戶的吸引值
f=zeros（ST）;%
h=zeros（SSM）;??
ww?=?zeros（MST）;
alpha?=?10;?????????%信道衰落因子
fsel=[1?2?3?4];
cq=[1?1.1?1.2?1.3];
ps=0.4;%parament
%%%%%%%%%%%%%%%%%????初???始???化??????%%%%%%%%%%%%%%%%%%%%%%%%%

for??x=1:S
????for?y=1:M???
????????f（1:21）=fsel（1）;
????????f（3:41）=fsel（2）;
????????f（5:61）=fsel（3）;
????????f（7:81）=fsel（4）;
????????w（:1）=3.0;%？？？？幾，任何信道下，用戶的w都初始化1
????????ww（1:1）=1.0;
????????ww（2:1）=1.5;%在信道2上，任何SU的效益是1.5
????????ww（3:1）=2.0;
????????ww（4:1）=2.5;
????????u（:1）=1.0;%?????自己設定幾
????end
end
pos_tx=position（rS）;??????%創建發射節點圖產生5個發射節點
pos_rv=position（rS）;??????%創建接收節點圖，產生5個接收節點
dis=distance（pos_txpos_rvS）;??????????%創建各發射到接收的距離矩陣
H=loss（disalphaS）;?%H變化區間大
H1=1./（（exp（-pi*（（H）.^2）））*10+10）*10;%H1元素的值比較接近??二維的
for?ii=1:M
????h（::ii）=H1.*cq（ii）;%4維的
end?????%各信道下各用戶間衰落，信道增益Hij
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%????畫圖?????????????%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
figure（）?????%畫網絡圖
for?ii=1:S
????x1=pos_tx（ii1）;
????y1=pos_tx（ii2）;
????x2=pos_rv（ii1）;
????y2=pos_rv（ii2）;???
????
xx=[x1x2];
yy=[y1y2];

plot（xxyy‘--bo‘‘LineWidth‘2...
????????????????‘MarkerEdgeColor‘‘k‘...
????????????????‘MarkerFaceColor‘‘w‘...
????????????????‘MarkerSize‘5）
????????????hold?on
??????????????????legend（‘認知用戶‘）
end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%????算???法????%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
temp=0;
%sir_node_sum=0;
flag=-1000000000;
for?t=2:T
????for?x=1:S%次級用戶數
??????for?y=1:M%信道數
???????if?y==f（xt-1）%!!!!w與信道有關的吧
???????????ww（yxt）=（1-ps）*w（xt-1）+ps*u（xt-1）;
???????else
???????????ww（yxt）=ww（yxt-1）;
???????end
???????if?ww（yxt）>flag??%?find?max
???????????flag=ww（yxt）;
???????????wmax=ww（yxt）;%?????寫維數嗎ww是暫時存儲量
???????????cha=y;
???????end
????end
??????f（xt）=cha;
??????w（xt）=wmax;%用戶1經歷4個信道?選出4個ww中最大的給w
???%[chawmax]=sir（ftMSh）;%1!!!!!!
???flag=-1000000000;?%完成在時刻t某個用戶遍歷4個信道后的選擇
????end
????
??for?x=1:S
????for?a=1:S
??????????if?a~=x??????%a不是用戶m
?????????????if?f（at）==f（xt）?????%a選的信道和當前的策略一致
?????????????????temp=temp+p*h（axf（at））;%h與*無關，是個數???????????????????????