function zeta=ff(x,a,b,theta)%theta需要给定，你没说是怎么来的
n=length(x);%样本容量即24
P=zeros(n,1);%初始化概率向量
for i=1:n
P(i)=1/(1+exp(-1.7*a(i)*(theta-b(i))));
end
T=mean(P)*ones(n,1);%T中每个元素都相等是P的均值
zeta=(P-T)'*(P-x)/sqrt(P'*diag((P-T).^2)*(1-P));%我的理解是p+q=1
希望对你有帮助！
function zeta=ff(x,a,b,theta)
[m,n]=size(x);
zeta=zeros(m,1);
P=zeros(n,1);%初始化概率向量
for k=1:m
for i=1:n
P(i)=1/(1+exp(-1.7*a(i)*(theta(k)-b(i))));
end
T=mean(P)*ones(n,1);%T中每个元素都相等是P的均值
zeta(k)=(P-T)'*(P-x(k,:)')/sqrt(P'*diag((P-T).^2)*(1-P));%我的理解是p+q=1
end
先建立一个函数，存储上述代码
然后在命令窗口调用zeta=ff(x,a,b,theta)就可以了（x,a,b,theta这些值都已经存到matlab里了）
扫码加好友，拉您进群
全部版块
我的主页
收藏
