如何用方程式写春联?

如何在春节坐在电脑前也很有过年气氛?尤其是面对着家里还是 Windows XP 的桌电。我翻到了一篇文章,介绍 GrafEq 这套付费但可无限试用且只支持到 Windows XP 和 OS 9.2 的古董级数学绘图软件。

我们这就拿它来做一点实验。

首先打入 x^2+y^2

宽客网,量化投资,宽客俱乐部

若想要这个圆方一点,可以打 x^6+y^6

宽客网,量化投资,宽客俱乐部

x^4+y^4

宽客网,量化投资,宽客俱乐部

现在我想在它身上挖个洞,于是想到要在不等式左边减掉一个在原点值很大、但是出原点就迅速变小的函数。

例如 1/(x^4+y^4) 就是这样的函数,它在 x, y 很小的时候有很小的分母、所以函数值很大,但是 x, y 变大之后四次方会让分母变大得很快,所以函数值就变得微不足道地小。

现在 x^4+y^4+1/(x^4+y^4)

宽客网,量化投资,宽客俱乐部

能不能让洞大一些呢?当然可以,把分子的 1 改成 100 就好

x^4+y^4+100/(x^4+y^4)

宽客网,量化投资,宽客俱乐部

能不能让洞只出现在上半部呢?当然也可以,将分母的 y 用 2y-2 取代

现在的式子是:x^4+y^4+100/(x^4+ (2y-2)^4)

宽客网,量化投资,宽客俱乐部

不过看起来有点勉强啊,毕竟 100/(x^4+ (2y-2)^4) 在左上角和右上角的地方降得不够快,和正在变大的 x^4+y^4 相加之后,函数值比 25 小的地方显得有点细(快不见得右上角和左上角)。我们可以从下面这张图看到在该函数在 x=y 这条线上的行为,因为我们要看沿着斜 45 度角这条线上的笔画粗细。

宽客网,量化投资,宽客俱乐部

我们现在考虑的是 f (x,y)= x^4+y^4+100/(x^4+ (2y-2)^4) 这个函数,要把它画出来需要三维空间,这里却只有两维,不得已只好取函数在 x=y 这条线上面的样子,于是上图的纵坐标代表 f(x,y)= x^4+y^4+100/(x^4+ (2y-2)^4) 函数的值,横坐标则是 x=y 这条线,往左是原图的左下角往右是原图的右上角,至于浅浅那条横线代表不等号右边的 25,于是在该线下方的函数值对应到被涂黑的部分,该线上方的函数值对应到留白的部分。

当然,缩放过了,你知道如何用 GrafEq 画出这张图吗?

我的答案是

(-5y+2(x/4)^4+100/((x/4)^4+ (2(x/4)-2)^4))(-5y+25)=0

为了让那个地方粗一些,我们决定在分母动手脚,让它在左上角和右上角不要增加太快,方法是让 y=1 附近的分母变大,例如把分母 x^4+ (2y-2)^4 再加上 (2y-1)^2,这会让 x=y 上的行为变成

(-5y+2(x/4)^4+100/((x/4)^4+ (2(x/4)-2)^4+ (2(x/4)-1)^2))(-5y+25)=0

宽客网,量化投资,宽客俱乐部

而 x^4+y^4+100/(x^4+ (2y-2)^4+ (2y-1)^2)

宽客网,量化投资,宽客俱乐部

可爱多了,忍不住想再戳第二个洞

x^4+y^4+100/(x^4+ (2y-2)^4+ (2y-1)^2) +100/(x^4+ (2y+2)^4+ (2y+1)^2)

一个「日字」就完成了

宽客网,量化投资,宽客俱乐部

可以再画其他笔画吗?没有问题,只要调整 x, y 前面的系数与常数,我们可以在任何位置画上任意长度

例如 -1/((y+3)^4+ (x/15)^4)-1/((y+4)^4+ (x/15)^4)-1/((y+5)^4+ (x/15)^4)

(知道这三条分别对应到哪一项吗?)

宽客网,量化投资,宽客俱乐部

将上式的左边放到还没挖洞的函数里

x^4+y^4-1/((y+3)^4+ (x/15)^4)-1/((y+4)^4+ (x/15)^4)-1/((y+5)^4+ (x/15)^4)

会发现这三条的粗细和长度都变了,离原点越远也就是越下面的笔画就会越细。这是因为离原点越远的话x^4+y^4 就越大,所以扣掉 XX 分之一之后就越容易超出 25(不等式右式给的边界),笔画就缩水了。

宽客网,量化投资,宽客俱乐部

话又说回来,这像不像旭日东升呢?

还是说象形字看起来比较有感觉?

x^4+y^4+100/(x^4+ (2y-2)^4+ (2y-1)^2) +100/(x^4+ (2y+2)^4+ (2y+1)^2)-1/((y+3)^4+ (x/15)^4)-1/((y+4)^4+ (x/15)^4)-1/((y+5)^4+ (x/15)^4)

宽客网,量化投资,宽客俱乐部

抑或是 x=0 的截面比较有启发性?

