| 查看: 243 | 回复: 3 | |||
| 当前主题已经存档。 | |||
quartzbj金虫 (正式写手)
|
[交流]
【讨论】矩阵部分数赋值问题？
|
||
|
一个矩阵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 ] |
回复此楼» 猜你喜欢
Sci. Bull. 悲剧经验
已经有5人回复
-
国自然评分
已经有4人回复
-
今年审到国自然15份，谈谈感受
已经有19人回复
-
提交了我也来说说感想
已经有4人回复
-
面上本子正文33页，违规吗？会被低分嘛？
已经有11人回复
-
评审有感
已经有22人回复
-
反应很差，大量原料没有反应
已经有6人回复
-
如果工作了想读博，可以边工作边读全日制嘛？
已经有5人回复
-
如何实现卤原子转化
已经有8人回复
-
河北省自然科学基金
已经有3人回复
zjliu1265
至尊木虫 (正式写手)
- 应助: 18 (小学生)
- 金币: 12147.2
- 散金: 200
- 帖子: 642
- 在线: 679.9小时
- 虫号: 585320
- 注册: 2008-08-08
- 性别: GG
- 专业: 光学
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
kuhailangyu(金币+2,VIP+0):谢谢参与 11-18 11:11
小木虫(金币+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 ] |
2楼2009-11-18 11:02:57
quartzbj
金虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 1811.5
- 散金: 250
- 帖子: 444
- 在线: 301.2小时
- 虫号: 312177
- 注册: 2006-12-30
- 专业: 应用数学方法
3楼2009-11-18 18:46:46
hitzhang
木虫 (正式写手)
- 仿真EPI: 1
- 应助: 0 (幼儿园)
- 贵宾: 2.15
- 金币: 1376.7
- 散金: 969
- 红花: 8
- 帖子: 862
- 在线: 226.4小时
- 虫号: 390575
- 注册: 2007-06-02
- 性别: GG
- 专业: 无机非金属类电介质与电解
4楼2009-11-18 21:06:50
5