首页 / 更多教程 / 【Matlab 064期】【物理应用11】计算晶体结构的x射线衍射图谱matlab源码

【Matlab 064期】【物理应用11】计算晶体结构的x射线衍射图谱matlab源码

内容导读

互联网集市收集整理的这篇技术教程文章主要介绍了【Matlab 064期】【物理应用11】计算晶体结构的x射线衍射图谱matlab源码，小编现在分享给大家，供广大互联网技能从业者学习和参考。文章包含11362字，纯文字阅读大概需要17分钟。

内容图文

% function xresult = refrac(varargin)
% REFRAC  Material property function interacted with X-rays 
%    r = refrac('Formula String',Energy,MassDensity,'PlotStyle') returns
%    disperion absorption, critical angle and attenuation length
%    correponding to the enery input. Also plot if energy is a list.
%
%    Input: 
%       1) Chemical formula is case sensitive (e.g. CO for Carbon Monoxide
%           vs Co for Cobalt);
%       2) Energy (0.03KeV~30KeV). Can be a single energy or an energy list;
%       3) Mass density (g/cm^3);
%       4) PlotStyle can be 'logxy','linear','logx' or 'logy'. Default
%           is 'linear'.
%    
%    Output structure:
%       1) Chemical formula;
%       2) Molecular weight;
%       3) Number of electrons per molecule;
%       4) Mass density (g/cm^3);
%       5) Electron density (1/m^3);
%       6) X-ray energy (KeV);
%       7) Corresponding X-ray wavelength (m);
%       8) Dispersion;
%       9) Absorption;
%       10) Critical angle (degree);
%       11) Attenuation length (m).
%       
%    Plot:
%       Only when energy is a list, plot dispersion, absorption, critical
%       angle and attenuation length.
%
%    Example 1: >> result=refrac('Si3N4',8.04778,3.44)
%           gives an output with structure below:
%               result = 
%                     formula: 'Si3N4'
%             molecularWeight: 140.28
%           numberOfElectrons: 70
%                 massDensity: 3.44
%             electronDensity: 1.4767e+028
%                      energy: 8.0478
%                  wavelength: 1.5406e-010
%                  dispersion: 1.1164e-005
%                  absorption: 1.6477e-007
%               criticalAngle: 0.27074
%                   attLength: 7.4403e-005
%    
%    Example 2: >>result=refrac('Si3N4',8:0.5:10,3.44)
%           gives an output with structure below together with a plot in
%           linear scale.
%               result =
%                     formula: 'Si3N4'
%             molecularWeight: 140.28
%           numberOfElectrons: 70
%                 massDensity: 3.44
%             electronDensity: 1.4767e+028
%                      energy: [5x1 double]
%                  wavelength: [5x1 double]
%                  dispersion: [5x1 double]
%                  absorption: [5x1 double]
%               criticalAngle: [5x1 double]
%                   attLength: [5x1 double]
%           result.energy, result.dispersion, etc. show the lists of this
%           values. >>result=refrac('Si3N4',8:0.5:10,3.44,'logxy') gives a
%           plot in log-log plot.
%
%    For more information about X-ray interactionw with matter, go to
%       http://www.cxro.lbl.gov
%       http://www.nist.gov/
%   
%    Atomic scattering factor table is taken from the above two websites.
%
%    Copyright 2004 Zhang Jiang
%    $Revision: 1.0 $  $Date: 2004/10/10 18:52:19 $
% refrac('Si3N4',8:0.5:10,3.44)
clc;
clear all;
% if nargin ~= 3 & nargin ~= 4
%     error('Incorrect inputment.');
%     return;
% end
formulaStr = 'Si3N4';%varargin{1};
energy = 8:0.5:10;%varargin{2};
massDensity =3.44;% varargin{3};
if ~ischar(formulaStr)
    error('Invalid chemical formula.');
    return;
end
if ~isnumeric(energy) | isempty(energy)
    error('Invalid x-ray energy.');
    return;
end
if min(energy)<0.03 | max(energy)>30
    error('Energy is out of range 0.03KeV~30KeV.');
    return;
end
if ~isnumeric(massDensity) | length(massDensity) > 1
    error('Invalid mass density.');
    return;
