Matlab bilenlerin fikrini almak isterim

[zeusmaster]

Matlab'da çalıştırdığımız bir kod bulunmakta.Subrutinlerin bazılarında değişiklik yapmak istiyorum.Nedeni ise yaptığımız hesaplamaların şuan tek malzeme içn yapılabiliyor olması.Material subrutini içinde değişlik yaptığım taktirde malzeme farklı bir malzeme olacak.Ben bazı değişiklikler yapıyorum ama varsa matlab bilgisi iyi olan onunda fikrini almak isterim.
material.m;

Kod Seç Genişlet

function a=gaasmaterial(x,prop)

%GAASMATERIAL material database for GaAs/AlGaAs mertial system
%
%a=gaasmaterial(x,prop)
%
%takes the x-content in field x and returns a material property
%according to the string prop.
%This database is somewhat rudimentary and incomplete and does not always
%give correct results. It should be ok for computations with Gamma electrons.
%X- and L-electrons are probably not described very well.
%Valid property strings are:
%
%all energies in eV, all masses in M0
%
%E_G6G8    energy between Gamma6 and Gamma8 point = normal gap
%E_G6G7   Gamma6-Gamma7 = split-off
%E_G8X6    Gamma8-X6 = indirect gap for x>.45
%E_G8L6    Gamma8-L6 = indirect gap f x>.45     
%
%M_eG      electron mass at Gamma point = DOS-mass = cond masse = t = l
%M_eXd     electron dos mass at X point
%M_eXl     electron longitudinal X mass
%M_eXt     electron transversal X mass
%M_eX      electron X mass
%M_eLd     electron dos mass at L point
%M_eLl     electron longitudinal L mass
%M_eLt     electron transversal L mass
%M_eL      electron L mass
%
%M_hh001   heavy hole 001 mass
%M_hh110   heavy hole 110 mass
%M_hh111   heavy hole 111 mass
%M_lh001   light hole 001 mass
%M_lh110   light hole 110 mass
%M_lh111   light hole 111 mass
%M_hhd     heavy hole DOS mass
%M_lhd     light hole DOS mass
%M_so      split-off hole mass
%
%Gamma1    Luttinger parameter
%Gamma2    Luttinger parameter
%Gamma3    Luttinger parameter
%
%The following energies give absolute values of the corresponding quantities.
%These energies are normalized such, that the GaAs CB edge at T=0 has energy 0.
%E_c       conduction band edge = min(E_g,E_x,E_l)
%E_v       valence band edge
%E_x       X band edge
%E_l       L band edge
%E_g       Gamma band edge
%E_so      split-off hole band edge
%E_hh      heavy hole band edge
%E_lh      light hole band edge
%E_d       donator level
%Delta_Ev  valence band offset
%
%eps       relative dielectric constant

%Copyright 1999 Martin Rother
%
%This file is part of AQUILA.
%
%AQUILA is free software; you can redistribute it and/or modify
%it under the terms of the BSD License as published by
%the Open Source Initiative according to the License Policy
%on MATLAB(R)CENTRAL.
%
%AQUILA is distributed in the hope that it will be useful,
%but WITHOUT ANY WARRANTY; without even the implied warranty of
%MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
%BSD License for more details.

global aquila_control aquila_material

%check for correct calling sequence
if bitget(aquila_control.progress_check,1)==0
   error('material: INITAQUILA must be run before accessing the materials database !')
end

a=zeros(size(x));
p=find(x>0.69);
if ~isempty(p)
   disp('material: x>0.69 not supported');
   return;
end

%some values are taken from tables
%the interpolation is computed only once and then used on subsequent calls.
if ~exist('aquila_material.alpha') %initialize values for interpolation, only done once
   aquila_material.e0=interp1([0 0.18 0.27 0.53 0.69],[0.505 0.590 0.644 0.750 0.828],[0:0.01:0.69]);
   aquila_material.alpha=interp1([0 0.18 0.27 0.53 0.69],[5.5 6.3 6.58 7.04 7.88]*1e-4,[0:0.01:0.69]);
   aquila_material.beta=interp1([0 0.18 0.27 0.53 0.69],[0.225 0.236 0.248 0.261 0.302],[0:0.01:0.69]);
