数学建模python实现基础编程练习1

数学建模python实现编程练习1

习题来自《python数学实验与建模》

1.利用sympy库进行画单叶双曲面

用pyplot画的不是很理想，就用sympy画了，利用plot3D函数

from sympy.plotting import plot3d
from sympy.abc import x,y
from sympy.functions import sqrt
plot3d((sqrt(((x**2)/8+(y**2)/10-1)*6),(x,-6,6),(y,-6,6)),(-sqrt(((x**2)/8+(y**2)/10-1)*6),(x,-6,6),(y,-6,6)))

这里范围取在了（-6,6）

2.利用scipy求解积分数值解

from numpy import inf
from scipy.integrate import quad
f = lambda x: np.sqrt(1+4*x)
print('Q1',quad(f,0,1))#quad() 一般取第一个数
f = lambda x: np.exp(-x)*np.sin(x)
print('Q2',quad(f,0,inf))

3. 利用sympy库求符号解

x,y = sp.symbols('x y')
print(sp.solve(x**3+6*x-4*x**2-8))

4.利用numpy组合不同矩阵并求行列式

import numpy.linalg as la
a1 = np.array([[1,3,5],[2,4,6]]).T
a2 = np.array([[1,2,3],[1,2,4]]).T
a3 = np.array([[2,6]])
a4 = np.array([[3,2]])
b1 = np.c_[a1,a2]
b2 = np.c_[a3,a4]
b3 = np.r_[b1,b2]
print('|b3|=',la.det(b3))

5.根据秩大小求解非齐次线性方程组

import sympy as sp
from numpy import linalg as la
a = np.array([[2,3,1],[1,-2,4],[3,8,-2],[4,-1,9]])
b = np.array([[4,-5,13,6]])
c = np.c_[a,b.T]
print('r(a)',sp.Matrix(a).rank())
print('r(a,b)',sp.Matrix(c).rank())
'''r(a)<r(a,b),原非齐次线性方程组无解，这里给出最小二乘解'''
print('最下二乘解=\n',la.pinv(a).dot(b.T))