Hi,
From my understanding, you are attempting to simulate the behaviour of a wave over time using a numerical solution to the wave equation in MATLAB. Your goal is to set the initial condition of the wave function u(x,0) such that it equals 1 within a specific range of the spatial domain (pi - 1 <= x <= pi + 1) and 0 elsewhere. However, you are not sure about how to implement this piecewise initial condition within the MATLAB function that you have written.
To implement the initial condition u(x,0)=1 for pi-1 <= x <= pi+1 and 0 otherwise, you need to modify the function f1(x) that represents the initial condition u(x,0) in your script. You can do this by defining ‘f1’ as an anonymous function that applies the given condition to the input x.
Here's how you can modify the ‘f1’ variable in your script:
f1 = @(x) double((x >= pi - 1) & (x <= pi + 1));
This function returns 1 when x is between pi - 1 and pi + 1, and 0 otherwise. The double function is used to ensure that the result is a floating-point number since logical expressions in MATLAB return logical true or false by default.
The final modified script will be as follows:
clc
clear all
% Define the initial condition function
f1 = @(x) double((x >= pi - 1) & (x <= pi + 1));
% Define the second initial condition derivative function
f2 = @(x) x.^2.*((sin(4.*x)).^2);
% Define the boundary condition functions
g0 = @(t) 0;
g1 = @(t) 0;
% Define the spatial and temporal domains
xspan = [0, 2*pi];
tspan = [0, 1];
% Define the number of spatial points and time points
nx = 20;
nt = 40;
% Define the wave speed squared
alpha = 4;
% Call the Wave function
[u, q] = Wave(f1, f2, g0, g1, xspan, tspan, nx, nt, alpha);
With this modification, the initial condition u(x,0) will be set to 1 within the specified range and 0 elsewhere. The Wave function will then use this initial condition to compute the evolution of u over time and space.
I hope this helps!
Sameer
