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

信息 灵感 创新

I? =Information,Inspiration,Innovation

 
 
 

日志

 
 
关于我

we are 5. Mathematics, Computation, Programming, Engineering, and Making fun of life.

网易考拉推荐

计算差商表函数  

2010-12-01 16:32:44|  分类: M&M |  标签: |举报 |字号 订阅

  下载LOFTER 我的照片书  |

用于插值函数求解中的差商表

%计算差商表函数
%输入参数:x:节点列向量,f:节点对于函数值的列向量
%本函数不做异常检查,请确保x和f均为n*1向量
%返回值:A:显示差商表,其主对角线上值即为f[x0,x1,...xi]
%注意A的第一列为0阶差商f(x)
%2010-12-1
function A=didiff(x,f)
    n=max(size(x));
    B=zeros(n,n+1);%B是包含了节点的差商表
    A=zeros(n,n);
    B(:,1)=x;
    B(:,2)=f;
    for i=3:n+1%列
        for j=i-1:n%行
            B(j,i)=(B(j,i-1)-B(j-1,i-1))/(B(j,1)-B(j-i+2,1));
        end
    end
    A=B(:,2:n+1);%获取差商
end

  评论这张
 
阅读(1920)| 评论(0)
推荐 转载

历史上的今天

评论

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

页脚

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