Effacer les filtres
Effacer les filtres

my code wont display figure

1 vue (au cours des 30 derniers jours)
Bijaya
Bijaya le 16 Oct 2023
Commenté : Torsten le 16 Oct 2023
% Parameters
L = 1; % Length of the string
c = 1; % Wave speed
num_segments = 200; % Number of segments
dx = L / num_segments; % Discretization step
dt = 0.005; % Time step
num_steps = 10 / dt; % Number of time steps
% Initial wave parameters
left_wave_length = 0.15 * L;
right_wave_length = 0.25 * L;
% Initialize the wave function
x = 0:dx:L;
u = sin(pi * x / left_wave_length) + sin(2 * pi * x / right_wave_length);
% Time integration
for step = 1:num_steps
% Compute second spatial derivative using finite differences
u_xx = (u(1:end-2) - 2*u(2:end-1) + u(3:end)) / dx^2;
% Update the wave equation using finite differences
u_new = 2*u(2:end-1) - u(2:end-1) + c^2 * dt^2 * u_xx;
% Update boundary conditions
u_new1 = [0, u_new, 0];
% Update the wave function
u = u_new1;
end
% Plot at t = 0
figure;
plot(x, u);
title('Wave Equation at t = 0');
xlabel('Position [L]');
ylabel('Amplitude');
xlim([0, L]);
% Plot at t = 10 [t]
figure;
plot(x, u);
title('Wave Equation at t = 10 [t]');
xlabel('Position [L]');
ylabel('Amplitude');
xlim([0, L]);
  1 commentaire
Torsten
Torsten le 16 Oct 2023
Choose another code from the internet - this one is completely wrong.

Connectez-vous pour commenter.

Réponses (0)

Catégories

En savoir plus sur Mathematics dans Help Center et File Exchange

Tags

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by