matlab中涉及的一些算法编程，必重谢

我的已知数据为x值，a值，b值，n=24，求拟合指标ζ值。求matlab编程，希望能用 %注释一下。
必当感激不急，以后你们有事，我也全力帮忙，如果赏的币不够，我会再赏。
备注：我是文科生，正在写论文，但是苦于对编程一窍不通，无能为力，所以才求救。并不是因为懒，希望你们能理解。

zeta =

    0.4563
    0.4563
    0.4563
    0.4563
    0.4563
    0.4563
    0.4563
    0.4563
   -0.1539
   -1.5875
    0.4563
   -0.1118
   -1.1677
    0.0564
    0.4563
    0.1884
   -0.4723
    0.6212
    0.0959
   -1.1677
    0.5055
    0.2120
   -0.4723
    0.7908
    0.2120
   -0.4723
    1.3881
    0.4563
    0.4563
   -1.4840
   -0.1539
   -1.5549
   -0.8179
   -1.0113
   -0.2330
   -0.9725
   -0.2883
   -0.9870
   -0.3768
   -0.9870
    2.7957
    0.4563
   -1.5875
   -1.4034
   -0.2298
   -0.5636
    0.3273
    1.1315
    0.4563
   -0.6471
   -1.0825
    0.1453
    0.1878
   -1.6112
    0.4259

自己顶一下吧

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
希望对你有帮助！

秋果55 发表于 2014-7-17 13:27
function zeta=ff(x,a,b,theta)%theta需要给定，你没说是怎么来的
n=length(x);%样本容量即24
P=zeros(n, ...

首先谢谢你的回复，但是当我输入到Matlab中，程序出错了，提示如下，请帮我分析一下，谢谢了

??? function zeta=ff(x,a,b,theta)%theta需要给定，你没说是怎么来的
    |
Error: Function definitions are not permitted in this context.

你是直接复制到命令窗口的吗，要建函数的，你提供一下x,a,b,theta 的数据

你好 数据如附件 请查收

秋果55 发表于 2014-7-20 21:40
你是直接复制到命令窗口的吗，要建函数的，你提供一下x,a,b,theta 的数据

谢谢你的回复，数据附件在下面我的回复里面。

貌似要用非线性拟合的函数啊

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里了）