Nylon_Min的博客

今天无意间发现了nympy这个库于是就想着能够使用它来实现以下牛顿迭代法：但是发现利用nympy如果写到函数里边,总是出现cannot unpack non-iterable Float object的错误调试无果后 一气之下 直接都不要函数了 直接写（生气.jpg牛顿迭代法流程图如下：代码如下：在这里插入代码片...

计算机数值分析：线性方程组迭代公式：高斯赛德尔迭代公式课本上高斯赛德尔流程图如下：这里直接就是使用了之前的雅可比迭代的文件：文件为：具体代码如下：#coding=gbkimport numpy as np;def GaoSiSaider(n,minner,datas,x): count=1; numVarity=datas.shape[0]; cols=...

计算机数值分析线性方程组的迭代法：雅可比迭代公式：书上的流程图如下：这里直接把输入输出弄到了文件之中具体代码如下：里面使用了python的numpy库#coding=gbkimport numpy as np;def Jacobi(n,minner,datas,x): count=1; numVarity=datas.shape[0]; cols=da...

拉格朗日插值公式流程图如下：python代码实现如下：# coding=gbk;def Data_in(list0,list1,list2): list0=input("请输入离散点的值：").split(" "); for i in range(0,len(list0)): if i%2==0: #如果是奇数 那么指的也就是x ...

曲线拟合最小二乘法其中使用了机器学习之中得正规化方程直接计算得到theta# coding=gbk;# 要输入的数据data=165 187 123 126 150 172 123 125 141 148import numpy as np;import matplotlib.pyplot as plt;def Data_in(list0,list1,list2): list0...

埃特金插值法这里要说明，我没有按照书中所给流程图来写。不过为了完整，还是将流程图贴上来实际上，我是按照下面这个表，相当于是递推获得代码如下#0.3 0.29850 0.4 0.39646 0.5 0.49311 0.6 0.58813 0.7 0.68122 #预测 0.462def Data_in(list1,list2): list0=input("请输入离散点的值...

这里并不能够说完全实现了复化辛普生公式，因为这里面涉及到具体的函数，而我们需要事先知道函数表达式，才能够求出来。故而这里是以书上的一道例题来写例题以及书上结果如下：这里我也并没有按照流程图来写（不过大同小异），而是按照表达式来写流程图和表达式如下代码如下：# coding=gbk;#因为使用复化辛普生公式会涉及到函数的具体形式# 所以这里我就暂且令 f(x)=sin(x)/x...

改进的欧拉公式流程图以及方程式如下（以下代码是按照方程式写就）由于使用改进的欧拉公式需要已知函数表达式的具体形式，故而这里用了书上的具体的习题习题如下：代码如下：# coding=gbk;#改进的欧拉方法需要函数的具体表达式 #这里就暂且用书上的题目所给公式了def compute_fx(x,y): return y-2*x/y;#进行数据初始化#需要输入初值...

计算机数值分析：牛顿插值法这里使用的是牛顿插值法的差商公式；说明： 牛顿插值法可以说是埃特金插值法的升级版，所以在很大的程度上，牛顿插值法与埃特金插值法有类似的地方，但是也有不相同的地方。而这个不相同的地方正是我们编程的时候应该要注意的。埃特金插值法代码实现：参考 埃特金Python实现首先来说明一下牛顿差商公式的定义：这里公式的说明参考：牛顿差商公式说明牛顿差商公式—...

计算机数值分析：线性方程组的直接解法：高斯消去法书本上的流程图如下：这里也是使用的文件录入数据文件如下：代码如下：#coding=gbk;import numpy as np;def GaussEliminate(datas,x): rows=datas.shape[0]; #这个行数其实同时也是变量的个数 cols=datas.shape[1]; ...