Thảo luận nhóm TMU báo cáo TỔNG hợp nhiệm vụ sử dụng python để giải các bài toán

BÔ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC THƯƠNG MẠI KHOA HỆ THỐNG THÔNG TIN KINH TẾ & THƯƠNG MẠI ĐIỆN TỬBẢN BÁO CÁO TỔNG HỢP HỌC PHẦN TOÁN CAO CẤP 2 Nhiệm vụ: “Sử dụng Python để giải các bà

Trang 1

BÔ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC THƯƠNG MẠI KHOA HỆ THỐNG THÔNG TIN KINH TẾ & THƯƠNG MẠI ĐIỆN TỬ

BẢN BÁO CÁO TỔNG HỢP

HỌC PHẦN TOÁN CAO CẤP 2

Nhiệm vụ: “Sử dụng Python để giải các bài toán”.

Lớp học phần: 2192FMAT0211 Nhóm thực hiện: 03

Giảng viên giảng dạy: Lê Văn Tuấn

HÀ NỘI – 2021

Trang 2

DANH SÁCH THÀNH VIÊN NHÓM 03 – LỚP HP: 2192FMAT0211

Nhóm trưởng: Lê Thị Giang

ST

21 Đinh Thùy Dương Làm nhiệm vụ 1 và kiểm tra các đáp án

22 Nguyễn Hưng Duy Làm nhiệm vụ 2

23 Lê Thị Giang Làm nhiệm vụ 6 và tổng hợp đáp án

24 Nguyễn Thị Hương Giang Làm nhiệm vụ 2

25 Nguyễn Thu Hà Làm nhiệm vụ 3 và word

27 Nguyễn Thị Hằng Làm nhiệm vụ 5

28 Nguyễn Thị Thúy Hằng Làm nhiệm vụ 4

29 Nguyễn Thúy Hằng Làm nhiệm vụ 1

30 Trần Thị Thu Hiền Làm nhiệm vụ 5

Trang 3

NHIỆM VỤ CỦA NHÓM

NHIỆM VỤ 1: Tính đạo hàm riêng cấp 1 & 2 của các hàm số trong bài 7.2 tại (3, 12) 1 NHIỆM VỤ 2: Tính các tích phân trong bài 8.1 với cận 3 đến 12 12 NHIỆM VỤ 3: Vẽ đồ thị nghiệm riêng của các PTVP trong bài 9.4 với điều kiện ban đầu y(3) = 2 trên miền [3, 12] 16 NHIỆM VỤ 4: Vẽ đồ thị nghiệm riêng của các PTVP trong bài 9.7 với điều kiện y(3) = 2, y’(3) = 1 trên miền [3, 12] 24 NHIỆM VỤ 5: Vẽ đồ thị nghiệm riêng của các PTSP trong bài 10.4 với điều kiện ban đầu x(0) = 3 trên miền [0, 9] 34 NHIỆM VỤ 6: Vẽ đồ thị nghiệm riêng của các PTSP trong bài 10.7 với điều kiện: x(0) = 3, x(1) = 4 trên miền [0, 9] 42

Trang 4

from sympy import *

x, y = symbols ('x y', real=True)

from sympy import *

Trang 5

from sympy import *

x, y = symbols('x y', real=True)

from sympy import *

z = log(x+sqrt(x**2 + y**2))

z_y = diff (z, y)

z ' y = 0.06312203124697226

Trang 6

z_y = lambdify ((x, y), z_y)

print("z ' y =",z_y(3,12))

import numpy as np

from sympy import *

Trang 7

from sympy import *

Trang 8

from sympy import *

z ' xx = 0.003170400327272324

Trang 9

from sympy import *

OverflowError: integer division result too

large for a float

Trang 10

from sympy import *

z= pow(x,(y**3))

z_y = diff (z, y)

z_y = lambdify ((x, y), z_y)

print("z ' y =",z_y(3,12))

large for a float

import numpy as np

from sympy import *

large for a float

import numpy as np

from sympy import *

large for a float

import numpy as np

from sympy import *

large for a float

Trang 11

from sympy import *

Trang 12

from sympy import *

Trang 13

from sympy import *

Trang 14

NHIỆM VỤ 2: Tính các tích phân trong bài 8.1 với cận 3 đến 12.

Trang 16

9.520275988605531e-from scipy import integrate

 0< x < 1

Trang 17

from scipy import integrate

f = lambda x: math.log(math.tan(x)) /(math.sin(x)*math.cos(x)))

result = integrate.quad(f, 3, 12)

print(result)

ValueError: math

domain error.

Trang 18

NHIỆM VỤ 3: Vẽ đồ thị nghiệm riêng của các PTVP trong bài 9.4 với điều kiện ban đầu

y(3) = 2 trên miền [3, 12].

y3 = 2

x = np.linspace(3,12)

y = odeint(model,y3,x)plt.plot(x,y)

plt.xlabel('Truc x')plt.ylabel('Truc y')plt.show()

Trang 19

y = odeint(model,y3,x)plt.plot(x,y)

3 y ' +2 xy=x e −x2

import math

import numpy as np

Trang 20

from scipy.integrate import odeint

Trang 21

from sympy import *

z = (x+y)**(sin(x)*sin(y))

z ' xy = -1.4366296429295926

Trang 22

from sympy import *

Trang 23

NHIỆM VỤ 2: Tính các tích phân trong bài 8.1 với cận 3 đến 12.

Trang 28

from scipy import integrate

f = lambda x: math.log(math.tan(x)) /(math.sin(x)*math.cos(x)))

ValueError: math

domain error.

Trang 29

14

Trang 30

NHIỆM VỤ 3: Vẽ đồ thị nghiệm riêng của các PTVP trong bài 9.4 với điều kiện ban đầu

y(3) = 2 trên miền [3, 12].

y3 = 2

y = odeint(model,y3,x) plt.plot(x,y)

plt.xlabel('Truc x') plt.ylabel('Truc y')

Trang 31

15

Trang 32

plt.xlabel('Truc x') plt.ylabel('Truc y') plt.show()

3 y ' +2 xy=x e −x2

import math

import numpy as np

Trang 33

from scipy.integrate import odeint

Trang 34

y ( ,y , ) plt.plot(x,y)

Trang 35

dydx = math.sin(x)*math.cos(x) + math.sin(x)*y return dydx

Trang 37

y3 = 2

ValueError: math domain error - lỗi vì tập xác định của arcsin(x) là x € [-1;1]

y3 = 2

Trang 38

19

Trang 39

y3 = 2

Trang 40

8. y ' √1−x2

+ y=arcsin x

Trang 41

ValueError: math domain error - lỗi vì tập xác định của arcsin(x) là x € [-1;1]

21

Trang 42

Trang 43

22

Trang 44

NHIỆM VỤ 4: Vẽ đồ thị nghiệm riêng của các PTVP trong bài 9.7 với điều kiện y(3) =

2, y’(3) = 1 trên miền [3, 12].

Trang 51

return [U[1], 2*U[1] + 3*U[0] + math.exp(4*x)]

Định dạng
Số trang	54
Dung lượng	1,26 MB