有关中继程序，出图，并附有对应文章

有关中继的程序，并有对应的文章说明，出图。关于中断概率的仿真图

clc
clear
a=0:0.005:0.05;
N=length（a）;
A=[];%线性不等式约束左边矩阵
b=[];%线性不等式约束右边向量
Aeq=[];%线性等式约束左边矩阵
beq=[];%线性等式约束右边向量
lb=zeros（13）; 
ub=[35;35;35];%自变量上限
x0=[1;1;1];%初始值

for i=1:N
      options=optimset（‘LargeScale‘‘off‘‘display‘‘iter‘）;
      [xfval]=fmincon（@trelayx0[][][][]lbub@（x）myconH（xa（i））options）;
      b1（i）=x（1）;
      b2（i）=fval;
    
      [xfval]=fmincon（@trelayx0[][][][]lbub@（x）myconAF（xa（i））options）;
      c1（i）=x（1）;
      c2（i）=fval;
    
      [xfval]=fmincon（@trelayx0[][][][]lbub@（x）myconDF（xa（i））options）;
      d1（i）=x（1）;
      d2（i）=fval;
    
end
figure（1）
plot（ab1‘gv-‘）;
hold on
plot（ac1‘ro:‘）;
plot（ad1‘ms-‘）;
xlabel（‘中断概率‘）;
ylabel（‘源节点功率‘）;
legend（‘ HDAF源节点功率‘‘AF源节点功率‘‘DF源节点功率‘1）;
axis（[00.051035]）;
set（gca‘XTick‘[0:0.01:0.05] ‘YTick‘[10:5:35]）;

figure（2）
plot（ab2‘gv-‘）;
hold on
plot（ac2‘ro:‘）;
plot（ad2‘ms-‘）;
xlabel（‘中断概率‘）;
ylabel（‘总功率‘）;
legend（‘ HDAF总功率‘‘AF总功率‘‘DF总功率‘1）;
axis（[0.0050.052060]）;
set（gca‘XTick‘[0.005:0.01:0.045]‘YTick‘[20:5:60]）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1193  2016-11-03 09:49  双中继\双中继仿真\main.m

     文件        148  2016-10-09 13:48  双中继\双中继仿真\myconAF.m

     文件        170  2016-10-17 09:28  双中继\双中继仿真\myconDF.m

     文件        145  2016-10-09 14:15  双中继\双中继仿真\myconH.m

     文件         39  2016-10-09 11:27  双中继\双中继仿真\trelay.m

     文件     542564  2016-09-10 15:37  双中继\双中继协作通信的功率分配改进方案_王丽元.pdf

     目录          0  2018-06-21 16:00  双中继\双中继仿真

     目录          0  2018-06-21 16:00  双中继

----------- ---------  ---------- -----  ----

               544259                    8