%%迭代功率选择，运用博弈论循环迭代功率分配策略，使系统性能达到最优
clear all
clc
t0=clock;
beta=[0.01,0.015];
m1=2;m2=1;
r1=0.2;r2=0.1;
alpha=1;
MAX_ITER=30;
beta=0.01;
syms x
b=[0,0];
pai=[0,0];
history_pai=pai'
history_b=b';
for j=1:MAX_ITER
for i=1:2
b(1)=2+(1/3)*b(1)+(1/3)*b(2);
b(2)=(1/3)*(b(1)+b(2));
pai(1)=(100/9)*(b(1)+b(2))*(4-b(1)+b(2));
pai(2)=(100/9)*(b(1)+b(2)-1)*(b(1)-b(2)-1);
end
history_b=[history_b b'];
history_pai=[history_pai pai'];
end
b1=history_b(1,:);
b2=history_b(2,:);
pai1=history_pai(1,:);
pai2=history_pai(2,:);
figure(1)
plot(0:MAX_ITER,b1,'-+g');hold on;
plot(0:MAX_ITER,b2,'-ob');hold off;
xlabel('重复博弈次数');
ylabel('信道价格');
figure(3)
plot(0:MAX_ITER,pai1,'-+g');hold on;
plot(0:MAX_ITER,pai2,'-ob');hold off;
xlabel('重复博弈次数');
ylabel('系统利润');
MAX_ITER=30;
beta=0.015;
syms x
b=[0,0];
pai=[0,0];
history_pai=pai'
history_b=b';
for j=1:MAX_ITER
for i=1:2
b(1)=3+0.5*b(2);
b(2)=0.5*b(1);
pai(1)=(100/3)*(1+0.5*b(2))*(4-b(1)+b(2));
pai(2)=(100/3)*(0.5*b(1)-1)*(b(1)-b(2)-1);
end
history_b=[history_b b'];
history_pai=[history_pai pai'];
end
b1=history_b(1,:);
b2=history_b(2,:);
pai1=history_pai(1,:);
pai2=history_pai(2,:);
figure(2)
plot(0:MAX_ITER,b1,'-+g');hold on;
plot(0:MAX_ITER,b2,'-ob');hold off;
xlabel('重复博弈次数');
ylabel('信道价格');
figure(4)
plot(0:MAX_ITER,pai1,'-+g');hold on;
plot(0:MAX_ITER,pai2,'-ob');hold off;
xlabel('重复博弈次数');
ylabel('系统利润');
time=etime(clock,t0)
