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

本文链接：https://blog.csdn.net/weixin_46039719/article/details/128680667

👨‍🎓个人主页：研学社的博客

💥💥💞💞欢迎来到本博客❤️❤️💥💥

🏆博主优势：🌞🌞🌞博客内容尽量做到思维缜密，逻辑清晰，为了方便读者。

⛳️座右铭：行百里者，半于九十。

📋📋📋本文目录如下：🎁🎁🎁

目录

💥1 概述

📚2 运行结果

🎉3 参考文献

🌈4 Matlab代码实现

💥1 概述

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

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

📚2 运行结果

部分代码：

%% 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)