Solar Symphony

KARUPPASAMYPANDIYAN M
/
Solar Symphony

on 7 Oct 2024

15
301
0
2
1553

Cite your audio source here (if applicable):

https://soundcloud.com/nasa/chorus-radio-waves-within-earths-atmosphere?in=nasa/sets/solar-system-beyond-sounds&si=c4fbe17dee9a4bbda687194ab3c04e91&utm_source=clipboard&utm_medium=text&utm_campaign=social_sharing

drawframe(1);

Write your drawframe function below

function drawframe(f)
    N = 600;
    [X0, Y0] = meshgrid(linspace(-12, 12, N)); 
    t = 2 * pi * f / 96;  % Time variable for orbit and rotation
    
    % Create Sun with glow
    sun_radius = 0.8; 
    sun_core = sqrt(X0.^2 + Y0.^2) <= sun_radius;  
    sun_glow_radius = 1.5;  
    sun_glow = (sqrt(X0.^2 + Y0.^2) > sun_radius) & (sqrt(X0.^2 + Y0.^2) <= sun_glow_radius);
    
    img = zeros(N, N, 3);
    sun_color = [1, 0.9, 0];  % Sun color (bright yellow)
    for i = 1:3
        img(:,:,i) = img(:,:,i) + sun_core * sun_color(i);
        img(:,:,i) = img(:,:,i) + sun_glow * sun_color(i) * 0.5;  % Sun glow
    end
    
    % Planet parameters (position, size, color)
    planets = {
        [0, 0], 0.5, [0, 0.4, 1];   % Earth
        [4, 0], 0.3, [1, 0.2, 0];   % Mars
        [5.5, 0], 0.6, [0.8, 0.5, 0.3]; % Jupiter
        [7, 0], 0.4, [0.6, 0.4, 0.9];   % Saturn
        [8.5, 0], 0.3, [0.5, 0.2, 0];   % Uranus
        [10, 0], 0.3, [0.6, 0.4, 0.3];  % Neptune
        [6.5, 1.5], 0.4, [0.9, 0.8, 0.2]; % Venus
        [4, -1.5], 0.2, [1, 1, 0]; % Mercury
        [9, 3], 0.4, [0.7, 0.7, 0.9]; % Pluto
    };
    
    % Speeds for orbital rotation and planet axis rotation
    orbit_speeds = [1, 0.5, 0.4, 0.3, 0.6, 0.4, 0.5, 0.7, 0.9]; 
    axis_speeds = [1, 0.8, 0.6, 0.7, 0.4, 0.3, 1, 1, 0.9];
    
    for j = 1:length(planets)
        % Orbit rotation
        orbit_radius = sqrt(planets{j, 1}(1)^2 + planets{j, 1}(2)^2);
        planet_x = orbit_radius * cos(t * orbit_speeds(j));
        planet_y = orbit_radius * sin(t * orbit_speeds(j));
        
        % Planet's self-rotation (spin)
        planet_theta = t * axis_speeds(j);
        X_rot = cos(planet_theta) * (X0 - planet_x) - sin(planet_theta) * (Y0 - planet_y);
        Y_rot = sin(planet_theta) * (X0 - planet_x) + cos(planet_theta) * (Y0 - planet_y);
        
        % Create planet's circular shape
        planet = sqrt(X_rot.^2 + Y_rot.^2) <= planets{j, 2};
        
        % Add planet's color
        for i = 1:3
            img(:,:,i) = img(:,:,i) + planet * planets{j, 3}(i);
        end
        
        % Saturn's ring
        if j == 4
            ring_mask = sqrt(X_rot.^2 + Y_rot.^2) <= (planets{j, 2} + 0.3) & sqrt(X_rot.^2 + Y_rot.^2) > planets{j, 2};
            img(:,:,1) = img(:,:,1) + ring_mask * 0.8;  % Ring color (whitish)
        end
    end
    
    % Starry background with brighter stars
    stars = rand(N, N) < 0.004;  % Increased density for more stars
    for i = 1:3
        img(:,:,i) = img(:,:,i) + stars * 1.0;  % Brighter stars
    end
    
    % Display the image
    imshow(img, 'InitialMagnification', 'fit');
    axis off;  % Hide axes for clean appearance
end

Solar Symphony

Movie

Audio

Remix Tree

Solar Symphony