mkpp

Делайте кусочный полином

Свернуть все на странице

Синтаксис

pp = mkpp(breaks,coefs)
pp = mkpp(breaks,coefs,d)

Описание

пример

pp = mkpp(breaks,coefs) создает кусочный полином pp от его пропусков и коэффициентов. Использовать ppval для оценки кусочного полинома в определенных точках, или unmkpp чтобы извлечь детали о кусочный полином.

pp = mkpp(breaks,coefs,d) задает, что кусочный полином является векторным, так что значение каждого из его коэффициентов является вектором длины d.

Примеры

свернуть все

Открыть Live Script

Создайте кусочный полином, который имеет кубический полином в интервале [0,4], квадратичный полином в интервале [4,10] и квартальный полином в интервале [10,15].

breaks = [0 4 10 15];
coefs = [0 1 -1 1 1; 0 0 1 -2 53; -1 6 1 4 77];
pp = mkpp(breaks,coefs)
pp = struct with fields:
      form: 'pp'
    breaks: [0 4 10 15]
     coefs: [3x5 double]
    pieces: 3
     order: 5
       dim: 1

Оцените кусочный полином во многих точках интервала [0,15] и постройте график результатов. Постройте график вертикальных штриховых линий в точках пропуска, где сходятся полиномы.

xq = 0:0.01:15;
plot(xq,ppval(pp,xq))
line([4 4],ylim,'LineStyle','--','Color','k')
line([10 10],ylim,'LineStyle','--','Color','k')

Figure contains an axes. The axes contains 3 objects of type line.

Открыть Live Script

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

Первые два подграфиков показывают квадратичный полином, и его отрицание сместилось на интервалы [-8, -4] и [-4,0]. Полином следующий

1-(x2-1)2=-x24+x.

Третий подграфик показывает кусочный полином, построенную путем чередования этих двух квадратичных частей на четырех интервалах. Вертикальные линии добавляются, чтобы показать точки, где полиномы встречаются.

subplot(2,2,1)
cc = [-1/4 1 0]; 
pp1 = mkpp([-8 -4],cc);
xx1 = -8:0.1:-4; 
plot(xx1,ppval(pp1,xx1),'k-')

subplot(2,2,2)
pp2 = mkpp([-4 0],-cc);
xx2 = -4:0.1:0; 
plot(xx2,ppval(pp2,xx2),'k-')

subplot(2,1,2)
pp = mkpp([-8 -4 0 4 8],[cc;-cc;cc;-cc]);
xx = -8:0.1:8;
plot(xx,ppval(pp,xx),'k-')
hold on
line([-4 -4],ylim,'LineStyle','--')
line([0 0],ylim,'LineStyle','--')
line([4 4],ylim,'LineStyle','--')
hold off

Figure contains 3 axes. Axes 1 contains an object of type line. Axes 2 contains an object of type line. Axes 3 contains 4 objects of type line.

Входные параметры

свернуть все

Пропуски, заданные как вектор длины L+1 со строго увеличивающимися элементами, которые представляют начало и конец каждого из L интервалы.

Типы данных: single | double

Полиномиальные коэффициенты, заданные как L-by- k матрица с i-й строкой coefs(i,:) содержит локальные коэффициенты порядка k полином i-го интервала, [breaks(i), breaks(i+1)]. Другими словами, полином coefs(i,1)*(X-breaks(i))^(k-1) + coefs(i,2)*(X-breaks(i))^(k-2) + ... + coefs(i,k-1)*(X-breaks(i)) + coefs(i,k).

Типы данных: single | double

Размерность, заданная как скаляр или вектор целых чисел. Задайте d чтобы обозначить, что кусочный полином имеет значения коэффициентов размера d.

Типы данных: single | double

Выходные аргументы

свернуть все

Кусочный полином, возвращается как структура. Используйте эту структуру с ppval функция для вычисления кусочного полинома в одной или нескольких точках запроса. Структура имеет эти поля.

ОбластьОписание
form

'pp' для кусочного полинома

breaks

Вектор длины L+1 со строго увеличивающимися элементами, которые представляют начало и конец каждого из L интервалы

coefs

L-by- k матрица с каждой строкой coefs(i,:) содержит локальные коэффициенты порядка k полином i-го интервала, [breaks(i),breaks(i+1)]

pieces

Количество штук, L

order

Порядок полиномов

dim

Размерность цели

Поскольку полиномиальные коэффициенты в coefs являются локальными коэффициентами для каждого интервала, необходимо вычесть нижнюю конечную точку соответствующего узлового интервала, чтобы использовать коэффициенты в обычном полиномиальном уравнении. Другими словами, для коэффициентов [a,b,c,d] по интервалу [x1,x2], соответствующий полином является

f(x)=a(x−x1)3+b(x−x1)2+c(x−x1)+d .

Расширенные возможности

.

См. также

| | |

Представлено до R2006a

Документация MATLAB

Поддержка

Для просмотра документации необходимо авторизоваться на сайте