| 查看: 1805 | 回复: 13 | ||
| 本帖产生 1 个 数学EPI ，点击这里进行查看 | ||
hui470新虫 (小有名气)
|
[求助]
MATLAB编程计算一个简单的数学式子。 已有1人参与
|
|
|
计算式子如图，其中输入参数为一个n*m的实矩阵，u(ij)表示矩阵第i行第j列的元素。 恳求各位大神帮忙编写一个MATLAB程序。在这里先谢谢各位了！！！  NN@$3[4}R)WN6[89WNCL285.png |
回复此楼» 猜你喜欢
336求调剂
已经有3人回复
-
286求调剂
已经有3人回复
-
一志愿华中科技大学071000，求调剂
已经有3人回复
-
279分求调剂 一志愿211
已经有14人回复
-
332求调剂
已经有3人回复
-
0856材料专硕353求调剂
已经有3人回复
-
材料学学硕080502 337求调剂-一志愿华中科技大学
已经有5人回复
-
085601调剂 358分
已经有4人回复
-
0856 271求调剂
已经有3人回复
-
能源材料化学课题组招收硕士研究生8-10名
已经有15人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- python之蟒开启理论计算化学的新时代
已经有22人回复
- 【分享】统计学学习顺序-统计学科技树
已经有6人回复
bluesine
铁杆木虫 (职业作家)
科苑小木虫
- 数学EPI: 5
- 应助: 132 (高中生)
- 贵宾: 1.991
- 金币: 9834.6
- 散金: 89
- 红花: 19
- 帖子: 3619
- 在线: 377.5小时
- 虫号: 869544
- 注册: 2009-10-12
- 性别: GG
- 专业: 数学物理
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与，应助指数 +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
bluesine
铁杆木虫 (职业作家)
科苑小木虫
- 数学EPI: 5
- 应助: 132 (高中生)
- 贵宾: 1.991
- 金币: 9834.6
- 散金: 89
- 红花: 19
- 帖子: 3619
- 在线: 377.5小时
- 虫号: 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 |
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
bluesine
铁杆木虫 (职业作家)
科苑小木虫
- 数学EPI: 5
- 应助: 132 (高中生)
- 贵宾: 1.991
- 金币: 9834.6
- 散金: 89
- 红花: 19
- 帖子: 3619
- 在线: 377.5小时
- 虫号: 869544
- 注册: 2009-10-12
- 性别: GG
- 专业: 数学物理
|
不好意思，连乘符号看漏了。。。下面补上： 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
- 金币: 9834.6
- 散金: 89
- 红花: 19
- 帖子: 3619
- 在线: 377.5小时
- 虫号: 869544
- 注册: 2009-10-12
- 性别: GG
- 专业: 数学物理
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
