báo cáo thí nghiệm lý thuyết mạch 2

ĐẠI HỌC BÁCH KHOA HÀ NỘITrường Điện- Điện tửKhoa Điện Báo cáo Thí nghiệm: Lý thuyết mạch 2Họ và tên: Nguyễn Văn AMSSV: 20XXXXXXLớp: Mã lớp : Nhóm: Hà Nội, tháng 3 năm 2022... NỘI DUNG TH

ĐẠI HỌC BÁCH KHOA HÀ NỘITrường Điện- Điện tử

Khoa Điện

Báo cáo Thí nghiệm: Lý thuyết mạch 2Họ và tên: Nguyễn Văn A

MSSV: 20XXXXXXLớp:

Mã lớp : Nhóm:

Hà Nội, tháng 3 năm 2022

trigơ và ổn áp sắt từII NỘI DUNG THÍ NGHIỆM

1 Đặc tính U(I) của cuộn dây lõi thép

- Tăng U từ 0V =>80VU(V

8 0.022 0.033U(V

I(A) 0.04

9 0.071 0.16 0.278 0.400 0.518 0.631 0.745- Giảm U từ 80V => 0V

I(A) 0.74

5 0.636 0.545 0.404 0.292 0.143 0.076 0.044 0.031U(V

0.336- Giảm U từ 80V => 0V

I(A) 0.336 0.322 0.296 0.273 0.253

0.201U(V

I(A)

II HƯỚNG DẪN SỬ DỤNG PCMATLAB để tính chế độ xác lập dừng trong mạch điện phi tuyến.

I2(U)=2.10 U (1)4 2

và đường thẳngI1(U) = (U - U)/R01(2)

0.3.20.2

Y: 0.05120

05 1 152030354045U

(V)H.2

Giao điểm của hai đường I và I sẽ cho các trị số của U và I (xem12

H.2) Chương trình tính viết bằng ngôn ngữ MATLAB như sau:

Ta ký hiệu y bằng x và viết phương trình (3) dưới dạng x = φ(x) (4)Ta sẽ tính lặp liên tiếp các trị số của x theo

công thức x… = φ(x) với k = 0, 1, 2, 3, …(5)cho đến khi đạt yêu cầu về độ chính xác |x…- x…| ≤ ε

Ở đây ε là một lượng nhỏ đã cho, ví dụ ε = 0.0001 với k = 0 ta có thể chọn xn là một trị số tùy ý, ví dụ xn = 0.Dễ thực hiện việc tính lặp bằng MATLAB ta dùng lệnh sau:

Ở Sau lệnh này sẽ có câu hỏi yêu cầu điền các giá trị với U0 = 20 ↲

R1=500↲Xn = 0 ↲ε = 0.0001 ↲

ở đây itera_1 là tên một chương trình lặp viết bằng ngôn ngữ MATLABđể giải bài toán đã cho.

Để đọc các trị số liên tiếp của xk từ màn hình ta dùng lệnh X (1:k)III NHIỆM VỤ.

Giải bằng phương pháp đồ thị và phương pháp tính lặp một mạch điện phi tuyến cho bởi giáo viên.

1 Phương pháp đồ thị

• Chương trình file mR1 = 500; U0 = 41; deltaU = 1;U = 0 : deltaU : U0;

I1=U0/R1*ones(1, length(U))-U/R1; plot(U,I1,U,I2);

grid I1 I2 U

• Kết quả

epsilon = input('Nhap gia tri cho Epsilon: ');

Xn = input('Nhap gia tri cho Xn: '); %Nhap gia tri cho rang gan voi nghiem chinh xac

str = 1; Max = 1000; % So lan lap toi da Tranh lap vo han khi epsilon qua nho.

while str == 1

x1 = g(Xn); cou = 1; %Lan lap dau tien

u(1,1) = Xn; u(1,2) = x1;u(1,3) = abs(x1 - Xn);

while (abs(x1 - Xn) >= epsilon) & (cou <= Max)cou = cou +1;

Xn = x1;x1 = g(Xn);

u(cou,1) = Xn; u(cou,2) = x1; u(cou,3)=abs(x1 - Xn);

%plot(x1, g(x1))

%hold on

if f(x1) == 0

disp 'Nghiem chinh xac: '

function [res_f] = f(x)

res_f = x- (20-100*sqrt(x/2))/500 ;

function [res_g] = g(x)

So lan lap:cou =

Nghiem xap xi:ans =

0.020017715498197