Если у вас возникли сложности с курсовой, контрольной, дипломной, рефератом, отчетом по практике, научно-исследовательской и любой другой работой - мы готовы помочь.
Предоплата всего
Подписываем
ЗАВДАННЯ № 1.
1. Розрахувати розподіл потенціалу в міжелектродному просторі польового транзистора
(варіант конструкції вибирається по передостанній цифрі номера залікової книжки) з точністю до
0,01 В. Номер варіанту вибирається по останній цифрі номера залікової книжки.
2. Побудувати картину поля з допомогою еквіпотенціалей і векторів напруженості
електричного поля. Побудувати косокутну проекцію потенціального рельєфу U(х,у) = –eV(х,у).
|
l1 = l2 = l3/2 = l4 = l5 = 1,0 мкм, d = 1,0 мкм, D = 3,0 мкм 1 – витік, 2 – заслін, 3 – стік –Vзи=0,3 B Vси=2 В |
КОД:
close all;clear all;clc;
m=90;
n=45;
hx=6/m;
hy=3/n;
sv=2
zv=-0.3
eps=0.01
v=zeros(m,n); V=0;
k=1;
e=-1.6e-19;
while k<1000
V=v;
k=k+1;
v(1:1:m,1)=sv;
v(1:1:m,n)=0;
v(m/6:1:2*m/6,2*n/3)=zv;
v(4*m/6:1:5*m/6,2*n/3)=zv;
v(1,1)=(v(1,2)+v(2,1))/2;
v(1,n)=(v(1,n-1)+v(2,n))/2;
v(m,1)=(v(m-1,1)+v(m,2))/2;
v(m,n)=(v(m-1,n)+v(m,n-1))/2;
for i=2:1:(m-1)
for j=2:1:(2*n/3-1)
v(i,j)=0.25*(v(i-1,j)+v(i+1,j)+v(i,j+1)+v(i,j-1));
end
end
j=(2*n/3);
for i=2:1:(m/6-1)
v(i,j)=0.25*(v(i-1,j)+v(i+1,j)+v(i,j+1)+v(i,j-1));
end
for i=(2*m/6+1):1:(4*m/6-1)
v(i,j)=0.25*(v(i-1,j)+v(i+1,j)+v(i,j+1)+v(i,j-1));
end
for i=(5*m/6+1):1:(m-1)
v(i,j)=(v(i-1,j)+v(i+1,j)+v(i,j+1)+v(i,j-1))/4;
end
for i=2:1:(m-1)
for j=(2*n/3+1):1:(n-1)
v(i,j)=0.25*(v(i-1,j)+v(i+1,j)+v(i,j+1)+v(i,j-1));
end
end
for i=2:1:m-1
v(i,n)=(v(i-1,n)+v(i+1,n)+v(i,n-1))/3;
end
for j=2:1:(n-1)
v(1,j)=(v(1,j-1)+v(1,j+1)+v(2,j))/3;
v(m,j)=(v(m,j-1)+v(m,j+1)+v(m-1,j))/3;
end
end
if min(abs(v-V))<eps
%Побудова сигналів
y1=1*hy:1*hy:n*hy;
x1=1*hx:1*hx:m*hx;
figure; mesh(x1,y1,v'/e)
xlabel('x, mkm'); ylabel('y, mkm'); zlabel('U, V');
title('Косокутнa проекція потенціального рельєфу');
Ex=zeros(m,n);
Ey=zeros(m,n);
for i=1:1:m-1
for j=1:1:n
Ex(i,j)=-(v(i+1,j)-v(i,j))/hx;
end
end
for i=1:1:m
for j=1:1:n-1
Ey(i,j)=-(v(i,j+1)-v(i,j))/hy;
end
end
figure
[X,Y] = meshgrid(1*hx:1*hx:m*hx,1*hy:1*hy:n*hy);
hold on
quiver(X,Y,Ex',Ey')
colormap hsv
xlabel('x, mkm'); ylabel('y, mkm'); zlabel('U, V');
title('Картина поля з векторами напруженості електричного поля')
hold off
figure('Color','w')
[C,h]=contour(x1,y1,v',25);
clabel(C,h);
xlabel('x, mkm'); ylabel('y, mkm'); zlabel('U, V');
title('Еквіпотенціальна картина поля')
break
end
end
Висновки: підчас виконання цієї роботи було виявлено що поле «стікає» в точки затвору. Також при збільшенні кількості ітерацій точність підвищується і як результат спад потенціалу стає більш пологим.
Навколо точок затвору видно замкнуті еквіпотенціальні лінії. Еквіпотенціальні лінії спочатку паралельні осі х але з наближенням до затвору вони викривляються і між точками затвору більше випирають по осі у. Також видно що спад до затвору більш різкий зі сторони стоку, де більші значення потенціалу поля.
З картини поля побудованої за допомогою векторів напруженості електричного поля замітно що поле направлина від стоку до витоку а також що вектори в районі затвору поступово міняють напрямок в його сторону.
Міністерство освіти і науки України
Національний технічний університет України
“Київський політехнічний інститут”
Факультет електроніки
Кафедра фізичної та біомедичної електроніки
Розрахунково-графічна робота
з дисципліни «Теорія поля»
Виконав: ст. 3 курсу
гр.. ДП-92
Сабов Т. М.
Перевірив:
ас. Федяй А. В.
Київ – 2011