хочу построить график барьерной емкости в чем моя ошибка?

@Skillet102 · Регистрация: 17.10.2012

Студворк — интернет-сервис помощи студентам

привет всем
хочу построить график барьерной емкости
в чем моя ошибка?

Matlab M
1
2
3
4
5
6
7
8
9
e0=9e-14;
e=12;
l0=8e-8;
f0=0.82;
for U=0:0.1:1;
C=((e0*e/l0)*(sqrt(f0/U)))
end;
plot(C,U);
grid;

@R2D2 · 06.05.2013, 15:51

Здравствуйте! Ошибка в том, что Вы не накапливаете точки при построении графика. Так происходит потому, что plot стоит у Вас вне цикла и рисует он только точку рассчитанную на последней итерации цикла for. Да и при таком количестве точек (у Вас их 11) сплошного графика Вы не получите. Нужно или увеличивать количество точек или отображать точки маркерами. Ну и напоследок стоит заметить, что мисье тонкий извращенец, раз любит считать такие вещи в цикле

. В общем я тут подправил кое что и свой вариант предложил.

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
%% Подправленый вариант
clear; close all; clc;
e0=9e-14;
e=12;
l0=8e-8;
f0=0.82;
k = 1;
for U=0:0.1:1;
C(k)=((e0*e/l0)*(sqrt(f0/U)));
plot(C(k),U, 'or'); hold on; grid on;
k = k+1;
end
%% Мой вариант
clear; close all; clc;
e0=9e-14; e=12; l0=8e-8; f0=0.82; U = 0:.1:1;
C = (e0*e/l0)*sqrt(f0./U);
plot(C, U, '-or'); grid on;

@Skillet102 · 06.05.2013, 18:51 **[ТС]**

спасибо,ошибку понял)
толк знаю

Добавлено через 36 минут

Matlab M
1
2
3
4
5
%Диффузионая емкость
Ie=1e-3;w=6e-8;q=1.6e-19;Nd=4e20;D=1e11; U = 0:.1:1;
Cd=((sqrt((e0*e)/(2*q*Nd)))*((Ie*w)/(D.*sqrt(U))))
figure(2);
plot(Cd,U,'-');grid on;

выдает ошибку
Matrix dimensions must agree.

определено же все
или нужно как-то иначе переопределить?
я отошел от циклов

Добавлено через 8 минут
хух,исправил
во всем виновата лишняя скоба)

@R2D2 · 06.05.2013, 18:58

Ну тут все дело в волшебных точках

.
Вот здесь:

Matlab M
1
D.*sqrt(U)

ее ставить не обязательно, т.к. в - это скаляр и матлаб воспримет запись без точки как умножение вектора sqrt(U) на вектор такой же длины, в котором все элементы равны D.
А вот при делении точку надо бы поставить:
Cd=((sqrt((e0*e)/(2*q*Nd)))*((Ie*w)./(D*sqrt(U))))

@Skillet102 · 06.05.2013, 22:03 **[ТС]**

упс
тяга у меня к ошибкам и недоделкам

Matlab M
1
2
3
4
%Распределение концентраций
Nis=5e16;X=2e-6;t=0:20:144;
n1=Nd*exp(-X^2./4*D*t);n2=Nis;
figure(4);semilogy(n1,'g',n2,'b');grid on;

хочу построить график распределения концентраций,но выдает ошибку
Not enough input arguments.

Добавлено через 1 час 8 минут
вот тут тоже строится распределение концентраций,но в моем случает отсутвует бор,а значит что Na у меня нет,как и Xk,Sk & DB
исходные данные здесь другие)

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
% Исходные данные
Nis=1e16;
Na=1e21;
Nd=4e18;
Xe=1.5e-4;
Xk=4e-4;
D=3.5e-13;
Se=50e-6;
Sk=100e-6;
q = 1.6e-19;
 
fprintf('Концентрации, в см^-3\n Np = %d\n Nn = %d\n No = %d\n',Na,Nd,Nis);
fprintf('Глубина залегания p-n переходов, в см\n Xa = %d\n Xd = %d\n',Xe,Xk);
DP = 3e-13;
DB = 1.3e-13;
fprintf('Коэффициенты диффузии фосфора и бора при T=1400 К, в см/с\n D_P = %d\n D_B = %d\n',DP,DB);
 
% Расчет времени диффузии фосфора
r  = 2e27;
td = 0;
  while (r>1e+013)
    td = td+1;
    N2 = Nd*exp(-(Xk^2)/(4*DP*td));
    r = abs(Nis-N2);
  end
fprintf('\nВремя диффузии фосфора, в с');
fprintf('\n td = %6.2f с',td);
td1 = fix(td/3600);
td2 = fix(rem(td,3600)/60);
td3 = fix(rem(rem(td,3600),60));
fprintf('  = %6.2f ч %6.2f мин %6.2f с',td1,td2,td3);
 
% Расчет времени диффузии бора 
r  = 2e27;
ta = 0;
    while (r>4.7e+015)
     ta = ta+1;
     N1 = Na*exp(-(Xe^2)/(4*DB*ta));
     N2 = Nd*exp(-(Xe^2)/(4*DP*td));
     r  = abs(N1-N2);        
   end
