?clear;clc;close?all;
xx=[-1??0??3??5??8??9?];
yy=[1??7??-4??0??7??3??];????%給定6個點



sp=spline（xxyy）;????%?spline（xy）??曲線插值????

x=xx（1）:0.01:xx（length（xx））;??????%通常取x（1）至x（n）?
y=ppval（xsp）;????%?ppval（??）???分段多項式的值??y的值變為數組y（）

subplot（311）
for?i=1:length（xx）
???plot（xx（i）yy（i）‘b*‘）?
???hold?on
end
plot（xy‘r‘）;
grid?on
title（‘B樣條插值曲線‘）
??
???for?i=1:（length（x）-1）
???????dx（i）=x（i+1）-x（i）;
???????dy（i）=y（i+1）-y（i）;
???????dddy（i）=?dy（i）/dx（i）;
???end
?????
?
???for?i=1:（length（x）-2）
???????ddx（i）=dx（i+1）-dx（i）;
???????ddy（i）=dy（i+1）-dy（i）;
???????K（i）=（dx（i）*ddy（i）-dy（i）*ddx（i））/（（dx（i）*dx（i）+dy（i）*dy（i））^1.5）;
???end???
??subplot（312）?
?????ud=linspace（xx（1）xx（length（xx））（length（x）-1））;
?????plot（uddddy）
?????grid?on
??title（‘曲線斜率變化圖‘）
??
??subplot（313）
??uu=linspace（xx（1）xx（length（xx））（length（x）-2））;
??plot（uuK）
??grid?on
???title（‘曲線曲率變化圖‘）

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件?????????962??2020-11-17?17:17??B-Spline_qulv.m