MATLAB---csape斜率拟合

function mycsap1()
%指定斜率csape显式拟合

%把测量得到的数据放入程序
 px= [0	0.6283	1.2566	1.8850	2.5133	3.1416	3.7699	4.3982	5.0265	5.6549	6.2832];
 py= [1.0000	0.8090	0.3090	-0.3090	-0.8090	-1.0000	-0.8090	-0.3090	0.3090	0.8090	1.0000];
 valconds=[0,0];%这是首末端点的斜率
 %把测量得到的数据放入程序
pp = csape(px,py,'complete',valconds);
x=0:6.2832/60:6.2832;
%规定密化点的横坐标
y=fnval(pp,x);%计算密化点的纵坐标
plot(x,y,'b.','MarkerSize',8);%绘制密化点便于观察
hold on
plot(x,y,'g');%把密化点连接成绿颜色的曲线
plot(px,py,'k.','MarkerSize',15);
%绘制型值点
hold off
axis equal
%%%%%%%%%%%%%%%%%%%%%%

function mycsap2()
%指定斜率csape显式拟合

%把测量得到的数据放入程序
 px= [-1.000	-0.8068	-0.3078 0.3078	0.8068	1.000];
 py= [0.0063 0.5908	0.9514	0.9514	0.5908	0.0063];
 valconds=[159.1529, -159.1529];%这是首末端点的斜率
 %把测量得到的数据放入程序
pp = csape(px,py,'complete',valconds);

x=1.0:-2.0/50:-1.0;%注意起始点的位置和点列的走向
%规定密化点的横坐标
y=fnval(pp,x);%计算密化点的纵坐标
plot(x,y,'b.','MarkerSize',8);%绘制密化点便于观察
hold on
plot(x,y,'g');%把密化点连接成绿颜色的曲线
plot(px,py,'k.','MarkerSize',15);
%绘制型值点
hold off
axis equal
%%%%%%%%%%%%%%%%%%%%%%

function mycsap3()
%指定斜率csape拟合
%参数方程拟合
 px= [-1.000	-0.8068	-0.3078 0.3078	0.8068	1.000];
 py= [0.0063