| 查看: 1807 | 回复: 13 | |||
| 本帖产生 1 个 数学EPI ，点击这里进行查看 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
hui470新虫 (小有名气)
|
[求助]
MATLAB编程计算一个简单的数学式子。 已有1人参与
|
||
|
计算式子如图，其中输入参数为一个n*m的实矩阵，u(ij)表示矩阵第i行第j列的元素。 恳求各位大神帮忙编写一个MATLAB程序。在这里先谢谢各位了！！！  NN@$3[4}R)WN6[89WNCL285.png |
回复此楼» 猜你喜欢
求助
已经有4人回复
-
0805材料320求调剂
已经有3人回复
-
297求调剂
已经有11人回复
-
268求调剂
已经有9人回复
-
一志愿西安交通大学材料工程专业 282分求调剂
已经有10人回复
-
278求调剂
已经有7人回复
-
26考研一志愿中国石油大学(华东)305分求调剂
已经有5人回复
-
280求调剂
已经有9人回复
-
279求调剂
已经有5人回复
-
286求调剂
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- python之蟒开启理论计算化学的新时代
已经有22人回复
- 【分享】统计学学习顺序-统计学科技树
已经有6人回复
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
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
