Определите точки образца, xи соответствующие выборочные значения, v.

x = 0:pi/4:2*pi; 
v = sin(x);

Определите точки запроса для более точной выборки в диапазоне x.

xq = 0:pi/16:2*pi;

Выполните интерполяцию функции в точках запроса и постройте график результата.

figure
vq1 = interp1(x,v,xq);
plot(x,v,'o',xq,vq1,':.');
xlim([0 2*pi]);
title('(Default) Linear Interpolation');

Теперь оцените v в тех же точках, используя 'spline' способ.

figure
vq2 = interp1(x,v,xq,'spline');
plot(x,v,'o',xq,vq2,':.');
xlim([0 2*pi]);
title('Spline Interpolation');

Определите набор значений функций.

v = [0  1.41  2  1.41  0  -1.41  -2  -1.41 0];

Определите набор точек запроса, которые попадают между точками по умолчанию. 1:9. В этом случае точками по умолчанию являются 1:9 потому что v содержит 9 значения.

xq = 1.5:8.5;

Оценить v в xq.

vq = interp1(v,xq);

Постройте график результата.

figure
plot((1:9),v,'o',xq,vq,'*');
legend('v','vq');

Определите набор точек образца.

x = 1:10;

Определите значения функции v $$(x)\mathrm{}={}^{\mathrm{5x}}$$+ x2i в точках выборки.

v = (5*x)+(x.^2*1i);

Определите точки запроса для более точной выборки в диапазоне x.

xq = 1:0.25:10;

Интерполировать v в точках запроса.

vq = interp1(x,v,xq);

Постройте график реальной части результата красным цветом, а воображаемой - синим.

figure
plot(x,real(v),'*r',xq,real(vq),'-r');
hold on
plot(x,imag(v),'*b',xq,imag(vq),'-b');

Интерполяция точек данных с временной меткой.

Рассмотрим набор данных, содержащий показания температуры, которые измеряются каждые четыре часа. Создайте таблицу с данными за один день и постройте график данных.

x = (datetime(2016,1,1):hours(4):datetime(2016,1,2))';
x.Format = 'MMM dd, HH:mm';
T = [31 25 24 41 43 33 31]';
WeatherData = table(x,T,'VariableNames',{'Time','Temperature'})

WeatherData=7×2 table
        Time         Temperature
    _____________    ___________

    Jan 01, 00:00        31     
    Jan 01, 04:00        25     
    Jan 01, 08:00        24     
    Jan 01, 12:00        41     
    Jan 01, 16:00        43     
    Jan 01, 20:00        33     
    Jan 02, 00:00        31     

plot(WeatherData.Time, WeatherData.Temperature, 'o')

Интерполяция набора данных для прогнозирования показаний температуры в течение каждой минуты дня. Поскольку данные являются периодическими, используйте 'spline' способ интерполяции.

xq = (datetime(2016,1,1):minutes(1):datetime(2016,1,2))';
V = interp1(WeatherData.Time, WeatherData.Temperature, xq, 'spline');

Постройте график интерполированных точек.

hold on
plot(xq,V,'r')

Определите точки образца, xи соответствующие выборочные значения, v.

x = [1 2 3 4 5];
v = [12 16 31 10 6];

Укажите точки запроса, xq, которые выходят за пределы области x.

xq = [0 0.5 1.5 5.5 6];

Оценить v в xq с использованием 'pchip' способ.

vq1 = interp1(x,v,xq,'pchip')

vq1 = 1×5

   19.3684   13.6316   13.2105    7.4800   12.5600

Далее оцените v в xq с использованием 'linear' способ.

vq2 = interp1(x,v,xq,'linear')

vq2 = 1×5

   NaN   NaN    14   NaN   NaN

Теперь используйте 'linear' с помощью 'extrap' вариант.

vq3 = interp1(x,v,xq,'linear','extrap')

vq3 = 1×5

     8    10    14     4     2

'pchip' экстраполяция по умолчанию, но 'linear' не делает.

Определите точки образца, xи соответствующие выборочные значения, v.

x = [-3 -2 -1 0 1 2 3];
v = 3*x.^2;

Укажите точки запроса, xq, которые выходят за пределы области x.

xq = [-4 -2.5 -0.5 0.5 2.5 4];

Теперь оцените v в xq с использованием 'pchip' метод и назначить любые значения за пределами области x к значению, 27.

vq = interp1(x,v,xq,'pchip',27)

vq = 1×6

   27.0000   18.6562    0.9375    0.9375   18.6562   27.0000

Определите точки образца.

x = (-5:5)';

Пример трех различных параболических функций в точках, определенных в x.

v1 = x.^2;
v2 = 2*x.^2 + 2;
v3 = 3*x.^2 + 4;

Создать матрицу v, столбцы которых являются векторами, v1, v2, и v3.

v = [v1 v2 v3];

Определение набора точек запроса, xq, чтобы быть более тонкой выборкой в диапазоне x.

xq = -5:0.1:5;

Оценить все три функции в xq и постройте график результатов.

vq = interp1(x,v,xq,'pchip');
figure
plot(x,v,'o',xq,vq);

h = gca;
h.XTick = -5:5;

Круги на графике представляют v, и сплошные линии представляют vq.