fun2.m

来自「随着蜂窝移动通信市场的巨大发展」· M 代码 · 共 18 行

function r=fun2(a)                             %求信干比
x=10;                                               %初始化信干比  

P=0.15898;
%x=fsolve('(40*x+1)*exp(-1/2*x)-1',x);           
%求解纳什均衡点的信干比

y(1)=fsolve('40*x*exp(-0.5*x)*((1-exp(-0.5*x)))^(79)-(1-exp(-0.5*x))^80-2.9564e+010*0.00025',x);%x是信干比
y(2)=fsolve('40*x*exp(-0.5*x)*((1-exp(-0.5*x)))^(79)-(1-exp(-0.5*x))^80-8.0052e+009*0.00305',x);
y(3)=fsolve('40*x*exp(-0.5*x)*((1-exp(-0.5*x)))^(79)-(1-exp(-0.5*x))^80-3.6996e+009*0.01311',x);
y(4)=fsolve('40*x*exp(-0.5*x)*((1-exp(-0.5*x)))^(79)-(1-exp(-0.5*x))^80-2.1824e+009*0.03693',x);
y(5)=fsolve('40*x*exp(-0.5*x)*((1-exp(-0.5*x)))^(79)-(1-exp(-0.5*x))^80-1.4457e+009*0.08247',x);
y(6)=fsolve('40*x*exp(-0.5*x)*((1-exp(-0.5*x)))^(79)-(1-exp(-0.5*x))^80-1.0441e+009*0.15898',x);
y(7)=fsolve('40*x*exp(-0.5*x)*((1-exp(-0.5*x)))^(79)-(1-exp(-0.5*x))^80-7.7475e+008*0.27692',x);
y(8)=fsolve('40*x*exp(-0.5*x)*((1-exp(-0.5*x)))^(79)-(1-exp(-0.5*x))^80-6.11e+008*0.44784',x);
y(9)=fsolve('40*x*exp(-0.5*x)*((1-exp(-0.5*x)))^(79)-(1-exp(-0.5*x))^80-4.8899e+008*0.68434',x); 

r=y;                                                %返回值