登陆注册
19349200000051

第51章 比例切割 算法初成

陈东风和李爷爷从董明竹家回来后已经快10点了。

进屋后,李爷爷对他说:“这段时间你婶婶不容易,她的事情你要多帮帮她。”

“当然,我是他侄子,不帮他帮谁?爷爷你放心。”陈东风肯定的说。

“那就好,你也知道,这次你文山叔叔去RB治病花了不少钱,光靠两个人在单位里的几十块工资是不够的。她想经商的想法很好,我很支持她。”

“现在的万元户,哪个不是做小生意的?婶婶只要敢闯,现在的华国是遍地黄金。”陈东风对董明竹和李文山以后从商也是很看好的。

“那就好,我老了,不过还是可以为她跑跑关系的。好了,你去洗漱吧,早点休息。”

陈东风现在哪里还能轻易睡得这么早,好不容易抓住灵感的尾巴,说什么也不能把他丢了。

其实陈东风在航模的设计阶段,不怕各种数学模型和理论计算。他最讨厌的是各种叶轮的工程设计图纸的绘制。

工程中表达叶轮表面的方法是投影图法,通常情况下使用的圆柱坐标系中的投影方法是旋转投影以及轴面投影、平面投影。叶轮的轴面投影图反映了叶轮的总体尺寸和特征,而平面投影相当于三视图里面的俯视图。需要叶片表面方程:Θ=(r,z)来一步步画出空间曲线。就算你画的再好,到了加工现场加工的时候也是晦涩难懂,需要专门的技术人员分析出一步步的运动指令,效率十分低下。

陈东风看到了贝塞尔曲线是光滑离散化成一段段的微小的直线段的实质,如果可以根据贝塞尔曲线的特征方程,设计出一套切实可行的算法——把各种曲率的曲线通过命令来生成,那么既可以简化设计中制图的复杂程度,也可以根据曲线各个离散化的坐标点反推出数控加工的路径。所以如果陈东风可以设计的出来那将是一举两得的。

70年代现在国际上的流行算法是(以一次方贝塞尔曲线为例)需要在两个定点A、B之间,在选定特定u的情况下在曲线上找到点C(u)。一个简单的方法是把u插到每一个基函数上,计算每个其与基函数的乘积以及其相应的控制顶点,最后将其相加。

虽然这种方法很好,但是缺乏数值稳定性,尤其是在计算伯恩斯坦多项式的时候可能引进数值误差。当然算法都是各个CAD软件的核心,是不会轻易示人的。

还有就是一款好的CAD软件不管是算法重要,它的控制核心也是非常重要,当然现在陈东风也没考虑到这些,只是想先把算法设计出来。

陈东风也是发了狠,不设计出这个算法,连门都懒的出了,过年前的卫生也忘记打扫了。好在李爷爷看他钻研辛亏,毫不计较,尽心尽力的为他做好了后勤工作。

终于,陈东风在大年三十前把这套他自己命名的比例切割算法给设计出来了。灵感来自于华国历史上著名的数学家祖冲之的割圆术,割圆术目的是等分圆之后取得正多边形,而他的比例切割是取一个特定的比例来逼近一个特殊的曲线。

简单来讲比例切割算法(以一次贝塞尔曲线为例)的基本观点是选择在AB中的一个点C,C将AB分为u:1-u(A到C的距离与AB之间的距离之比是u),让我们找到决定C在哪里的方法。

从A到B的向量是B-A。因为u是在0和1之间的比率,点C位于u(B-A)。将A的位置加以考虑,点C为A+u(B-A)=(1-u)A+uB。因此,对于给定的u,(1-u)A+uB是在A和B之间的点C,将AB分为u:1-u的两段。

更加具有普遍性的比例切割算法的想法如下是假设我们想要找到C(u),u在[0,1]中。由第一个多段线P0-P1-P2-P3...-Pn开始,利用上面的法则找到在线段上的点P1i,P1i在P0i到P0(i+1)的连线上并且将这段线分为u:1-u的两部分。依次地,我们可以得到n个点10,11,12,...,1(n-1),他们定义了一个新的多段线,一共有n-1段。

新点由1i进行标记,再次利用上面的规则我们可以得到第二个多段线,具有n-1个点(20,21,...,2(n-2))和n-2条边。从这个多段线开始,进行第三次,得到新的多段线,由n-2个点30,31,...,3(n-3)和n-3条边组成。重复这个过程n次得到一个点n0。

