Showing Page:
1/4
Отчет по учебной практике на тему:
« ОСНОВЫ РАБОТЫ С MATLAB »
Конечные разности для волнового уравнения Реализация простых
численных схем для волнового уравнения.
Применяя центрированные разности второго порядка для аппроксимации производных,,
Можно организовать получение
Showing Page:
2/4
Proboem Setup
N = 101;
Showing Page:
3/4
Initial Condition
Boundary condition
L = 4*pi; x =
linspace(0,L,N);
% It has three data
set; 1: past, 2:
current, 3: future.
u = zeros(N,3); s =
axis([0,L,-2,2]);
xlabel('x'); ylabel('u');
title('Wave equation');
% Dirichet Boundary conditions
u(1,:) = 0;
u(end,:) = 0;
% Gaussian Pulse y =
2*exp(-(x-L/2).^2);
u(:,1) = y;
u(:,2) = y;
% Plot the initial condition. handle_line
Showing Page:
4/4
Simulation
filename = 'wave.gif';
for ii=1:100
disp(['at ii= ', num2str(ii)]); u(2:end-
1,3) = s*(u(3:end,2)+u(1:end-2,2)) ...
+ 2*(1-s)*u(2:end-1,2) ...
- u(2:end-1,1);
u(:,1) = u(:,2);
u(:,2) = u(:,3);
handle_line.YData = u(:,2);
drawnow;
frame = getframe(gcf);
im = frame2im(frame);
imwrite(A,map,filename,'gif','WriteMode','append','DelayTime',0.05);
end
end

Unformatted Attachment Preview

Отчет по учебной практике на тему: « ОСНОВЫ РАБОТЫ С MATLAB » Конечные разности для волнового уравнения Реализацияпростых численных схем для волнового уравнения. Применяя центрированные разности второго порядка для аппроксимации производных,, Можно организовать получение Proboem Setup N = 101; L = 4*pi; x = linspace(0,L,N); % It has three data set; 1: past, 2: current, 3: future. u = zeros(N,3); s = Initial Condition axis([0,L,-2,2]); xlabel('x'); ylabel('u'); title('Wave equation'); Boundary condition % Dirichet Boundary conditions u(1,:) = 0; u(end,:) = 0; % Gaussian Pulse y = 2*exp(-(x-L/2).^2); u(:,1) = y; u(:,2) = y; % Plot the initial condition. handle_line Simulation filename = 'wave.gif'; for ii=1:100 disp(['at ii= ', num2str(ii)]); u(2:end- 1,3) = s*(u(3:end,2)+u(1:end-2,2)) ... + 2*(1-s)*u(2:end-1,2) ... - u(2:end-1,1); u(:,1) = u(:,2); u(:,2) = u(:,3); handle_line.YData = u(:,2); drawnow; frame = getframe(gcf); im = frame2im(frame); imwrite(A,map,filename,'gif','WriteMode','append','DelayTime',0.05); end end Name: Description: ...
User generated content is uploaded by users for the purposes of learning and should be used following Studypool's honor code & terms of service.
Studypool
4.7
Trustpilot
4.5
Sitejabber
4.4