Как я могу построить производную функцию кубической сплайновой интерполяции?

Question

Я хочу спросить, как я могу построить производную функцию кубической сплайновой интерполяции? Я показываю, как я это сделал в своем коде. Ошибка:

Ошибка с помощью сюжета недопустимого аргумента вторых данных Ошибки в function_MTU4000_Real (строку 90) Plot1 = участок (x1, SPEED1, 'б');

%calculation of lifting of intake valve (approximation spline function) 

x1=0.0:0.1:202.1; 
y1=xlsread('Steuerzeiten_Schrittweite_MTU4000.xlsx',1,'D2:D2023'); 
Lifting1=spline(x1,y1); 
x2=202.1:0.1:701.9; 
Lifting2=0*x2; 
x3=702.0:0.1:720.0; 
y3=xlsread('Steuerzeiten_Schrittweite_MTU4000.xlsx',1,'D7022:D7202'); 
Lifting3=spline(x3,y3); 

%calculation and plot of speed intake 

figure(2);hold on; grid on; 
Speed1=fnder(Lifting1); 
plot1=plot(x1,Speed1,'b'); 
Speed2=Lifting2; 
plot2=plot(x2,Speed2,'b'); 
Speed3=fnder(Lifting3); 
plot3=plot(x3,Speed3,'b'); 
hold off 
legend([plot1,plot2,plot3],'Intake') 
set(gca,'XTickLabel',{'OT','90','UT','270','ZOT','450','UT','630','OT'}); 
title('Intake Valve Speed') 
xlabel('Crank Angle [°]') 
ylabel('Speed [m/°]') 

Answer 1

fnder возвращает структуру, а не массив для построения графика. Вы должны использовать ppval или fnval оценить его ...

plot1 = plot(x1, ppval(Speed1,x1)) 

Docs:

https://lost-contact.mit.edu/afs/cs.stanford.edu/pkg/matlab-r2015b/matlab/r2015b/help/curvefit/examples/cubic-spline-interpolation.html

Для использования csapi см документации: https://uk.mathworks.com/help/curvefit/csapi.html

% docs example 
bcs = csapi({x,y}, z); 
fnplt(bcs)