{"id":2483,"date":"2024-05-15T12:43:14","date_gmt":"2024-05-15T16:43:14","guid":{"rendered":"https:\/\/blogs.mathworks.com\/matlab\/?p=2483"},"modified":"2024-05-20T15:14:41","modified_gmt":"2024-05-20T19:14:41","slug":"celebrating-spring-matlab-tulip","status":"publish","type":"post","link":"https:\/\/blogs.mathworks.com\/matlab\/2024\/05\/15\/celebrating-spring-matlab-tulip\/","title":{"rendered":"Celebrating Spring: MATLAB Tulip"},"content":{"rendered":"
\r\n
\"\"<\/a><\/div>\r\n
Joining us again is Eric Ludlam, development manager of MATLAB\u2019s charting team. Discover more about Eric on our contributors bio page<\/span>.<\/a> Last time Eric was here, Daffodils were on his mind<\/a>. Now, he focuses on tulips.<\/div>\r\n
\"MWSignTulips.jpg\"<\/div>\r\n
Spring is here in Natick and the tulips are blooming! While tulips appear only briefly here in Massachusetts, they provide a lot of bright and diverse colors and shapes. To celebrate this cheerful flower, here's some code to create your own tulip!<\/div>\r\n
This script uses some fun tricks, like extracting a subset of the 'hot' colormap for a red-to-yellow colormap, and specifying several material properties to give it a more sunny day look. There's also the sinpi and cospi functions that Mike likes so much<\/a>!<\/div>\r\n
Explore the code and create some new breeds of tulips for your MATLAB.<\/div>\r\n
\r\n
\r\n
%% MATLAB Tulip<\/span><\/span><\/div>\r\n<\/div>\r\n
\r\n
np = 3; % number of petals<\/span><\/span><\/div>\r\n<\/div>\r\n
\r\n
n=np*20+1; % Theta resolution<\/span><\/span><\/div>\r\n<\/div>\r\n
\r\n
theta=linspace(0, np*(1\/np)*2,n);<\/span><\/div>\r\n<\/div>\r\n
\r\n
r=linspace(0,1,80)';<\/span><\/div>\r\n<\/div>\r\n
\r\n
<\/div>\r\n<\/div>\r\n
\r\n
newplot<\/span><\/div>\r\n<\/div>\r\n
\r\n
for <\/span>k=[0 1\/np] % 2 layers of petals<\/span><\/span><\/div>\r\n<\/div>\r\n
\r\n
\u00a0\u00a0\u00a0\u00a0x=1-(.8*(1-mod(np*(theta+k),2)).^3-.05).^2\/2;<\/span><\/div>\r\n<\/div>\r\n
\r\n
\u00a0\u00a0\u00a0\u00a0Z1=((x.*r).^6);<\/span><\/div>\r\n<\/div>\r\n
\r\n
\u00a0\u00a0\u00a0\u00a0R2=x.*r*(1-k*.38);<\/span><\/div>\r\n<\/div>\r\n
\r\n
\u00a0\u00a0\u00a0\u00a0R3=x.*(r.^8)*.4.*Z1;<\/span><\/div>\r\n<\/div>\r\n
\r\n
\u00a0\u00a0\u00a0\u00a0X=R2.*cospi(theta)-R3.*cospi(theta);<\/span><\/div>\r\n<\/div>\r\n
\r\n
\u00a0\u00a0\u00a0\u00a0Y=R2.*sinpi(theta)-R3.*sinpi(theta);<\/span><\/div>\r\n<\/div>\r\n
\r\n
\u00a0\u00a0\u00a0\u00a0C=repmat(r,1,n).^25;<\/span><\/div>\r\n<\/div>\r\n
\r\n
%% Petals<\/span><\/span><\/div>\r\n<\/div>\r\n
\r\n
\u00a0\u00a0\u00a0\u00a0surface(X,Y,Z1*2,C,FaceColor=\"interp\"<\/span>,EdgeColor=\"none\"<\/span>,FaceLighting=\"gouraud\"<\/span>);<\/span><\/div>\r\n<\/div>\r\n
\r\n
\u00a0\u00a0\u00a0\u00a0line(X(end,:),Y(end,:),Z1(end,:)*2,Color=\"#cc0\"<\/span>,LineWidth=1);<\/span><\/div>\r\n<\/div>\r\n
\r\n
end<\/span><\/span><\/div>\r\n<\/div>\r\n
\r\n
%% Stem<\/span><\/span><\/div>\r\n<\/div>\r\n
\r\n
[CX,CY,CZ]=cylinder;<\/span><\/div>\r\n<\/div>\r\n
\r\n
surface(CX*.1,CY*.1,(CZ-1)*2,[],FaceColor=\"#9b8\"<\/span>,EdgeColor=\"none\"<\/span>);<\/span><\/div>\r\n<\/div>\r\n
\r\n
%% Decorate!<\/span><\/span><\/div>\r\n<\/div>\r\n
\r\n
cmap=hot(200); % Colors of tulip<\/span><\/span><\/div>\r\n<\/div>\r\n
\r\n
colormap(cmap(80:150,:)) % Extract subset of hot colormap.<\/span><\/span><\/div>\r\n<\/div>\r\n
\r\n
axis off equal<\/span><\/span><\/div>\r\n<\/div>\r\n
\r\n
view([0,10]);<\/span><\/div>\r\n<\/div>\r\n
\r\n
material([.6 .9 .3 2 .5])<\/span><\/div>\r\n<\/div>\r\n
\r\n
camlight(40,40)<\/span><\/div>\r\n<\/div>\r\n
\r\n
\r\n
\r\n\r\n\r\n
<\/div>\r\n
<\/div>\r\n
<\/div>\r\n
<\/div>\r\n<\/div>\r\n<\/div>\r\n<\/div>\r\n<\/div>\r\n \r\n\r\n<\/div>\r\n