资源描述:
《三次样条插值,数值计算,上机报告,c》由会员上传分享,免费在线阅读,更多相关内容在应用文档-天天文库。
1、为了适应公司新战略的发展,保障停车场安保新项目的正常、顺利开展,特制定安保从业人员的业务技能及个人素质的培训计划三次样条插值,数值计算,上机报告,c 1、题目 x964 y 以上数据点的插值可以得到平方根函数的近似,在区间[0,64]上作图 用三次样条程序求S(x) 2、解题方法 在TC中建立一维数组,用以记录X,Y,h,μ,λ,d各值,建立二维数组用以记录三对角矩阵。将数据排列在三对角矩阵中,根据高斯消去法的原理解方程,算出各M值。通过解出的M值和h值计算三次样条的方程。 3、计算结果 区间[0,1],[1,4],[4,9],[9,16],[
2、16,25],[25,36],[36,49],[49,64]上的三次样条函数分别是 *(1-x)*(1-x)*(1-x)-*x*x**(1-x)+*x -*(4-x)*(4-x)*(4-x)+*(x-1)*(x-1)*(x-1)+*(4-x)-*(x-1)*(9-x)*(9-x)*(9-x)-*(x-4)*(x-4)*(x-4)-*(9-x)+*(x-4)目的-通过该培训员工可对保安行业有初步了解,并感受到安保行业的发展的巨大潜力,可提升其的专业水平,并确保其在这个行业的安全感。为了适应公司新战略的发展,保障停车场安保新项目的正常、顺利开展,特制定安保从业人
3、员的业务技能及个人素质的培训计划 -*(16-x)*(16-x)*(16-x)+*(x-9)*(x-9)*(x-9)+*(16-x)+*(x-9)*(25-x)*(25-x)*(25-x)-*(x-16)*(x-16)*(x-16)+*(25-x)+*(x-16)-*(36-x)*(36-x)*(36-x)+*(x-25)*(x-25)*(x-25)+*(36-x)+*(x-25) *(49-x)*(49-x)*(49-x)+*(x-36)*(x-36)*(x-36)+*(49-x)+*(x-36)*(64-x)*(64-x)*(64-x)-*(x-49)*
4、(x-49)*(x-49)+*(64-x)+*(x-49)C程运行结果截图 Matlab画图结果截图 4、结果分析 运算过程中f' 的值不存在,当x趋向于0时,fx '趋近于正无穷,因 此f' 取值不准确,在计算中,f' 取了100,因此前几个区间上的三 次样条误差较大,而随着自变量的增大,误差减小。与拉格朗日插值 相比三次样条插值对于此题的估计更为准确。如下图所示,绿色线条 表示的是八次拉格朗日插值,尽管在x较小时拟合比较准确,但随着目的-通过该培训员工可对保安行业有初步了解,并感受到安保行业的发展的巨大潜力,可提升其的专业水平,并确保其
5、在这个行业的安全感。为了适应公司新战略的发展,保障停车场安保新项目的正常、顺利开展,特制定安保从业人员的业务技能及个人素质的培训计划 自变量的增大,误差明显变大,稳定性较差。 5、源程序 #include main() { doublex[9]={0,1,4,9,16,25,36,49,64}; doubley[9]={0,1,2,3,4,5,6,7,8}; doublet,s,matrix[9][10],h[8],u[9],v[8],xishu[8][4]; inti,j; for(j=0;j=0;j--){ s=matrix[j][j+
6、1]; matrix[j][i]=matrix[j][i]-matrix[j+1][i]*s; } } for(j=0;j>%定义一个语句函数g(x)=exp(x^10*exp(x-1)) I= 例秦九韶算法 a0=3, ak=2ak-1+3, Pn(x)=anx^n+an-1x^(n-1)+...+a1x+a0目的-通过该培训员工可对保安行业有初步了解,并感受到安保行业的发展的巨大潜力,可提升其的专业水平,并确保其在这个行业的安全感。为了适应公司新战略的发展,保障停车场安保新项目的正常、顺利开展,特制定安保从业人员的业务技能及个人素质的培训计
7、划 求I1=P100(),I2=P150(13) >>x=input('x='); n=input('n='); a=3; fori=1:n a=2*a+3; end s=z; b=(a-3)/2; form=1:100 s=x*s+b; b=(b-3)/2; end disp(s); >>x= n=100 >>x=3 n=100 结论:目的-通过该培训员工可对保安行业有初步了解,并感受到安保行业的发展的巨大潜力,可提升其的专业水平,并确保其在这个行业的安全感。为了适应公司新战略的发展,保障停车场安保新项目的正常、顺利开展,
8、特制定安保从业人员的业务