已阅读5页,还剩2页未读, 继续免费阅读
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
Lagrange插值MATLAB源程序function y=lagrange(x0,y0,x); % x0自变量取值向量已知 y0为已知对应x0的函数取值,x为要求插值点坐标n=length(x0); m=length(x); for i=1:m z=x(i); s=0.0; for k=1:n p=1.0; for j=1:n if j=k p=p*(z-x0(j)/(x0(k)-x0(j); %插值基函数 end end s=p*y0(k)+s; %lagrange插值多项式 end y(i)=s;end 测试:x0=0:2 y0=2 3 5 x=0.5Lagrange(x0,y0,x)x0=0:2 y0=2 3 5 x=0:0.01:2Lagrange(x0,y0,x) Newton插值MATLAB源程序function f=Newton(x,y,x0) syms t; n = length(x); c(1:n) = 0.0; f = y(1); y1 = 0; l = 1; for(i=1:n-1) for(j=i+1:n) y1(j) = (y(j)-y(i)/(x(j)-x(i); end c(i) = y1(i+1); l = l*(t-x(i); f = f + c(i)*l; simplify(f); y = y1; if(i=n-1) if(nargin = 3) f = subs(f,t,x0); else f = collect(f); end end endtest: x=1 -1 2 y=0 -3 4 x0=-1:0.1:2Newton(x,y,x0)分段插值MATLAB源程序function y = div_linear(x0,y0,x,n) for i = 1:n-1 if (x = x0(i) & (x =j t(i,j)=(t(i,j-1)-t(i-1,j-1)/(x(i)-x(i-j+1); end end end for k=1:n s=1; m=1; for j=1:k if jnewtonNewton插值MATLAB源程序function yi=newton(x,y,xi) m=length(x);n=length(y); if m=n error(x and y must same);end f=zeros(n+1,1); k=2; f(1)=y(1) while k=n+1 f(1)=y(k);k,x(k) for i=1:k-1 if i=k-1 f(i+1)=(f(i)-y(i)/(x(k)-x(i); end end cs(i)=f(i+1); y(k)=f(k); k=k+1;end cfwh=0; for i=1:n-2 w=1; for j=1:i w=w*(xi-x(j); end cfwh=cfwh+cs(i)*w; end yi=y(1)+cfwh;x=0:2 y=2 3 5 xi=0.5newton(x,y,xi)Newton2)functionf=Newton(x,y,x0)symst;if(length(x)=length(y) n=length(x);c(1:n)=0.0;else disp(x和y的维数不相等!);return;endf=y(1);y1=0;l=1;for(i=1:n-1)for(j=i+1:n) y1(j)=(y(j)-y(i)/(x(j)-x(i);end c(i)=y1(i+1);l=l*(t-x(i);f=f+c(i)*l;simplify(f);y=y1;if(i=n-1) if(nargin=3) f=subs(f,t,x0);else f=collect(f);%将插值多项式展开f=vpa(f,6);endend endx=0:2 y=2 3 5 x0=0.5Newton(x,y,x0) function=Newton(x,y,v)x=input(X数组=:);y=input(Y数组=:);v=input(插值点数值=:);n=length(x);t=zeros(n,n);u=0;fori=1:nt(i,1)=y(i);endforj=2:nfori=2:nifi=jt(i,j)=(t(i,j-1)-t(i-1,j-1)/(x(i)-x(i-j+1);endendendfork=1:ns=1;m=1;forj=1:kifjks=s*(v-x(j);endendm=s*t(k,k);u=u+m;enddisp
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 2024年企业劳动合同电子版(4篇)
- 2024年法律服务合同格式范本(四篇)
- 2024年个人汽车抵押借款合同样本(2篇)
- 2024年上海市农村房屋买卖合同(二篇)
- 2024年工厂简单转让合同(二篇)
- 机械制图与技术测量(第二版)课题三投影基本知识
- 易顺煤矿安全生产责任制
- 购车合同简易版4篇
- 小学武术教育的发展前景分析研究 教育教学专业
- 肠系膜结核的护理查房
- 转轮除湿设计方案
- 蝴蝶兰组培苗的工厂化生产成本分析
- 危险化学品储罐全面检验报告表格(共20页)
- kg舞台升降液压系统设计说明书
- 《英语语法》课程标准
- 音乐八年级下册欣赏课《山丹丹开花红艳艳》
- 桩板墙施工技术交底
- 汽车荷载冲击系数计算(自编)
- 张公中学中考考务工作方案
- V303新会计准则年审 会计报表审计调整及分析系统1
- 供应商PPAP审核检查表--示例
评论
0/150
提交评论