Lê Hải Đăng
Điện 5 – K54
MSSV : 20090710
Bàn 10

Báo cáo thí nghiệm

BÀI THÍ NGHIỆM TRƯỜNG ĐIỆN TỪ SỐ 1
Giải phương trình Poisson và phương trình Laplace
dạng sai phân bằng Matlab

Nội dung :
Giả sử ta cần tính sự phân bố của điện thế V của điện trường tĩnh trong một
miền phẳng hình chữ nhật ABCD với các kích thước a=15m, b=16m, nếu đã cho
điện thế V=0 trên các cạnh AB, BC, CD, DA và mật độ điện tích khối tự do phân bố
đều ở miền bên trong các biên ρ = 2 ε
Để giải phương trình Possion dạng sai phân bằng Matlab, ta thực hiện các
lệnh sau :

% Định nghĩa các thông số đã cho của bài toán

n = 16; m = 17; h = 1; rotd = 2; delta = 0.01;

% ác định các đi u i n i n

i = 1; for j = 1: m, V0(i,j) = 0 ; end;
i = n; for j = 1: m, V0(i,j) = 0 ; end;
j = 1; for j = 1: n, V0(i,j) = 0 ; end;
j = m; for j = 1: n, V0(i,j) = 0 ; end;

% Thực hi n l nh


7

8

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16

0.000
0.000
0.000
0.000
0.000
0.000
0.000
0.000

13.778
15.199
16.105
16.548
16.553
16.121
15.224
13.807
11.779
8.988
5.193
0.000

0.000
6.671
11.766
15.618
18.468
20.484
21.775
22.409
22.417
21.798
20.518
18.510
15.660
11.802
6.693
0.000


15.366
8.554
0.000

0.000
9.040
16.316
22.042
26.401
29.544
31.582
32.589
32.602
31.617
29.598
26.466
22.108
16.372
9.074
0.000

0.000
9.340
16.898
22.872
27.437
30.737
32.882
33.942
33.955

27.784
31.138
33.318
34.396
34.409
33.355
31.194
27.852
23.219
17.150
9.475
0.000

0.000
9.346
16.909
22.888
27.458
30.761
32.906
33.968
33.980
32.942
30.815
27.523
22.955
16.965
9.380
0.000


15.390
8.565
0.000

0.000
7.770
13.864
18.561
22.083
24.594
26.211
27.006
27.014
26.235
24.631
22.127
18.606
13.902
7.793
0.000

0.000
6.688
11.799
15.664
18.526
20.550
21.846
22.481
22.488

8.489
8.959
9.188
9.190
8.965
8.499
7.759
6.692
5.201
3.101
0.000

0.000
0.000
0.000
0.000
0.000
0.000
0.000
0.000
0.000
0.000
0.000
0.000
0.000
0.000
0.000
0.000

___________________________________________

Hình 1 :
Hình 2 :
left :
-0.65
left :
-0.75
bottom :
-0.65
bottom :
-0.75
width :
1.3
width :
1.5
heigh :
1.3
heigh :
1.5

___________________________________________
Hải Đăng – Điện 5 K54

*3*


Ta được như hình sau :

+ Trong cửa sổ Set Fomula, gõ lệnh: R2- R1 để xác định miền cần khảo sát trường
điện từ
+ Vào menu Options/Application/Electrostatics. Để thiết lập điều kiện bờ cho

sử dụng hình biểu diễn 3 chiều bằng cách tick vào ô Height (3-D plot). Hình sau thu
được khi chọn Height (3-D plot) và Deformed Mesh.

___________________________________________
Hải Đăng – Điện 5 K54

*6*