| 查看: 2058 | 回复: 13 | |||
| 本帖产生 1 个 数学EPI ，点击这里进行查看 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
hui470新虫 (小有名气)
|
[求助]
MATLAB编程计算一个简单的数学式子。 已有1人参与
|
||
|
计算式子如图，其中输入参数为一个n*m的实矩阵，u(ij)表示矩阵第i行第j列的元素。 恳求各位大神帮忙编写一个MATLAB程序。在这里先谢谢各位了！！！  NN@$3[4}R)WN6[89WNCL285.png |
回复此楼» 猜你喜欢
澳大利亚 Murdoch University 全奖博士招生（3个名额）地质化工冶金领域
已经有24人回复
-
阿姨
已经有4人回复
-
一作与独作在应聘高校教师时区别大吗
已经有4人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有9人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有18人回复
-
静等基金结果
已经有21人回复
-
FileCode能看出啥？
已经有28人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有6人回复
-
filecode
已经有16人回复
-
重要来源：本周末出结果
已经有10人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- python之蟒开启理论计算化学的新时代
已经有22人回复
- 【分享】统计学学习顺序-统计学科技树
已经有6人回复
bluesine
至尊木虫 (职业作家)
科苑小木虫
- 数学EPI: 5
- 应助: 132 (高中生)
- 贵宾: 1.991
- 金币: 10420.7
- 散金: 89
- 红花: 19
- 帖子: 3802
- 在线: 382.1小时
- 虫号: 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
- 金币: 10420.7
- 散金: 89
- 红花: 19
- 帖子: 3802
- 在线: 382.1小时
- 虫号: 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
