BÁO CÁO THỰC HÀNH MATLAB pdf

Bài 3.1 a Hãy thiết lập một ma trận Acoi là ma trận gốc với các phần tử bất kỳ và với số dòng bằng số kí tự của họ và số cột bằng số kí tự của tên người giải.b Ma trận B có số cột gấp đô

BÁO CÁO THỰC HÀNH

MATLAB

yi=7.4882e-002

Bài 1.2 Một mạch điện xoay chiều có các đại lượng điện áp U=Ur +jUx (V) và giá trị tổng trở Z=R+jX (om);(dữ liệu lấy trong bảng theo chữ cái họ tên và tên đệm người giải)

a) Thực hiện phép các thao tác cần thiết để hiển thị tên người giải trên màn hình.b) Sử dụng các lệnh Matlab để xác định giá trị dòng điện chạy trong mạch

c) Biểu thị kết quả tính toán dòng điện dưới dạng Imjα, với Im là modul và α là góc pha

a/ >> ten= 'Nguyen Quang Tung'

ten =

Nguyen Quang Tung

b/ >>U=116.2+110.4i; Z=6.73+4.67i; I=U/Z

Bài 1.3 Điện áp định mức của mạng điện là U kV, điện thành phần trở tác dụng là

Ω, công suất truyền tải trên đường dây là S=P+j*Q Hãy áp dụng các lệnh Matlab

để xác định tổn thất điện năng và áp dụng lệnh fprintf để biểu thị kết quả trên màn

hình với n chữ số sau dấu phẩy tĩnh/động(t/đ), biết thời gian tổn thất công suất cực đại trong năm là τ h Các dữ liệu tính toán lấy trong bảng 1.2 với họ tên người giải

>> P=45.33; Q= 36.2;U= 10; R=5.87; to=3890;

deltaP= (P^2+Q^2)*R/(U^2);

fprintf('deltaA=%10.4e kWh\n',deltaP*to)

deltaA=7.6843e+005

Bài 1.4 Thực hiện phép tính toán (dữ liệu lấy theo số thứ tự ghi trong đầu bài của

người giải theo bảng 1.3):

a)Số Z1 cho dưới dạng đại số(hoặc dạng mũ), hay chuyển về dạng mũ(hoặc dạng đại số), kiểm tra và ghi lại kết quả

b) Số Z2 cho dưới dạng đại số(hoặc dạng mũ), hay chuyển về dạng mũ(hoặc dạng đại số), kiểm tra và ghi lại kết quả

c) Hãy tính giá trị của biểu thức và ghi lại kết quả dưới dạng mũ, đối số(argument) biểu thị trong khoảng từ (-pi:pi)

Bài 2.2 Hãy xây dựng hàm xác định điện trở của mạch gồm n nhánh song song và

áp dụng hàm vừa xây dựng để tính điện trở tương đương của mạch gồm 5 nhánh song song với các điện trở lấy theo chữ cái đầu của họ và tên người giải trong bảng 2.1:với

R1=4,78; R2=5,35; R3=5,34 ; R5=6.54

Bài 2.3 Hãy hiển thị hàm f1(x), cho trong bảng 2.2 với chữ cái đầu của tên người

giải, xác định giá trị của hàm này với x=1.3

Bài 2.4 Biểu thị hàm véc tơ F=[f2 f3], cho trong bảng 2.2 với chữ cái đầu của họ

và tên người giải, dưới dạng inline objects và xác định giá trị của hàm ứng với x1

và x2

Bài 3.1 a) Hãy thiết lập một ma trận A(coi là ma trận gốc) với các phần tử bất kỳ

và với số dòng bằng số kí tự của họ và số cột bằng số kí tự của tên người giải.b) Ma trận B có số cột gấp đôi số cột của ma trận gốc;

c) Hãy thiết lập ma trận C gồm các phần tử từ dòng thứ 1 dến dòng thứ 2 và từ cột thứ 1 dến cột thứ 3 của ma trận gốc

d)Hãy thiết lập ma trận D cùng kích thước như ma trận A

e) Hãy gọi lại các phần tử của ma trận gốc:

Bài 4.1 Cho ma trận A như bài toán 3.1, hãy xác định :

a) Giá trị cực đại của mỗi cột (C1);

b) Giá trị cực tiểu của mỗi cột (C2);

c) Giá trị cực đại của mỗi dòng (C3);

d) Giá trị cực tiểu của mỗi dòng (C4);

e) Giá trị cực đại của cột và chỉ số của chúng (C1, I);

a/ >> C1=max(A)

C1 =

5 9 9 6

Bài 4.2 Số liệu thống kê về các đại lượng x và y cho trong bảng 4.1, lấy theo chữ

cái đầu của tên người giải:

1) Hãy sử dụng các hàm trong Matlab để xác định các đại lượng:

