9楼: Originally posted by bluesine at 2015-11-19 08:44:16
完整代码：
U = ;
= size(U);

y = (13/12)^m;

for i=1:n
    tmpZ =1;
    for k=1:m
        tmp = abs(U(i,k)-1/2);
        tmpZ =tmpZ*(1+1/2*tmp-1/2*tmp^2);
    end
    y = y-2/n*tmpZ;
en ...

9楼: Originally posted by bluesine at 2015-11-19 08:44:16
完整代码：
U = ;
= size(U);

y = (13/12)^m;

for i=1:n
    tmpZ =1;
    for k=1:m
        tmp = abs(U(i,k)-1/2);
        tmpZ =tmpZ*(1+1/2*tmp-1/2*tmp^2);
    end
    y = y-2/n*tmpZ;
en ...

7楼: Originally posted by bluesine at 2015-11-19 08:35:29
不好意思，连乘符号看漏了。。。下面补上：

U = ;
= size(U);
y = (13/12)^m;
for i=1:n
    for k=1:m
        tmp = abs(U(i,k)-1/2);
        y = y -2/n*(1+1/2*tmp- 1/2*tmp^2);
    end
end

...