{"id":8344,"date":"2022-02-14T12:59:53","date_gmt":"2022-02-14T17:59:53","guid":{"rendered":"https:\/\/blogs.mathworks.com\/community\/?p=8344"},"modified":"2022-02-14T13:01:34","modified_gmt":"2022-02-14T18:01:34","slug":"happy-valentines-day-2","status":"publish","type":"post","link":"https:\/\/blogs.mathworks.com\/community\/2022\/02\/14\/happy-valentines-day-2\/","title":{"rendered":"Happy Valentines Day!"},"content":{"rendered":"
\n
I found this lovely algorithmic heart in the comments of an old Walking Randomly post<\/a>. Research software superstar Mike Croucher<\/a> gets credit for the blog, but the comment and the associated code are courtesy of Mariano Beguerisse<\/a> (with some minor edits from me).<\/div>\n
Thanks Mariano!<\/div>\n
\n
\n
i = sqrt(-1);<\/span><\/div>\n<\/div>\n
\n
<\/div>\n<\/div>\n
\n
% Fuzzy heart-shaped Mandelbrot fractal<\/span><\/span><\/div>\n<\/div>\n
\n
% Mariano Beguerisse<\/span><\/span><\/div>\n<\/div>\n
\n
% October 2010<\/span><\/span><\/div>\n<\/div>\n
\n
<\/div>\n<\/div>\n
\n
% iterations<\/span><\/span><\/div>\n<\/div>\n
\n
n = 300;<\/span><\/div>\n<\/div>\n
\n
% resolution<\/span><\/span><\/div>\n<\/div>\n
\n
N = 200;<\/span><\/div>\n<\/div>\n
\n
% Create grid<\/span><\/span><\/div>\n<\/div>\n
\n
x = linspace(-1, 1, N);<\/span><\/div>\n<\/div>\n
\n
y = linspace(-1.4, .6, N);<\/span><\/div>\n<\/div>\n
\n
[X,Y] = meshgrid(x, y);<\/span><\/div>\n<\/div>\n
\n
<\/div>\n<\/div>\n
\n
Z = X + i*Y;<\/span><\/div>\n<\/div>\n
\n
Zn = Z;<\/span><\/div>\n<\/div>\n
\n
<\/div>\n<\/div>\n
\n
for <\/span>j=1:n<\/span><\/div>\n<\/div>\n
\n
% Mandelbrot map with random noise<\/span><\/span><\/div>\n<\/div>\n
\n
Zn = -i*(Zn).^2 + (rand(N,N).^(1\/5)).*Z;<\/span><\/div>\n<\/div>\n
\n
M = abs(Zn);<\/span><\/div>\n<\/div>\n
\n
ind1 = find(M<2);<\/span><\/div>\n<\/div>\n
\n
ind2 = find(M>=2);<\/span><\/div>\n<\/div>\n
\n
M(ind1) = 0;<\/span><\/div>\n<\/div>\n
\n
M(ind2) = -1;<\/span><\/div>\n<\/div>\n
\n
end<\/span><\/span><\/div>\n<\/div>\n
\n
<\/div>\n<\/div>\n
\n
imagesc(x, y, M)<\/span><\/div>\n<\/div>\n
\n
colormap([1 1 1; 1 0 0])<\/span><\/div>\n<\/div>\n
\n
set(gca, YDir=\"normal\"<\/span>)<\/span><\/div>\n<\/div>\n
\n
title(\"Love or whatever\"<\/span>, FontSize=16)<\/span><\/div>\n
\n
\n
<\/div>\n
<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n","protected":false},"excerpt":{"rendered":"
<\/div>\n

\nI found this lovely algorithmic heart in the comments of an old Walking Randomly post. Research software superstar Mike Croucher gets credit for the blog, but the comment and the associated code are... read more >><\/a><\/p>\n","protected":false},"author":69,"featured_media":0,"comment_status":"open","ping_status":"closed","sticky":false,"template":"","format":"standard","meta":[],"categories":[1],"tags":[],"_links":{"self":[{"href":"https:\/\/blogs.mathworks.com\/community\/wp-json\/wp\/v2\/posts\/8344"}],"collection":[{"href":"https:\/\/blogs.mathworks.com\/community\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/blogs.mathworks.com\/community\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/blogs.mathworks.com\/community\/wp-json\/wp\/v2\/users\/69"}],"replies":[{"embeddable":true,"href":"https:\/\/blogs.mathworks.com\/community\/wp-json\/wp\/v2\/comments?post=8344"}],"version-history":[{"count":3,"href":"https:\/\/blogs.mathworks.com\/community\/wp-json\/wp\/v2\/posts\/8344\/revisions"}],"predecessor-version":[{"id":8353,"href":"https:\/\/blogs.mathworks.com\/community\/wp-json\/wp\/v2\/posts\/8344\/revisions\/8353"}],"wp:attachment":[{"href":"https:\/\/blogs.mathworks.com\/community\/wp-json\/wp\/v2\/media?parent=8344"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/blogs.mathworks.com\/community\/wp-json\/wp\/v2\/categories?post=8344"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/blogs.mathworks.com\/community\/wp-json\/wp\/v2\/tags?post=8344"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}