1.1) giá trị trung bình của các đại lượng x1,x2 và y: Xtb1,Xtb2,Y_tb;

1.2) Phương sai của các đại lượng Cx1, Cx2, Cy:

1.3) Độ lệch chuẩn của các đại lượng x1, x2 và y

1.4 ) Độ lệch trung bình bình phương của các đại lượng x1, x2 và y

1.5) Hệ số tương quan của các giá trị quan sát giữa các đại lượng Rx1_y, Rx2_y:2) Hãy xây dựng hàm hồi quy thực nghiệm và đánh giá sai số của các hàm:2.1) Bậc hai: y=f(x1); y=f(x2)

2.2) Hàm tuyến tính: y=f(x1); y=f(x2)

2.3) Hàm hội quy bội: y=f(x1,x2)

1.1)

>> x1=[2.83 5.64 7.85 10.94 13.54 16.33 19.24 21.45 26.52 31.12];

>> x2=[5.65 7.36 8.77 16.24 19.33 22.63 24.22 26.82 30.32 37.43];

3.0930e-014

Bài 4.3 Hãy xác định các giá trị hàm nội suy của ma trận y ứng với các giá trị của

x(dữ liệu lấy theo chữ cái đầu của tên người giải trong bảng 4.2):

Bài 4.4 Xây dựng hàm y trong khoảng biến thiên x(theo bảng 4.2) và các hàm nội

suy: bậc thang, tuyến tính hóa, lập phương mô phỏng và hàm lập phương và đánh giá sai số so với hàm gốc

'tuyen tinh hoa',

'lap phuong mo phong',

Bài 5.1 Vẽ đồ thị hàm số y trong khoảng biến thiên x=a dến b bới bước tính h(lấy

theo chữ cái đầu của tên người giải trong bảng 5.1)bằng lệnh:

grid;

xlabel('truc x'); ylabel('truc y');

title(' Ve ham so y=f(x) voi lenh plot')

Bài 5.2 Vẽ đồ thị hàm số y trong khoảng biến thiên của x (xem dữ liệu bảng 5.2

Dữ liệu được lấy ứng với chữ cái đầu theo họ của người giải):

x=0:2*pi;

y=2.5.*x.^2+3.7.*x+6;

a) Với các trục tọa độ chia theo tỷ lệ log cơ số 10;

b) Với các trục x chia theo tỷ lệ log cơ số 10;

c) Với các trục y chia theo tỷ lệ log cơ số 10;

a)

>> x=0:2*pi;

y=2.5.*x.^2+3.7.*x+6;

Bài 5.3 Xây dựng đồ thị hàm d trong khoảng biến thiên của phi(xem dữ liệu bảng

5.1) trong hệ tọa độ cực( dữ liệu được lấy ứng với chữ cái đầu theo tên đệm của người giải)

Bài 5.4 Vẽ đồ thị của hàm Z trong khoảng biến thiên x và y (xem dữ liệu bảng 5.2)

với bề mặt đặc; bề mặt các đường mức; bề mặt trong suốt và bề mặt thảng đứng (thác nước) lần lượt trên 4 phần của một trang giấy.( Dữ liệu được lấy ứng với chữ cái đầu theo tên người giải)

Norm of First-order Trust-region

Iteration Func-count f(x) step optimality radius

1 3 1830.39 2.1e+003 1

2 6 158.656 1 258 1

3 9 4.81873 0.864875 25.8 2.5

4 12 0.129974 0.179451 2.87 2.5

5 15 0.000376465 0.046744 0.131 2.5

6 18 5.44277e-009 0.00298961 0.000432 2.5

7 21 1.35665e-018 1.1771e-005 5.67e-009 2.5

Optimization terminated successfully:

First-order optimality is less than options.TolFun

r =

0 -3 31

c)

function addpoly(f,g)

nf= length(f); % do dai cua da thuc f

ng= length(g); % do dai cua da thuc g

Chương 9 : Giải mạch điện

Bài tập 9.1

Cho mạch điện một chiều (hình vẽ) với suất điện động Ei và các điện trở:

Ri (Ω) cho trong bảng 9.1.(dữ liệu lấy theo chữ cái đầu của tên người giải) Hãy áp dụng chương trình Matlab xác định các giá trị dòng điện chạy trong mạch

Số liệu bảng 9.1 như sau:

R1=2.2 Ω; R2=2.05 Ω; R3=1.7 Ω; R4=3.5 Ω; R5=2.8 Ω;

E1=100 V; E2=96 V; E3=0 V ; E4=0 V ; E5=0 V;

>> disp('giai mach dien theo phuong phap dong dien nhanh');

R=[2.2 2.05 1.7 3.5 2.8]; % ma tran dien tro nhanh

E=[100 96 0 0 0]; %ma tran suat dien dong nguon

