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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 数学 » 计算数学 » MATLAB编程计算一个简单的数学式子。
141/2返回列表12下一页
查看: 2065  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 时间戳又变了8-15 +6 archvillain 2026-08-15 10/500 2026-08-15 14:29 by 枫叶之巅
[基金申请] 有时候，自然基金真的不能太认真 （我的申报经验） +7 majunge000 2026-08-11 9/450 2026-08-15 13:34 by kiz2011
[基金申请] 关于Filecode分析方法 +10 majunge000 2026-08-10 13/650 2026-08-15 12:26 by hanpeng972
[考博] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 HFw0lei2R37i 2026-08-14 5/250 2026-08-15 11:41 by KxMI1BYBxWX1
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜？ +6 Ldrop2023 2026-08-13 6/300 2026-08-15 10:09 by tuanggou
[教师之家] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +5 7lpolszZVXgi 2026-08-14 6/300 2026-08-15 09:28 by sunzitan
[考博] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +4 k0dTPqJtl0jt 2026-08-14 5/250 2026-08-15 04:45 by 4wMiSEwB6436
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 k0dTPqJtl0jt 2026-08-14 4/200 2026-08-15 01:28 by 4wMiSEwB6436
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 22/1100 2026-08-14 23:45 by Noways
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +23 Lanmanbaby 2026-08-09 37/1850 2026-08-14 13:40 by Equinoxhua
[基金申请] filecode +15 documentary 2026-08-10 17/850 2026-08-14 10:08 by kissu88
[基金申请] 我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人 +10 家与远方 2026-08-10 15/750 2026-08-14 02:08 by 绵羊哥哥
[基金申请] 静等基金结果 +8 gjjjzhong 2026-08-10 21/1050 2026-08-13 17:56 by 且听虎啸
[基金申请] 重要来源：本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] 2019年青年基金涵评意见，大家看看几个A，几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 确定了，国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] 关于filecode，很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
信息提示
关闭
请填处理意见
关闭