end
T=aquila_control.T;

%return the requested material property
switch prop     
case {'E_G6G8'} % Gamma6-Gamma8 = normal gap
   index=find(x>=0);
   x2=x(index);
   s=size(x2);
   x2=x2(:);
   e0=interp1([0:0.01:0.69],aquila_material.e0,x2,'nearest');
   alpha=interp1([0:0.01:0.69],aquila_material.alpha,x2,'nearest');
   beta=interp1([0:0.01:0.69],aquila_material.beta,x2,'nearest');
   a2=e0-T*T*alpha./(beta+T);
   size(a2);
   a(index)=a2;
   index=find(x<0); %formulas for InGaAs, implemented extremely rudimentary
   x2=x(index);
   s=size(x2);
   x2=-x2(:);
   a2=0.505-5.5e-4*T^2/(T+200)+(0.42-2.76e-4*T^2/(T+83)-0.505+5.5e-4*T^2/(T+200))*x2-0.475*x2.*(1-x2);
   a(index)=a2;
case {'E_G6G7'} % G6-G7 = split-off
   a=1.786+1.028*x+0.407*x.*x;
case {'E_G8X6'} % G8-X6 = ind Gap for x>.45
   a=1.900+1.25*x+0.143*x.*x;
case {'E_G8L6'} % G8-L6 = ind Gap for x>.45     
   a=1.708+1.642*x;
   
case {'M_eG'} % El-G-Mass (DOS-Mass=cond-Mass=t=l)
   index=find(x>=0);
   a(index)=0.067+0.083*x(index);
   index=find(x<0);
   a(index)=0.023*ones(size(x(index)));
case {'M_eXd'} % El-X-Mass (DOS)
   a=(0.85-x*0.14);
case {'M_eXl'} % El-X-Mass long
   a=(1.3-x*(1.3-1.1));
case {'M_eXt'} % El-X-Mass trans
   a=(0.23-x.*(0.23-0.19));
case {'M_eX'} % El-X-Mass
   ext=gaasmaterial(x,'M_eXt');
   exl=gaasmaterial(x,'M_eXl');
   a=3.0*(ext.*exl./(ext+2.0*exl));
   
case {'M_eLd'} % El-L-Mass (DOS)
   a=(0.56+x*0.22);
case {'M_eLl'} % El-L-Mass long
   a=(1.9-x*(1.9-1.32));
case {'M_eLt'} % El-L-Mass trans
   a=(0.0754+x*(0.15-0.0754));
case {'M_eL'} % El-L-Mass
   elt=gaasmaterial(x,'M_eLt');
   ell=gaasmaterial(x,'M_eLl');
   a=3.0*(elt.*ell./(elt+2.0*ell));
   
case {'M_hh001'} % HH-Mass 001
   g1=gaasmaterial(x,'Gamma1');
   g2=gaasmaterial(x,'Gamma2');
   a=1./(g1-2*g2);
case {'M_hh110'} % HH-Mass 110
   g1=gaasmaterial(x,'Gamma1');
   g2=gaasmaterial(x,'Gamma2');
   g3=gaasmaterial(x,'Gamma3');
   a=2./(2*g1-g2-3*g3);
case {'M_hh111'} % HH-Mass 111
   g1=gaasmaterial(x,'Gamma1');
   g3=gaasmaterial(x,'Gamma3');
   a=1./(g1-2*g3);
case {'M_lh001'} % LH-Mass 001
   g1=gaasmaterial(x,'Gamma1');
   g2=gaasmaterial(x,'Gamma2');
   a=1./(g1+2*g2);
case {'M_lh111'} % LH-Mass 111
   g1=gaasmaterial(x,'Gamma1');
   g3=gaasmaterial(x,'Gamma3');
   a=1./(g1+2*g3);
case {'M_lh110'} % LH-Mass 110
   g1=gaasmaterial(x,'Gamma1');
   g2=gaasmaterial(x,'Gamma2');
   g3=gaasmaterial(x,'Gamma3');
   a=2./(2*g1+g2+3*g3);
