MATLAB: Need help on difficulty 3D Plot. – Problem resolved

3d plotsMATLAB

I got below question. Need help on coding for below issue.

Best Answer

Here's some code to get you started

% Data
T=300;
k=1.38064852e-23;
Ea=7.2*10^(-21);
Ed=7.2*10^(-21);
Nv=2.5*10^(25)*((0.59*T/300)^(3/2))*10^(-6);
Nc=2.5*10^(25)*((1.08*T/300)^(3/2))*10^(-6);
Ec=0.5*(1.166-0.000473*T*T/(636+T))*1.6*10^(-19);
Ev=0;
% Collect data to be passed to function
data = [T, k , Ea, Ed, Nv, Nc, Ec, Ev];
% Dopant concentrations
Nd = 10^17;
Na = 10^5;
 
Ef0 = 10^-21;   % Initial guess at fermi energy
% Use fzero to find fermi energy, i.e. the value of Ef that makes 
% function Efn return zero
Ef = fzero(@Efn, Ef0,[],data,Nd,Na);  
disp(Ef)
function F = Efn(Ef,data, Nd, Na)
T = data(1);
k = data(2);
Ea = data(3);
Ed = data(4);
Nv = data(5);
Nc = data(6);
Ec = data(7);
Ev = data(8);
kT = k*T;
F = Nc*exp(-(Ec-Ef)/kT) + Na/(1+4*exp(-(Ef-Ea)/kT)) - Nv*exp(-(Ef-Ev)/kT) - Nd/(1+2*exp(-(Ed-Ef)/kT));
end

The above will calculate the Fermi level for one pair of dopant concentrations. See if you can take it from here.

Related Solutions

MATLAB: Mesh plot – Z must be a matrix – problem resolved

% Data
T=300;
k=1.38064852e-23;
Ea=7.2*10^(-21);
Ed=7.2*10^(-21);
Nv=2.5*10^(25)*((0.59*T/300)^(3/2))*10^(-6);
Nc=2.5*10^(25)*((1.08*T/300)^(3/2))*10^(-6);
Ec=0.5*(1.166-0.000473*T*T/(636+T))*1.6*10^(-19);
Ev=0;
% Scale energy levels
kT = k*T;
Ea = Ea/kT; Ed = Ed/kT; Ec = Ec/kT; Ev = Ev/kT;
% Collect data to be passed to function
data = [Ea, Ed, Nv, Nc, Ec, Ev];
% Dopant concentrations
Nd = [4E13, 6E13, 7E18, 1E17, 1E15, 1E15, 1E13, 1E15, 1E15, 2E16, 2.5E15];
Na = [1E15, 1E13, 1E15, 1E15, 1E17, 1E15, 8E14, 5E13, 6.E15, 6.4E13, 6.3E17];
n = numel(Nd) ; m = numel(Na) ; 
Ef = zeros(m,n); 
[Na,Nd] = meshgrid(Na,Nd) ; 
Ef0 = 10^(-21)/kT;   % Initial guess at scaled fermi energy
% Use fzero to find fermi energy, i.e. the value of Ef that makes 
% function Efn return zero
for i = 1:n
  for j = 1:m
    Ef(i,j) = fzero(@Efn, Ef0,[],data,Nd(i,j),Na(i,j));  
  end
end
Ef = Ef*kT;  % Rescale  
fprintf('%g\n',Ef)
pcolor(Na,Nd,Ef)
%patch(Na,Nd,Ef)
%plot3(Na,Nd,Ef)
grid on

MATLAB: How to load multiple files and do the computation for each of them and write the output in consequent rows

Your loop will want to start with grabbing all files in a folder. uigetfile with the MultiSelect option will work wonders here:

[Name,Path]=uigetfile({'*.xlsx*'},'MultiSelect','on');
entirefile=fullfile(Paths,Names);

Then start your loop. It wants to run a number of times equal to the number of elements in your entirefile cell, as you also have in your own code. Within this loop, do all your loading, calculation and writing. To avoid cluttering the workspace, you might want to look into a structure and dot indexing.

for i=1:numel(entirefile);
    %your entire code including
    enginedata(i).raw=importdata(entirefile{i});
    enginedata(i).temperature = your temperature calculation
    enginedata(i).oilpressure = oil pressure calculation
    %and so on
    %finally, to write in consecutive fields the results, simply keep indexing your things:
    celltarget=['B',num2str(i)]; %B1 first run, B2 2nd run etc. 
    writecell(mean_engine_temperature,filename,'Sheet','Temperatures','Range',celltarget);
end

Best Answer

Related Solutions

MATLAB: Mesh plot – Z must be a matrix – problem resolved

MATLAB: How to load multiple files and do the computation for each of them and write the output in consequent rows

Related Question