以上想法只是给定了比例切割想法的几何解释,而实际计算需要一个具体的计算方法。

首先,对于每一对临近的控制点,可以画出一条右上方和右下方的箭头(类似于杨辉三角),并且在两个箭头的交点处写下一个新点。例如相邻的两个点分别为ij 和i(j+1),新点是(i+1)j,右下方(相对应的左下方)的箭头表示将其尾数ij(相对应的为i(j+1))乘以1-u(相对应的乘以u),新的点是两个的和。

因此,从初始的第0列开始,我们计算第1列。之后从第1列得到第2列。最终,在n次计算之后我们最终到达了一个单个的点n0并且这个点就是在曲线上的点。下面的算法总结了上面我们讨论的内容,输入的是具有n+1个点的数列P和在0到1之间的u,最终得到在贝塞尔曲线上的点C(u)。

这个计算过程可以用递归的方法表示,对于j=0,1,...,n用P0,j表示Pj,也就是P0,j是第0列的第j项元素,在第i列计算第j项如下:P(i,j)=(1-u)P(i-1,j)+uP(i-1,j+1),(i=1,2....,n;j=0,1,2...,n-i)

元素Pi,j是(1-u)Pi-1,j(左上方元素)和 uPi-1,j+1(左下方元素)的和,最终的结果(在曲线上的点)是Pn,0.在这种想法的基础上,通过编程就可以得到基本的算法程序。

在这个基本算法的基础上,陈东风还需要对螺旋线、球面螺旋线、双弧外摆线和星行线、心脏线、圆内螺旋线、正弦曲线、太阳线和费马曲线等等几百种曲线给出需要选定的控制点数量和控制比例u。这个工作如果没有计算机的帮助的话,估计他这辈子都得耗在这上面了。

“好在,通用的算法已经计算出来了,可能有迭代算法效率的不高的问题,但是计算机应该可以克服。”陈东风一边自言自语,一边站了起来,往窗外一看天快黑了,正好肚子有点饿,把桌子上堆成小山的草稿纸整理下后,就出了房间。

同类推荐
  • 千面首富

    千面首富

    于哲从小便和做古董商的父亲一起生活。直到有一天发现了家里的秘密,父亲乃是祖传做人皮面具的,从此于哲便开始了自己波澜壮扩的一生。做家族富少,继承那五十亿遗产。做偶像明星,跟女明星谈恋爱,闹绯闻。做一线导演,亲手颠覆娱乐圈规则。戴上面具,就成为了别人,再也没有自我!
  • 银鳞河畔

    银鳞河畔

    超能力者陈拾,他和普通人唯一不同的地方在,他除了正常生活外还要维护其他人的正常生活。
  • 都市极品水神

    都市极品水神

    一个平凡的高中生,机缘巧合之下踏上修仙之路,从此机缘不断,踏上了人生的巅峰!PS:新建交流群号876968864
  • 辉煌从菜园子开始

    辉煌从菜园子开始

    告别996的生活,吴飞回到了乡下老家。偶得神奇空间的他,养养鱼,种种地,偶尔开开同学会,……简介无力,起名也不会,所以大家就随便看看好了
  • 一路向医

    一路向医

    科学的尽头是哲学;哲学的尽头是神学;那,医学的尽头是什么?王铮,一步步的向着医学的尽头,艰难前行。686864088,群号。
