2024年10月15日

习题2.12

摘要: 1.代码实现 点击查看代码 import numpy as np from sympy import Matrix, symbols # 定义一个矩阵 A = np.array([[-1,1,0],[-4,3,0],[1,0,2]]) # 使用numpy.linalg.eig求解特征值和特征向量 e 阅读全文

posted @ 2024-10-15 16:58 克卜勒星球 阅读(16) 评论(0) 推荐(0)

习题2.8

摘要: 1.代码实现 点击查看代码 import numpy as np n=1000 A=np.zeros((n,n)) b=np.arange(1,n+1) np.fill_diagonal(A[1:,:-1],1) np.fill_diagonal(A[:-1,1:],1) np.fill_diago 阅读全文

posted @ 2024-10-15 16:56 克卜勒星球 阅读(28) 评论(0) 推荐(0)

习题2.7(2)

摘要: 1.代码实现 点击查看代码 import numpy as np a = np.array([[2, 3, 1], [1, -2, 4], [3, 8, -2], [4, -1, 9]]) b = np.array([4, -5, 13, -6]) # 最小二乘解 x3 = np.linalg.pi 阅读全文

posted @ 2024-10-15 16:54 克卜勒星球 阅读(19) 评论(0) 推荐(0)

习题2.7(1)

摘要: 1.代码实现 点击查看代码 import numpy as np a=np.array([[4,2,-1],[3,-1,2],[11,2,0]]) b=np.array([2,10,8]) #唯一解 x1=np.linalg.inv(a)@b x2=np.linalg.solve(a,b) prin 阅读全文

posted @ 2024-10-15 16:53 克卜勒星球 阅读(22) 评论(0) 推荐(0)

习题2.5(2)

摘要: 1.代码实现 点击查看代码 import numpy as np import matplotlib.pyplot as plt from mpl_toolkits.mplot3d import Axes3D # 设置椭圆抛物面的参数a和b a = 2 b = 1 # 定义x和y的范围 x = np 阅读全文

posted @ 2024-10-15 16:51 克卜勒星球 阅读(20) 评论(0) 推荐(0)

习题2.5(1)

摘要: 1.代码实现 点击查看代码 import numpy as np import matplotlib.pyplot as plt #横纵坐标 x=np.linspace(-5,5,100) y=np.linspace(-5,5,100) #网格生成 X,Y=np.meshgrid(x,y) #写法一 阅读全文

posted @ 2024-10-15 16:49 克卜勒星球 阅读(38) 评论(0) 推荐(0)

习题2.4

摘要: 1.代码实现 点击查看代码 import numpy as np import matplotlib.pyplot as plt # 定义x的范围 x = np.linspace(-10, 10, 400) # 创建一个2行3列的子图网格 fig, axs = plt.subplots(2, 3, 阅读全文

posted @ 2024-10-15 16:48 克卜勒星球 阅读(27) 评论(0) 推荐(0)

习题2.3

摘要: 1.代码实现 点击查看代码 import numpy as np import matplotlib.pyplot as plt # 定义x的范围 x = np.linspace(-10, 10, 400) # 生成从-10到10的400个点 # 初始化图形 plt.figure(figsize=( 阅读全文

posted @ 2024-10-15 16:47 克卜勒星球 阅读(33) 评论(0) 推荐(0)

习题2.2

摘要: 1.代码实现 点击查看代码 import numpy as np import matplotlib.pyplot as plt from scipy.special import gamma # 定义x的范围 x = np.linspace(0.1, 10, 400) # 注意从0.1开始,因为G 阅读全文

posted @ 2024-10-15 16:45 克卜勒星球 阅读(26) 评论(0) 推荐(0)

习题2.1

摘要: 1.代码实现 点击查看代码 import matplotlib as plt from cProfile import label import matplotlib.pyplot as plt import numpy as np from matplotlib import font_manag 阅读全文

posted @ 2024-10-15 16:44 克卜勒星球 阅读(30) 评论(0) 推荐(0)

ex6-10

摘要: 1.代码实现 点击查看代码 import numpy as np import networkx as nx L=[(1,2,20),(1,5,15),(2,3,20),(2,4,60),(2,5,25),(3,4,30),(3,5,18),(5,6,15)] G=nx.Graph() G.add_ 阅读全文

posted @ 2024-10-15 16:36 克卜勒星球 阅读(8) 评论(0) 推荐(0)

ex6-11

摘要: 1.代码实现 点击查看代码 import numpy as np import cvxpy as cp L=[(1,2,18),(1,5,15),(2,3,20),(2,4,60),(2,5,12),(3,4,30),(3,5,18),(4,6,10),(5,6,15)] a=np.ones((6, 阅读全文

posted @ 2024-10-15 14:05 克卜勒星球 阅读(9) 评论(0) 推荐(0)

习题3.2

摘要: 1.代码实现 点击查看代码 def generate_sequence(n): # 初始化序列列表,包含前两个初始项 sequence = [-2, -2] # 使用递推关系式填充序列的其余部分 for i in range(2, n): # xn+2 = xn+1 + 2*xn next_term 阅读全文

posted @ 2024-10-15 00:28 克卜勒星球 阅读(30) 评论(0) 推荐(0)

习题3.3

摘要: 1.代码实现 点击查看代码 import numpy as np from scipy.sparse.linalg import eigs import pylab as plt L=[(1,2),(1,3),(1,4),(1,5),(1,6),(2,4),(2,5),(2,6),(3,1),(3, 阅读全文

posted @ 2024-10-15 00:27 克卜勒星球 阅读(28) 评论(0) 推荐(0)

习题4.3

摘要: 1.代码实现 点击查看代码 import matplotlib.pyplot as plt import numpy as np import cvxpy as cp x=cp.Variable(6,pos=True) obj=cp.Minimize(x[5]) a1=np.array([0.025 阅读全文

posted @ 2024-10-15 00:26 克卜勒星球 阅读(21) 评论(0) 推荐(0)

习题5.3

摘要: 1.代码实现 点击查看代码 from scipy.optimize import minimize import numpy as np c1=np.array([1,1,3,4,2]) c2=np.array([-8,-2,-3,-1,-2]) A=np.array([[1,1,1,1,1],[1 阅读全文

posted @ 2024-10-15 00:25 克卜勒星球 阅读(13) 评论(0) 推荐(0)

习题5.5

摘要: 1.代码实现 点击查看代码 import numpy as np from scipy.optimize import minimize def objective(x): x1, x2, x3 = x return -(2 * x1 + 3 * x1 ** 2 + 3 * x2 + x2 ** 2 阅读全文

posted @ 2024-10-15 00:18 克卜勒星球 阅读(19) 评论(0) 推荐(0)

习题5.4

摘要: 代码实现 点击查看代码 import numpy as np from scipy.optimize import minimize def objective(x): return -np.sqrt(x[0]) def constraint1(x): return x[0] - 10 def co 阅读全文

posted @ 2024-10-15 00:16 克卜勒星球 阅读(36) 评论(0) 推荐(0)

导航