【答案】应助回帖

» 猜你喜欢

请问哪里可以有青B申请的本子可以借鉴一下。已经有4人回复
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗已经有6人回复
孩子确诊有中度注意力缺陷已经有14人回复
三甲基碘化亚砜的氧化反应已经有4人回复
请问下大家为什么这个铃木偶联几乎不反应呢已经有5人回复
请问有评职称，把科研教学业绩算分排序的高校吗已经有5人回复
2025冷门绝学什么时候出结果已经有3人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有4人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复
AI论文写作工具：是科研加速器还是学术作弊器？已经有3人回复

» 本主题相关价值贴推荐，对您同样有帮助:

python之蟒开启理论计算化学的新时代已经有22人回复
【分享】统计学学习顺序-统计学科技树已经有6人回复

1楼2015-11-17 13:32:46

回帖置顶 ( 共有1个 )

bluesine

铁杆木虫 (职业作家)

科苑小木虫

数学EPI: 5
应助: 132 (高中生)
贵宾: 1.991
金币: 9496.3
散金: 89
红花: 19
帖子: 3534
在线: 363.8小时
虫号: 869544
注册: 2009-10-12
性别: GG
专业: 数学物理

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与，应助指数 +1
hui470: 金币+100, ★★★★★最佳答案 2015-11-18 18:06:08
Edstrayer: 数学EPI+1 2016-01-02 00:28:32

代码如下：下面例子以U= [1 2 3; 4 5 6; 7 8 9; 10 11 12]为例子计算得到的结果，其他值的时候只需要改版U即可
U = [1 2 3; 4 5 6; 7 8 9; 10 11 12];
[n,m] = 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
for i=1:n
for j=1:n
      for k=1:m
         tmp1 = U(i,k)-1/2;
         tmp2 = U(j,k)-1/2;
         tmp3 = tmp1 -tmp2;
         y = y + 1/n^2*(1+1/2*abs(tmp1)+1/2*abs(tmp2) -1/2*abs(tmp3));
      end
end
end
y

此例子结果
y =

  136.3964

赞一下(2人)

板凳要做十年冷文章不发一个字

4楼2015-11-18 10:13:16

回帖支持 ( 显示支持度最高的前 50 名 )

bluesine

铁杆木虫 (职业作家)

科苑小木虫

数学EPI: 5
应助: 132 (高中生)
贵宾: 1.991
金币: 9496.3
散金: 89
红花: 19
帖子: 3534
在线: 363.8小时
虫号: 869544
注册: 2009-10-12
性别: GG
专业: 数学物理

完整代码：
U = [0.16666667       0.83333333       0.50000000
0.50000000       0.16666667       0.83333333
0.16666667       0.16666667       0.16666667
0.50000000       0.50000000       0.50000000
0.83333333       0.16666667       0.50000000
0.16666667       0.50000000       0.83333333
0.50000000       0.83333333       0.16666667
0.83333333       0.50000000       0.16666667
0.83333333       0.83333333       0.83333333 ];
[n,m] = 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;
end

for i=1:n
for j=1:n
      tmpZ =1;
      for k=1:m
         tmp1 = 1/2*(U(i,k)-1/2);
         tmp2 = 1/2*(U(j,k)-1/2);
         tmp3 = (tmp1 -tmp2);
         tmpZ = tmpZ * (1+abs(tmp1)+abs(tmp2)-abs(tmp3));
      end
   y = y + 1/n^2*tmpZ;
end
end
y
结果为：
y =

0.0330

赞一下(3人)

板凳要做十年冷文章不发一个字

9楼2015-11-19 08:44:16

普通回帖

hui470

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 1607.5
散金: 100
红花: 1
帖子: 234
在线: 57.7小时
虫号: 3540436
注册: 2014-11-16
专业: 数理统计

输出结果是cd的平方就可以了

2楼2015-11-17 13:40:00

hui470

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 1607.5
散金: 100
红花: 1
帖子: 234
在线: 57.7小时
虫号: 3540436
注册: 2014-11-16
专业: 数理统计

自己顶下，求大神来解救！！！！！！

3楼2015-11-18 00:49:27

hui470

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 1607.5
散金: 100
红花: 1
帖子: 234
在线: 57.7小时
虫号: 3540436
注册: 2014-11-16
专业: 数理统计

引用回帖:

4楼: Originally posted by bluesine at 2015-11-18 10:13:16
代码如下：下面例子以U= 为例子计算得到的结果，其他值的时候只需要改版U即可
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- ...

非常感谢！！！

5楼2015-11-18 18:07:00

hui470

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 1607.5
散金: 100
红花: 1
帖子: 234
在线: 57.7小时
虫号: 3540436
注册: 2014-11-16
专业: 数理统计

引用回帖:

亲这个程序哪里不对。。。。例如
0.16666667 0.83333333 0.50000000
0.50000000 0.16666667 0.83333333
0.16666667 0.16666667 0.16666667
0.50000000 0.50000000 0.50000000
0.83333333 0.16666667 0.50000000
0.16666667 0.50000000 0.83333333
0.50000000 0.83333333 0.16666667
0.83333333 0.50000000 0.16666667
0.83333333 0.83333333 0.83333333
这个矩阵的输出应该是：0.033034
这个程序输出的是一个负数。。。。。。。。。。。。。。。。。

6楼2015-11-18 19:38:01

bluesine

铁杆木虫 (职业作家)

科苑小木虫

数学EPI: 5
应助: 132 (高中生)
贵宾: 1.991
金币: 9496.3
散金: 89
红花: 19
帖子: 3534
在线: 363.8小时
虫号: 869544
注册: 2009-10-12
性别: GG
专业: 数学物理

引用回帖:

6楼: Originally posted by hui470 at 2015-11-18 19:38:01
亲这个程序哪里不对。。。。例如
0.16666667 0.83333333 0.50000000
0.50000000 0.16666667 0.83333333
0.16666667 0.16666667 0.16666667
0.50000000 0.50000000 0.50000000
0.83333333 0.166 ...

不好意思，连乘符号看漏了。。。下面补上：

U = [1 2 3; 4 5 6; 7 8 9; 10 11 12];
[n,m] = 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

for i=1:n
for j=1:n
   tmpZ =1;
      for k=1:m
         tmp1 = 1/2*(U(i,k)-1/2);
         tmp2 = 1/2*(U(j,k)-1/2);
         tmp3 = 1/2*(tmp1 -tmp2);
         tmpZ = tmpZ * (1+abs(tmp1)+abs(tmp2)-abs(tmp3));
      end
   y = y + 1/n^2*tmpZ;
end
end
y

板凳要做十年冷文章不发一个字

7楼2015-11-19 08:35:29

bluesine

铁杆木虫 (职业作家)

科苑小木虫

数学EPI: 5
应助: 132 (高中生)
贵宾: 1.991
金币: 9496.3
散金: 89
红花: 19
帖子: 3534
在线: 363.8小时
虫号: 869544
注册: 2009-10-12
性别: GG
专业: 数学物理

引用回帖:

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

...

tmp3 = 1/2*(tmp1 -tmp2);
改成tmp3 = (tmp1 -tmp2);

板凳要做十年冷文章不发一个字

8楼2015-11-19 08:38:38

getengqing

木虫 (正式写手)

数学EPI: 1
应助: 104 (高中生)
金币: 3759.8
散金: 30
红花: 5
帖子: 592
在线: 110.5小时
虫号: 2986702
注册: 2014-02-23
性别: GG
专业: 自然语言理解与机器翻译

引用回帖:

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 ...

大神啊

一起交流学习/分享优秀资源

10楼2015-11-19 09:10:44