解析式大赛的获奖作品代码和公式公布啦！！！
2016/6/7 超级数学建模
上方超级数学建模可加关注
传播数学干货，学会理性的方式去思考问题
大家期待已久的运行代码终于
新鲜出炉了!!!
抱歉让各位粉丝久等了
接下来
就是揭秘奇迹的时刻
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、兔子
解析式：无
代码：
运行软件：不详
解析式揭秘告此一段落了
大家是否还意犹未尽？
小编告诉你们个秘密哦
在不久的将来
下一波比赛即将来临!!!
敬请期待...
http://www.duyihua.cn
返回 超级数学建模 返回首页 返回百拇医药