解析式大赛的获奖作品代码和公式公布啦!

上方超级数学建模可加关注

传播数学干货，学会理性的方式去思考问题

大家期待已久的运行代码终于

新鲜出炉了！！！

抱歉让各位粉丝久等了

接下来

就是揭秘奇迹的时刻

1、emoji

解析式：无

代码：

hold on

axis([-30 25 -25 23])

axis equal

grid off

h=ezplot('x.^2+y.^2=22.5^2+0.000001*((y+7.5)*(y+12.6)).^0.5',[-30,30],[-30,30])

set(h,'Color','k')

h=ezplot('y=0.04*x.^2-17',[-12.4,12.4],[-20,-10])%lower lip

set(h,'Color','k')

h=ezplot('y=0.013*x.^2-10',[-11.8,11.8],[-12,-8])%teeth bot

set(h,'Color','k')

h=ezplot('y=0.012*x.^2-6',[-11.25,11.25],[-6,-4])%upper lip bot

set(h,'Color','k')

h=ezplot('y=0.012*x.^2-5',[-11.1,11.1],[-5,-3])%upper lip top

set(h,'Color','k')

h=ezplot('y=0.65*x+8',[-20,-10.5],[-5,1.5])%left tear upper line

set(h,'Color','k')

h=ezplot('y=-0.65*x+8',[10.5,20],[-5,1.5])%right tear upper line

set(h,'Color','k')

h=ezplot('y=7.45*x+79.4',[-12,-10.5],[-10,1.2])%left tear lower line

set(h,'Color','k')

h=ezplot('y=-7.45*x+79.4',[10.5,12],[-10,1.2])%right tear lower line

set(h,'Color','k')

h=ezplot('(x+8.5).^2/36+(y+2.5).^2/20=1', [-14.2709,-2.81], [-2,2])%left eye upper curve

set(h,'Color','k')

h=ezplot('(x-8.5).^2/36+(y+2.5).^2/20=1', [2.81,14.2709], [-2,2])%right eye upper curve

set(h,'Color','k')

h=ezplot('-16.5+23.14^.5*cos(t)','-8.3+23.14^.5*sin(t)',[2.3855,5.921981])%left tear curve

set(h,'Color','k')

h=ezplot('16.5+23.14^.5*cos(t)','-8.3+23.14^.5*sin(t)',[3.502797,7.0393])%left tear curve

set(h,'Color','k')

h=ezplot('(x+8.5).^2/50+(y+6).^2/20=1', [-10.9,-2.7], [-2.4498,0])%left eye lower curve

set(h,'Color','k')

h=ezplot('(x-8.5).^2/50+(y+6).^2/20=1', [2.7,10.9], [-2.4498,0])%right eye lower curve

set(h,'Color','k')

h=ezplot('-3.78+1.05*cos(t)','-1.49+1.05*sin(t)',[4.3009,6.69])%left eye connection

set(h,'Color','k')

h=ezplot('3.78+1.05*cos(t)','-1.49+1.05*sin(t)',[2.7348,5.12388])%left eye connection

set(h,'Color','k')

h=ezplot('(x+8).^2/110+(y+2).^2/90=1', [-19,-5], [0.4,6.902])%left eyebrow upper curve

set(h,'Color','k')

h=ezplot('(x-8).^2/110+(y+2).^2/90=1', [5,19], [0.4,6.902])%right eyebrow upper curve

set(h,'Color','k')

h=ezplot('(x+6).^2/150+(y+6).^2/140=1', [-18,-8], [-1,5.2])%left eyebrow lower curve

set(h,'Color','k')

h=ezplot('(x-6).^2/150+(y+6).^2/140=1', [8,18], [-1,5.2])%right、

运行软件：MALTAB

2、缤纷五环

解析式：

代码：

ParametricPlot3D[{

{4 + (3 + Cos[v]) Sin[u], 4 + (3 + Cos[v]) Cos[u], 4 + Sin[v]},

{8 + (3 + Cos[v]) Cos[u], 3 + Sin[v], 4 + (3 + Cos[v]) Sin[u]},

{12 + (3 + Cos[v]) Sin[u], 4 + (3 + Cos[v]) Cos[u], 4 + Sin[v]},

{16 + (3 + Cos[v]) Cos[u], 3 + Sin[v], 4 + (3 + Cos[v]) Sin[u]},

{20 + (3 + Cos[v]) Sin[u], 4 + (3 + Cos[v]) Cos[u], 4 + Sin[v]}

},

{u, 0, 2 Pi},

{v, 0, 2 Pi}]

运行软件：mathematica 10.0

3、蓝宝石帘

解析式：

代码：

clear

x=-10:0.01:10;

y=-10:0.01:10;

[X,Y]=meshgrid(x,y);

Z=abs(sin(X.^2-Y.^2)-exp(cos(X.*Y)));

surf(X,Y,Z)

shading flat

view(0,90)

shading flat

xlabel('X'),ylabel('Y'),zlabel('Z')

title('神奇的蓝宝石帘’)

运行软件：MATLAB

4、Math&Life

解析式：

代码：

无

运行软件： GrafEq 2.13

5、水晶之恋

解析式：

代码：

x=-20*pi:pi/20:2*pi;%角度范围

y=(cos(2*x-pi/2))-2.*sin(4.*(2*x-pi/2))+(cos((2*x-pi/2)./12)).^5;

%半径取值

polar(x,y,'k')%极坐标画图

运行软件：MATLAB

6、兔子

运行软件：不详

解析式揭秘告此一段落了

大家是否还意犹未尽？

小编告诉你们个秘密哦

在不久的将来

下一波比赛即将来临！！！

敬请期待...