matlab_D2D中双中继仿真程序，双中继协作通信的功率改进方案

clc clear a=0:0.005:0.05; N=length(a); A=[];%线性不等式约束左边矩阵 b=[];%线性不等式约束右边向量 Aeq=[];%线性等式约束左边矩阵 beq=[];%线性等式约束右边向量 lb=zeros(1,3); ub=[35;35;35];%自变量上限 x0=[1;1;1];%初始值 for i=1:N options=optimset('LargeScale','off','display','iter'); [x,fval]=fmincon(@trelay,x0,[],[],[],[],lb,ub,@(x)myconH(x,a(i)),options); b1(i)=x(1); b2(i)=fval; [x,fval]=fmincon(@trelay,x0,[],[],[],[],lb,ub,@(x)myconAF(x,a(i)),options); c1(i)=x(1); c2(i)=fval; [x,fval]=fmincon(@trelay,x0,[],[],[],[],lb,ub,@(x)myconDF(x,a(i)),options); d1(i)=x(1); d2(i)=fval; end figure(1) plot(a,b1,'gv-'); hold on plot(a,c1,'ro:'); plot(a,d1,'ms-'); xlabel('中断概率'); ylabel('源节点功率'); legend(' HDAF源节点功率','AF源节点功率','DF源节点功率',1); axis([0,0.05,10,35]); set(gca,'XTick',[0:0.01:0.05], 'YTick',[10:5:35]); figure(2) plot(a,b2,'gv-'); hold on plot(a,c2,'ro:'); plot(a,d2,'ms-'); xlabel('中断概率'); ylabel('总功率'); legend(' HDAF总功率','AF总功率','DF总功率',1); axis([0.005,0.05,20,60]); set(gca,'XTick',[0.005:0.01:0.045],'YTick',[20:5:60]);