注册 登录  
 加关注
   显示下一条  |  关闭
温馨提示!由于新浪微博认证机制调整,您的新浪微博帐号绑定已过期,请重新绑定!立即重新绑定新浪微博》  |  关闭

Mr.Right

不顾一切的去想,于是我们有了梦想。脚踏实地的去做,于是梦想成了现实。

 
 
 

日志

 
 
关于我

人生一年又一年,只要每年都有所积累,有所成长,都有那么一次自己认为满意的花开时刻就好。即使一时不顺,也要敞开胸怀。生命的荣枯并不是简单的重复,一时的得失不是成败的尺度。花开不是荣耀,而是一个美丽的结束,花谢也不是耻辱,而是一个低调的开始。

网易考拉推荐

Neville插值MATLAB源码  

2016-12-04 23:22:31|  分类: 编程 |  标签: |举报 |字号 订阅

  下载LOFTER 我的照片书  |
Neville插值是一种线性逐次插值,在精密轨道的插值中应用。
Neville插值MATLAB源码 - 阿英 - Mr.Right
 
xn=[-1 -0.5 0 0.5]
f=[0.86199480 0.95802009 1.0986123 1.2943767]
x0 = 0.25;
[y, Table] = neville(xn,f,x0)
[y0, ylist ]= nev(x0, xn, f)

------------------------------------------------------------------

function [y, Table]=neville(xn, yn, x)
% Example:
% xn=[-1 -0.5 0 0.5]
% f=[0.86199480 0.95802009 1.0986123 1.2943767]
% [y, Table] = neville(xn,f,0.25)
n=length(xn)-1; % n = order of interpolation (n+1 = # of points)
Table(:, 1) = xn;
Table(:, 2) = yn;
for i=2:n+1;
    for j=3:i+1;
        Table(i,j)=((x-xn(i))*Table(i-1,j-1)-...
            (x-xn(i-j+2))*Table(i,j-1))/(xn(i-j+2)-xn(i));
    end;
end;
y = Table(n+1,n+2);





--------------------------------------------------------------------

function [y0, y] = nev(x0, x, y)
% Neville's algorithm as a function (save as "nev.m")
%
% inputs:
%    x(1),...,x(n+1)    x coords
%    y(1),...,y(n+1)    y coords
%    xx=evaluation point for interpolating polynomial p
%
% output:  p(x0)
%  Example:
% xn=[-1 -0.5 0 0.5]
% f=[0.86199480 0.95802009 1.0986123 1.2943767]
% x0 = 0.25;
% [y0, ylist ]= nev(x0, xn, f)
n=length(x)-1; % n = order of interpolation (n+1 = # of points)
for i = n:-1:1
   for j = 1:i
      y(j) = (x0-x(j))*y(j+1) - (x0-x(j+n+1-i))*y(j);
      y(j) = y(j)/(x(j+n+1-i)-x(j));
   end
end
y0 = y(1);




-----------------------------------------
https://en.wikipedia.org/wiki/Neville's_algorithm
  评论这张
 
阅读(154)| 评论(0)
推荐 转载

历史上的今天

在LOFTER的更多文章

评论

<#--最新日志,群博日志--> <#--推荐日志--> <#--引用记录--> <#--博主推荐--> <#--随机阅读--> <#--首页推荐--> <#--历史上的今天--> <#--被推荐日志--> <#--上一篇,下一篇--> <#-- 热度 --> <#-- 网易新闻广告 --> <#--右边模块结构--> <#--评论模块结构--> <#--引用模块结构--> <#--博主发起的投票-->
 
 
 
 
 
 
 
 
 
 
 
 
 
 

页脚

网易公司版权所有 ©1997-2017