Hello everybody i have been trying to code a 1D wave problem using the Leap Frog Scheme but i am running into an error that i am not familiar with. Any suggestions?freq = 1; % Hz
omega = 2*pi*freq; % angular speed m/s
c = 10; % m/s
amp = 1;% Space
dx = 0.001;xmax = 1; % domain of unit length
lambda = 0.05; % size of wave
k = 2*pi/lambda; % wavenumber
x = 0:dx:xmax;Vx = length(x); % matrix of propagation in the x direction
u_0 = zeros(size(x)); % defining wave initial conditions
for i = 1:Vx; if x(i) < lambda; u_0(i) = (1 - cos(k*x(i)))*0.5; else break endend %Time
dt = 0.01; t = 0:dt:1; Vt = length(t);figure(1)plot(x,u_0,'.-b');axis([0 lambda 0 1])phi = zeros(Vx,Vt);for n = 1:Vx; phi(1,i) = u_0(i); phi(2,i) = u_0(i); end for i = 2:Vt for n = 2:Vx-1 phi(i,n+1) = phi(i,n-1) + (u_0*(dt/dx))*(phi(i+1,n) - phi(i-1,n)/phi(i,n-1)); end end
MATLAB: Subscripted assignment dimension mismatch – Help
time
Related Question
- Index in position 1 exceeds array bounds (must not exceed 301) and Error in IMR_2psv (line 34) u(i,:)=(1/16)*(-u(i-2,:)+4*u(i-1,:)+10*u(i,:)+4*u(i+1,:)-u(i+2,:));
- In-Situ Combustion Model – MATLAB code
- Exceeds the number of array
- How to ‘crop’ a graph
- How to plot a solution from a PDE within a unit square
Best Answer