Z=[1,1,1,1,1;R(1),-R(2),0,0,0;0,R(2),-R(3),0,0;0,0,R(3),-R(4),0;0,0,0,R(4),-R(5)];

% ma tran dien tro mach vong

U=[0;E(1)-E(2);E(2)-E(3);E(3)-E(4);E(4)-E(5)]; %ma tran dien ap

I=Z\U; %dong dien tren cac nhanh

disp('ket qua la:')

disp(' I1 I2 I3 I4 I5')

fprintf('%g'),disp([I(1) I(2) I(3) I(4) I(5)])

giai mach dien theo phuong phap dong dien nhanh

ket qua la:

I1 I2 I3 I4 I5

26.1547 26.1172 -24.9763 -12.1314 -15.1642

Bài tập 9.2

Giải bài toán ví dụ 9.1 với các tham số phức bảng 9.2

Số liệu như sau:

Z1=130+j*120Ω; Z2=0 Ω ; Z3=125+j*94 Ω ; Z4=128+j*112 Ω ; Z5=0 Ω;

E1=4+j*6.3 V; E2=6.2+j*7.2 V; E3=6.2+j*8.5 V; E4=5.4+j*6.0 V; E5=0 V;

disp('giai mach dien theo phuong phap dong dien nhanh');

Z=[4+j*6.3;6.2+j*7.2;6.2+j*8.5;5.4+j*6.0;0]; % ma tran dien tro nhanh

E=[130+j*120;0;125+j*94;128+j*112;0]; %ma tran suat dien dong nguonZ=[1,1,1,1,1;Z(1),-Z(2),0,0,0;0,Z(2),-Z(3),0,0;0,0,Z(3),-Z(4),0;0,0,0,Z(4),-Z(5)];

% ma tran dien tro mach vong

U=[0;E(1)-E(2);E(2)-E(3);E(3)-E(4);E(4)-E(5)]; %ma tran dien ap

I=Z\U; %dong dien tren cac nhanh

disp('ket qua la:')

disp(' I1 I2 I3 I4 I5')

fprintf('%g'),disp([I(1) I(2) I(3) I(4) I(5)])

giai mach dien theo phuong phap dong dien nhanh

ket qua la:

I1 I2 I3 I4 I5

22.9126 - 6.0873i 0.0000 - 0.0000i 14.2199 - 4.3337i 20.9208 - 2.5046i -58.0533 +12.9256i

Bài tập 9.3

Hãy xác định dòng điện và công suất của mạch điện(hình vẽ) biết điện áp

u(i)=Um*cos(ωt+thetau) và phụ tải Z=Zm∠thetaz (Dữ liệu trong bảng 9.2 theo chữ cái của họ người giải) Vẽ đồ thị các đại lượng

Số liệu bảng 9.2 như sau:

Um=95V; thetau=220; Zm=3.25Ω; thetaz=820;

>> disp('Xac dinh dong dien va cong suat Ve do thi cac dai luong');

Um=95;

thetau=22;

Bài tập 9.4

Nguồn điện xoay chiều(hình vẽ) có điện áp u(t)=Um*sin(ωt), tần số 50Hz được nối với điện trở R và cuộn dây L, Điện trở sun Rsh (dữ liệu bảng 9.2 theo chữ cái của tên đệm của người giải) Hãy xác định giá trị hiệu dụng của điện áp URL, các thành phần UR, UL Vẽ đồ thị các đại lượng này

i

Số liệu bảng 9.2 như sau:

Số liệu bảng 9.2 như sau:

Um = 120,V; R =31.5Ω; C =120μF;

R

i C

>> disp('Gia tri hieu dung dien ap va dong dien chay trong mach');

title('Do thi dien ap: Sinh vien:Nguyen Quang Tung-Lop D3H3 ')

Gia tri hieu dung dien ap va dong dien chay trong mach

I,A Ur,V Uc,V

2.0605 64.9054 54.6561

Bài tập 9.6

Hãy xác định dòng điện và công suất tác dụng trên đầu vào mạng hai cực (hình vẽ), biết điện áp là U, các giá trị điện dẫn tác dụng và phản kháng tương ứng là G và B (dữ liệu lấy theo bảng 9.3 theo chữ cái đầu của tên người giải)

Số liệu bảng 9.3 như sau:

∠ Hãy xác định các giá trị hiệu dụng của dòng điện chạy trong mạch, biết C,

R, L (lấy trong bảng dữ liệu 9.3 theo chữ cái đầu của tên người giải)

Số liệu bảng 9.3 như sau:

>> disp('Cac Gia tri hieu dung cua dong dien chay trong mach dien');

áp nguồn 1 từ -300:+300.

Số liệu bảng 9.4 như sau:

e1=125; e2=122; a1=-20; a2=1,7;Z=0.8+j*1.6;

e1 Z = R + j*X e2