(-50x+y^4+100/((2y-2)^4+ (2y-1)^2) +100/((2y+2)^4+

(2y+1)^2)-1/((y+3)^4)-1/((y+4)^4)-1/((y+5)^4))(-50x+25)=0

宽客网,量化投资,宽客俱乐部

下一题:要怎么画出斜的笔画?毕竟调常数只会平移图形,调整系数顶多让矮胖的变瘦高的,那要怎么调整「角度」呢?

这就要用到所谓「线性变换」的概念了,例如将

-1/(4(x+4)^4+ ((y+1)/5)^4)

宽客网,量化投资,宽客俱乐部

套用变换「 x ? x+y 、 y ? x-y 」会得到

-1/(4(x+y+4)^4+ ((x-y+1)/5)^4)

宽客网,量化投资,宽客俱乐部

那么再下一个问题:要怎么画出弯的笔画?

这可是再多线性变换都做不到的事,穷途末路的我们回想起圆弧是弯的 (x-19)^2+ (y+12)^2=400 ,所以我们只要设定「到圆弧的距离够小」,或者说「到圆心的距离界在某两个数中间」就好了吧!((x-19)^2+ (y+12)^2-400)^2

不过等等,这笔一画下去可就是整个圆了,能不能只取其中一段呢?

宽客网,量化投资,宽客俱乐部

我想只靠一个圆大概不够,那你有没有试第二个?

((x-4)^2+ (y+5)^2-13)((x-19)^2+ (y+12)^2-400)=0

宽客网,量化投资,宽客俱乐部

也许设定「到两个圆的距离都够小」会有意外收获

((x-4)^2+ (y+5)^2-13)^2+ ((x-19)^2+ (y+12)^2-400)^2

宽客网,量化投资,宽客俱乐部

唉呀,为了把这笔画整合到本来的图形里,我们需要把「够靠近零」的部分放到分母,变成一个「够大」的条件

-100/(((x-4)^2+ (y+5)^2-13)^2+ ((x-19)^2+ (y+12)^2-400)^2)

(需要看看它的样子吗?我猜是不用。)

不过它其实还是有点脆弱,根据努力尝试的结果,我建议狠下心来将分子分母同时平方四次,这样该大的才会大 ,该小的才会小。

-100^16/(((x-4)^2+ (y+5)^2-13)^2+ ((x-19)^2+ (y+12)^2-400)^2)^16

宽客网,量化投资,宽客俱乐部

(这个图形和前一个一样吗?)(如果不,那它为什么和前一个那么像?)

最后将所有这些东西通通加在一起

x^4+y^4+100/(x^4+ (2y-2)^4+ (2y-1)^2) +100/(x^4+ (2y+2)^4+ (2y+1)^2)-1/((y+3)^4+ (x/15)^4)-1/((y+4)^4+ (x/15)^4)-1/((y+5)^4+ (x/15)^4)-1/(4(x+y+4)^4+ ((x-y+1)/5)^4)-100^16/(((x-4)^2+ (y+5)^2-13)^2+ ((x-19)^2+ (y+12)^2-400)^2)^16

就得到了:

宽客网,量化投资,宽客俱乐部

至于红色的春联纸怎么画呢?我希望字出现在正中间,所以应该要有些平移,我希望它够尖,所以指数部分应该够大,我希望它是 45 度的,所以需要线性变换,就决定是你了 (x+y+2)^100+ (x-y-2)^100

宽客网,量化投资,宽客俱乐部

这就是结果啦

宽客网,量化投资,宽客俱乐部

最后,就来道规规矩矩的练习题让大家试试吧:

宽客网,量化投资,宽客俱乐部

祝福各位新年快乐!

编按:如果以上这串你有看没有懂没关系,为了犒赏你能够阅读这篇文章到底,把以下式子贴入 google 搜索,会有意想不到的惊喜哦!

x^4+y^4+100/(x^4+ (2y-2)^4+ (2y-1)^2) +100/(x^4+ (2y+2)^4+ (2y+1)^2)-1/((y+3)^4+ (x/15)^4)-1/((y+4)^4+ (x/15)^4)-1/((y+5)^4+ (x/15)^4)-1/(4(x+y+4)^4+ ((x-y+1)/5)^4)-100^16/(((x-4)^2+ (y+5)^2-13)^2+ ((x-19)^2+ (y+12)^2-400)^2)^16-25

关于GrafEq

古董级函数作图软件—— GrafEq 。这个软件只有 2M 大小,它的功能就只有一个:作出形如 x2 + y2 = 1 的二元等式或者不等式的图像。令人惊叹的是,这个软件的图像绘制能力异常强大, Mathematica 等大型专业数学软件完全不是它的对手。
金融工程, 数学算法

风险提示及免责条款

市场有风险,投资需谨慎。本文不构成个人投资建议,也未考虑到个别用户特殊的投资目标、财务状况或需要。用户应考虑本文中的任何意见、观点或结论是否符合其特定状况。据此投资,责任自负。

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场,不承担相关法律责任。如若转载,请注明出处。 如若内容造成侵权/违法违规/事实不符,请点击下方“内容举报”进行投诉反馈!
立即
投稿

微信公众账号

微信扫一扫加关注

返回
顶部