Thảo luận nhóm TMU 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

plt.xlabel'Truc x'plt.ylabel'Truc y'plt.showplt.xlabel'Truc x'plt.ylabel'Truc y'plt.show ValueError: math domain error - lỗi vì tập xác định của arcsinx là x € -1;1 TIEU LUAN MOI downloa

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

T

Hà Nội, ngày 5/10/2021Nhóm trưởng

Trang 3

NHIỆM VỤ CỦA NHÓM

NV1 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).

NV2 Tính các tích phân trong bài 8.1 với cận 3 đến 12.

NV3 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

Trang 4

from sympy import *

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

from sympy import *

Trang 5

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

from sympy import *

Trang 6

TIEU LUAN MOI download : skknchat123@gmail.com

import numpy as np

from sympy import *

Trang 7

z_x = lambdify ((x, y), z_x)

print (z_x(3,12)

import numpy as np

from sympy import *

Trang 8

from sympy import *

Trang 9

z_xx = lambdify((x, y), z_xx)

print(z_xx(3,12))

import numpy as np

from sympy import *

Trang 10

from sympy import *

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

z_y = diff(z, y)

z_yy = (diff(z_y,y))

z_yy = lambdify((x, y), z_yy)

OverflowError: integer division result toolarge for a float

z ' y(3,12) = 312

Trang 11

Trang 12

from sympy import *

Trang 13

import numpy as np Màn hình hiển thị:

from sympy import *

Trang 14

from sympy import * -1.4366296429295926

Trang 15

from scipy import integrate

f = lambda x: math.sqrt(3*x+5) / x

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

print(result)

6.7758332587582695,1.1351502655320377e-10

1−x2 > 0-1 < x < 1

Giá trị tích phân là:

0.0817030676982187

Trang 16

4−x2 > 0-2 < x < 2

Trang 17

Màn hình hiện kết quả là: Điều kiện bài

Trang 18

0 0

Trang 19

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()

import math

Trang 20

Trang 21

return dydxy3 = 2

Trang 22

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

Trang 23

y3 = 2

Trang 24

def model(y,x):

dydx = (2*x**4 + 2*y)/xreturn dydx

y3 = 2

Trang 25

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

Trang 26

0 0

Trang 27

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].

import mathimport numpy as npfrom scipy.integrate import odeintimport matplotlib.pyplot as plt

import math

Trang 28

import numpy as npfrom scipy.integrate import odeintimport matplotlib.pyplot as plt

0 0

Trang 29

def model(U, x):

return [U[1], 4*U[1] - 12*x**2 - 6*x - 4]

U0 = [2, 1]

xs = np.linspace(3, 12)

Trang 30

Us = odeint(model, U0, xs)

ys = Us[:,0]

plt.xlabel("Truc x")plt.ylabel("Truc y")plt.title("PTVP")plt.plot(xs,ys)plt.show()

Trang 31

Trang 32

7. y ' ' − y =2 sin x−4 cos x

0 0

Trang 33

8. y ' ' −4 y ' =e x [(−4 x +4 ) cos x−(2 x +6 )sin x ]

Trang 34

9. y ' ' + y =cos x+cos 2 x

Trang 35

10. y ' ' − y =2 e x −x2

Trang 36

0 0

Trang 37

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].

for n in range(1, N):

x[n] = n*x[n-1]+math.factorial(n)plt.plot(x)

2 x(n+1) + nx(n) = n !

Trang 38

x[n] = math.factorial(n-1)-(n-1)*x[n-1]

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

3 x(n+1) – 9n x(n) = 3n2 +

2 n import mathimport numpy as np

import matplotlib.pyplot as plt

Trang 39

N = 9

x = np.zeros(N, float)x[0] = 3

x[n] = math.pow(9,n-1)*x[n-1]+math.pow(3,(n-1)**2+2*(n-1))plt.plot(x)

Trang 40

x[n] = math.factorial(n-1)/math.pow(3,n-1)-(n-1)/3*x[n-1]

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

x[n] =(n-1)/2*x[n-1]+math.factorial(n-1)/math.pow(2,n)plt.plot(x)

0 0

Trang 41

x[n] = math.pow(4,-(n-1))*x[n-1]+(n-1)*math.pow(2,-(n-1)**2-(n-1))plt.plot(x)

Trang 42

def model(y,x):

dydx = x*math.log(x) + y/(x*math.log(x))

return dydx y3 = 2

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

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

Trang 45

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].

import math import numpy as np from scipy.integrate import odeint import matplotlib.pyplot as plt

Trang 46

import math

import numpy as np from scipy.integrate import odeint import matplotlib.pyplot as plt

def model(U, x):

TIEU LUAN returnMOI[U[1],download9*U[0]+ math:skknchat123@gmail.exp(3*x)*math.cos(x)].com

Trang 47

import math import numpy as np

from scipy.integrate import odeint

import matplotlib.pyplot as plt

Trang 48

Trang 49

4. y ' ' −4 y ' =−12 x 2−6 x −4

def model(U, x):

return [U[1], 4*U[1] - 12*x**2 - 6*x - 4]

U0 = [2, 1]

xs = np.linspace(3, 12)

Trang 50

Us = odeint(model, U0, xs)

ys = Us[:,0]

plt.xlabel("Truc x") plt.ylabel("Truc y") plt.title("PTVP") plt.plot(xs,ys) plt.show()

def model(U, x):

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

Trang 52

plt.plot(xs,ys) plt.show()

Trang 54

import mathimport numpy as np from scipy.integrate import odeint import matplotlib.pyplot as plt

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