24小时热门版块排行榜    

查看: 1947  |  回复: 4

张益坤

新虫 (初入文坛)

[求助] Matlab技术在某个圆内的点的个数 已有3人参与

请看图片,谢谢大神教教我,怎么做,或者我要怎么去学习。新手很迷茫

Matlab技术在某个圆内的点的个数
QF2Z1W]BYPAYQR6WIGG~}J8.png
回复此楼

» 收录本帖的淘帖专辑推荐

规范

» 猜你喜欢

科研有我,我在科研
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
计算出每个点到圆心的距离,小于半径的视为在园内,反之则园外。
2楼2015-11-23 20:38:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

onelightwang

金虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
2楼给出了基本思路,判断各点是否在圆内,然后每个点作为圆心循环一遍即可。
不过这样循环计算量略大,若数据点数为N,基本要判断N*(N-1)次
当N很大时,计算会比较费时间。

楼主可以通过向量化来提高计算效率。
假设从数据库中提取的圆心坐标序列分别为X,Y。(X,Y都为列向量,假设长度为n)
1.构建一个矩阵【 (x_i-x_j)^2 +(y_i-y_j)^2】_nxn,通过与 r^2的差值,就可以轻易的判断某个点在半径r内的点数。

2,目标矩阵的构造可简化成【(x_i-x_j)^2 】_nxn 矩阵的构建,因为后半部分类似

3.【(x_i-x_j)^2 】_nxn的矩阵可分解为[(x_i)^2  +  (x_j)^2   -2 *(x_i)*(x_j)]_nxn

4.继续分解成3个n x n矩阵,并分别实现
           [(x_i)^2 ]_nxn 的实现:X.^2*ones(size(X,2),size(X,1))
           [(x_j)^2 ]_nxn 的实现:(X.^2*ones(size(X,2),size(X,1)))',把上面的转置一下即可
         【 -2 *(x_i)*(x_j)]_nxn的实现:-2*X*X'
ok,任务完成
楼主按照4到1的顺序编程,就可以看到符合要求的点分布了,最后再用一个sum函数,就可以得到需要的结果饿了。
避免过多循环,可以提高运行速度!
3楼2015-11-23 23:15:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

快乐时光lk

新虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
还有一个方法,跟楼上各位都不同,就是把这数据转换成两张栅格矩阵,第一张圆变成一个面状的矩阵,圆内部设为1,外部设为0。第二张图,点直接转换成1,其他区域转换成0,把两个矩阵相点乘.* 得到的结果图,我们只需统计其中1的个数即为点数

发自小木虫IOS客户端
4楼2015-11-24 01:38:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

快乐时光lk

新虫 (小有名气)

选哪个都要根据数据的特点来选,如果圆是由公式确定的,用上面的的方法很好,但若是不知道圆心或半径或者是shp文件可以考虑我这个方法

发自小木虫IOS客户端
5楼2015-11-24 01:42:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 张益坤 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085600材料与化工 +7 安全上岸! 2026-03-16 7/350 2026-03-20 20:37 by zhukairuo
[考研] 一志愿南昌大学,327分,材料与化工085600 +7 Ncdx123456 2026-03-19 7/350 2026-03-20 20:02 by zhukairuo
[考研] 求调剂 +3 eation27 2026-03-20 3/150 2026-03-20 19:32 by JourneyLucky
[考研] 材料与化工专硕调剂 +7 heming3743 2026-03-16 7/350 2026-03-20 19:31 by zhukairuo
[考研] 能源材料化学课题组招收硕士研究生8-10名 +5 脱颖而出 2026-03-16 14/700 2026-03-20 09:30 by kkcoco25
[考研] 307求调剂 +9 冷笙123 2026-03-17 9/450 2026-03-19 22:44 by 学员8dgXkO
[考研] 288求调剂 +15 于海海海海 2026-03-19 15/750 2026-03-19 22:41 by 学员8dgXkO
[考研] 304求调剂 +6 司空. 2026-03-18 6/300 2026-03-18 23:03 by 星空星月
[考研] 085601专硕,总分342求调剂,地区不限 +5 share_joy 2026-03-16 5/250 2026-03-18 14:48 by haxia
[考研] 收复试调剂生 +4 雨后秋荷 2026-03-18 4/200 2026-03-18 14:16 by elevennnne
[考研] 0854,计算机类招收调剂 +3 胡辣汤放糖 2026-03-15 6/300 2026-03-18 12:09 by 上岸上岸……..
[考研] 0703化学336分求调剂 +6 zbzihdhd 2026-03-15 7/350 2026-03-18 09:53 by zhukairuo
[考研] 有没有道铁/土木的想调剂南林,给自己招师弟中~ +3 TqlXswl 2026-03-16 7/350 2026-03-17 15:23 by TqlXswl
[考博] 26申博 +4 八6八68 2026-03-16 4/200 2026-03-17 13:00 by 轻松不少随
[考研] 11408 一志愿西电,277分求调剂 +3 zhouzhen654 2026-03-16 3/150 2026-03-17 07:03 by laoshidan
[基金申请] 今年的国基金是打分制吗? 50+3 zhanghaozhu 2026-03-14 3/150 2026-03-16 17:07 by 北京莱茵润色
[考研] 277材料科学与工程080500求调剂 +3 自由煎饼果子 2026-03-16 3/150 2026-03-16 14:10 by 运气yunqi
[考研] 327求调剂 +6 拾光任染 2026-03-15 11/550 2026-03-15 22:47 by 拾光任染
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +3 嘉年新程 2026-03-15 3/150 2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 289求调剂 +4 这么名字咋样 2026-03-14 6/300 2026-03-14 18:58 by userper
信息提示
请填处理意见