- 博客(12)
- 收藏
- 关注
RSS订阅原创 快速弦截法解方程
/*****快速弦截法*** 公式：x(k+1) = x(k) - ( x(k) - x(k-1) ) * f( x(k) ) / f( x(k) - f(x(k-1) ) 属性：两步迭代法 描述：计算x(k+1)时需要利用前两步信息x(k),x(k-1).免去了Newton法中需要求解一阶导函数的繁琐 《数值分析简明教程》-2 Editon -高等教育出版社 -page 140 -无
2005-06-14 21:19:00
4279
原创 牛顿迭代法求方程根
/*****牛顿迭代法求方程根*** 牛顿公式：x(k+1) = x(k) - f(x(k)) / f(x(k)) 迭代函数：Ф(x) = x - f(x) / f(x) 属性：方程求根迭代法 此时的迭代函数必须保证X(k)有极限，即迭代收敛。 《数值分析简明教程》-2 Editon -高等教育出版社 -page 136 -算法流程图 代码维护：2005.6.14 Drag
2005-06-14 21:08:00
5385
1
原创 四阶Runge-Kutta法解常微分方程
/*****四阶Runge-Kutta法*** 经典格式： y(n+1) = y(n) + h/6 ( K1 + 2*K2 + 2*K3 + K4 ) K1 = f( x(n) , y(n) ) K2 = f( x(n+1/2) , y(n) + h/2*K1 ) K3 = f( x(n+1/2) , y(n) + h/2*K2 ) K4 = f( x(n+1) , y(n) + h*K
2005-06-14 21:00:00
8451
1
原创 改进Eular方法解常微分方程
/*****改进Eular方法*** 预报：*Y(n+1)=Y(n) + h * f( x(n) , y(n) ); 改进：Y(n+1)=Y(n) + h/2 *［ f( x(n) , y(n) ) + f( x(n+1) , *Y(n+1) ) ］ h=x(n+1)-x(n) 平均化形式： Yp= Yn + hf( X(n),Y(n) ) Yc= Yn + hf( X(n+1),
2005-06-14 20:50:00
1685
原创 Eular方法解常微分方程
/*****Eular方法*** Eular公式：Y(n+1)=Y(n) + h * f( x(n) , y(n) ) h=x(n+1)-x(n) 精度(误差)：I= Y(x(n+1)) - Y(n+1) = h*h/2 Y(ε) 局部截断误差 由泰勒公式可知Eular方法仅为一阶方法，是很粗糙的计算方法 属性：差分方法
2005-06-14 20:45:00
2369
原创 变步长梯形求解微分方程
/*****变步长梯形求解微分方程*** 获取初值： T1 = h/2［ f(x(k)) + f(x(k+1)) ］ n-1 变步长梯形公式： T2n = 1/2*Tn + h/2 * ∑ f ( x(k+1/2) ) k=
2005-06-14 16:36:00
2751
1
原创 定步长复化Simpson解积分方程
/*****定步长复化Simpson解积分方程*** n-1 定步长复化Simpson公式：Sn=h/6｛f(a)-f(b)+ ∑ ［ 4f(x(k-1/2)) + 2f(x(k)) ］｝
2005-06-14 16:26:00
2465
原创 定步长复化梯形求解积分方程
/*****定步长复化梯形求解积分方程*** n-1 n-1 定步长复化梯形公式： Tn= ∑ h/2〔 f(a) + f(b) + 2 ∑ f(x(k)) 〕
2005-06-14 16:17:00
1972
原创 高斯消元法解线性方程组（选主元）
/*****高斯消元法解线性方程组*** 《数值分析简明教程》-2 Editon -高等教育出版社- page 177 消元回代算法流程图 - page 180 选主元素算法流程图 描述：
2005-06-14 16:04:00
3308
原创 拉格朗日（lagrange）插值计算
/**拉格朗日插值计算拉格朗日插值公式： n n Pn(x(i))= ∑〔 ∏ （x-x(j)）/（x(k)-x(j)） 〕y(k) k=0 j=0 j≠k 属性：插值计算法
2005-06-12 16:27:00
4937
原创 分段线性插值
/*****分段线性插值*** 计算公式：S1(x)=Φ0（x-x(i)/h(i))）y(i)+Φ1（x-x(i)/h(i)）y(i+1) 步长：h(i)=x(i+1)-x(i) 领域：x(i)≤x≤x(i+1) 属性：插值计算法 精度：| f(x) - S1(x) | = h(i)*h(i)/8 * max | f(x) |
2005-06-12 16:16:00
3184
原创 方程求根二分法
/*****方程求根二分法*** 在实际计算允许的误差范围(ε)内 对所求根区间(a,b)不断缩小直得到所期望精度的数据 属性：数值逼近法 误差（精度）：| x-x|≤( b-a)/2 = ( b -a )/pow(2,k+1)=ε. x为精确值，x为第k次二分后数据
2005-06-12 16:05:00
1704
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人