case {'M_hhd'} % HH-DOS-Mass
   a=(0.55+0.26*x);
case {'M_lhd'} % LH-DOS-Mass
   a=(0.08+0.08*x);
case {'M_so'} % SO-Mass
   a=(0.165+0.135*x);     
   
case {'Gamma1'} % Lutt Param
   a=1.0./((1.0-x)./7.10+x./3.76);
case {'Gamma2'} % Lutt Param
   a=1.0./((1.0-x)./2.02+x./0.90);
case {'Gamma3'} % Lutt Param
   a=1.0./((1.0-x)./2.91+x./1.42);
   
case {'E_c'} % CB-Edge
   p=find(x<0.45);
   eg=gaasmaterial(x,'E_g');
   a(p)=eg(p);
   p=find(x>=0.45);
   ex=gaasmaterial(x,'E_x');
   a(p)=ex(p);
case {'E_v'} % VB-Edge
   a=gaasmaterial(x,'E_g')-gaasmaterial(x,'E_G6G8');
case {'E_x'} % X-Point
   a=gaasmaterial(x,'E_v')+gaasmaterial(x,'E_G8X6');
case {'E_l'}    % L-Point
   a=gaasmaterial(x,'E_v')+gaasmaterial(x,'E_G8L6');
case {'E_g'} % G-Point
   a=gaasmaterial(x,'E_G6G8')-0.505-gaasmaterial(x,'Delta_Ev'); % therefore Eg=0 for x=0
case {'E_so'} % SO-Edge
   a=gaasmaterial(x,'E_g')-gaasmaterial(x,'E_G6G7');
   
case {'Delta_Ev'} % VB-Offset
   index=find(x>=0);
   a(index)=0.51*x(index);
   index=find(x<0);
   p=1.5+(4.88-1.5)/.9*(-x(index)-0.1);
   a(index)=-(gaasmaterial(zeros(size(x(index))),'E_G6G8')-gaasmaterial(x(index),'E_G6G8'))./(1+p);     
case {'eps'} % Epsilon
   index=find(x>=0);
   a(index)=13.18-3.12*x(index);
   index=find(x<0);
   a(index)=12;
   
case {'E_d'} % Donatorlevel
   ec=gaasmaterial(x,'E_c');
   p=find(x<0.22);
   a(p)=ec(p)-.006;
   p=find((x>=0.22)&(x<0.35));
   a(p)=ec(p)-(0.006+(x(p)-0.22)/(.35-.22)*(0.019-0.006));
   p=find(x>=0.35);   
   a(p)=ec(p)-0.019;
otherwise
   error('material: property undefined!');
end

[hyperThread]

 Burada bir hesaplama sistemin var, material.m fonksiyonu da bir malzeme için çalışıyor. Bunun üzerinde değişiklik yapcan ama o değişikliğin neler olmasını istediğini de yazarsan daha sağlıklı bir biçimde yardımcı olabilirim

[zeusmaster]

Kod Seç Genişlet

aquila_material.e0=interp1([0 0.18 0.27 0.53 0.69],[0.505 0.590 0.644 0.750 0.828],[0:0.01:0.69]);

Burada ki [0 0.18 0.27 0.53 0.69] malzemenin içindeki Al oranı,
[0.505 0.590 0.644 0.750 0.828] ise her x değerine karşın enerji değeri.İşte ben bu örneğin 0.505 i değiştirdiğimde aşağıda neleri daha değiştireceğim şimdi biraz değişti ama simetrik değil yani olmadı


Mesaj tekrarı yüzünden mesajınız birleştirildi. Bu mesajın gönderim tarihi : 09 Mart 2012 - 16:09:21
Hatta ben ekran çıktısı üzerinde anlatayım
Değişiklik yapmadan önce

Burada en sağ üstte 2 tane kuyu gibi birşey var.İşte bunlar kuantum kuyuları :)Kuyunun altındaki kırmızı çıkıntının ona yaklaşmasını istiyorum(Band gap in küçülmesi).Değişikliği yapınca yaklaşma oldu ama simetrik değil.Yukarısı ne kadar aşağı iniyorsa alt tarafta okadar çıkmalı aslında.