fprintf('\nВремя диффузии бора, в с');
fprintf('\n ta = %6.2f с',ta);
ta1 = fix(ta/3600);
ta2 = fix(rem(ta,3600)/60);
ta3 = fix(rem(rem(ta,3600),60));
fprintf('  = %6.2f ч %6.2f мин %6.2f с\n',ta1,ta2,ta3);
 
% Расчет концентраций примесей
i = 1; p1=0; n2=0;
 
   for x = 0:2e-7:6e-4;
     N1 = Na*exp(-(x^2)/(4*DB*ta));
     N2 = Nd*exp(-(x^2)/(4*DP*td));
     N3 = Nis;
 
     if (x <= Xe)
       N = N1-N2+N3;
       p1=p1+N;  % *0.1e-7
       if (x<Xe)
         N_ae=N/(Xe-x)*1e6 ;      
       end
     end
     if (x > Xe & x <= Xk)
         N = N2-N3-N1;
    if (x<=3.8e-4)
           n2=n2+N; %*0.1e-7   
           N_de=N/(x-Xe)*1e6;   
         end
         if (x>3.8e-4 & x<Xk)
           N_dk=N/(Xk-x)*1e6;  
         end
     end
 
      if (x > Xk)
         N = N1-N2+N3;
         N_ak=N/(x-Xk)*1e6;  
      end
 
  xx(i)  = x;
  yy1(i) = N1;
  yy2(i) = N2;
  yy3(i) = N3;
  yy(i)  = N;
  i = i+1;
 end
 
figure(1); title('Концентрации примесей в полупроводнике.'); 
semilogy(xx,yy1,'r',xx,yy2,'g',xx,yy3,'b'); grid;
figure(2); title('Результирующая концентрация примесей\n в полупроводнике.'); 
semilogy(xx,yy,'g'); grid;

Добавлено через 53 минуты
фигня какая-то получается

@R2D2 · 06.05.2013, 22:57

В последней стоке пишите так:

Matlab M
1
semilogy(n1, n2,'ob')

Вот только по оси Y у Вас ничего не меняется... Впрочем Вам решать, правильно ли это...

@Skillet102 · 07.05.2013, 08:09 **[ТС]**

спасибо,буду думать)

Зосима · 07.05.2013, 11:51

Может ты хотел так:

Matlab M
1
semilogy(t,n1,'g', t,n2,'b')

@Skillet102 · 07.05.2013, 12:49 **[ТС]**

я хочу как тут,зависит от х-расстояния
тольок вот не знаю как решить проблемку с отсутствием бора

Добавлено через 10 минут

Matlab M
1
2
3
4
5
%Распределение концентраций
Nis=5e16;t=0:20:144;x = 0:2e-7:6e-4;
n1=Nd*exp(-x^2./4*D*t);n2=Nis;
figure(4);%semilogy(n1,n2,'-ob');
semilogy(n1,'g', n2,'b');grid on;

попробовал так,выдает ошибку
Matrix must be square.

Зосима · 07.05.2013, 12:55

Сообщение от Skillet102

выдает ошибку

Забыл точку в n1=Nd*exp(-x.^2./(4*D*t) );
И еще у тебя t и x разных размерностей...

что ты пытаешься изобрести?

@Skillet102 · 07.05.2013, 13:35 **[ТС]**

t это время,как его изменишь?

@Skillet102 · 07.05.2013, 14:34 **[ТС]**

хочу построить график барьерной емкости в чем моя ошибка?

по идее у меня должен получиться вот такой график,но только без красной линии и с моими данными
код у этого графика такой

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
format long;
 
% Исходные данные
Nis=1e16;
Na=1e21;
Nd=4e18;
Xe=1.5e-4;
Xk=4e-4;
D=3.5e-13;
Se=50e-6;
Sk=100e-6;
q = 1.6e-19;
 
fprintf('Концентрации, в см^-3\n Np = %d\n Nn = %d\n No = %d\n',Na,Nd,Nis);
fprintf('Глубина залегания p-n переходов, в см\n Xa = %d\n Xd = %d\n',Xe,Xk);
DP = 3e-13;
DB = 1.3e-13;
fprintf('Коэффициенты диффузии фосфора и бора при T=1400 К, в см/с\n D_P = %d\n D_B = %d\n',DP,DB);
 
% Расчет времени диффузии фосфора
r  = 2e27;
td = 0;
  while (r>1e+013)
    td = td+1;
    N2 = Nd*exp(-(Xk^2)/(4*DP*td));
    r = abs(Nis-N2);
  end
fprintf('\nВремя диффузии фосфора, в с');
fprintf('\n td = %6.2f с',td);
td1 = fix(td/3600);
td2 = fix(rem(td,3600)/60);
td3 = fix(rem(rem(td,3600),60));
fprintf('  = %6.2f ч %6.2f мин %6.2f с',td1,td2,td3);
 
