使用 Scherrer 方程计算晶粒尺寸 [英] Using Scherrer equation for calculating the grain size
本文介绍了使用 Scherrer 方程计算晶粒尺寸的处理方法,对大家解决问题具有一定的参考价值,需要的朋友们下面随着小编来一起学习吧!
问题描述
我正在尝试通过 Scherrer 方程计算晶粒尺寸,但我一直坚持 FWHM.
I'm trying to calculate the grain size by Scherrer equation but I have stuck in FWHM.
import numpy as np
#import math
k = 0.94
wave_length = 1.5406e-10
data = np.genfromtxt("G3.txt")
indice = np.argmax(data[:,1])
peak = (data[indice, :])
#D = (k*wave_length) / (beta*cos((math.radian(theta))
信息:Scherrer 方程、半高全宽,相关问题
推荐答案
我不知道如何在 python 中解决这个问题(至少目前是这样).所以我在Matlab中制作了它.
I don't know how to solve this problem in python (at least for this moment). So I made it in Matlab.
clear all
clc
A = dlmread('YOUR DATAS'); %Firstly add to path
plot(A(:,1),A(:,2)) %Plotting the graph
hold on
min_peak = input('Just write a value that is higher than minimum peak values: ');
%This value must be between requested peaks and non-requested peaks (you can see this in graph)
[yval, yval_i] = findpeaks(A(:,2),'MinPeakHeight',min_peak); %Finding peaks
scatter(A(yval_i,1), yval); %Showing peaks
Beta = [];
xval = [];
for k = 1:size(yval_i,1) %Finding x values corresponding to y (peak) values
xval1 = A(yval_i(k),1);
xval = [xval xval1];
end
Theta = xval / 2;
for i = 1:size(yval,1) %Finding half of max. peak values
yval_i1 = yval_i(i,1);
while (yval(i,1))/2 < A(yval_i1+1,2)
yval_i1 = yval_i1+1;
end
yval_i2 = yval_i(i,1);
while (yval(i,1))/2 < A(yval_i2-1,2)
yval_i2 = yval_i2-1;
end
plot(A(yval_i2,1)*ones(size(A(:,2))), A(:,2));
plot(A(yval_i1,1)*ones(size(A(:,2))), A(:,2));
% hold on
% scatter(A(yval_i1,1),A(yval_i1,2))
% scatter(A(yval_i2,1),A(yval_i2,2))
B = abs(A(yval_i1,1)-A(yval_i2,1));
Beta = [Beta B];
end
Beta
K = 0.94;
Lambda = 1.5406e-10;
To = [];
for j = 1:size(Beta,2)
To1 = (K*Lambda)/(Beta(j)*cos(Theta(j)));
To = [To To1];
end
To = abs(To)
这篇关于使用 Scherrer 方程计算晶粒尺寸的文章就介绍到这了,希望我们推荐的答案对大家有所帮助,也希望大家多多支持IT屋!
查看全文