end
if nargin == 4
    plotStyle = varargin{4};
    if ~strcmp(plotStyle,'logxy') &...
            ~strcmp(plotStyle,'linear') &...
            ~strcmp(plotStyle,'logx') &...
            ~strcmp(plotStyle,'logy')
        error('Invalid PlotStyle.');
        return;
    end
else
    plotStyle = 'linear';
end
%===============================================================
% some constants
%===============================================================
THOMPSON = 2.81794092e-15;          % m
SPEEDOFLIGHT = 299792458;           % m/sec
PLANCK = 6.626068e-34;              % m^2*kg/sec
ELEMENTCHARGE = 1.60217646e-19;     % Coulombs
AVOGADRO = 6.02214199e23;           % mole^-1
%===============================================================
% sort energy list and convert to wavelength
%===============================================================
energy = sort(energy);  % sort energy from min to max
wavelength = (SPEEDOFLIGHT*PLANCK/ELEMENTCHARGE)./(energy'*1000.0);
%===============================================================
% determine elements and number of atoms in the formula
%===============================================================
nElements = 0;
formula.elements = {};
formula.nAtoms = {};
try
    for iFormulaStr = 1:length(formulaStr)
        formulaChar = formulaStr(iFormulaStr);
        if formulaChar <= 'Z' &  formulaChar >= 'A'
            nElements = nElements + 1;
            formula.elements{nElements} = formulaChar;
            formula.nAtoms{nElements} = '0';
        elseif formulaChar <= 'z' &  formulaChar >= 'a'...
                & ((formulaStr(iFormulaStr-1) <= 'Z'& formulaStr(iFormulaStr-1) >= 'A')...
                | (formulaStr(iFormulaStr-1) <= 'z'& formulaStr(iFormulaStr-1) >= 'a'))
            formula.elements{nElements} = ...
                [formula.elements{nElements},formulaChar];
        elseif (formulaChar <= '9' &  formulaChar >= '0') | formulaChar == '.'
            formula.nAtoms{nElements} = ...
                [formula.nAtoms{nElements},formulaChar];
        else
            error('Invalid chemical formula.');
            return;
        end
    end
catch
    error('Invalid chemical formula.');
    return;
end
for iElements = 1:nElements
    formula.nAtoms{iElements} = str2num(formula.nAtoms{iElements});
    if formula.nAtoms{iElements} == 0
        formula.nAtoms{iElements} = 1;
    end
end
%===============================================================
% read f1 and f2 from tables
%===============================================================
formula.f1f2Table = {};
folder='.\AtomicScatteringFactor\';
file1= dir(fullfile(folder,'*.nff'));
DD=file1.name;
for iElements = 1:nElements
   
%     file = fullfile(matlabroot,'toolbox','xrayrefraction',...
%         'AtomicScatteringFactor',lower(formula.elements{iElements}));
    TT=formula.elements{iElements};
    file = [folder,file1(iElements).name];
    try 
        fid = fopen(file);
        fgetl(fid);
    catch
        
        error(['Element ''',formula.elements{iElements},...
            ''' is NOT in the table list.']);
        return;
    end
    table = cell2mat(textscan(fid,'%f %f %f'));
    table(find(table(:,1)<29),:) = [];
    formula.f1f2Table{iElements} = table;
    fclose(fid);
end
%===============================================================
% determine the atomic number and atomic weight
%===============================================================
formula.atomicNumber = {};
formula.atomicWeight = {};
% file = fullfile(matlabroot,'toolbox','xrayrefraction','periodictable.mat');
file='.\periodictable.mat';
load(file);
for iElements = 1:nElements
    for iAtomicnum =1:length(elementAbbr)
        if strcmp(elementAbbr{iAtomicnum},formula.elements{iElements})
            formula.atomicNumber{iElements} = iAtomicnum;
            formula.atomicWeight{iElements} = atomicWeight(iAtomicnum);
            break;
        end
    end
end
%===============================================================
% determine molecular weight and number of electrons
%===============================================================
formula.molecularWeight = 0;
formula.numberOfElectrons = 0;
for iElements = 1:nElements
    formula.molecularWeight = formula.molecularWeight + ...
        formula.nAtoms{iElements}*formula.atomicWeight{iElements};
    formula.numberOfElectrons = formula.numberOfElectrons +...
        formula.atomicNumber{iElements}*formula.nAtoms{iElements};
end
%===============================================================
% interpolate to get f1 and f2 for given energies
%===============================================================
formula.interpf1 = {};
formula.interpf2 = {};
for iElements = 1:nElements
%     formula.interpf1{iElements} = exp(interp1(...
%         formula.f1f2Table{iElements}(:,1),...
%         log(formula.f1f2Table{iElements}(:,2)),...
%         energy'*1000,'cubic'));
%     formula.interpf2{iElements} = exp(interp1(...
%         formula.f1f2Table{iElements}(:,1),...
%         log(formula.f1f2Table{iElements}(:,3)),...
%         energy'*1000,'cubic'));
    formula.interpf1{iElements} = interp1(...
        formula.f1f2Table{iElements}(:,1),...
        formula.f1f2Table{iElements}(:,2),...
        energy'*1000,'cubic');
    formula.interpf2{iElements} = interp1(...
        formula.f1f2Table{iElements}(:,1),...
        formula.f1f2Table{iElements}(:,3),...
        energy'*1000,'cubic');
end
%===============================================================
% calculate dispersion, absorption, critical angle and attenuation length
%===============================================================
% initialize
xresult.formula = formulaStr;
xresult.molecularWeight = formula.molecularWeight;
xresult.numberOfElectrons = formula.numberOfElectrons;
xresult.massDensity = massDensity;
xresult.electronDensity = 1e6*massDensity/formula.molecularWeight*AVOGADRO...
    *formula.numberOfElectrons;
xresult.energy = energy';
xresult.wavelength = wavelength;
xresult.dispersion = zeros(length(energy),1);
xresult.absorption = zeros(length(energy),1);
xresult.criticalAngle  = zeros(length(energy),1);
xresult.attLength = zeros(length(energy),1);
for iElements = 1:nElements
    xresult.dispersion = xresult.dispersion + ...
        wavelength.^2/(2*pi)*THOMPSON*AVOGADRO*massDensity*1e6 / ...
        formula.molecularWeight * ...
        formula.nAtoms{iElements} .* ...
        formula.interpf1{iElements};
    xresult.absorption = xresult.absorption + ...
        wavelength.^2/(2*pi)*THOMPSON*AVOGADRO*massDensity*1e6 / ...
        formula.molecularWeight * ...
        formula.nAtoms{iElements} .* ...
        formula.interpf2{iElements};
end
xresult.criticalAngle = sqrt((2*xresult.dispersion))*180/pi;
xresult.attLength = xresult.wavelength./xresult.absorption/(4*pi);
% assignin('base','formula',formula);
%===============================================================
% plot
%===============================================================
% plot only when energy is a list
if length(energy) == 1
    return;
end
switch plotStyle
    case 'logxy'
        xscale = 'log';
        yscale = 'log';
    case 'logx'
        xscale = 'log';
        yscale = 'linear';
    case 'logy'
        xscale = 'linear';
        yscale = 'log';
    case 'linear'
        xscale = 'linear';
        yscale = 'linear';
end
xresultFig = figure(...
    'Name',['X-Ray Interactions with ',formulaStr],...
    'NumberTitle','off',...
    'PaperOrientation','landscape');
subplot(2,2,1);
% plot(xresult.energy,xresult.dispersion,'b-o',...
%     'MarkerEdge','r','MarkerSize',0.1);
plot(xresult.energy,xresult.dispersion);
xlabel('Energy (KeV)');
ylabel('Dispersion');
set(gca,'XScale',xscale,'YScale',yscale);
grid on;
box on;
title('Dispersion');
subplot(2,2,2);
% plot(xresult.energy,xresult.absorption,'b-o',...
%     'MarkerEdge','r','MarkerSize',0.1);
plot(xresult.energy,xresult.absorption);
xlabel('Energy (KeV)');
ylabel('Absorption');
set(gca,'XScale',xscale,'YScale',yscale);
grid on;
box on;
title('Absorption');
subplot(2,2,3);
% plot(xresult.energy,xresult.criticalAngle,'b-o',...
%     'MarkerEdge','r','MarkerSize',0.1);
plot(xresult.energy,xresult.criticalAngle);
xlabel('Energy (KeV)');
ylabel('Critical Angle (degree)');
set(gca,'XScale',xscale,'YScale',yscale);
grid on;
box on;
title('Critical Angle');
subplot(2,2,4);
% plot(xresult.energy,xresult.attLength,'b-o',...
%     'MarkerEdge','r','MarkerSize',0.1);
plot(xresult.energy,xresult.attLength);
xlabel('Energy (KeV)');
ylabel('Attenuation Length (m)');
set(gca,'XScale',xscale,'YScale',yscale);
grid on;
box on;
title('Attenuation Length');
suptitle(['X-Ray Interactions with ',formulaStr,...
    ', \rho_{mass}=',num2str(xresult.massDensity),'g/cm^3']);

【Matlab 064期】【物理应用11】计算晶体结构的x射线衍射图谱matlab源码 - 文章图片

内容总结

以上是互联网集市为您收集整理的【Matlab 064期】【物理应用11】计算晶体结构的x射线衍射图谱matlab源码全部内容，希望文章能够帮你解决【Matlab 064期】【物理应用11】计算晶体结构的x射线衍射图谱matlab源码所遇到的程序开发问题。如果觉得互联网集市技术教程内容还不错，欢迎将互联网集市网站推荐给程序员好友。

内容备注

版权声明：本文内容由互联网用户自发贡献，该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务，不拥有所有权，不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容，请发送邮件至 gblab@vip.qq.com 举报，一经查实，本站将立刻删除。

内容手机端

扫描二维码推送至手机访问。

本文链接：https://qyyshop.com/info/612462.html

来源：【匿名】

【上一篇】node 深入浅出【下一篇】关于IE的RegExp.exec的问题

更多 ►

【【Matlab 064期】【物理应用11】计算晶体结构的x射线衍射图谱matlab源码】教程文章相关的互联网学习教程文章

【Matlab 064期】【物理应用11】计算晶体结构的x射线衍射图谱matlab源码【代码】

% function xresult = refrac(varargin) % REFRAC Material property function interacted with X-rays % r = refrac('Formula String',Energy,MassDensity,'PlotStyle') returns % disperion absorption, critical angle and attenuation length % correponding to the enery input. Also plot if energy is a list. % % Input: % 1) Chemical formula is case sensitive (e.g. CO for Carbon Monoxide % ...

跨应用的访问 contentprovider【图】

首先需要明确一点就是插入的位置还是这个应用的sqllite的数据库650) this.width=650;" src="/upload/getfiles/default/2022/11/14/20221114082637500.jpg" title="快照16.jpg" />650) this.width=650;" src="/upload/getfiles/default/2022/11/14/20221114082637877.jpg" title="快照17.jpg" />思路:provider 提供一个可以访问的uri，然后外部直接来访问，要是和内部提供的uri匹配上了。就直接返回某一个操作。650) this.width=650...

padding应用-导航

<!DOCTYPE html> <html lang="en"> <head> <meta charset="UTF-8"> <meta http-equiv="X-UA-Compatible" content="IE=edge"> <meta name="viewport" content="width=device-width, initial-scale=1.0"> <title>Document</title> <style> .nav { height: 41px; border-top: 3px solid #ff8500; border-bottom: 1px solid #edeef0; background-color...

【计网实验】实验一：应用协议与数据包分析实验(使用Wireshark)【图】

实验目的通过本实验，熟练掌握Wireshark的操作和使用，学习对HTTP协议进行分析。实验步骤 1、下载安装Wireshark； 2、在cmd中输入命令“ipconfig”，得到本机IPv4 地址；（该地址用于帮助分析） 3、开启抓包，访问 www.china.com，页面加载完毕之后停止抓包。抓包结果如下：（记以及要把抓包结果保存下来，也许验收的时候助教要看吧） 4、利用显示过滤器进行筛选，选择有用信息进行分析。报文分析 1、综合分析截获的报文，查看...

栈应用：最小栈（第二题）【代码】

题目：设计一个支持 push，pop，top 操作，并能在常数时间内检索到最小元素的栈。push(x) -- 将元素 x 推入栈中。pop() -- 删除栈顶的元素。top() -- 获取栈顶元素。getMin() -- 检索栈中的最小元素分析：1.每次压入两个元素，一个是元素本身，另外一个进行判断，是否为最小，即存放当前栈最小值，代码实现如下 1/** initialize your data structure here. */ 2public MinStack() {3 stack = new Stack<Integer>();4}5 6public...

PX4/Pixhawk---基于NSH调试的uORB第一个应用测试【图】

PX4/Pixhawk—基于NSH调试的uORB第一个应用测试1 NSH连接测试??(1) 测试前准备：PX4FMU（已刷好固件）USB线安装好PX4 Toolchain??(2) Pixhawk通过USB连接电脑，并安装好了驱动： ????(3) 打开Tera Term软件（PX4 Toolchain->TeraTerm）： ??2 开始第一个应用??http://www.pixhawk.com/start?id=zh/dev/px4_simple_app/*‘).addClass(‘pre-numbering‘).hide();$(this).addClass(‘has-numbering‘).parent().append($numbering)...

数字证书在WEB应用中登录

1数字证书登录认证的优点作为企业信息系统的第一道大门，身份认证是确保企业信息资源只能被合法用户所访问的重要保障。传统的口令认证方式虽然简单，但是由于其易受到窃听、重放等攻击的安全缺陷，使其已无法满足当前复杂网络环境下的安全认证需求。传统账号+口令登录的弊端：1. 口令易被猜测。由于有的用户为了方便记忆，使用非常简单的口令，比如“1234”、“abcd”等这些非常容易被猜测的口令。2. 口令易被窃听...

PLSQL导入Excel亲测各种问题-应用程序与驱动程序不匹配【图】

天才第一步:打开PLSQL-TOOLS-ODBC Importer问题1: 解决方法: 首先，确定自己PLSQL，Oracle，Office分别是32位还是64位必须全部一致(现在一般都是64位) 然后添加自己配置的数据源 64位：C:\Windows\System32\odbcad32.exe 32位：C:\Windows\SysWOW64\odbcad32.exe点击添加点击确定如果这里没有excel驱动,只有Oracle和SQL Server两个驱动的话, 建议查看Office32位还是64位附Office下载链接：https://msdn.itellyou.cn/ 如图然后...

[cobra]-强大的CLI应用程序库【代码】【图】

关于我我的博客|文章首发 cobra介绍 Cobra是一个用于创建强大的现代CLI应用程序的库，也是一个用于生成应用程序和命令文件的程序。 Cobra用于许多Go项目，如Kubernetes、Hugo和Github CLI等。此列表包含了使用Cobra的更广泛的项目列表。概述 Cobra是一个库，提供了一个简单的界面来创建功能强大的现代CLI界面，类似于git＆go工具。 Cobra还是一个应用程序，它将生成您的应用程序框架以快速开发基于Cobra的应用程序。眼镜蛇提供：...

钉钉应用云开发实战营第5次打卡任务【代码】【图】

钉钉应用云开发实战营第5次打卡任务任务目标：根据今天的课程内容，布置一个会议室预定系统上线钉钉应用，需要用户将手机号码读取成功，显示在会议室预定的首页，然后截图提交即可。会议室预定教程：https://developers.dingtalk.com/document/tutorial/quickly-build-h5-microapplications-meeting-room-reservation 1、创建钉钉应用 1.1、登录钉钉开发者后台地址：https://open-dev.dingtalk.com/1.2、创建应用1.3、创建成功...

点酷听书手机客户端应用源码

点酷听书，让你不再总是疲劳看书了，真人发音准确非常，还有海量电子书等你来看哦。客户端代码包括如下几个方面本地阅读器这个比较多我就不明说了，我选的也不是最好的，我选的是最简单的。语音播放是之前论坛首页的foc离线语音朗读库，这个效果其实非常好，真心的。但是好像目前到期了，不过在我的努力下，延长了使用时间网络下载部分我选择是LNreader这个开源项目，当然由于这个项目采用了jsoup库，主要做html的解析显示，...

凡客传说4-DelayQueue使用方式，源码分析,以及应用场景【图】

DelayQueue 顾名思义,它是一个延时队列使用方式 :假设我们生产者提交一个任务，消费者5秒钟之后才可以执行，那么我们可以把任务定义为如下格式，并实现Delayed接口，其中data是任务存储的信息。/*** 具体的任务* @author wangshixiang*/ public class Task implements Delayed {/*** 数据*/private final String data;/*** 任务执行时间*/private final long time;public Task(String data,TimeUnit timeUnit,long time){this.dat...

POI解析Excel之应用反射等技术实现动态读取【代码】【图】

目录树背景技术选型问题分析技术要点及难点分析源码分析测试用例背景Tip：因为产品提的需求我都开发完了，进行了项目提测；前天老大走过来说：你用spring-boot开发一个解析Excel的jar包.....详细对话如下：A：世生，你用spring-boot开发一个解析Excel的jar包。B：为什么不在原来的项目上进行开发呢？(很纳闷，这个东西不是一般用于前端上传数据的嘛，这是后端层，咋搞这个)A：因为xxxx这个需求有比较多的数据需要初始化，这个jar...

统一 – 如何使启动器进度条与我的应用程序一起使用？【代码】

背景研究我正在尝试更新Unity启动器中的进度条,以便使用Quickly调用test创建一个简单的python / Gtk应用程序;但是,按照this video中的说明,我无法成功更新Unity启动器中的进度条.在Unity Integration视频中,没有使用Quickly,因此应用程序的结构方式略有不同,如果没有在默认的Quickly ubuntu应用程序模板应用程序中进行修改,视频中使用的代码似乎无法正常运行. 截图以下是Unity Launcher中当前显示的应用程序图标的屏幕截图. 以下...

如何应用执行的SQL的执行计划【代码】

通过explain ，我们可以获取特定SQL 的执行计划。但对于同一条SQL，不同的变量、不同的系统负荷，其执行计划可能不同。我们要如何取得SQL执行时间点的执行计划？KingbaseES 提供了 auto_explain 扩展插件，可以自动跟踪SQL执行计划。1、设置参数设置参数：shared_preload_libraries = ‘auto_explain‘auto_explain.log_min_duration = 1000 auto_explain.log_min_duration 是最短语句执行时间（以毫秒为单位），将此设置为0 将记...

首页 / 更多教程 / 【Matlab 064期】【物理应用11】计算晶体结构的x射线衍射图谱matlab源码

【Matlab 064期】【物理应用11】计算晶体结构的x射线衍射图谱matlab源码

内容导读

内容图文

内容总结

内容备注

内容手机端

【【Matlab 064期】【物理应用11】计算晶体结构的x射线衍射图谱matlab源码】教程文章相关的互联网学习教程文章

【Matlab 064期】【物理应用11】计算晶体结构的x射线衍射图谱matlab源码【代码】

跨应用的访问 contentprovider【图】

padding应用-导航

【计网实验】实验一：应用协议与数据包分析实验(使用Wireshark)【图】

栈应用：最小栈（第二题）【代码】

PX4/Pixhawk---基于NSH调试的uORB第一个应用测试【图】

数字证书在WEB应用中登录

PLSQL导入Excel亲测各种问题-应用程序与驱动程序不匹配【图】

[cobra]-强大的CLI应用程序库【代码】【图】

钉钉应用云开发实战营第5次打卡任务【代码】【图】

点酷听书手机客户端应用源码

凡客传说4-DelayQueue使用方式，源码分析,以及应用场景【图】

POI解析Excel之应用反射等技术实现动态读取【代码】【图】

统一 – 如何使启动器进度条与我的应用程序一起使用？【代码】

如何应用执行的SQL的执行计划【代码】

更多教程 - 最新教程

更多教程 - 最热教程

首页 / 更多教程 / 【Matlab 064期】【物理应用11】 计算晶体结构的x射线衍射图谱matlab源码

【Matlab 064期】【物理应用11】 计算晶体结构的x射线衍射图谱matlab源码

内容导读

内容图文

内容总结

内容备注

内容手机端

【【Matlab 064期】【物理应用11】 计算晶体结构的x射线衍射图谱matlab源码】教程文章相关的互联网学习教程文章

更多教程 - 最新教程

更多教程 - 最热教程

首页 / 更多教程 / 【Matlab 064期】【物理应用11】计算晶体结构的x射线衍射图谱matlab源码

【Matlab 064期】【物理应用11】计算晶体结构的x射线衍射图谱matlab源码

【【Matlab 064期】【物理应用11】计算晶体结构的x射线衍射图谱matlab源码】教程文章相关的互联网学习教程文章