用MATLAB实现拉格朗日插值和分段线性插值讲课稿.doc

用MATLAB实现拉格朗日插值和分段线性插值讲课稿.doc

ID:60808714

大小:24.50 KB

页数:7页

时间:2020-12-20

用MATLAB实现拉格朗日插值和分段线性插值讲课稿.doc_第1页
用MATLAB实现拉格朗日插值和分段线性插值讲课稿.doc_第2页
用MATLAB实现拉格朗日插值和分段线性插值讲课稿.doc_第3页
用MATLAB实现拉格朗日插值和分段线性插值讲课稿.doc_第4页
用MATLAB实现拉格朗日插值和分段线性插值讲课稿.doc_第5页
资源描述:

《用MATLAB实现拉格朗日插值和分段线性插值讲课稿.doc》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库

1、精品好文档,推荐学习交流用MATLAB实现拉格朗日插值和分段线性插值1、实验内容:     用MATLAB实现拉格朗日插值和分段线性插值。2、实验目的:1)        学会使用MATLAB软件;2)        会使用MATLAB软件进行拉格朗日插值算法和分段线性差值算法;3、实验原理:利用拉格朗日插值方法进行多项式插值,并将图形显式出来。4、实验步骤及运行结果(1)实现lagrange插值1)定义函数:f=1/(x^2+1) 将其保存在f.m文件中,具体程序如下:functiony=f1(x)y=1./(x.^2+1);

2、2) 定义拉格朗日插值函数:将其保存在lagrange.m文件中,具体实现程序编程如下:functiony=lagrange(x0,y0,x)m=length(x);    /区间长度/n=length(x0);仅供学习与交流,如有侵权请联系网站删除谢谢7精品好文档,推荐学习交流fori=1:n l(i)=1;endfori=1:mforj=1:nfork=1:n ifj==k    continue; endl(j)=(x(i)-x0(k))/(x0(j)-x0(k))*l(j);end endend y=0;fori=1:n

3、 y=y0(i)*l(i)+y;    end3) 建立测试程序,保存在text.m文件中,实现画图:x=-5:0.001:5;y=(1+x.^2).^-1;p=polyfit(x,y,n);py=vpa(poly2sym(p),10)plot_x=-5:0.001:5;f1=polyval(p,plot_x);figureplot(x,y,‘r',plot_x,f1)输入n=6,出现下面的图形:通过上图可以看到当n=6是没有很好的模拟。仅供学习与交流,如有侵权请联系网站删除谢谢7精品好文档,推荐学习交流于是重新运行text.M

4、并选择n=11由此可见n=11时的图像是可以很好的实现模拟(2)分段线性插值:   建立div_linear.m文件。具体编程如下/*分段线性插值函数:div_linear.m文件*/functiony=div_linear(x0,y0,x,n)%forj=1:length(x)fori=1:n-1   if(x>=x0(i))&&(x<=x0(i+1))      y =(x-x0(i+1))/(x0(i)-x0(i+1))*y0(i)+(x-x0(i))/(x0(i+1)-x0(i))*y0(i+1);   else    

5、 continue;   endend%end  测试程序(text2.m):   n=input(‘输入n=:’);x0=linspace(-5,5,n);forx=-5:0.01:5    y=div_linear(x0,f(x0),x,n);   holdon; plot(x,y,'r');   plot(x,f(x),'b');end仅供学习与交流,如有侵权请联系网站删除谢谢7精品好文档,推荐学习交流2)运行测试程序,这是会出现:输入n=:2)输入n=6,并按Enter键,出现:4)关掉图形界面后,重新运行程序,输入n=

6、11,并按enter键后出现:5)再次关掉图形界面,输入n=100,并按enter键,出现:此时。图形将于原函数图形基本吻合,说明分割区间越多,图像接近真实的图像。(3)用lagrange插值观察y=

7、sin(k*π*x)

8、的误差分析:1)编写函数文件,保存在f2.m中x=0:0.01:1;k=input('输入k:')n=input('输入n:');y=abs(sin(k*pi*x));p=polyfit(x,y,n-1);py=vpa(poly2sym(p),8);plot_x=0:0.01:1;f1=polyval(p,p

9、lot_x);plot(x,y,plot_x,f1);2)运行该程序:输入k=:1仅供学习与交流,如有侵权请联系网站删除谢谢7精品好文档,推荐学习交流输入n=:2出现如下图形界面:关掉图形界面后重新运行f2.m,输入k=:1,n=:3出现如下界面: 再次关掉图形界面,输入k=:1,n=:6后出现:此时图形基本吻合。类推,输入k=2,n=3后出现:                 k=2,n=11,出现如下图形:                          k=2,n=15,这时出现:k=2,n=19,出现:当k=2,n=2

10、1时,图形如下:此时基本吻合。5、实验总结:仅供学习与交流,如有侵权请联系网站删除谢谢7精品好文档,推荐学习交流   通过本次课程设计,我初步掌握了MATLAB运用,加深了对于各种线性插值的理解;培养了独立工作能力和创造力;综合运用专业及基础知识,解决实际数学问

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。