24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

fluenter

银虫 (正式写手)

应助: 4 (幼儿园)
金币: 131.9
散金: 222
红花: 2
帖子: 399
在线: 172.2小时
虫号: 1754779
注册: 2012-04-14
性别: MM
专业: 工程热力学

[求助] 查找矩阵中符合条件的元素，并赋值给另外一个矩阵。

[a,b]=meshgrid(1:10,-5:10)
c=a+b
d=a+b+c
查找d>10的元素，
并将相应的a，b，c复制给e，f，g。

百度上有个答案，但是我验证一下，还是不对。
http://zhidao.baidu.com/question/337920323.html
必须用线性索引，行列下标索引的话有问题。举一个你容易理解的例子x(1:3,1:3)取的是x第1到3行，第1到3列的3*3矩阵，而不是x(1,1),x(2,2),x(3,3)这三个值。因此你上面X(m,n)是一个范围内的矩阵，而不是你要的大于3小于12的X值。这时候需要用sub2din把[m,n]转成线性索引:

X=[3,2,6;3,7,2;8,9,13]
ind = find(X>3&X<12)
[m,n] = ind2sub(size(X),ind)
X(m,n)

n =

   1
   2
   2
   3

ans =

   8    9    9 13
   3    7    7    2
   8    9    9 13
   3    2    2    6

[ Last edited by fluenter on 2012-5-11 at 15:42 ]

回复此楼

» 猜你喜欢

2025年遐想已经有4人回复
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球已经有8人回复
自然科学基金委宣布启动申请书“瘦身提质”行动已经有4人回复
求个博导看看已经有18人回复

» 本主题相关价值贴推荐，对您同样有帮助:

fortran能求带字母矩阵的本征值吗？（矩阵元素带有字母，能得到形式解吗）已经有5人回复
请高手指点一下怎么确定一个符号矩阵中的等值元素已经有6人回复
一个矩阵问题的推导已经有7人回复
一个矩阵的正交补矩阵怎么求啊已经有9人回复
对称正定矩阵中绝对值大的元素在逆矩阵中相应元素的绝对值也大？已经有5人回复
请教一个矩阵求逆问题已经有4人回复
【讨论】求助一个矩阵公式是否正确已经有6人回复
【整理】今天整理硬盘上的资料，偶然发现的好玩的小东西就直接发上来吧已经有105人回复
【交流】请教一下各位高手关于量子力学中的矩阵力学的一个问题已经有4人回复
【转载】p(素数)个元素的有限域上的n阶可逆矩阵的个数是多少?如何考虑? 已经有4人回复

每天坚持一件事情

1楼 2012-05-11 15:10:25

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

xiegangmai

版主 (职业作家)

我没头衔

仿真EPI: 3
应助: 157 (高中生)
贵宾: 5.735
金币: 40609
散金: 2394
红花: 63
沙发: 4
帖子: 4046
在线: 1330.2小时
虫号: 653607
注册: 2008-11-13
性别: GG
专业: 机械动力学
管辖: 仿真模拟

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与，应助指数 +1
fluenter: 金币+5, ★★★★★最佳答案, 谢谢 2012-05-14 09:33:54
fluenter: 回帖置顶 2012-05-29 16:42:43

不是很理解你最终要得到什么，（1）e, f, g是与a, b, c相同维度的矩阵呢，（2）还是仅符合d>10条件元素组成的矩阵？

（1）给e,f,g初始化为0

CODE:

clear

[a,b] = meshgrid(1:10,-5:10);

c=a+b ;

d=a+b+c;

[ m, n ] = find( d> 10 );

e = zeros( size( a ) );

f = zeros( size( a ) );

g = zeros( size( a ) );

for i = 1 : length( m )

    e( m(i), n(i) ) = a( m(i), n(i) );

    f( m(i), n(i) ) = b( m(i), n(i) );

    g( m(i), n(i) ) = c( m(i), n(i) );

end

（2）e,f,g只能为一维数组（一行或一列）

CODE:

clear

[a,b] = meshgrid(1:10,-5:10);

c=a+b ;

d=a+b+c;

[ m, n ] = size( a );

ind = find( d > 10 );

a = reshape(a, m * n, 1 );

b = reshape(b, m * n, 1 );

c = reshape(c, m * n, 1 );

e = a( ind );

f = b( ind );

g = c( ind );

赞一下(1人)

回复此楼

明德厚学、求是创新

4楼2012-05-14 00:03:10

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主 fluenter 的主题更新

返回列表