вилина
u2 = analitic (x, t, 5), (2);
mesh (u2);
%
%% Побудова графіків при обліку 100 гармонік (n = 100)
% Для дрібної сітки: = 0:0.01:1;% 1 - один метр
t = 0:0.01:1;% 1 - одна хвилина = analitic (x, t, 5), (3); (u3);
Додаток 3
file В«pointВ»
clc
clear all
%% Побудова графіків при обліку 5 гармонік (n = 5)
% Для дрібної сітки: = 0:0.01:1;% 1 - один метр
t = 0:0.01:1;% 1 - одна хвилина = analitic (x, t, 5); (1); (u1); ('x'); ('t') ; ('u (x, t)');
title ('Висновок заданої точки x = 0.81 (м), t = 45 (сек) на графік аналітичного рішення задачі')
% = 0.81;% 81 сантиметр = 45/60;% 45 секунд = analitic (x, t, 5);
x = 82;% номер вузла по осі x = 76;% номер вузла по осі t
hold on (x, t, u1, 'k *'); off
%
Додаток 4
file В«numericalВ»
function u = numerical (l, t, a, h);
h = 0.2; = 0.1; = 1; = 6;% по осі x = 11;% по осі t
% = zeros (m, n);
% j = 1: n (1, j) = 0;
% i = 1: m (i, 1) = 0; (i, 2) = 0;
% (m, 2) = 0.01;
% j = 3: ni = 2: m-1 (i, j) = 0.25 * (u (i +1, j-1) + u (i-1, j-1) ) +1.5 * u (i, j-1) - u (i, j-2); (m, j) = u (m-1, j) +0.02 * j; = u '
Додаток 5
file В«n_reshenieВ»
u = numerical (1,1,1,0.2); (1); (u); ('x'); ('t'); ('u (x, t)' );
title ('Чисельне рішення задачі для великої сітки')
u
Додаток 6
file В«small_numericalВ»
function u = small_numerical (l, t, a, h);
a = 1; = h ^ 2/(2 * a ^ 2);% умова стійкості = round (1 + l/h);% по осі x
n = round (1 + t/tau);% по осі t
% = zeros (m, n);
% j = 1: n (1, j) = 0;
% i = 1: m (i, 1) = 0; (i, 2) = 0;
% (m, 2) = 0.01;
% j = 3: ni = 2: m-1 (i, j) = 0.25 * (u (i +1, j-1) + u (i-1, j-1) ) +1.5 * u (i, j-1) - u (i, j-2); (m, j) = u (m-1, j) +0.02 * j; = u ';
Додаток 7
file В«n_reshenie_smВ»
clcall = small_numerical (1,1,1,0.05); (1); (u); (...