$excrt$ $\text{是用来解决当模数不互质的情况下的普通}crt\text{情况的}$ 具体见 "P4777" 事实上和$crt$没有什么关系 主要思想是:我们不断的合并两个同余方程,最后合并到只剩一个。 对于方程组 $\begin{cases} x \equiv b_1\ ({\rm m ...
分类:
其他好文 时间:
2018-12-31 19:15:42
阅读次数:
223
我开始以为矩阵是为了把线性方程组的系数抽取出来,方便方程组化简和求解,后来发现矩阵的用处不止如此,不然就不会写一本书了。 矩阵可以方便的用来表示线性空间,一个简单的二维数阵,就可以表示成n维线性空间。 一个毫无意义的有序数阵,我们赋予它意义,他就可以表示成一个空间。那为什么要这么做呢?这是因为矩阵的 ...
分类:
其他好文 时间:
2018-12-30 16:21:14
阅读次数:
223
在物理学和数学中,在动力系统领域,双摆是一个摆锤,另一个摆锤连接在其末端,是一个简单的物理系统,具有丰富的动态特性,对初始条件具有很强的敏感性。双摆的运动由一组耦合的常微分方程控制并且是混沌的。 由于双摆运动涉及到二阶微分方程组,在matlab中对双摆的仿真需要用到ode45求其数值解。 (代码来自 ...
分类:
其他好文 时间:
2018-12-29 23:06:26
阅读次数:
587
#include #include #define ll long long using namespace std; ll k; ll a1,r1; ll a2,r2; ll x,y; ll g; void init() { scanf("%lld",&k); } void exgcd(ll a,... ...
分类:
其他好文 时间:
2018-12-14 22:53:35
阅读次数:
199
? 书中第六章部分程序,加上自己补充的代码,包括高斯消元法求解线性方程组,高斯 - 约旦消元法求解线性方程组 ● 高斯消元法求解线性方程组,将原方程转化为上三角矩阵,然后从最后一个方程开始求解 ● 高斯 - 约旦消元法求解线性方程组,将原方程转化为约旦标准型,无解时产生一个 yA == 0,yb ! ...
分类:
编程语言 时间:
2018-12-13 20:23:13
阅读次数:
131
利用首次积分法(First Integral)求解对称形式的常微分方程组:\[\frac{{\rm\,d}x}{-x+y+z}=\frac{{\rm\,d}y}{x-y+z}=\frac{{\rm\,d}z}{x+y-z}\] \[\frac{{\rm\,d}x}{-x^2+y^2+z^2}=\fr ...
分类:
其他好文 时间:
2018-12-03 10:26:11
阅读次数:
339
资源链接:https://pan.baidu.com/s/10MWPlZS60D08HLADR_O5Ow高清中文PDF版,290页,为完整版。本书使用向量的概念对国内高校工科“线性代数”的课程内容进行了较全面的几何分析。从向量的几何意义开始,分别讲述了向量组、向量空间、行列式、矩阵、线性方程组和二次型的几何意义或几何解释,其中不乏重要概念的物理意义的解释。这本书就像一串项梁,把上百个概念和定理的几
分类:
其他好文 时间:
2018-11-27 22:52:46
阅读次数:
3494
1. 最小二乘 $Ax=b$ 经常会没有解,当方程个数大于未知数个数,也即 $m n$ 时,列空间并不是 $R^m$ 空间的全部,因此 $b$ 可能不在列空间中,这时候方程组就无解,但我们不应该就此而停止。 也就是误差 $e = b Ax$ 并不总是能得到 0,这时候,如果误差 $e$ 的长度尽可能 ...
分类:
其他好文 时间:
2018-11-26 00:21:09
阅读次数:
269
1. 分别用Gauss消去法、列主元Gauss消去法、三角分解方法求解方程组 程序: (1)Guess消去法: function x=GaussXQByOrder(A,b) %Gauss消去法 N = size(A); n = N(1); x = zeros(n,1); for i=1:(n-1) ...
分类:
其他好文 时间:
2018-11-20 21:12:40
阅读次数:
227
题目描述: 有一棵n 个点的以 1 为根的树, 以及 n 个整数变量xi。树上 i 的父亲是 fi, 每条边(i,fi)有一个权值wi,表示一个方程 xi + xfi = wi,这 n-1个方程构成了一个方程组。 现在给出q 个操作,有两种类型: 1 u v s,表示询问加上 xu + xv = s ...
分类:
其他好文 时间:
2018-11-07 22:13:34
阅读次数:
279