cass10样条拟合快捷命令

cass10样条拟合快捷命令表示如下：

1、DD通用绘图命令。

2、A画弧ARC。

3、RR符号重新生成。

4、FF绘制多点房屋。

5、PE复合线编辑PEDIT。

如下：

function

fit_tulun 

%多项式拟合方式

x=linspace(0,4*pi,6) 

y_jingque=sin(x)

%决定稀疏样本点数据

p_poly=polyfit(x,y_jingque,5) 

x_poly_fit=linspace(0,4*pi,100) 

y_poly_fit=polyval(p_poly,x_poly_fit) 

%三次样条拟合方式

sp=csapi(x,y_jingque) 

%求三次样条函数的导数. 

s_diff=fnder(sp,1) 

plot(x_poly_fit,y_poly_fit,'ko',x_poly_fit,y_poly_fit,'b:') 

%plot(x_poly_fit,y_poly_fit,'b:') 

hold

on 

fnplt(sp,'r') 

fnplt(s_diff,'c') 

x1=linspace(0,4*pi,200) 

plot(x1,sin(x1),'m','linewidth',1.8) 

legend('多项式拟合样本点','多项式拟合曲线','三次样条拟合曲线','三次样条导数曲线','正弦曲线精确图形').

答案是可以得到具体表达式。

先运行sp=spaps(x,y,tol)；

得到结果后运行

>>sp.coefs

得到每个多项式的系数向量，其实也就是多项式函数的具体表达式，相信你能看得出来。

---