转自mathwork,Sudharsana Iyengar
关注微信公众号:数学建模BOOM,回复“情人节”,获取绘图代码文件;同时有国赛与美赛获奖论文可供下载。
代码:
n=800;
p=pi;
[R,T]=ndgrid(linspace(0,1,n),linspace(-2,20*p,n));
x=1-(.5)*((5/4)*(1-mod(3.6*T,2*p)/p).^2-.25).^2;
U=2*exp(-T/(8*p));
L=sin(U);
J=cos(U);
y=1.99*(R.^2).*(1.2*R-1).^2.*L;
K=x.*(R.*L+y.*J);
X=K.*sin(T);
Y=K.*cos(T);
Z=x.*(R.*J-y.*L);
surf(X,Y,Z,'LineStyle','none')
grid,axis off;
colormap(jet)
关注微信公众号:数学建模BOOM,
回复 “群”,获得交流群群号,b站课件、数模资料在群文件
回复“课程”,查看入门级数学建模精品课