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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 数学 » 计算数学 » MATLAB编程计算一个简单的数学式子。
141/2返回列表12下一页
查看: 1806  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 301求调剂 +10 yy要上岸呀 2026-03-17 10/500 2026-03-21 03:14 by JourneyLucky
[考研] 一志愿武理材料305分求调剂 +6 想上岸的鲤鱼 2026-03-18 7/350 2026-03-21 01:03 by JourneyLucky
[考研] 304求调剂 +7 司空. 2026-03-18 7/350 2026-03-20 23:08 by JourneyLucky
[考研] 考研调剂求学校推荐 +3 伯乐29 2026-03-18 5/250 2026-03-20 22:59 by JourneyLucky
[考研] 317求调剂 +5 申子申申 2026-03-19 9/450 2026-03-20 22:26 by JourneyLucky
[考研] 329求调剂 +9 想上学吖吖 2026-03-19 9/450 2026-03-20 22:01 by luoyongfeng
[考博] 招收博士1-2人 +3 QGZDSYS 2026-03-18 3/150 2026-03-20 11:58 by 呱呱呱呱叫
[考研] 320求调剂0856 +3 不想起名字112 2026-03-19 3/150 2026-03-19 22:53 by 学员8dgXkO
[考研] 0703化学调剂 +4 18889395102 2026-03-18 4/200 2026-03-19 16:13 by 30660438
[考研] 328求调剂，英语六级551，有科研经历 +4 生物工程调剂 2026-03-16 12/600 2026-03-19 11:10 by 生物工程调剂
[考研] 311求调剂 +6 26研0 2026-03-15 6/300 2026-03-18 14:43 by haxia
[考研] 312求调剂 +8 陌宸希 2026-03-16 9/450 2026-03-18 12:39 by Linda Hu
[考研] 生物学071000 329分求调剂 +3 我爱生物生物爱� 2026-03-17 3/150 2026-03-18 10:12 by macy2011
[考博] 26博士申请 +3 1042136743 2026-03-17 3/150 2026-03-17 23:30 by 轻松不少随
[考研] 326求调剂 +5 上岸的小葡 2026-03-15 6/300 2026-03-17 17:26 by ruiyingmiao
[考研] 考研化学学硕调剂，一志愿985 +4 张vvvv 2026-03-15 6/300 2026-03-17 17:15 by ruiyingmiao
[考研] 085601求调剂 +4 Du.11 2026-03-16 4/200 2026-03-17 17:08 by ruiyingmiao
[考研] 275求调剂 +4 太阳花天天开心 2026-03-16 4/200 2026-03-17 10:53 by 功夫疯狂
[考研] 考研调剂 +3 淇ya_~ 2026-03-17 5/250 2026-03-17 09:25 by Winj1e
[考研] 0703化学调剂 290分有科研经历，论文在投 +7 腻腻gk 2026-03-14 7/350 2026-03-16 10:12 by houyaoxu
信息提示
关闭
请填处理意见
关闭