地域密着型エリア広告配信リクルートの地域密着型広告ドコイク?アドネットワーク Ads by ドコイク?[無料でホームページを作成] [通報・削除依頼]
[無料でホームページを作成] [通報・削除依頼]

Igor macro(Gaussian broadening)

Macro gaussian_borodening(x,y,FWHM,xmin,xmax,estep)
string x="eV";
string y="Si1";
variable FWHM=0.2;
variable xmin=-13, xmax=10,estep=0.02720;
silent 1; pauseupdate;
variable ndata, nstep;
variable i,j;
variable sigma;
sigma=FWHM/(2*sqrt(2*ln(2))); print "sigma=",sigma, ", sigma^2=", sigma^2;
ndata=(DimSize($x,0));
nstep=round((xmax-xmin)/estep); //print nstep;
make/O/N=(nstep) xgb,ygb;
i=0;
i=0;
do
  xgb[i] = xmin + estep * i;
  ygb[i] = 0.0
  i=i+1;
while(i<=nstep)
variable A=(1/(sqrt(2*3.1415926)*sigma))
variable B=(2*sigma^2)
variable four_sigma = 4*sigma
i=0;
do
  j=0;
  do
    if(xgb[i]+four_sigma>=$x[j] && $x[j]>=xgb[i]-four_sigma)
      ygb[i] = ygb[i] + $y[j]*A*exp(-(xgb[i]-$x[j])^2/B)
    endif
    j=j+1;
  while(j<ndata)
  i=i+1
while(i<nstep)
//
Display $y vs $x
ModifyGraph rgb=(0,0,0),mode($y)=1
//
string new_xw, new_yw
new_xw = "xgb_" + y
new_yw = "ygb_" + y
make/O/N=(nstep) $new_xw,$new_yw;
$new_xw=xgb;
$new_yw=ygb;
killwaves xgb, ygb;
//
AppendToGraph $new_yw vs $new_xw
//AppendToGraph ygb vs xgb
ModifyGraph tick=2,mirror=1,fSize=18,standoff=0,font="Arial";DelayUpdate
Label left "DOS, \\f02N\\f00\\M(\\F'Symbol'e\\F'Arial')";DelayUpdate
Label bottom "Energy, \\F'Symbol'\\f02e-e\\f00\\B\\F'Arial'f\\M / eV"
ModifyGraph width=340,height=227,gFont="Arial"
ModifyGraph rgb=(0,0,0)
ModifyGraph lsize=2
//Legend/C/N=text0/J/F=0/B=1/A=MC "\\s('total_DOSx0.5') Band calc. (FLPAW-GGA)"
ModifyGraph gfSize=18,gmSize=18
SetAxis bottom xmin,xmax
SetAxis left 0.0,5
EndMacro
QRコード
携帯用QRコード
アクセス数
ページビュー数