数值分析
小盼你最萌哒
这个作者很懒,什么都没留下…
展开
-
Jacobi和Gauss-Seidel迭代法求解方程组
迭代法简介迭代法也称辗转法,是一种不断用变量的旧值递推新值的过程,跟迭代法相对应的是直接法(或者称为一次解法),即一次性解决问题。迭代算法是用计算机解决问题的一种基本方法,它利用计算机运算速度快、适合做重复性操作的特点,让计算机对一组指令(或一定步骤)进行重复执行,在每次执行这组指令(或这些步骤)时,都从变量的原值推出它的一个新值,迭代法又分为精确迭代和近似迭代。比较典型的迭代法如“二分法”和"...原创 2018-10-18 21:35:32 · 8128 阅读 · 2 评论 -
最小二乘法拟合
import numpy as np#P145 例题2x=np.array([-1,-0.75,-0.5,-0.25,0,0.25,0.5,0.75,1])y=np.array([-0.2209,0.3295,0.8826,1.4392,2.0003,2.5645,3.1334,3.7601,4.2836])#习题六 2x=np.array([1,3,4,5,6,7,8,9,10]...原创 2018-12-04 20:46:55 · 851 阅读 · 0 评论 -
复化梯形公式和复化Simpon公式计算积分的近似值
习题七 10.1import numpy as npk=np.arange(0,9)xi=np.arange(0,1.1,0.125)yi=xiprint(k)print(xi)for x in range(0,9): yi[x]=xi[x]/(1+xi[x]**2)print(yi)ans=0for x in range(1,8): ans=ans+yi[...原创 2018-12-04 21:59:43 · 9368 阅读 · 0 评论 -
对分区间法求解方程组
import numpy as npimport matha=1b=2for x in range(0,20): #Xn=(a+b)/2 Xn=round((a+b)/2,8) fxn=round(Xn**3+4*Xn**2-10,8) Xn=(a+b)/2 fxn=Xn**3+4*Xn**2-10 #fxn=Xn**3+4*Xn**2...原创 2018-12-05 14:55:48 · 2239 阅读 · 0 评论 -
Newton法解非线性方程
习题八 10.(2)#P218 例4# def f(x):# return (x**3+10*x-20) #'''定义 f(x) = (x-3)^3''' # def fd(x):# return (3*x**2+10) #'''定义 f'(x) = 3*((x-3)^2)'''# xi=1.5# for i in range(1,10):#...原创 2018-12-06 16:06:55 · 469 阅读 · 0 评论 -
Newton法解非线性方程组
#P223 例7import numpy as npdef f1(x1,x2): return (4*x1**2+3*x2**2-1)def f2(x1,x2): return (x1**3-8*x2**3-1)def fd11(x1): return 8*x1def fd12(x2): return 6*x2def fd21(x1): r...原创 2018-12-06 16:54:51 · 2552 阅读 · 0 评论 -
欧拉方法
#P233 例题1import numpy as npx=0y=0yn=0for i in range(0,11): x1=round(x+0.1,6) y1=round(y+0.1*(x-y),6) yn=round(x+np.e**(-x)-1,6) ans=round(yn-y,6) print(i,x,y,yn,ans) x=x1...原创 2018-12-06 19:34:07 · 2448 阅读 · 0 评论 -
改进欧拉法
#P236 例2import numpy as npx=0y=0yn=1h=0.1accuracy=6for i in range(0,10): yp=round(y+h*(x-y),accuracy) yq=round(y+h*(x-yp+h),accuracy) #yn=round(x*(3+x**2)/3/(1+x**2),accuracy)...原创 2018-12-06 20:53:08 · 4697 阅读 · 0 评论 -
Runge-Kutta法
#P240 例3import numpy as npaccuracy=6def Fx(x): return x-1+np.e**(-x)def fx(x,y): return x-yx=0y=0h=0.2for i in range(0,5): k1=fx(x,y) k2=fx(x+h/2,y+h/2*k1) k3=fx(x+h/2...原创 2018-12-06 22:11:00 · 938 阅读 · 0 评论