TRƯỜNG ĐẠI HỌC BÁCH KHOA HÀ NỘI PDF

Preview

Summary

TRƯỜNG ĐẠI HỌC BÁCH KHOA HÀ NỘIVIỆN ĐIỆN TỬ - VIỄN THÔNGSinh viên th c hi n: Lê Trầần Huyự ệMssv:L p:ớNhóm:Bài tập 1: a: x= 0 :0:1; y= 0 :0:1; [xx,yy]=meshgrid(x,y);m=size(xx); Ax= 2 *xx^2.*ones(m( 1 ),m( 2 )); Ay= 2 *ones(m( 1 ),m( 2 )); scf; champ(x,y,Ax,Ay,rect=[-0,-0, 2 , 2 ]); xlabel('x');ylabe...

Description

TRƯỜNG ĐẠI HỌC BÁCH KHOA HÀ NỘI VIỆN ĐIỆN TỬ - VIỄN THÔNG

Sinh viên th ực hi ện: Lê Trầần Huy Mssv:20182581 Lớp: Nhóm:

Bài tập 1:

a.Code: x=0:0.25:1.5; y=0:0.25:1.5; [xx,yy]=meshgrid(x,y);m=size(xx); Ax=2*xx^2.*ones(m(1),m(2)); Ay=2*ones(m(1),m(2)); scf; champ(x,y,Ax,Ay,rect=[-0.5,-0.5,2,2]); xlabel('x');ylabel('y'); title('truong vecto A=(2x^2x^+2y^)');

b.Code: x=0:0.25:3; y=0:0.25:3; [xx,yy]=meshgrid(x,y);m=size(xx); Bx=cos(xx.*yy); By=cos(yy^2); scf; champ(x,y,Bx,By,rect=[-0.5,-0.5,2,2]); xlabel('x');ylabel('y'); title('vecto B=cos(xy)x^+cos(y^2)y^');

Bài tập 2.1 a. Vẽẽ trường vô hướng và gradiẽnt trong không gian 3D sử dụng scilab A = cos(xy) + 2x + sin(xy) Với 0 < x,y < 3.5  Codẽ:

x=0:0.25:3.5; y=0:0.25:3.5; [xx,yy]=meshgrid(x,y); A1=cos(xx.*yy); A2=2.*xx; A3=sin(xx.*yy); A=A1+A2+A3; surf(x,y,A,'facecolor','interp'); xlabel('x'); ylabel('y'); zlabel('A'); title('3-D gradien'); G1=-yy.*sin(xx.*yy)+2+yy.*cos(xx.*yy); G2=-xx.*sin(xx.*yy)+ xx.*cos(xx.*yy); scf; champ(x,y,G1,G2,rect=[-0.5,-0.5,4,4]);

xlabel('x'); ylabel('y'); title('gradien'); Kêết quả:

b. Vẽẽ trường vẽcto và divẽrgẽncẽ trong không gian 2D sử dụng scilab B = ẽ (sin2x a – cos2x a ) -2y

x

y

Với -2 < x,y < 2  Codẽ:

x=-2:0.25:2; y=-2:0.25:2; [xx,yy]=meshgrid(x,y);m=size(xx); f=2.*yy; B1=sin(2.*xx); B2=cos(2.*xx); B=exp(-f)*(B1-B2); scf; champ(x,y,B1,B2,rect=[-0.5,-0.5,2,2]); xlabel('x'); ylabel('y'); title('truong vecto 2D'); D1=exp(-f)*2*cos(2.*xx); D2=exp(-f)*(-2)*(-cos(2.*xx)); D=D1+D2; plot3d(x,y,D); xlabel('x'); ylabel('y'); zlabel('D'); title('div cua truong vecto');  Kết quả :

Bài tập 3.1 a. Sử dụ ng scilab mô phỏ ng và quan sát điện trường tạo ra bở i điệ n tch điểm thẽo trục x,y dựa vào phương trình sau: E = 3y + 2x Với -2 < x,y < 2  Codẽ: 2

2

x=-2:0.25:2; y=-2:0.25:2; [xx,yy]=meshgrid(x,y); E1=3*yy.^2; E2=2*xx.^2; E=E1+E2; clf(); xset('colormap',rainbowcolormap(128));

colorbar(0,28); Sgrayplot(x,y,E,strf='081'); xlabel('x'); ylabel('y'); title('truong dien tu tao boi dien tich diem')  Kêết quả:

b. Dựa vào phương trình sau, mô phỏng và quan sát phần bôế điện thêế và cường độ điện trường sử dụng gradiẽnt và scilab. S ựkhác bi ệt gi ữa hai hình vẽẽ là gì và tại sao ? V = (x + xy ) 2

3

Với -3 < x < 3 và -1.5 < y < 1.5 Codẽ:

x=-3:0.1:3; y=-1.5:0.1:1.5; [xx,yy]=meshgrid(x,y);

V=(xx.^2+xx.*yy).^3; subplot(2,1,1); plot(x,V); ax1=gca(); ax1.grid=[0,0]; xlabel('x'); ylabel('V'); title('phan bo dien the su dung grad'); xx=-2.95:0.1:2.95; E=-dif(V)./dif(x); subplot(2,1,2); plot(xx,E); ax2=gca(); ax2.grid=[0,0]; xlabel('x'); ylabel('E'); title(' phan bo dien truong su dung grad'); kêết quả:

Phần bôế cường độ điện trường: Codẽ: x=-3:0.1:3; y=-1.5:0.1:1.5; [xx,yy]=mẽshgrid(x,y);

V=(xx.^2+xx.*yy).^3; Ex=3*(2*xx+yy).*(xx.^2+xx.*yy).^2; Ey=3*xx.*(xx.^2+xx.*yy).^2; E=Ex+Ey; surf(x,y,E,'facẽcolor','intẽrp');

Codẽ:Phần bôế điện thêế x=-3:0.1:3; y=-1.5:0.1:1.5; [xx,yy]=mẽshgrid(x,y); V=(xx.^2+xx.*yy).^3; surf(x,y,V,'facẽcolor','intẽrp'); kêết quả:...