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