我们把 0 阶差商到 n 阶差商都放在一个一维数组里面就会发现其中的规律,也可以根据这些规律来指导我们编程。
0 阶差商 n 个,一阶差商 n-1个 .... 直到 n 阶差商 1个。
求差商代码为:
bool CShuzhijisuanDlg::Diff()
{
int cur;
cur=m_length*(m_length+1)/2+1;
m_diff=new double[cur];
for(int i=0;i<m_length;i++){//先存放0阶差商
m_diff[i]=parmy[i];
m_diff[i];
}
//存放二阶以上差商
cur=m_length;
int temp;
for(int internal=1;internal<m_length;internal++){
for(i=internal;i<m_length;i++){
temp=cur-m_length+internal-1;
m_diff[cur++]=(m_diff[temp]-m_diff[temp+1])/(parmx[i-internal]-parmx[i]);
}
}
m_isdiff=true;
return true;
}
差商求完我们开始求牛顿插值了:
double CShuzhijisuanDlg::Newton(double xx)
{
if(m_isdiff==false)
Diff();
short int cur=m_length;
double y,z=parmy[0];
for(int i=1;i<m_length;i++){
y=1;
for(int j=0;j<i;j++){
y=y*(xx-parmx[j]);
}
z+=m_diff[cur]*y;
cur=cur+m_length-i;
}
return z;
}
责任编辑:小草