бчислювальної (інструментальної) ? в похибок цілочисельного алгоритму обчислення функції для поточних значень аргументу x з інтервалу апроксимації.
В
7. Опис і лістинг програми
У даному розділі наводиться MatLab-програма, що реалізує розглянуту вище модель для експериментального аналізу повної ? , методичної ? м та обчислювальної (інструментальної) ? в похибок цілочисельного алгоритму обчислення функції.
Лістинг програми виглядає наступним чином:
clear
% Обчислення табличних значень аргументу і функції з кроком h = 2 ^ -3;% крок таблиці = 0.5: h: 1;% таблиця значень аргументу = sqrt (xs);% таблиця значень функції = 2 ^ 7; = 2 ^ 7; = round (xs * Mx);% таблиця цілочисельних значень аргументу = round (fs * Mf);% таблиця цілочисельних значень функції
% Обчислення поточних значень функції з кроком h/4 = 0.5: h/4:1; = sqrt (x); (5,1,1); plot (x, fx, 'g' )% графік функцііon; 'Graphic sqrt (x)'; = round (x * Mx); = floor ((x-0.5)/h) +1;% значення індексу
% обчислення не масштабованих значень функції
% за методом лінійної інтерполяцііj = 1:1:17 (i (j) <5) (j) = fs (i (j)) + (x (j)-xs (i ( j))) * ((fs (i (j) +1)-fs (i (j)))/h); (j) = fs (i (j));;
% Обчислення значень функції з целочисленному алгоритмом (i (j) <5) (j) = FS (i (j)) + round (2 ^ -4 * (X (j)-XS (i (j))) * (FS (i (j) +1)-FS (i (j)))); (j) = FS (i (j));;
% Демаcштабірованіе значень функції = F/Mf; (5,1,2); plot (x, fdem, 'g'); on; 'Graph f demash';
% Повна похибка = sqrt (x) - fdem; (5,1,3); plot (x, eps, 'b'); on; 'Graph polnoi pogreshnosti';
% Методична погрешность_m = sqrt (x) - f; (5,1,4); plot (x, e_m, 'r'); on; 'Graphic methodich pogreshnosti';
% Інструментальна погрешность_instr = fF/Mf; (5,1,5); plot (x, e_instr, '-k'); on; 'Graphic instrum pogreshnosti';
. Результати моделювання
Нижче зображені графіки відповідно функції, повної ? , методичної ? м та обчислювальної (інструментальної) ? в похибок, отримані в рез...
