| 查看: 2252 | 回复: 13 | ||
| 本帖产生 1 个 数学EPI ，点击这里进行查看 | ||
hui470新虫 (小有名气)
|
[求助]
MATLAB编程计算一个简单的数学式子。 已有1人参与
|
|
|
计算式子如图，其中输入参数为一个n*m的实矩阵，u(ij)表示矩阵第i行第j列的元素。 恳求各位大神帮忙编写一个MATLAB程序。在这里先谢谢各位了！！！  NN@$3[4}R)WN6[89WNCL285.png |
回复此楼» 猜你喜欢
考博求助
已经有3人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有5人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有9人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有10人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有6人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有9人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- python之蟒开启理论计算化学的新时代
已经有22人回复
- 【分享】统计学学习顺序-统计学科技树
已经有6人回复
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与，应助指数 +1
hui470: 金币+100, ★★★★★最佳答案 2015-11-18 18:06:08
Edstrayer: 数学EPI+1 2016-01-02 00:28:32
感谢参与，应助指数 +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 |
4楼2015-11-18 10:13:16
|
完整代码： 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 |
9楼2015-11-19 08:44:16
hui470
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1607.5
- 散金: 100
- 红花: 1
- 帖子: 234
- 在线: 57.7小时
- 虫号: 3540436
- 注册: 2014-11-16
- 专业: 数理统计
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
- 专业: 数理统计
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
|
不好意思，连乘符号看漏了。。。下面补上： 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
8楼2015-11-19 08:38:38
getengqing
木虫 (正式写手)
- 数学EPI: 1
- 应助: 104 (高中生)
- 金币: 3759.8
- 散金: 30
- 红花: 5
- 帖子: 592
- 在线: 110.5小时
- 虫号: 2986702
- 注册: 2014-02-23
- 性别: GG
- 专业: 自然语言理解与机器翻译
10楼2015-11-19 09:10:44
