赞
踩
如图,我的右侧纵坐标轴刻度重复了,研究半天,调不出来,求大神指导哈。
代码如下:
clc
clear
syms a;
VA=5;
V=VA+1;
epsilon=0.05;
alpha=0.2;
beta=0.95;
RR=[];
BB=[];
LL=[];
for L=0:1:100
T=10^(-alpha*L/10);
kai=epsilon+(1-T)/T;
b=a*sqrt(T/(2*(kai+1)));
F=V^2*(1-2*T)+2*T+T^2*(V+kai)^2;
D=T*(V*kai+1);
x=(T*(V+kai)+1)/(T*(kai+1)+1);
y1=sqrt(0.5*(F+sqrt(F^2-4*D^2)));
y2=sqrt(0.5*(F-sqrt(F^2-4*D^2)));
[IAB,Sb]=het_b(x,y1,y2);
lambda=(T*(V*kai+1)+V)/(T*(V+kai)+1);
Sa=entropy_a(lambda);
R=beta*IAB-Sb+Sa;
f=matlabFunction((1./(sqrt(2*pi*VA)))*exp(-a.^2./(2*VA))*0.5*erfc(b));
Copyright © 2003-2013 www.wpsshop.cn 版权所有,并保留所有权利。