h**********c
发帖数: 4120 | 1 写了一篇小文章，computational geometry方面的，
哪位牛人肯给看看，point to ellipse distance: a binary search approach,
有matlab的code
我老板不是搞这个方向的，这纯属于个人利用假期搞的，忙了一春假。
请牛人看看够不够发篇文章，latex 一共２１页。
如果能发在权威刊物最好，不行发哪个的大学的学报上也行。我现在很需要发文。
节选若干。
愿意帮忙请留言或站内方式联系。 |
g****t
发帖数: 31659 | 2 你的方法比牛顿法好？
写了一篇小文章，computational geometry方面的，
哪位牛人肯给看看，point to ellipse distance: a binary search approach,
有matlab的code
我老板不是搞这个方向的，这纯属于个人利用假期搞的，忙了一春假。
请牛人看看够不够发篇文章，latex 一共２１页。
如果能发在权威刊物最好，不行发哪个的大学的学报上也行。我现在很需要发文。
节选若干。
愿意帮忙请留言或站内方式联系。
【在 h**********c 的大作中提到】
: 写了一篇小文章，computational geometry方面的，
: 哪位牛人肯给看看，point to ellipse distance: a binary search approach,
: 有matlab的code
: 我老板不是搞这个方向的，这纯属于个人利用假期搞的，忙了一春假。
: 请牛人看看够不够发篇文章，latex 一共２１页。
: 如果能发在权威刊物最好，不行发哪个的大学的学报上也行。我现在很需要发文。
: 节选若干。
: 愿意帮忙请留言或站内方式联系。
|
h**********c
发帖数: 4120 | 3 感谢回帖。
首先说，牛顿法本人没有试过，对于四阶多项式求导比较麻烦，可不可导也是个问题。
在哪里收敛哪里不收敛。
求出来根，虚根怎么办，四阶可能会有四个根，取哪一个。
我的方法就是全是简单的乘积加和，开方，比较麻烦的地方就是每次iteration要求一
次一个3x3矩阵的determinant，不过这还是简单的乘积加和运算. 讨论degnerate case
要更简单一些，程序上code比较好写。
对于数学上的方法比方构造一个polynomial quadratic convergent 超出本人能力。从
计算几何讲，离散化，近似，叠代，看起来是可行的。
【在 g****t 的大作中提到】
: 你的方法比牛顿法好？
:
: 写了一篇小文章，computational geometry方面的，
: 哪位牛人肯给看看，point to ellipse distance: a binary search approach,
: 有matlab的code
: 我老板不是搞这个方向的，这纯属于个人利用假期搞的，忙了一春假。
: 请牛人看看够不够发篇文章，latex 一共２１页。
: 如果能发在权威刊物最好，不行发哪个的大学的学报上也行。我现在很需要发文。
: 节选若干。
: 愿意帮忙请留言或站内方式联系。
|
g****t
发帖数: 31659 | 4 我不是做你这行的,
但是从计算的经验上来讲,我觉得这种东西,首先要试的肯定就是牛顿法.
你查查文献看吧.多半别人都做过了.另外四个根什么的不成问题.
你查查一元多项式牛顿法的文献看.收敛到哪个根是初始点决定的.
初始点的选择多元主要靠猜,一元则有定理可以用.
另外四次一元多项式,其实是有解析解的,剩下就是开平方,开立方这类东西的
快速算法,这些也都有文献.为什么it is desired 不解方程求解这个问题?
所以,你肯定你这个题目是个合适的研究题目吗?
感谢回帖。
首先说，牛顿法本人没有试过，对于四阶多项式求导比较麻烦，可不可导也是个问题。
在哪里收敛哪里不收敛。
求出来根，虚根怎么办，四阶可能会有四个根，取哪一个。
我的方法就是全是简单的乘积加和，开方，比较麻烦的地方就是每次iteration要求一
次一个3x3矩阵的determinant，不过这还是简单的乘积加和运算. 讨论degnerate case
要更简单一些，程序上code比较好写。
对于数学上的方法比方构造一个polynomial quadratic converge
【在 h**********c 的大作中提到】
: 感谢回帖。
: 首先说，牛顿法本人没有试过，对于四阶多项式求导比较麻烦，可不可导也是个问题。
: 在哪里收敛哪里不收敛。
: 求出来根，虚根怎么办，四阶可能会有四个根，取哪一个。
: 我的方法就是全是简单的乘积加和，开方，比较麻烦的地方就是每次iteration要求一
: 次一个3x3矩阵的determinant，不过这还是简单的乘积加和运算. 讨论degnerate case
: 要更简单一些，程序上code比较好写。
: 对于数学上的方法比方构造一个polynomial quadratic convergent 超出本人能力。从
: 计算几何讲，离散化，近似，叠代，看起来是可行的。
|
B********e
发帖数: 10014 | 5 test your code, compare your results with old ones
confident? then just go ahead to submit
show 2pages here on bbs, who the hell can make decision for you?
【在 h**********c 的大作中提到】
: 写了一篇小文章，computational geometry方面的，
: 哪位牛人肯给看看，point to ellipse distance: a binary search approach,
: 有matlab的code
: 我老板不是搞这个方向的，这纯属于个人利用假期搞的，忙了一春假。
: 请牛人看看够不够发篇文章，latex 一共２１页。
: 如果能发在权威刊物最好，不行发哪个的大学的学报上也行。我现在很需要发文。
: 节选若干。
: 愿意帮忙请留言或站内方式联系。
|
g****t
发帖数: 31659 | 6 nod,如果比文献上的其他的方法好,
那就去投稿就可以了.
实践是第一位的.
test your code, compare your results with old ones
confident? then just go ahead to submit
show 2pages here on bbs, who the hell can make decision for you?
【在 B********e 的大作中提到】
: test your code, compare your results with old ones
: confident? then just go ahead to submit
: show 2pages here on bbs, who the hell can make decision for you?
|
h**********c
发帖数: 4120 | 7 牛顿法也许很好，但是但是用来解方程的。
它不能用来生成一个椭圆或者类似的convex shape，也就是不是一个表达椭圆的办法。
在计算几何里,convex hull是不连续的，没法求导的。而convex hull是可以用2维矢量
空间的一个generating set 来近似。
当然本贴只是想引起相关方向上的人士的兴趣。如果愿意的话我可以把２０页的pdf发
过去。
my point may be niave. 但是，您必须证明您是serious的.
严格讲，计算几何都是numerical guy这点我们是应该承认的。加减乘除越简单越好。
另外，我没有找到求椭圆到点距离的现成文章或程序例题。再找找，如果哪位有，请提
供出处。
我没有发过任何文章，不太清楚过程，还望赐教。
【在 g****t 的大作中提到】
: nod,如果比文献上的其他的方法好,
: 那就去投稿就可以了.
: 实践是第一位的.
:
: test your code, compare your results with old ones
: confident? then just go ahead to submit
: show 2pages here on bbs, who the hell can make decision for you?
|
c*******h
发帖数: 1096 | 8 The distance from a point (x,y) to an ellipse (x/a)^2 + (y/b)^2 = 1 is
ya^2*sin(t) + xb^2*cos(t)
- -----------------------------
a^2*sin^2(t) + b^2*cos^2(t)
where t satisfies
x^2+y^2-a^2-b^2
sin(2t+s) = ---------------------------------------
sqrt[ 4x^2y^2 + (b^2-y^2-a^2+x^2)^2 ]
and s satisfies
2xy
cos(s) = ---------------------------------------
sqrt[ 4x^2y^2 + (b^2-y^2-a^2+x^2)^2 ]
b^2-y^2-a^2+
【在 h**********c 的大作中提到】
: 写了一篇小文章，computational geometry方面的，
: 哪位牛人肯给看看，point to ellipse distance: a binary search approach,
: 有matlab的code
: 我老板不是搞这个方向的，这纯属于个人利用假期搞的，忙了一春假。
: 请牛人看看够不够发篇文章，latex 一共２１页。
: 如果能发在权威刊物最好，不行发哪个的大学的学报上也行。我现在很需要发文。
: 节选若干。
: 愿意帮忙请留言或站内方式联系。
|
h**********c
发帖数: 4120 | 9 请给出空间复杂度
【在 c*******h 的大作中提到】
: The distance from a point (x,y) to an ellipse (x/a)^2 + (y/b)^2 = 1 is
: ya^2*sin(t) + xb^2*cos(t)
: - -----------------------------
: a^2*sin^2(t) + b^2*cos^2(t)
: where t satisfies
: x^2+y^2-a^2-b^2
: sin(2t+s) = ---------------------------------------
: sqrt[ 4x^2y^2 + (b^2-y^2-a^2+x^2)^2 ]
: and s satisfies
: 2xy
|
g****t
发帖数: 31659 | 10 现实应用的牛顿法大部分都不显示求导。
另外:
给定一个点集，用最小二乘法辨识出来椭圆的a,b,
然后求解析解。这样也肯定比你的方法快。
【在 h**********c 的大作中提到】
: 牛顿法也许很好，但是但是用来解方程的。
: 它不能用来生成一个椭圆或者类似的convex shape，也就是不是一个表达椭圆的办法。
: 在计算几何里,convex hull是不连续的，没法求导的。而convex hull是可以用2维矢量
: 空间的一个generating set 来近似。
: 当然本贴只是想引起相关方向上的人士的兴趣。如果愿意的话我可以把２０页的pdf发
: 过去。
: my point may be niave. 但是，您必须证明您是serious的.
: 严格讲，计算几何都是numerical guy这点我们是应该承认的。加减乘除越简单越好。
: 另外，我没有找到求椭圆到点距离的现成文章或程序例题。再找找，如果哪位有，请提
: 供出处。
|
h**********c
发帖数: 4120 | 11 我写这个程序的本意就是能让计算机本科二三年级的学生就理解，并且不需要任何特别
package或算法就可以用c或java把matlab的script重写出来。
这个版在数学以及算法上的高手是很多的，办法肯定有更好的。 由于时间关系， 我没
法试验了。牛顿法以前写过，挺长时间了。
如果说数学上更好的办法，这个点是一个fixed point,能找到一个quadratic
convergence，那么这是最快的。
我的稿还是打算投一投试试，主要希望一些定理，概念能应用到离散的convex hull 中。
大家对我的观点还是比较宽容的，如果您愿意看，过些日子可以发给您。
我打算先试再说，如果您也是相关刊物的editor，我现在就可以给您发。
手上一大堆数据要处理，这个问题需要放一放。
我相信数学方法应该是去构造一个|f''(x)|=0的函数算这个fixed point，quadratic
convergence 估计是不能再快了， 解一个系统用 1e-7左右的精度，通常就三到四个
iteration.
【在 g****t 的大作中提到】
: 现实应用的牛顿法大部分都不显示求导。
: 另外:
: 给定一个点集，用最小二乘法辨识出来椭圆的a,b,
: 然后求解析解。这样也肯定比你的方法快。
|
g****t
发帖数: 31659 | 12 我是工程师。很多年没写过文章了。
我写这个程序的本意就是能让计算机本科二三年级的学生就理解，并且不需要任何特别
package或算法就可以用c或java把matlab的script重写出来。
这个版在数学以及算法上的高手是很多的，办法肯定有更好的。 由于时间关系， 我没
法试验了。牛顿法以前写过，挺长时间了。
如果说数学上更好的办法，这个点是一个fixed point,能找到一个quadratic
convergence，那么这是最快的。
我的稿还是打算投一投试试，主要希望一些定理，概念能应用到离散的convex hull 中。
大家对我的观点还是比较宽容的，如果您愿意看，过些日子可以发给您。
我打算先试再说，如果您也是相关刊物的editor，我现在就可以给您发。
手上一大堆数据要处理，这个问题需要放一放。
我相信数学方法应该是去构造一个|f''(x)|=0的函数算这个fixed point，quadratic
convergence 估计是不能再快了， 解一个系统用 1e-7左右的精度，通常就三到四个
iteration.
【在 h**********c 的大作中提到】
: 我写这个程序的本意就是能让计算机本科二三年级的学生就理解，并且不需要任何特别
: package或算法就可以用c或java把matlab的script重写出来。
: 这个版在数学以及算法上的高手是很多的，办法肯定有更好的。 由于时间关系， 我没
: 法试验了。牛顿法以前写过，挺长时间了。
: 如果说数学上更好的办法，这个点是一个fixed point,能找到一个quadratic
: convergence，那么这是最快的。
: 我的稿还是打算投一投试试，主要希望一些定理，概念能应用到离散的convex hull 中。
: 大家对我的观点还是比较宽容的，如果您愿意看，过些日子可以发给您。
: 我打算先试再说，如果您也是相关刊物的editor，我现在就可以给您发。
: 手上一大堆数据要处理，这个问题需要放一放。
|
h**********c
发帖数: 4120 | 13 这东西就得碰，有投十个发一，两个就行，今年有干有稀地先编个十篇八篇的。
万事开头难。
中。
【在 g****t 的大作中提到】
: 我是工程师。很多年没写过文章了。
:
: 我写这个程序的本意就是能让计算机本科二三年级的学生就理解，并且不需要任何特别
: package或算法就可以用c或java把matlab的script重写出来。
: 这个版在数学以及算法上的高手是很多的，办法肯定有更好的。 由于时间关系， 我没
: 法试验了。牛顿法以前写过，挺长时间了。
: 如果说数学上更好的办法，这个点是一个fixed point,能找到一个quadratic
: convergence，那么这是最快的。
: 我的稿还是打算投一投试试，主要希望一些定理，概念能应用到离散的convex hull 中。
: 大家对我的观点还是比较宽容的，如果您愿意看，过些日子可以发给您。
|
