- 博客(12)
- 收藏
- 关注
原创 快速弦截法解方程
/*****快速弦截法*** 公式: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 4291
原创 牛顿迭代法求方程根
/*****牛顿迭代法求方程根*** 牛顿公式: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 5392 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 8467 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 1691
原创 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 2376
原创 变步长梯形求解微分方程
/*****变步长梯形求解微分方程*** 获取初值: 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 2764 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 2473
原创 定步长复化梯形求解积分方程
/*****定步长复化梯形求解积分方程*** n-1 n-1 定步长复化梯形公式: Tn= ∑ h/2〔 f(a) + f(b) + 2 ∑ f(x(k)) 〕
2005-06-14 16:17:00 1977
原创 高斯消元法解线性方程组(选主元)
/*****高斯消元法解线性方程组*** 《数值分析简明教程》-2 Editon -高等教育出版社- page 177 消元回代算法流程图 - page 180 选主元素算法流程图 描述:
2005-06-14 16:04:00 3317
原创 拉格朗日(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 4963
原创 分段线性插值
/*****分段线性插值*** 计算公式: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 3196
原创 方程求根二分法
/*****方程求根二分法*** 在实际计算允许的误差范围(ε)内 对所求根区间(a,b)不断缩小直得到所期望精度的数据 属性:数值逼近法 误差(精度):| x-x|≤( b-a)/2 = ( b -a )/pow(2,k+1)=ε. x为精确值,x为第k次二分后数据
2005-06-12 16:05:00 1708
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人