% Расчет времени диффузии бора 
r  = 2e27;
ta = 0;
    while (r>4.7e+015)
     ta = ta+1;
     N1 = Na*exp(-(Xe^2)/(4*DB*ta));
     N2 = Nd*exp(-(Xe^2)/(4*DP*td));
     r  = abs(N1-N2);        
   end
fprintf('\nВремя диффузии бора, в с');
fprintf('\n ta = %6.2f с',ta);
ta1 = fix(ta/3600);
ta2 = fix(rem(ta,3600)/60);
ta3 = fix(rem(rem(ta,3600),60));
fprintf('  = %6.2f ч %6.2f мин %6.2f с\n',ta1,ta2,ta3);
 
% Расчет концентраций примесей
i = 1; p1=0; n2=0;
 
   for x = 0:2e-7:6e-4;
     N1 = Na*exp(-(x^2)/(4*DB*ta));
     N2 = Nd*exp(-(x^2)/(4*DP*td));
     N3 = Nis;
 
     if (x <= Xe)
       N = N1-N2+N3;
       p1=p1+N;  % *0.1e-7
       if (x<Xe)
         N_ae=N/(Xe-x)*1e6 ;      
       end
     end
     if (x > Xe & x <= Xk)
         N = N2-N3-N1;
    if (x<=3.8e-4)
           n2=n2+N; %*0.1e-7   
           N_de=N/(x-Xe)*1e6;   
         end
         if (x>3.8e-4 & x<Xk)
           N_dk=N/(Xk-x)*1e6;  
         end
     end
 
      if (x > Xk)
         N = N1-N2+N3;
         N_ak=N/(x-Xk)*1e6;  
      end
 
  xx(i)  = x;
  yy1(i) = N1;
  yy2(i) = N2;
  yy3(i) = N3;
  yy(i)  = N;
  i = i+1;
 end
 
figure(1); title('Концентрации примесей в полупроводнике.'); 
semilogy(xx,yy1,'r',xx,yy2,'g',xx,yy3,'b'); grid;
figure(2); title('Результирующая концентрация примесей\n в полупроводнике.'); 
semilogy(xx,yy,'g'); grid;

@Skillet102 · 07.05.2013, 14:36 **[ТС]**

я смог получить вот такой график

с таким кодом

Matlab M
1
2
3
4
5
%Распределение концентраций
Nis=5e16;t=144;x = 0:2e-7:6e-4;
n1=Nd*exp(-x.^2./4*D*t);n2=Nis;
figure(4);%semilogy(n1,n2,'-ob');
semilogy(x,n1,'g', x,n2,'b');grid on;

Новые блоги и статьи Все статьи Все блоги /
Символические и жёсткие ссылки в Linux. algri14 15.03.2026 Существует два типа ссылок — символические и жёсткие. Ссылка в Linux — это дополнительная запись в каталоге, которая может указывать либо на inode «файла-ИСТОЧНИКА», тогда это будет «жёсткая. . .	[Owen Logic] Поддержание уровня воды в резервуаре количеством включённых насосов: моделирование и выбор регулятора ФедосеевПавел 14.03.2026 Поддержание уровня воды в резервуаре количеством включённых насосов: моделирование и выбор регулятора ВВЕДЕНИЕ Выполняя задание на управление насосной группой заполнения резервуара,. . .	делаю науч статью по влиянию грибов на сукцессию anaschu 13.03.2026 прикрепляю статью	SDL3 для Desktop (MinGW): Создаём пустое окно с нуля для 2D-графики на SDL3, Си и C++ 8Observer8 10.03.2026 Содержание блога Финальные проекты на Си и на C++: hello-sdl3-c. zip hello-sdl3-cpp. zip Результат:
Установка CMake и MinGW 13.1 для сборки С и C++ приложений из консоли и из Qt Creator в EXE 8Observer8 10.03.2026 Содержание блога MinGW - это коллекция инструментов для сборки приложений в EXE. CMake - это система сборки приложений. Здесь описаны базовые шаги для старта программирования с помощью CMake и. . .	Как дизайн сайта влияет на конверсию: 7 решений, которые реально повышают заявки Neotwalker 08.03.2026 Многие до сих пор воспринимают дизайн сайта как “красивую оболочку”. На практике всё иначе: дизайн напрямую влияет на то, оставит человек заявку или уйдёт через несколько секунд. Даже если у вас. . .	Модульная разработка через nuget packages DevAlt 07.03.2026 Сложившийся в . Net-среде способ разработки чаще всего предполагает монорепозиторий в котором находятся все исходники. При создании нового решения, мы просто добавляем нужные проекты и имеем. . .	Модульный подход на примере F# DevAlt 06.03.2026 В блоге дяди Боба наткнулся на такое определение: В этой книге («Подход, основанный на вариантах использования») Ивар утверждает, что архитектура программного обеспечения — это структуры,. . .

@Skillet102 2 / 2 / 0 Регистрация: 17.10.2012 Сообщений: 88
	07.05.2013, 08:09 [ТС]
	спасибо,буду думать) 0

@Skillet102 2 / 2 / 0 Регистрация: 17.10.2012 Сообщений: 88
	07.05.2013, 13:35 [ТС]
	t это время,как его изменишь? 0