![](https://img-blog.csdnimg.cn/20201014180756927.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
计算方法
油醋三椒
StuTan
展开
-
Euler 格式和改进Euler 格式求解初值问题及代码
介绍关于初值问题的解实例代码#include<iostream>using namespace std;double h;double x[100],y[100];double f(double x0, double y0){ return -2*x0/y0 + y0; }void Eluer(){ for(int i=1;i<10;i++...原创 2019-12-24 15:16:22 · 5519 阅读 · 1 评论 -
杜利特尔 (Doolittle)矩阵分解法求线性方程组的解
简介 若方阵 A 可以分解为一个下三角矩阵 L 和一个上三角矩阵 U的乘积,即 A = LU,则这种分解称为 A 的一种三角分解或 LU分解。如果 L 为单位下三角矩阵,则称为杜利特尔 (Doolittle)。 以四阶矩阵为例,可分解为以下形式:实例代码#include<iostream>using namespace std;int n;double a[1...原创 2019-12-24 15:00:19 · 37970 阅读 · 0 评论 -
高斯消元法直接求解线性方程组
简介实例代码#include<iostream>using namespace std;int n;double a[100][100],b[100],x[100];void input(){ for(int i=0;i<n;i++) { for(int j=0;j<n;j++) cin>>a[i][j]; cin>...原创 2019-12-24 14:51:57 · 454 阅读 · 0 评论 -
Gauss-Seidel 算法求线性方程组的解及其代码
简介则高斯赛德尔迭代式为:实例代码#include<iostream>#include<math.h>using namespace std;int n;double a[100][100],b[100],x[100][100];double e;void input(){ for(int i=0;i<n;i++) { for...原创 2019-12-24 14:40:21 · 1653 阅读 · 0 评论 -
雅克比 (Jacobi) 迭代法求线性方程组的解及其代码
简介实例可以看到,当迭代次数 k 增大时,迭代值会越来越接近解 x1 = x2 = x3 = 1。代码#include<iostream>#include<math.h>using namespace std;int n;double a[100][100],b[100],x[100][100];double e;void input(){...原创 2019-12-24 14:33:59 · 12203 阅读 · 2 评论 -
Newton 牛顿迭代法求根公式及其改进算法
简介原创 2019-12-24 14:22:56 · 3950 阅读 · 0 评论 -
Gauss 求积公式及代码
简介Newton-Cotes 公式在构造时,限定积分区间的等分点作为求积节点,这样做在简化处理的同时也限制了代数精度。如果求积节点也可自由选择,即机械求积公式: ∫abf(x)dx≈(b−a)∑i=0nλif(xi)\ \int_a^b{f\left( x \right) dx\approx \left( b-a \right) \sum_{i=0}^n{\lambda _if\...原创 2019-12-24 14:07:14 · 6504 阅读 · 0 评论 -
Newton-Cotes 公式及代码
简介Newton-Cotes 公式是指在等距节点下使用 Lagrange 插值多项式建立的数值求积公式。设将求积区间 [a, b] 划分为 n 等分,选取等分点xi=a+ih,h=b−an,i=0,1,2,⋅⋅⋅,nx_i = a + ih , h=\frac{b-a}{n}, i = 0, 1, 2, · · · , nxi=a+ih,h=nb−a,i=0,1,2,⋅⋅⋅,n作为...原创 2019-12-24 13:54:45 · 9297 阅读 · 0 评论 -
埃尔米特插值及其代码
简介Lagrange 插值、Newton 插值虽然构造比较简单,但仅要求他们的函数值与被插函数的函数值在节点上相同,但没有考虑到导数值。Hermite 插值 (切触插值) 是要构造一个插值函数,不但在给定的节点上取已知函数值,而且取已知导数值,使插值函数与被插函数的密合程度更好。关于两点三次 Hermite 插值即,Hermite插值通式为:实例代码#include&...原创 2019-12-24 13:36:32 · 4285 阅读 · 2 评论 -
牛顿插值法及其代码
简介Lagrange 插值方法虽然易算,但若要增加一个节点时,全部基函数 φ(x) 都需重新算过。本节介绍另外一种插值方法-牛顿插值法。牛顿插值法是一种重要的插值方法,与拉格朗日插值法在同阶时产生的多项式在化简以后是一样的,余项也是一样的,因此两者只是写法形式不同而已。在针对变化的数据时,其各自独特的写法形式给其适合的应用场景带来了区别。牛顿插值法适用于被插的点不断增多,而且可以很好的复...原创 2019-12-24 10:34:30 · 9356 阅读 · 0 评论 -
拉格朗日插值法及代码
简介在数值分析中,拉格朗日插值法是以法国十八世纪数学家约瑟夫·拉格朗日命名的一种多项式插值方法。许多实际问题中都用函数来表示某种内在联系或规律,而不少函数都只能通过实验和观测来了解。如对实践中的某个物理量进行观测,在若干个不同的地方得到相应的观测值,拉格朗日插值法可以找到一个多项式,其恰好在各个观测的点取到观测到的值。这样的多项式称为拉格朗日 (插值) 多项式。数学上来说,拉格朗日插值...原创 2019-12-24 09:52:35 · 10984 阅读 · 2 评论 -
切比雪夫插值
简介使用平均分布的点作为插值多项式的基点Xi很普遍。在很多情况下,用于插值的数据点仅以这种形式存在,例如当数据由相同时间间隔分布的一起读取的数据所组成时。在其他情况下,我们可以在认为合适的地方自由选择基点。事实证明,基点间距选取的方式对于插值误差有很大的影响,切比雪夫插值是一种特定最优的点间距选取方式。理论分析切比雪夫理论切比雪夫插值的动机是在插值区间上,提高对如下插值误差(如:牛顿差商公...原创 2019-12-23 21:46:24 · 17424 阅读 · 3 评论 -
有关向量的叉乘
1.当两个向量均为二维向量时:a→ = (x1 , y1)\overrightarrow{a}\ =\ \left( x_1\ ,\ y_1 \right) a = (x1 , y1)b→ = (x2 , y2)\overrightarrow{b}\ =\ \left( x...原创 2019-12-10 21:00:48 · 1089 阅读 · 0 评论