Bài tập matlab - Chương trình mô phỏng chuyển động của electron
Bạn đang xem bản rút gọn của tài liệu. Xem và tải ngay bản đầy đủ của tài liệu tại đây (25.71 KB, 2 trang )
clc;
fprintf('Chuong trinh mo phong chuyen dong cua electron\n')
dentat=10^(-10);
q=1.6*10^(-19);
me=9.109*10^(-31);
%Nhap chieu dai ong
d=input('Moi ban nhap vao chieu dai cua ong (cm)\n');
while d>200
fprintf('Do dai ban nhap qua lon!\n')
fprintf('Xin hay nhap chieu dai trong khoang tu 100cm den 200cm\n')
d=input('Moi ban nhap vao chieu dai cua ong (cm)\n');
end
while d<100
fprintf('Do dai ban nhap qua nho!\n')
fprintf('Xin hay nhap chieu dai trong khoang tu 100cm den 200cm\n')
d=input('Moi ban nhap vao chieu dai cua ong (cm)\n');
end
d=d/100;%doi don vi ra met
%Nhap ban kinh ong
r=input('Moi ban nhap vao ban kinh cua ong (cm)\n');
while r>20
fprintf('Ban kinh ban nhap qua lon!\n')
fprintf('Xin hay nhap ban kinh trong khoang tu 10cm den 20cm\n')
r=input('Moi ban nhap vao ban kinh cua ong (cm)\n');
end
while r<10
fprintf('Ban kinh ban nhap qua nho!\n')
fprintf('Xin hay nhap ban kinh trong khoang tu 10cm den 20cm\n')
r=input('Moi ban nhap vao ban kinh cua ong (cm)\n');
end
r=r/100;%doi don vi ra met
V=input('Moi ban nhap vao dien the cua ong (v)\n');
while V>500
fprintf('Dien the ban nhap qua lon!\n')
fprintf('Xin hay nhap dien the trong khoang tu 100v den 500v\n')
V=input('Moi ban nhap vao dien the cua ong (v)\n');
end
while V<100
fprintf('Dien the ban nhap qua nho!\n')
fprintf('Xin hay nhap dien the trong khoang tu 100v den 500v\n')
V=input('Moi ban nhap vao dien the cua ong (v)\n');
end
E=V/d;
F=q*E;
x(1)=0;
vantoc(1)=0;
chieudai(1)=0;
n=2;
tam=chieudai(1);
lap=0;
while lap<=10
y=randint(1,1,[0,2*r*100])/100;
y1=randint(1,1,[0,2*r*100])/100;
y2=randint(1,1,[0,2*r*100])/100;
while tam<=d
vantoc(n)=(dentat*F)/me+vantoc(n-1);
x(n)=dentat*vantoc(n)+x(n-1);
chieudai(n)=chieudai(n-1)+x(n);
tam=chieudai(n);
plot(chieudai(n),y,'r .')
plot(chieudai(n),y1,'b .')
plot(chieudai(n),y2,'g .')
axis([0 d 0 2*r])
title('Chuyen dong cua dien tu trong ong co dien truong')
xlabel('chieu dai (m)'),ylabel('duong kinh ong (m)')
pause(0.001)
n=n+1;
end
lap=lap+1;
x(1)=0;
vantoc(1)=0;
chieudai(1)=0;
n=2;
tam=chieudai(1);
end
