应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 专业学科区 » 数学 » 计算数学 » MATLAB编程计算一个简单的数学式子。
141/2返回列表12下一页
查看: 1863  |  回复: 13
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
本帖产生 1 个 数学EPI ，点击这里进行查看

hui470

新虫 (小有名气)

[求助] MATLAB编程计算一个简单的数学式子。 已有1人参与

计算式子如图，其中输入参数为一个n*m的实矩阵，u(ij)表示矩阵第i行第j列的元素。
恳求各位大神帮忙编写一个MATLAB程序。在这里先谢谢各位了！！！

MATLAB编程计算一个简单的数学式子。
NN@$3[4}R)WN6[89WNCL285.png
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2015-11-17 13:32:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

bluesine

铁杆木虫 (职业作家)

科苑小木虫

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与，应助指数 +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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bluesine

铁杆木虫 (职业作家)

科苑小木虫
完整代码：
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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

hui470

新虫 (小有名气)
输出结果是cd的平方就可以了
赞一下 回复此楼
2楼2015-11-17 13:40:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hui470

新虫 (小有名气)
自己顶下，求大神来解救！！！！！！
赞一下 回复此楼
3楼2015-11-18 00:49:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hui470

新虫 (小有名气)
引用回帖:
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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hui470

新虫 (小有名气)
引用回帖:
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-  ...

亲 这个程序哪里不对。。。。例如
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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bluesine

铁杆木虫 (职业作家)

科苑小木虫
引用回帖:
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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bluesine

铁杆木虫 (职业作家)

科苑小木虫
引用回帖:
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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

getengqing

木虫 (正式写手)
引用回帖:
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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 hui470 的主题更新
141/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 279学硕食品专业求调剂院校 20+4 孤独的狼爱吃羊 2026-04-12 18/900 2026-04-13 01:17 by 幸免 ..
[考研] 290调剂生物0860 +28 哇哈哈，。 2026-04-11 31/1550 2026-04-13 01:16 by 幸免 ..
[考研] 291求调剂 +8 关忆北. 2026-04-11 8/400 2026-04-12 09:32 by 逆水乘风
[考研] 材料与化工300求调剂 +39 肖开文 2026-04-09 43/2150 2026-04-12 01:30 by 秋豆菜芽
[考研] 求调剂，一志愿大连理工大学354分 +5 雨声余生 2026-04-11 6/300 2026-04-11 16:12 by 雨声余生
[考研] 农学0904 312求调剂 +6 Say Never 2026-04-10 6/300 2026-04-11 10:33 by wwj2530616
[考研] 22408 352分求调剂0854类 +4 努力的夏末 2026-04-09 4/200 2026-04-11 09:57 by zhq0425
[考研] 284求调剂 +12 archer.. 2026-04-10 13/650 2026-04-11 08:44 by zhq0425
[考研] 一志愿华南理工大学331分材料求调剂 +9 天下ww 2026-04-09 9/450 2026-04-10 22:58 by Ftglcn90
[考研] 求调剂 +5 不会飞的鱼@ 2026-04-10 5/250 2026-04-10 19:07 by chemisry
[考研] 071000生物学调剂求助 +17 zzzzwww 2026-04-09 20/1000 2026-04-10 15:55 by 求调剂zz
[考研] 考研二轮调剂 +8 故人?? 2026-04-09 8/400 2026-04-10 09:44 by 青梅duoduo
[考研] 材料化工总分334求调剂 +16 Riot2025 2026-04-08 17/850 2026-04-09 20:19 by maddjdld
[考研] 284求调剂 +7 让我上岸吧阿西 2026-04-09 7/350 2026-04-09 18:59 by haironglove
[考研] 复试调剂，一志愿郑州大学材料与化工289分 +31 硕星赴 2026-04-08 31/1550 2026-04-09 16:54 by Delta2012
[考研] 085400电子信息类（川大控制工程）求调剂可跨专业 求老师联系 +3 626776879 2026-04-08 3/150 2026-04-09 16:05 by 猪会飞
[考研] 一志愿西南大学生物学学硕344 求生物学相关调剂/生物与医药 +7 超人不会飞@ 2026-04-08 7/350 2026-04-09 09:35 by gong120082
[考研] 331求调剂 +5 luoxin0706. 2026-04-08 5/250 2026-04-08 22:15 by zhouyuwinner
[考研] 一志愿吉大化学327求调剂 +12 王王白石 2026-04-06 13/650 2026-04-08 16:05 by luoyongfeng
[考研] 085100建筑学 寻求跨专业调剂 一志愿南大294分 校级省级国家级奖项若干 踏实肯干 +3 1021075758 2026-04-06 4/200 2026-04-07 09:23 by 蓝云思雨
信息提示
关闭
请填处理意见
关闭