有时候我们要数值求解超越方程的多个根，但是数值方法都要给定一个初值。 matlab有内建函数fsolve求解非线性方程（组），但是只能求一组给定初值的解，如果要求多个根（如频率方程），可以先用mathematica画图，然后观察得到零点个数和大概位置，然后调用fsolve求解。 这里提供了一个 ...

详细实验指导见上一篇，此处只写内容啦 实验内容： 求解如下4元线性方程组的近似解。 • Jacobi迭代过程 运行结果： •Gauss-Seidel迭代 运行 ...

Doolittle分解复习 LDU分解 定理：矩阵A的各阶顺序主子式≠0，A可唯一分解为：A=LDU 其中，L单位下三角矩阵，U单位上三角矩阵，D非奇 ...

1.追赶法--先对矩阵A作Crount分解A=LU 例：用追赶法求解Ax=b 2.向量的范数（绝对者的和，绝对值的最大值，平方和开方） 1.绝对值的和： 2.绝对值的最大值： 3.平方和开方： 3.矩阵的范数（最大列 ...

一、由于任何偏微分方程初边值问题通过半离散之后都会转化为两点边值问题，以下仅以两点边值问题开始简要讨论Galerkin-Legendre 谱方法的数值实现过程 考虑如下分数阶两点边值问题 \begin{align}&-(-\Delta)^{\frac{\alpha}{2}}u(x)=f ...

　　在乘法逆元里我们对于仅满足b，m互质的情况，我们需要求解的是一个同余方程：b*x≡1（mod m），那么接下来我们就讨论一下类似的线性同余方程的求解。 线性同余方程： 　　给定整数a，b，m，求一个整数满足：a*x≡b（mod m），或给出无解。 　　因为未知数的次数为1，所以我 ...

简介 van der Pol 方程 code ...

1）求解线性不定方程 　　ax + by = c 　　先求出一组解， 然后考虑如何表示通解， 设d = gcd(a, b), 假设c不是d的倍数， 则左边是d的倍数而右边不是， 则方程无解， 所以方程有解当且仅当d | c. 　　设c = c' * d, 我们先考虑方程 ax ...