最新三次样条插值的Matlab实现(自然边界和第一边界条件).doc

《最新三次样条插值的Matlab实现(自然边界和第一边界条件).doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、__________________________________________________（第一边界条件）源代码：functiony=yt1(x0,y0,f_0,f_n,x)_____________(1)%第一类边界条件下三次样条插值；%xi所求点；%yi所求点函数值；%x已知插值点；%y已知插值点函数值；%f_0左端点一次导数值；%f_n右端点一次导数值；n=length(x0);z=length(y0);h=zeros(n-1,1);k=zeros(n-2,1);l=zeros(n-2,1);S=2*eye(n);fori=1

2、:n-1h(i)=x0(i+1)-x0(i);endfori=1:n-2k(i)=h(i+1)/(h(i+1)+h(i));l(i)=1-k(i);end%对于第一种边界条件：k=[1;k];_______________________(2)l=[l;1];_______________________(3)%构建系数矩阵S：fori=1:n-1S(i,i+1)=k(i);S(i+1,i)=l(i);end%建立均差表：F=zeros(n-1,2);fori=1:n-1F(i,1)=(y0(i+1)-y0(i))/(x0(i+1)-x0(i

3、));endD=zeros(n-2,1);收集于网络，如有侵权请联系管理员删除__________________________________________________fori=1:n-2F(i,2)=(F(i+1,1)-F(i,1))/(x0(i+2)-x0(i));D(i,1)=6*F(i,2);end%构建函数D：d0=6*(F(1,2)-f_0)/h(1);___________(4)dn=6*(f_n-F(n-1,2))/h(n-1);___________(5)D=[d0;D;dn];______________(6)m

4、=SD;%寻找x所在位置，并求出对应插值：fori=1:length(x)forj=1:n-1if(x(i)<=x0(j+1))&(x(i)>=x0(j))y(i)=(m(j)*(x0(j+1)-x(i))^3)/(6*h(j))+...(m(j+1)*(x(i)-x0(j))^3)/(6*h(j))+...(y0(j)-(m(j)*h(j)^2)/6)*(x0(j+1)-x(i))/h(j)+...(y0(j+1)-(m(j+1)*h(j)^2)/6)*(x(i)-x0(j))/h(j);break;elsecontinue;endend

5、end（2）（自然边界条件）源代码：仅仅需要对上面部分标注的位置做如下修改:__(1):functiony=yt2(x0,y0,x)__(2):k=[0;k]__(3):l=[l;0]__(4)+(5):删除—(6):D=[0:D:0]流程图、PFMEA、控制计划的关系过程编号和描述变差来源过程流程图产品特性过程特性收集于网络，如有侵权请联系管理员删除__________________________________________________项目功能潜在失效模式潜在失效后果严重度S级别潜在失效起因/机理频度O现行设计控制预防现行设计控

6、制探测探测度DRPN建议措施责任及目标完成日期措施结果采取的措施SODRPN零件/过程编号过程名称/操作描述生产设备材料特性特殊特性分类方法反应计划和纠正措施编号产品过程产品/过程规范/公差评价/测量技术样本控制方法容量频率收集于网络，如有侵权请联系管理员删除