![](https://img-blog.csdnimg.cn/20190928154355986.jpg?x-oss-process=image/resize,m_fixed,h_224,w_224)
数值分析实验
数值分析实验作业
xfdelovely
这个作者很懒,什么都没留下…
展开
-
数值分析实验 实验5-1 数值微分的几个方法 python3实现
题目代码import math# 求解区间a, b = list(map(int, input("请输入求解区间,格式为 a b :").split()))# 初值x0 = float(input("请输入初值:"))h = 0.1# yn的精确值def yn(x): return math.sqrt(1 + 2 * x)# fxdef fx(x, y...原创 2019-10-12 21:33:04 · 495 阅读 · 0 评论 -
数值分析实验 实验4-2 选主元高斯消去法 python3实现
题目代码import numpy as np# 行列式的维度n = int(input("请输入行列式的维度:"))# 要求解的行列式x = np.empty([n, n])print("请输入行列式:")for i in range(n): x[i] = np.array(list(map(float, input().split())))# 行列式的符号re...原创 2019-10-12 21:29:50 · 765 阅读 · 0 评论 -
数值分析实验 实验4-1 高斯-塞德尔迭代法 python3实现
题目代码import numpy as np# 系数矩阵A的行数m = int(input("请输入系数矩阵A的行数:"))# 系数矩阵A的列数n = int(input("请输入系数矩阵A的列数:"))# 系数矩阵Aa = np.empty([m, n])print("请输入系数矩阵A:")for i in range(m): a[i] = np.array(l...原创 2019-10-12 21:27:04 · 435 阅读 · 0 评论 -
数值分析实验 实验3-1 牛顿下山法 python3实现
题目代码# 存放各步的迭代值result = []# 存放各步的下山因子all_r = []# 判断是否奇异,True为奇异,False为非奇异def strange(xk): return True if (3*xk*xk-1) == 0 else False# 原函数def fx(xk): return xk * xk * xk - xk - 1...原创 2019-10-12 21:23:11 · 1170 阅读 · 0 评论 -
数值分析实验 实验1-2 牛顿插值公式 python3实现
题目:代码:#数据点集的大小n = int(input("请输入数据点集的大小:"))#数据点集dic = {}for i in range(n): a, b = list(map(float, input("请输入数据点对,格式为 x y :").split())) dic[a] = bpoints = list(dic.keys())#预测点pre = fl...原创 2019-09-28 15:56:54 · 738 阅读 · 0 评论 -
数值分析实验 实验1-1 拉格朗日插值公式 python3实现
题目:代码:#预测点pre_points = list(map(float, input("请输入预测点序列,格式为x1 x2 x3 :").split()))#构造插值点,为等差数列input_points = list(range(-5, 6))#预测点的值result = []def fx(x): return 1 / (1 + x * x)#拉格朗日差值...原创 2019-09-28 16:22:49 · 953 阅读 · 0 评论 -
数值分析实验 实验2-1 复化梯形公式 python3实现
题目:设计复化梯形公式求积算法代码:# 要积的函数,字符串形式fx_str = str(input("请输入要积的函数,格式为 x*x+1 (即变量用x表示,乘号用*表示):"))# 求积的区间[a, b]a, b = list(map(float, input("请输入求积区间[a, b],格式为 a b :").split()))# 区间划分的份数n = int(input("...原创 2019-09-28 17:25:56 · 1600 阅读 · 0 评论 -
数值分析实验 实验2-2 复化辛甫生公式 python3实现
题目:设计复化辛浦生求积算法代码:# 要积的函数,字符串形式fx_str = str(input("请输入要积的函数,格式为 x*x+1 (即变量用x表示,乘号用*表示):"))# 求积的区间[a, b]a, b = list(map(float, input("请输入求积区间[a, b],格式为 a b :").split()))# 区间划分的份数n = int(input("请...原创 2019-09-28 17:43:09 · 381 阅读 · 0 评论 -
数值分析实验 实验2-3 龙贝格公式 python3实现
题目:代码:import math# 求积区间a, b = 0, 1# 区间二分的最大次数n = 4# 分点x = []# 要积的函数def fx(xk): if xk == 0: return 1 return math.sin(xk) / xk# 梯形公式 x2 > x1def trapezium(x1, x2): ...原创 2019-09-28 19:48:21 · 837 阅读 · 0 评论