matlab模拟黑洞包含吸积盘和喷流，简单模拟

本文介绍

黑洞的简单实现和模拟

代码

% Black Hole Simulation in 3D

% Clear workspace and figures
clear; close all; clc;

% Create figure and set axis properties
figure;
axis([-10 10 -10 10 -10 10]);
hold on;
grid on;
view(3);

% Parameters for the black hole and accretion disk
numParticles = 1000;  % Number of particles in the accretion disk
diskRadius = 5;       % Radius of the accretion disk
diskHeight = 0.5;     % Height of the accretion disk
jetHeight = 10;       % Height of the jet
jetRadius = 1;        % Radius of the jet

% Create particles for the accretion disk
theta = 2*pi*rand(numParticles,1);
r = diskRadius * sqrt(rand(numParticles,1));
x = r .* cos(theta);
y = r .* sin(theta);
z = diskHeight * (rand(numParticles,1) - 0.5);

% Create particles for the jet
jetParticles = 500;
jetX = jetRadius * (rand(jetParticles,1) - 0.5);
jetY = jetRadius * (rand(jetParticles,1) - 0.5);
jetZ = linspace(0, jetHeight, jetParticles)';

% Plot the black hole
[xSphere, ySphere, zSphere] = sphere(20);
surf(0.5*xSphere, 0.5*ySphere, 0.5*zSphere, 'FaceColor', 'k', 'EdgeColor', 'none');

% Plot the accretion disk
hDisk = scatter3(x, y, z, 1, 'r');

% Plot the jet
hJet = scatter3(jetX, jetY, jetZ, 1, 'b');

% Animation loop
numFrames = 200;
for t = 1:numFrames
    % Update accretion disk particles (simple circular motion)
    theta = theta + 0.1;
    x = r .* cos(theta);
    y = r .* sin(theta);
    set(hDisk, 'XData', x, 'YData', y, 'ZData', z);
    
    % Update jet particles (move upwards)
    jetZ = jetZ + 0.1;
    set(hJet, 'XData', jetX, 'YData', jetY, 'ZData', jetZ);
    
    % Loop jet particles back to start
    jetZ(jetZ > jetHeight) = 0;
    
    % Pause to create animation effect
    pause(0.05);
end