热门推荐
  • 愿化一面镜子,常常照你笑

    愿化一面镜子,常常照你笑

    本书中精选了朱生豪先生的120封情书。这些情书无不展示了他在爱情中的浪漫情怀和真心实意,一字一句,都充满了生活的情趣和小哲理,更让我们重新认识了一个立体的、栩栩如生的朱生豪。细细咀嚼,满口生香。愿化一面镜子,常常照你笑。
  • 学妹不太乖

    学妹不太乖

    (1v1苏爽,女主超强,作者自我感觉不虐) 她和平常人不太一样……好吧,很不一样。她在地下长大,只偶尔出次任务。而最后一次任务,就葬送了她的一生。她最初只是游戏人间,直到看到他的那串琉璃手链。她想要“幻色”。她和他朝夕相处,学校、家里、老宅……但他们彼此眼中的景象,究竟是真实,还是只是幻境的底色?
  • 都市之超级医仙

    都市之超级医仙

    【十万追读,都市精品】一个被家族抛弃的弃子叶凡。一纸不公平的婚姻合约却将他和冷艳美女总裁绑在了一起。他,左手惊天医术救世人,右手绝世武功踩强敌!他,纵横都市,一路嚣张,所向无前!
  • 傲娇苟少:笨蛋,宠你没商量

    傲娇苟少:笨蛋,宠你没商量

    你好,我是苟晞,江城荀氏集团董事长兼总经理。事情是这样的,我正在会所庆生,突然灵魂脱窍附身到了一条哈士奇身上。现在身无分文,手机是哈士奇饲主的,我乱输入一个号码,就找到了你,实在是缘分。能借给我1000块钱做路费么?等我回到荀氏,我还你10000!还让你做执行总监!这是我的X付宝账号goldenxi@#.com,谢谢!……因为随口一句玩笑的赌咒,江城大少苟晞,灵魂脱窍,附身在了普通少女石榴家的哈士奇身上。苟晞:???(我是谁?我在哪?发生了什么?)石榴:你是谁!苟晞:……我是你老公。石榴:老公?苟晞:哎,脑婆。石榴:为什么我老公会变成我的狗?
  • 史上最强归来

    史上最强归来

    《道德经》:“天地不仁,以万物为刍狗。”通俗点:就是天地看待万物是一样的,不对谁特别好,也不对谁特别坏,一切随其自然发展!换句话说,不管万物变成什么样子,那是万物自己的行为(包括运气),与天地无关,天还是干天的事,地还是干地的事,一切犹如随风入夜,润物无声,天地最是自然不过的。书友群:442038609
  • 追妻无门:女boss不好惹

    追妻无门:女boss不好惹

    青涩蜕变,如今她是能独当一面的女boss,爱了冷泽聿七年,也同样花了七年时间去忘记他。以为是陌路,他突然向他表白,扬言要娶她,她只当他是脑子抽风,他的殷勤她也全都无视。他帮她查她父母的死因,赶走身边情敌,解释当初拒绝她的告别,和故意对她冷漠都是无奈之举。突然爆出她父母的死居然和冷家有丝毫联系,还莫名跳出个公爵未婚夫,扬言要与她履行婚约。峰回路转,破镜还能重圆吗? PS:我又开新文了,每逢假期必书荒,新文《有你的世界遇到爱》,喜欢我的文的朋友可以来看看,这是重生类现言,对这个题材感兴趣的一定要收藏起来。
  • 太古尊帝

    太古尊帝

    【热血爽文,最火爆,挑战你的脑洞!】绝代丹皇转世而来,魔武双修,碎裂惊天阴谋!这一世他绝不受任何欺辱,御龙灵凤,五行聚顶,敢诛武帝灭仙魔!
  • 极品都市仙尊

    极品都市仙尊

    【网站合作作品】一个穷的连房租都交不起的落魄小民工,却更衰的出了车祸。但是天有福将,在这场车祸中他竟然意外得到了天地玄黄无极体,成为天纵奇才。都市、仙界、神界、虚无界,任我逍遥。主角在得到异能之后,入修真界狂妄天下。八千横里,我自笑刀自傲又如何?天地万物法则惟我傲视!且看强大少年如何称霸天下
  • 罪孽道途

    罪孽道途

    大夏皇宫,金碧辉煌的宫殿在骄阳下熠熠生辉,汉白玉铺筑的石阶直达天听,华美到过分的宫廷建筑一座又一座,让人数都数不过来,极尽强盛的夏朝到达了一个朝代所能达到了极致,这是至高点,然而,在有心人眼中,却也是大厦倾塌的预兆。
  • 快穿病娇出门左转

    快穿病娇出门左转

    何为病娇?也许是,悄无声息中将你的一切毁灭,温柔拥抱,极尽诱惑:“别怕,你还有我呢。”————眼前的病态少年,手指染血,轻言细语:“阿白…不要让我沾上你的血。”凌白:冷漠脸jpg绝美少年眯眼轻笑,笑容的弧度安静美好,隐藏的危险难以察觉。“阿白,别想逃。”“从未想。”(现代世界她的不二臣不要订阅,不要不要不要!已坑。)