- 博客(9)
- 收藏
- 关注
RSS订阅原创 python用内在open读写txt文件（流程和一些可能遇到的问题）
print()函数遇到这个 \n 就会换行，在加上一次打印结束后 print()本身也会换行，等于换行了两次，所以就会出现两行之间空一行。上面的代码会循环文件夹里面的所有文件，但有时候这个文件夹可能里面会有一些无关文件，比如我们只想打开‘.txt’文件，但文件夹有一些其他文件，这时候可以用下面的代码来操作。这个错误的原因是，我这里的txt文件里面有中文字符，而open默认的编码方式是‘gbk’，为了能正确读取中文，需要更改为“UTF-8”,之类的，我们也不想每个文件分开写代码来读取，这时候可以用。
2023-05-06 23:50:10
1953
1
原创 利用 Python 和 Selenium 自动下载知网期刊文件
利用Selenium自动下载知网期刊文件1 下载 Chrome 的 Webdriver2 代码2.1 导入相关的包2.2 设置 WebDriver2.3 找到期刊名在网页上对应元素2.4 下载期刊安装 Python 和Selenium 的过程就不讲了，下面直接上代码。目标是下载知网某个期刊某一期的全部文件（PDF格式）。这里以《教育研究》这个期刊为例，下载其中2021年第5期的所有文献。1 下载 Chrome 的 Webdriver想要用 Selenium 操控浏览器，首先需要下载浏览器对应的WebD
2021-11-19 21:42:57
2665
原创 Mathematica斐波那契线搜索代码
(*参数说明：f为函数，a为区间下界，b为区间上界，e为精度*)(*输出说明：输出为一个列表{a,b,n}，其中[a,b]代表最后求得的解区间，n代表迭代次数*)goldenSearch[f_, a_, b_, e_] := Module[{n = 2, d = b - a}, While[d > e, d = Fibonacci[n]/Fibonacci[n + 1]*d; u = a + d; l = b - d; If[f[u] < f[l], a = l,
2020-08-07 22:52:51
407
原创 L^2准则下的多项式逼近，一个比较有意思的问题吧
问题描述实际上：f(a)=∫(g−∑0naixi)2=∫g2+∫(∑0naixi)2−2∫g∑0naixif(a)=\int(g-\sum_0^na_ix^i)^2=\int g^2+\int (\sum_0^na_ix^i)^2-2\int g\sum_0^na_ix^if(a)=∫(g−∑0naixi)2=∫g2+∫(∑0naixi)2−2∫g∑0naixi求偏导数：∂f(a)∂ak=2∫(∑0naixi)xk−2∫g∗xk=0\frac{\partial f (a)}{\partia
2020-07-09 00:06:39
351
原创 Mathmatica多项式带余除法代码
几乎没有调用内置函数，除了求多项式最高次数时用了一下 Exponent[](*解析多项式*)(*将f=a0+a1*x+...+an*x^n解析成{{a0,0},{a1,1},...,{an,n}}的形式*)polyCoefficients[f_] := Module[{ rules1 = { c_*base_^power_ -> {c, power}, base_^power_ -> {1, power}, c_*x_ -> {c, 1},
2020-07-01 13:17:48
1613
原创 求解n阶方阵零化多项式的mathematica代码
Remove["`*"];zeroPolynomial[A_] := Module[{n = Length[A], base, m, s}, base = NestList[Dot[A, #] &, IdentityMatrix[n], n^2]; m = Transpose[Flatten /@ base]; s = NullSpace[m][[1]]; Total[Table[x^i*s[[i]], {i, Length[s]}]] ]
2020-06-30 22:10:53
430
原创 若儿当标准型
1. 引言这里略去若尔当块和若尔当矩阵的定义，主要是想谈谈自己对若尔当标准型的理解。我们都知道的简单结论是，对于nnn阶方阵AAA，如果AAA有nnn个互异的特征值λ1,...,λn\lambda_1,...,\lambda_nλ1,...,λn，那么AAA是可以对角化的，并且向量空间V=ker(A−λ1I)+...+ker(A−λnI)V=ker(A-\lambda_1 I)+...+ker(A-\lambda_n I)V=ker(A−λ1I)+...+ker(A−λnI)。实际上对角化的过
2020-06-28 00:12:45
1355
原创 ProjectEuler Mathematica 答案
50f[1] = Prime[1];f[n_] := f[n] = f[n - 1] + Prime[n];i = 1;While[f[i] < 1000000, i++];maxIter = i - 1;primeList = Table[Prime[i], {i, maxIter}]; i = 1; While[ Not[Or @@ PrimeQ[MovingAve...
2020-04-21 22:22:24
189
原创 杜利特尔分解(LU分解）的JAVA代码
杜利特尔分解(LU分解）的JAVA代码数学基础：设A为n阶方阵，满足其各阶主子式： det (Subscript[A, p]) != 0 A=LU其代码如下： public static double[][] LUDecomposition(double a[][],String X)//X代表返回矩阵的种类，若为U，返回u矩阵；若为l,返回l矩阵。若为其他，返回a...
2018-08-15 17:46:58
2274
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人