前两天对学习了R里面计算的基本范围，以及一些求解方程的方法，今天来看看积分，其实上个学期学了数值分析，对这部分的算法是有所了解的，当时是用matlab写了一遍，已经忘了怎么实现的了，现在用R重新写一遍吧，算法有梯形积分法，辛普森积分法，自适应积分法。 梯形积分法 梯形积分法可以用下图很好的解释 就是 ...
分类：
编程语言 时间：
2018-06-07 21:49:26
阅读次数：
1248
有时候我们需要进行一些复杂的数学计算，比如求导， 求积分，解方程，还是用abcd字母代表变量的方程等，这就需要进行复杂的数学运算还需要具备良好的数学基础。不过现在有一个非常方便的在线工具，只需要几秒钟， 就能告诉我们所有的答案。 sage &ems;sage是一个 ...
分类：
其他好文 时间：
2018-06-04 14:28:16
阅读次数：
581
高斯消元法 可以用于求解线性方程组，即n元1次方程组。利用矩阵，大致思路与普通解方程方法类似。只是更具一般性。将系数与右侧的常数存成一个矩阵，然后每次用第i行消去下面每行的第i个系数，最后就会得到一个一元方程，然后从后到前依次代回即可。 然后就是精度的问题，因为计算机中没有分数，所以只能用doubl ...
分类：
其他好文 时间：
2018-05-19 15:40:12
阅读次数：
189
https://www.cnblogs.com/xinbaby829/p/7289431.html 我们每个人都会在我们的生活或者工作中遇到各种各样的最优化问题，比如每个企业和个人都要考虑的一个问题“在一定成本下，如何使利润最大化”等。最优化方法是一种数学方法，它是研究在给定约束之下如何寻求某些因素 ...
分类：
编程语言 时间：
2018-05-18 20:07:11
阅读次数：
233
题目：https://www.lydsy.com/JudgeOnline/problem.php?id=2337 异或就一位一位考虑； x为到n的概率，解方程组即可； 考虑了n就各种蜜汁错误，所以索性不管n了，这样的题好像不管n比较方便。 代码如下： ...
分类：
其他好文 时间：
2018-05-17 19:57:45
阅读次数：
187
逐次下降法的定义： 下面给出Python实现 # -*- coding: utf-8 -*-import numpy as np from numpy import *from common_libs import *import matplotlib.pyplot as plt #消元发求解方程组 ...
分类：
编程语言 时间：
2018-05-13 19:52:08
阅读次数：
295
如何解方程a*x≡b(mod m)呢？因为a*x-b|m, 故令a*x-b=-y*m，即a*x+m*y=b。根据Bezout定理，该方程有解当且仅当gcd(a,m)|b。我们把等式两边同乘以gcd(a,m)/b，得到a*x0+m*y0=gcd(a, m)。这个方程可以用扩展欧几里得算法求得得到x0。 ...
分类：
其他好文 时间：
2018-04-22 17:21:10
阅读次数：
176
SymPy是一个符号计算的Python库。它的目标是成为一个全功能的计算机代数系统，同时保持代码简 洁、易于理解和扩展。它完全由Python写成，不依赖于外部库。SymPy支持符号计算、高精度计算、模式匹配、绘图、解方程、微积分、组合数学、离散 数学、几何学、概率与统计、物理学等方面的功能 正如你看 ...
分类：
编程语言 时间：
2018-04-18 00:55:05
阅读次数：
522
二分法是一种求解方程近似根的方法。对于一个函数 f(x)f(x)，使用二分法求 f(x)f(x) 近似解的时候，我们先设定一个迭代区间（在这个题目上，我们之后给出了的两个初值决定的区间 [-20,20][?20,20]），区间两端自变量 xx 的值对应的 f(x)f(x) 值是异号的，之后我们会计算 ...
分类：
编程语言 时间：
2018-04-12 19:56:15
阅读次数：
511
代码： 1)扩展欧几里得算法求ax+by=c时 2）求解模线性方程 ax = b(mod n) 同余方程 ax≡b (mod n)对于未知数 x 有解，当且仅当 gcd(a,n) | b。且方程有解时，方程有 gcd(a,n) 个解。 求解方程 ax≡b (mod n) 相当于求解方程 ax+ ny ...
分类：
编程语言 时间：
2018-04-11 00:08:55
阅读次数：
243
