方程的艺术【Batman Equation？】，与大家共分享！

路人

人人上看到的，换了张超大清晰图【Powered by Google Pic search】，求逆向推测创作思路！

路人

呕心沥血花了两个小时的成果········



PS：已知此方程在NS环境下拖动显示区域时会出现如下问题（输入的方程没有问题）


CP Manager表示画整个方程式的话表示内存不够，但可以分式作图
CP报错，求高手


wolframalpha表示字符太长、TI-92+模拟器表示无法黏贴。
方程线性表示详见4L



tns下载：（第一也不是图像，请放心~）
注意！
这个方程在模拟器上的共花费了10s才显现出来，
经估计：模拟器上运行需要1s的文档，实体机至少要10~∞s
所以请做好ns"死机"的准备（可能也就等个2Mins吧~）

Batman.tns (4.92 KB, 下载次数: 82)

2011-7-31 23:53 上传

点击文件名下载附件

diameter

2# 路人

笔记本上没装Mathematica（我也不常用），
用wolframAlpha试了试，线性的写了一长串才完成方程的一点点……
再输要吐血了

路人

2# 路人  

笔记本上没装Mathematica（我也不常用），
用wolframAlpha试了试，线性的写了一长串才完成方程的一点点……
在输要吐血了
diameter 发表于 2011-7-31 21:52

1. f1(x)=zeros(((((x)/(7)))^(2)*√(((abs(abs(x)-3))/(abs(x)-3)))+(((y)/(3)))^(2)*√(((abs(y+((3*√(33))/(7))))/(y+((3*√(33))/(7)))))-1)*(abs(((x)/(2)))-((3*√(33)-7)/(112))*x^(2)-3+√(1-(abs(abs(x)-2)-1)^(2))-y)*(9*√(((abs((abs(x)-1)*(abs(x)-0.75)))/((1-abs(x))*(abs(x)-0.75))))-8*abs(x)-y)*(3*abs(x)+0.75*√(((abs((abs(x)-0.75)*(abs(x)-0.5)))/((0.75-abs(x))*(abs(x)-0.5))))-y)*(2.25*√(((abs((x-0.5)*(x+0.5)))/((0.5-x)*(0.5+x))))-y)*(((6*√(10))/(7))+(1.5-0.5*abs(x))*√(((abs(abs(x)-1))/(abs(x)-1)))-((6*√(10))/(14))*√(4-(abs(x)-1)^(2))-y),y)

复制代码

突然感觉用
TI-89ti & v200 & 92+
的同学是多么的强大！

zhaochenze

那么V5的方程是外星人写的么？

路人

那么V5的方程是外星人写的么？
zhaochenze 发表于 2011-8-1 06:59

我认为是几个F(X)=0的方程乘起来的结果，
至于方程是怎么算出来的，
我认为是NS之类的计算器根据图像枚举出的方程式

zhengran14

帅！！！！！！！！！！！！！！！！！！

diameter

wolframAlpha表示压力很大，方程输不完……

路人

wolframAlpha表示压力很大，方程输不完……
diameter 发表于 2011-8-1 09:40

恩，这个方程式实在是太长了（526个字符）！

~呼叫CP高手解决2L问题~

lch

果然很强大

zhj

1. ((x/7)^2sqrt(abs(abs(x)-3))/(abs(x)-3)+(y/3)^2*sqrt(abs(y+3sqrt(33)/7))/(y+3sqrt(33)/7)-1)(abs(x/2)-((3sqrt(33)-7)/112)x^2-3+sqrt(1-(abs(abs(x)-2)-1)^2)-y)(9sqrt(abs((abs(x)-1)(abs(x)-75))/((1-abs(x))(abs(x)-75)))-8abs(x)-y)(3abs(x)+0.75sqrt(abs((abs(x)-0.75*(abs(x)-0.5)))/(((0.75-abs(x))*(abs(x)-0.5))))-y)*(6sqrt(10)/7+(1.5-0.5abs(x))sqrt(abs(abs(x)-1)/(abs(x)-1))-6sqrt(10)/14*sqrt(4-(abs(x)-1)^2)-y)=0

复制代码

我用grafeq画不出来

imath

乘积就是多条的合成
根号估计是控制定义域的
其实大约的确是一堆二次型曲线

hcz

类似ES编程的方式，把多分支算式合并

hcz

= =Mathematica杯具中

hcz

1. p := Abs[x]; ContourPlot[{
   
2.   If[((6 Sqrt[10])/7 - y) (y + (3 Sqrt[33])/7) > 0,
   
3.     p^2/49 + y^2/9 - 1] == 0,
   
4.   p/2 - (3 Sqrt[33] - 7)/112 p^2 - 3 + Sqrt[1 - (Abs[p - 2] - 1)^2] ==
   
5.     y,
   
6.   If[(1 - p) (p - 3/4) > 0, 9 - 8 p - y] == 0,
   
7.   If[(3/4 - p) (p - 1/2) > 0, 3 p + 3/4 - y] == 0,
   
8.   If[(1/2 - p) (p + 1/2) > 0, 9/4 - y] == 0,
   
9.   If[p - 1 > 0, (6 Sqrt[10])/7 + (3 - p)/
   
10.      2 - (3 Sqrt[10])/7 Sqrt[4 - (p - 1)^2] - y] == 0
    
11.   }, {x, -8, 8}, {y, -4, 4}, PlotPoints -> 60]

复制代码

一种可以画出图的。。。不过实在是太慢了

Cyvre

乘积就是多条的合成
根号估计是控制定义域的
其实大约的确是一堆二次型曲线
imz 发表于 2011-8-3 12:21

我也是这样想的

imath

改天写个表达式生成器

路人

16# Cyvre

我也是这么想的，
没人去解决2LCP报错的问题么。

noivan

((x/7)^2sqrt(abs(abs(x)-3))/(abs(x)-3)+(y/3)^2*sqrt(abs(y+3sqrt(33)/7))/(y+3sqrt(33)/7)-1)(abs(x/2)-((3sqrt(33)-7)/112)x^2-3+sqrt(1-(abs(abs(x)-2)-1)^2)-y)(9sqrt(abs((abs(x)-1)(abs(x)-75))/((1-abs(x)) ...
zhj 发表于 2011-8-3 12:01

这个软件挺不错，就是只支持到xp啊

imath

nspire数学模式写好，复制到字符串里就是线性模式