
计算方法
肥仔甜甜
学习日记,欢迎评论。
展开
-
矩阵三角分解(LU+高斯+平方根+求解方法+求逆)
待更新#include <iostream>#include <vector>#include <math.h>using namespace std;typedef vector<vector<double>> Matrix;#define MAX 6void dispMatrix(Matrix m, int n); ...原创 2019-05-30 23:18:40 · 1586 阅读 · 1 评论 -
复化Cotes+复化Simpson+复化梯形+变步长梯形+Romberg求积
待完善。。Cotes求积公式Simpson求积公式变步长梯形这里以Simpson求积流程作Sample:需要注意这里定义函数模板,Solve 第一个参数为函数名。typedef double (*Func)(int, double, double);vector<double> Solve(Func func, double a, double b);#inc...原创 2019-06-05 17:11:02 · 1687 阅读 · 0 评论 -
线性方程组求解(高斯、高斯若尔当)+应用(求逆,求行列式)
高斯消去法#include <iostream>#include <vector>#include <math.h>using namespace std;typedef vector<vector<double>> Matrix;#define MAX 6void dispMatrix(Matrix m, int n); ...原创 2019-05-28 21:00:55 · 1492 阅读 · 0 评论 -
迭代法(雅可比+高斯塞尔德+SOR)
#include <iostream>#include <vector>#include <math.h>using namespace std;#define eps 0.00001#define maxn 200typedef vector<vector<double>> Matrix;void dispRes(vect...原创 2019-05-31 18:10:00 · 1965 阅读 · 0 评论 -
预报校正欧拉求积(改进欧拉)
先用欧拉法求得一个初步的近似值,称为预报值,然后用它替代梯形法右端的yi+1再直接计算fi+1,得到校正值yi+1,这样就可以得到预报校正公式。{yp=yn+hf(xn,yn)yc=yn+hf(xn+1,yp)yn+1=12(yp+yc)\left\{\begin{array}{l}{y_{p}=y_{n}+h f\left(x_{n}, y_{n}\right)} \\ {y_{c}=y_{n...原创 2019-06-06 00:03:20 · 2195 阅读 · 0 评论 -
埃特金插值 aitken
#include <iostream>#include<iomanip>#include <vector>using namespace std;typedef vector<double> List;double aitken(List x, List y, double data);int main(){ int n; ...原创 2019-06-04 21:28:02 · 1262 阅读 · 0 评论