24小时热门版块排行榜

返回列表

当前主题已经存档。

quartzbj

金虫 (正式写手)

应助: 1 (幼儿园)
金币: 1783.5
散金: 250
帖子: 444
在线: 299.7小时
虫号: 312177
注册: 2006-12-30
专业: 应用数学方法

[交流] 【讨论】矩阵部分数赋值问题？

一个矩阵A：
4 5    16    17 18 9
2 20 30 40 9    0
3 20 8    6 7    0
31  40  18 8 19 6
1 83 2    6 7 18

希望把矩阵A大于10的数先挑出来赋值为1，其它为0，

第二步，但是如果同一行连续有两个1，就把原来矩阵A中相应数比较小的再赋值为0，不连续的点保持不变。

第三步，如果某个1周围8个点都为0，把原来矩阵A中对应数比较大的再赋值为1.

老是循环使用不好，怎么办？

[ Last edited by quartzbj on 2009-11-18 at 08:54 ]

回复此楼

» 猜你喜欢

到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗已经有7人回复
申请2026年博士已经有5人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有5人回复
寻求一种能扛住强氧化性腐蚀性的容器密封件已经有6人回复
2025冷门绝学什么时候出结果已经有7人回复
请问有评职称，把科研教学业绩算分排序的高校吗已经有6人回复
Bioresource Technology期刊，第一次返修的时候被退回好几次了已经有7人回复
请问哪里可以有青B申请的本子可以借鉴一下。已经有4人回复
请问下大家为什么这个铃木偶联几乎不反应呢已经有5人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复

1楼2009-11-17 21:19:14

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

zjliu1265

至尊木虫 (正式写手)

应助: 18 (小学生)
金币: 11778.2
散金: 200
帖子: 642
在线: 678.3小时
虫号: 585320
注册: 2008-08-08
性别: GG
专业: 光学

★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
kuhailangyu(金币+2,VIP+0):谢谢参与 11-18 11:11

第一步不用循环容易实现，但是后面两步不用循环比较麻烦。
其实循环也不是不好，我个人认为写程序首要的是写出来准确
的程序，其次才是去优化程序。

顺便问一下，你的问题是和元胞自动机相关么？

我写了一段程序如下，仅供参考：

A=[4,5,16,17,18,9;...
      2,20,30,40,9,0;...
      3,20,8,6,7,0;...
      31,40,18,8,19,6;...
      1,83,2,6,7,18];

% step 1
Ao=A;
A(A<=10)=0;
A(A>10)=1;
A
% step 2
[x,y]=find([A(1:end-1,

+A(2:end,

]>1.5);
Ix1=sub2ind(size(A),x,y);
Ix2=sub2ind(size(A),x,y+1);
Ix1,Ix2
[xx,Ix]=min([Ao(Ix1),Ao(Ix2)],[],1);
Ao(Ix1(Ix<1.5))=0;
Ao(Ix2(Ix>1.5))=0;

% step 3
S=A(1:end-2,2:end-1)+A(3:end,2:end-1)+...
A(2:end-1,1:end-2)+A(2:end-1,3:end)+...
A(1:end-2,1:end-2)+A(3:end,1:end-2)+...
A(1:end-2,3:end)+A(3:end,3:end);
[x,y]=find(S<0.5);
dx=[-1,0,1,-1,1,-1,0,1];
dy=[-1,-1,-1,0,0,1,1,1];
for k=1:length(x);
xn=x(k)+dx;
yn=y(k)+dy;
Ix=sub2ind(size(Ao),xn,yn);
[xx,Ik]=max(Ao(Ix));
Ao(Ix(Ik))=1;
end

[ Last edited by zjliu1265 on 2009-11-18 at 11:04 ]