1、1实验 3 Matlab 编程实现 Lagrange 插值算法复习:1、 输出一个正整数,求该正整数的阶乘。函数参考:2、编写函数实现对任意输入一个向量的排序(向量里的元素从小到大)函数参考:2Lagrange 插值算法一、理论知识:1、线性插值 101010)( yxyxL2、二次插值 210)()()(lll, , )(20100xxl )(21011xxl )()(1202xxl3、 n 次 Lagrange 插值 nkknylylyllL010 )()()()( )()()()( 111 nkkkkokk xxxxl nk kkj jjn yxxL0)()(二、实验题目:1、 已知 ,
2、 , ,用线性和二次插值求 的近似值。124395线性插值你选择的节点是:4.9你的程序:(5-9)/(4-9)*2+(5-4)/(9-4)*3插值结果:2.20003二次插值你的程序:(5-4)*(5-9)/(1-4)*(1-9)*1+(5-1)*(5-9)/(4-1)*(4-9)*2+(5-1)*(5-4)/(9-1)*(9-4)*3插值结果:2.26672、 编写函数,实现拉格朗日插值多项式求近似值,并用你的程序验证 1 题。要求: 1、程序具有验证插值节点互异的功能。2、函数头如下:function yi=Lagrange(x,y,xi)%x 为向量,全部的插值节点%y 为向量,插值节
3、点处的函数值%xi 为被估计函数自变量,可以为单个值,也可以为向量。%yi 为 xi 处的函数估计值end3、可能会用的函数:length(x) % 向量 x 中分量的个数4、 的实现:=你的函数:function yi=Lagrange(x,y,xi)%x%y%xi%yixin=length(x);% x=hy(x);% for i=1:n-1% if xi=xi+1% disp()% return% end% end4for i=1:nfor j=1:nif x(i)=x(i+1)disp()returnendendends=0;for i=1:nl=1;for j=1:nif j=il=
4、l*(xi-x(j)/(x(i)-x(j);endends=s+l*y(i);endyi=s;end调用你的函数运行 1 题的结果: lagrange(x,y,xi)ans =2.2667附:分组名单星期二下午 5-6 节第 1 组组长:陈絮莹 缪 妃 何贵堂 刘 钰 马 倩第 2 组组长:李杰玉 黎筱惠 雷 霞 肖 娴 林碧珍 朱元正第 3 组组长:陈 静 苏小丽 李 郑 何淑楠 田冬秀 曾敬军第 4 组组长:杨 欣 王雪梅 徐莉萍 石小芳 雷 敏 唐 嘉第 5 组组长:杨佳悦 郭 滢 李媛媛 何 可 陈思露第 6 组组长:王钰琪 寇 玠 杨 丹 熊晨曦 周 丹第 7 组组长:姚 瑶 高倩倩
5、 金 杨 周海宁 杨 琴第 8 组组长:雷 芳 陈 艳 王玉娇 余 非 张 雪 王海燕星期三下午 5-6 节第 9 组组长: 刘超慧 王玉利 秦佳丽 张青梅 廖 婷 程思远第 10 组组长:杨 琴 冯康欣 黄宜纯 田晓东 郑美艳5第 11 组组长:黄 倩 肖雪梅 舒玉秀 杨 阳 黄 倩 宋亚超第 12 组组长:乔 欢 曹人月 万袁源 刘学勤 师小诚 沈金勇第 13 组组长:张全兴 程德超 冯 啸 魏 丹 李 茜 罗凤菊第 14 组组长:张 洋 何婷婷 刘云丹 彭英萍 马 静第 15 组组长:杨 丽 王书琪 袁 杰 宋慧玲 杨璐萍星期三下午 7-8 节第 16 组组长:李 欢 蒋书丽 康斯梦 王菊花 李芝琴第 17 组组长:杨 梅 郑雨来 李 维 刘玉兰 羊 玲第 18 组组长:左艳君 古 月 黄文凤 杨 娟 胡 洲 黄 川第 19 组组长:吴 星 谭 婷 张 欢 向 巧 钱 强 陈虹弟第 20 组组长:曾大超 胡 敏 马树述 罗玉婷 李琳玲第 21 组组长:石章波 拉吉石明 岳榆川 金小刚 张泽松第 22 组组长:贾孙鹏 袁 鹏 颜冬芹 陈 诚 张 博第 23 组组长:李自强 黄金辉 彭琦 岳 琪 李 宾 李 闯 第 24 组组长:王文媛 林小渝 刘 燕 严 英 何思敏 穆芦芸