使用x射线与大块材料中原子相互作用的基本理论模拟相对强度XRD图案(Matlab代码实现)

           目录

💥1 概述

📚2 运行结果

🎉3 参考文献

👨‍💻4 Matlab代码


💥1 概述

这一系列 MATLAB 代码旨在使用x射线与大块材料中原子相互作用的基本理论模拟相对强度XRD图案。

除非需要混合物成分规则,否则无需预定义参数。在这种情况下应编辑Scatteringdata.xlsx并将所需的组件添加到工作表的末尾。

要操作该程序,请使用以下过程:

  1. 运行XRD_Intensity_Simulation.m
  2. 选择 X 射线源

📚2 运行结果

 

 

🎉3 参考文献

[1]胡林彦,张庆军,沈毅.X射线衍射分析的实验方法及其应用[J].河北理工学院学报,2004(03):83-86+93.

👨‍💻4 Matlab代码

主函数部分代码:

%% Theoretical Peak Intensity Calcuclator
% This program takes user input about a specific crystal structure and
% determines the peak positions and relative intensities. Based on theory
% that can be found in Structure of Materials: An Introduction to
% Crystallography, Diffraction, and Symmetry by Marc De Gaef and Michael E.
% McHenry
%% 

%Extract the necessary simulation parameters from the user
[crystaltype,hklvec,lambda,atoms,lattice_spacing] = dialogue();

%Calcualte the d spacing of individual crystal planes. 
dspacevec=dspace(hklvec,crystaltype,lattice_spacing);

%Calcualte the theta and s values at different hkl peaks. 
[thetavec,svec]=theta_s_calc(lambda,dspacevec);

%Calculate the scattering factor for each atom involved. 
scatterfactorvec=scatterfactor(svec, crystaltype, atoms);

%Calculate the strucutre factor for each plane. 
structurefactors=structurefactor(scatterfactorvec,crystaltype);

%Calculate absolute intensity
intensityvec=intensitycalc(structurefactors,crystaltype,thetavec);

%Compare Absolute and Relative Intensities 
compareintensity(intensityvec,hklvec,thetavec,crystaltype,atoms,lattice_spacing)

 


本文来自互联网用户投稿,文章观点仅代表作者本人,不代表本站立场,不承担相关法律责任。如若转载,请注明出处。 如若内容造成侵权/违法违规/事实不符,请点击【内容举报】进行投诉反馈!

相关文章

立即
投稿

微信公众账号

微信扫一扫加关注

返回
顶部