Can anyone help me in getting this graph using the attached equation (16).

I want to plot the graph between mobility vs temp. for different donor concentration using eq. 16 . I am not getting how i'll write the code for this

[Walter Roberson]

Define the constants. Then define N1 and . Then define beta_cw__2 as an anonymous function of T . You can then use that to define an anonymous function that returns . Now with all the parts in place, you can create a vector of actual temperatures and pass it into the function to calculate the list of outputs; you can the plot() as appropriate.

[Nudrat Sufiyan]

I have tried to write code. But I am not getting how i'll define Nd and 1/mu

clear all;
a=2.61e-4;
b=2.90e-4;
c=1.70e-2;
kc=0.3;
teta=1065;
Nd=[10e16:10e17:10e18];
N1=(1+kc).*Nd;

[Walter Roberson]

Nd = [10e16, 10e17, 10e18];

and when you define the vector of temperatures, define it as a column vector.

[Nudrat Sufiyan]

Sir,

How i'll define vector of temperature as column vector as value of T we can take between 300k to 600k. Please help me in this

[Walter Roberson]

num_temperatures = 50;
t = linspace(300, 600, num_temperatures) .';

Answer 1

[Sam Chak]

Hi @Nudrat Sufiyan

Guess you need to type out this equation

and insert the parameters

Example only, not true equation. One of the simpler way without using loop.

% code that you have typed so far
T     = linspace(300, 600, 30001);
a     = 2.61e-4;
b     = 2.90e-4;
c     = 1.70e-2;
kc    = 0.3;
theta = 1065;
Nd    = 10e16;
Ni    = (1 + kc)*Nd;

% additional stuffs that you should type
beta  = 1;       % a bit lazy to type out, trying doing it yourself
iMu   = a*(Ni/1e17)*log(1 + beta.^2).*(T/300).^(-1.5) + b*(T/300).^(1.5) + c./(exp(theta./T) - 1);
plot(T, 1./iMu), hold on

% copy/paste and modify the value of Nd
Nd    = 10e17;
iMu   = a*(Ni/1e17)*log(1 + beta.^2).*(T/300).^(-1.5) + b*(T/300).^(1.5) + c./(exp(theta./T) - 1);
plot(T, 1./iMu)

Nd    = 10e18;
iMu   = a*(Ni/1e17)*log(1 + beta.^2).*(T/300).^(-1.5) + b*(T/300).^(1.5) + c./(exp(theta./T) - 1);
plot(T, 1./iMu), hold off, grid on, xlabel('T, [K]'), ylabel('\mu, [put the unit here]')

[Nudrat Sufiyan]

yes, I need help in this equation

[Sam Chak]

@Nudrat Sufiyan, See example above. Type out the equation correctly.

[Nudrat Sufiyan]

Thank you so much sir, I am really grateful to you for this and now i'll try to write code according to the example you have provided. Thanks alot sir.

[Sam Chak]

You are welcome, @Nudrat Sufiyan. If you find the example and the MATLAB code are helpful, please consider accepting ✔ and voting 👍 the Answer. Thanks!

[Nudrat Sufiyan]

sure sir. Thankyou