24小时热门版块排行榜    

查看: 2036  |  回复: 13
本帖产生 1 个 数学EPI ,点击这里进行查看

hui470

新虫 (小有名气)

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

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

MATLAB编程计算一个简单的数学式子。
NN@$3[4}R)WN6[89WNCL285.png
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖置顶 ( 共有1个 )

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
板凳要做十年冷文章不发一个字
4楼2015-11-18 10:13:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

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
板凳要做十年冷文章不发一个字
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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 【2027博士申请】纳米药物递送方向 20+3 13586093586 2026-08-03 5/250 2026-08-08 21:54 by 13586093586
[基金申请] 关于豆爷回答的JTJC与%2F数量 +5 yang182083 2026-08-06 7/350 2026-08-08 18:29 by zhanghaozhu
[基金申请] 国基金的申报应该改成非等额制,评价高的钱多评价低的钱少,但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 17:27 by oEVWOejN9taj
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 17:19 by oEVWOejN9taj
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 qTvzQBRAHjCi 2026-08-07 5/250 2026-08-08 17:02 by oEVWOejN9taj
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 qTvzQBRAHjCi 2026-08-07 5/250 2026-08-08 16:47 by oEVWOejN9taj
[基金申请] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 Vi50GxzrFcSG 2026-08-07 5/250 2026-08-08 16:27 by oEVWOejN9taj
[论文投稿] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 Vi50GxzrFcSG 2026-08-07 6/300 2026-08-08 15:07 by oEVWOejN9taj
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 WQyTGMbfH7Cx 2026-08-07 4/200 2026-08-08 15:02 by oEVWOejN9taj
[基金申请] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 4/200 2026-08-08 14:27 by oEVWOejN9taj
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 KXLV3nuBVBY7 2026-08-07 5/250 2026-08-08 14:07 by oEVWOejN9taj
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 娱乐 +6 Tide man 2026-08-03 6/300 2026-08-07 22:40 by 铁帽子农民
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] 听说今天filecode变了 +24 布布和一二 2026-08-06 47/2350 2026-08-07 16:02 by zhiyanjiang
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[基金申请] 影响面上的因素 +8 布布和一二 2026-08-05 11/550 2026-08-06 10:41 by 宝贝虫子
信息提示
请填处理意见