MATLAB | MATLAB版玫瑰祝伟大女性节日快乐！！

news2024/11/18 23:41:16

妇女节到了，这里祝全体伟大的女性，节日快乐，事业有成，万事胜意。

作为MATLAB爱好者，这里还是老传统画朵花叭，不过感觉大部分样式的花都画过了，这里将一段很古老的2012年的html玫瑰花代码转成MATLAB版本

html版本代码出自：

http://www.romancortes.com/blog/1k-rose/

网站有代码详细原理，可以去瞅一眼：

完整代码

MATLAB代码也不是很长：

function htmlRose
f = 500;
h = -250;
ax                 = gca;
ax.NextPlot        = 'add';
ax.DataAspectRatio = [1,1,1];
ax.YDir            = 'reverse';
ax.XColor          = 'none';
ax.YColor          = 'none';
function V = calc(a, b, c)
if c > 60 
    V.x = sin(a.*7).*(13+5./(.2+power(b.*4, 4)))-sin(b).*50;
    V.y = b.*f+50;
    V.z = 625+cos(a.*7).*(13+5./(.2+power(b.*4, 4)))+b.*400;
    V.r = a-b./2;
    V.g = a;
else
    A = a.*2-1;
    B = b.*2-1;
    ind = A.*A+B.*B<1;
    a = a(ind); b = b(ind);
    A = A(ind); B = B(ind);
    if c > 37
        j = double(bitand(int32(c),1));
        n = double(j.*6+~j.*4);
        o = .5./(a+.01)+cos(b.*125).*3-a.*300;
        w = b.*h;
        V.x = o.*cos(n)+w.*sin(n)+j.*610-390;
        V.y = o.*sin(n)-w.*cos(n)+550-j.*350;
        V.z = 1180+cos(B+A).*99-j.*300;
        V.r = .4-a.*.1+power(1-B.*B, -h.*6).*.15-a.*b.*.4+cos(a+b)./5+...
              power(cos((o.*(a+1)+(B>0).*w-(B<=0).*w)./25), 30).*.1.*(1-B.*B);
        V.g = o./1e3+.7-o.*w.*3e-6;
    elseif c > 32
        c = c.*1.16-.15;
        o = a.*45-20;
        w = b.*b.*h;
        V.z = o.*sin(c)+w.*cos(c)+620;
        V.x = o.*cos(c)-w.*sin(c);
        V.y = 28+cos(B.*.5).*99-b.*b.*b.*60-V.z./2-h;
        V.r = (b.*b.*.3+power((1-(A.*A)), 7).*.15+.3).*b;
        V.g = b.*.7;
    else
        o = A.*(2-b).*(80-c.*2);
        w = 99-cos(A).*120-cos(b).*(-h-c.*4.9)+cos(power(1-b, 7)).*50+c.*2;
        V.z = o.*sin(c)+w.*cos(c)+700;
        V.x = o.*cos(c)-w.*sin(c);
        V.y = B.*99-cos(power(b, 7)).*50-c./3-V.z./1.35+450;
        V.r = (1-b./1.2).*.9+a.*.1;
        V.g = power((1-b), 20)./4+.05;
    end
end
end
zBuffer = zeros(f,f);
for k = 1:25
    for i = 0:45
        V = calc(rand(1, 5e3), rand(1, 5e3), i./.74);
        z = round(V.z + 0.5);
        x = int32(V.x.*f./z-h + 0.5);
        y = int32(V.y.*f./z-h + 0.5);
        ind2 = y<f;
        x = x(ind2); y = y(ind2); z = z(ind2);
        V.r = V.r(ind2); V.g = V.g(ind2);
        zBufferIndex = y.*f+x;
        zBufferBool = (~zBuffer(zBufferIndex))|zBuffer(zBufferIndex)>z;
        zBuffer(zBufferIndex(zBufferBool)) = z(zBufferBool);
        RGB = [bitcmp(int32(V.r.*h));
               bitcmp(int32(V.g.*h));
               bitcmp(int32(-V.r.*V.r.*80))].';
        RGB(RGB<0) = 0;
        RGB(RGB>255) = 255;
        RGB = double(RGB)./255;
        scatter(x(zBufferBool),y(zBufferBool),5,'filled',...
            'CData',RGB(zBufferBool,:),'MarkerEdgeColor','none',...
            'MarkerFaceAlpha',.7)
    end
    drawnow
end
end