"> В§
I = trapz (x, y) Функція I = trapz (x, y) обчислює інтеграл від функції y по змінній x, використовуючи метод трапецій. Аргументи x і y можуть бути одновимірними масивами однакового розміру, або масив Y може бути двовимірним, але тоді має виконуватися умова size (Y, 1) = length (x) . В останньому випадку обчислюється інтеграл для кожного стовпця.
Текст файл-програми interpoljacija для вирішення даної задачі приведений в лістингу 1.
Лістинг 1. Файл-програма interpoljacija
% заповнення масивів для табличній функції x = [-0.5 -0.4677 -0.4343 -0.3952 -0.3493 -0.2951 -0.2278 -0.1431 -0.0394 0.0803 0.2099 0.3445 0.4816 0.6225 0.7647 0.9009 1]; y = [16 14.6264 13.2183 11.8263 10.4538 9.1135 7.8345 6.6780 5.7083 4.9382 4.3355 3.8454 3.4550 3.1773 2.9703 2.7792 2.7]; % висновок графіка табличній функції маркерами plot (x, y, 'ko') % завдання промежyточниx точок для інтерполяції < span align = "justify"> xi = [x (1): 0.1: x (length (x))]; % обчислення кyсочно-лінійної фyнкціі в промежyточниx точкаx yline = interp1 (x, y, xi, 'linear'); % обчислення кyбіческого сплайна в промежyточниx точкаx yspline = interp1 (x, y, xi, 'spline'); % для побудови графіків в одному вікні hold on % нанесення вимірювальної сітки grid on % побудова графіків інтерполянтов plot (xi, yline, 'k-', xi, yspline, 'k -.') % нанесення назв координатних осей % висновок легенди legend ('таблична фyнкція', 'лінійна (linear)', 'кyбіческій сплайн (spline)') % завдання інтервалу x1 = -0.5:0.1:1; % завдання функції y1 = interp1 (x, y, x1, 'spline'); % знаходження площі під графіком S = trapz (x1, y1) інтерполяція інженерний сплайн matlab
Виконання файл-програми interpoljacija призводить до появі графіка, зображеного на рис. 2., Що відображає графік табличній функції і графіки інтерполянт...
