origin6[1].0磁性拟合过程.doc
- 文档编号:7305747
- 上传时间:2023-05-11
- 格式:DOC
- 页数:8
- 大小:411KB
origin6[1].0磁性拟合过程.doc
《origin6[1].0磁性拟合过程.doc》由会员分享,可在线阅读,更多相关《origin6[1].0磁性拟合过程.doc(8页珍藏版)》请在冰点文库上搜索。
回应小木虫及晶体解析联合网站上部分虫友恳求,将origin软件的拟合过程做成视频文件,遗憾的是视频文件太大上传不了网站,故以截图的方式给各位同行粗略演示拟合过程。
1.先将数据导入origin软件,并建立温度(T),磁化率(Xm)和磁化率温度乘积(XmT)各列。
如先右击空白栏勾点addnewcolumn,出现空白栏,在空白栏的顶端右击勾点setcolumnvalues,出现对话框(如下),在“Col(X)=”右边填写:
col(温度栏)*相应式量/外场强度/样品质量,如col(b)*1589.68/5000/0.02641,――OK
2.挑选A及D栏(ctrl键选择)再点plot—scatter做图
3.建立拟合公式(该例子为双核锰(S1=S2=5/2)
(A)双核锰(S1=S2=5/2)的公式:
见MolecularMagnetism---OlivierHahn114页
Origin拟合过程:
对Mn2拟合cmT曲线
符号:
p1---g;p2---J/kBA=2*exp(2*p2/x)+10*exp(6*p2/x)+28*exp(12*p2/x)+60*exp(20*p2/x)+110*exp(30*p2/x);
B=1+3*exp(2*p2/x)+5*exp(6*p2/x)+7*exp(12*p2/x)+9*exp(20*p2/x)+11*exp(30*p2/x);
y=x*0.375*p1*p1*(A/B)
(B)考虑分子场近似:
Origin拟合过程:
对Mn2拟合cmT曲线
符号:
p1---g;p2---J/kB;p3---zJ/kB
A=2*exp(2*p2/x)+10*exp(6*p2/x)+28*exp(12*p2/x)+60*exp(20*p2/x)+110*exp(30*p2/x);
B=1+3*exp(2*p2/x)+5*exp(6*p2/x)+7*exp(12*p2/x)+9*exp(20*p2/x)+11*exp(30*p2/x);
s=0.375*p1*p1*(A/B);
y=s/(1-2*p3*s/(0.375*p1*p1))
(C)考虑顺磁杂质及温度无关项:
Origin拟合过程:
对Mn2拟合cmT曲线
符号:
p1---g;p2---J/kB;p3---zJ/kB;p4---cm’;p5---cTIP
A=2*exp(2*p2/x)+10*exp(6*p2/x)+28*exp(12*p2/x)+60*exp(20*p2/x)+110*exp(30*p2/x);
B=1+3*exp(2*p2/x)+5*exp(6*p2/x)+7*exp(12*p2/x)+9*exp(20*p2/x)+11*exp(30*p2/x);
s=0.375*p1*p1*(A/B)/x;
t=s/(1-2*p3*s/(0.375*p1*p1))*(1-p4)+0.125*2.5*3.5/x*p4+p5;
y=t*x
4.拷贝拟合式子(若同时考虑分子场近似、顺磁杂质及温度无关项)其式子如下:
A=2*exp(2*p2/x)+10*exp(6*p2/x)+28*exp(12*p2/x)+60*exp(20*p2/x)+110*exp(30*p2/x);
B=1+3*exp(2*p2/x)+5*exp(6*p2/x)+7*exp(12*p2/x)+9*exp(20*p2/x)+11*exp(30*p2/x);
s=0.375*p1*p1*(A/B)/x;
t=s/(1-2*p3*s/(0.375*p1*p1))*(1-p4)+0.125*2.5*3.5/x*p4+p5;
y=t*x
(注意每行后边“;”号)
5.打开XmT-T曲线origin软件界面,点analysis-nonlinearcurvefit
出来如下界面
6.先建立函数,(如果已建有该函数点More再选择)
若需创建新函数即点New按钮,
Name,Form(选Y-script),NumberofParam个数,Definition(此栏就将从别处拷来的拟合式子粘贴),如下填写
Save保存,再点BasicMode进入上级界面,
再点accept
7.选择数据列,点SelectDataset
预先记住温度及磁化率温度乘积两列的栏名(columnname),如A及D
再点对话框其中---->xindep?
?
?
[1…0]
再从availabledatasets中选后边是a的那一行,如a9129mt300k_a
再点Assign,则自变量x的值就变为A列了,同时把因变量y的值变为D列。
当然可以改动温度范围,改动<=Row<=两侧的值就行。
8.好,现在就开始拟合了,点StartFitting
记住P1、P2、P3、P4----各代表什么量,估计其值并填入value栏中,有些量不想参与拟合或定值的,可以在Value值改为0或某个定值,Vary处的勾去掉;
此处P1---g,先给个初值2.0;P2---J为二核锰簇内的交换作用从曲线知很弱,可以设个小的负值,如-1.0;P3---zJ’二核锰簇间的磁交换作用。
如果无法同时拟合就先固定一些量来拟合。
再点10lter,一直拟合直到拟合值及拟合诊断因子值(Chi^2)不变.如果Chi^2很小,同时拟合值也很合乎情理,估计就差不多了,如果二者任何一值不合理,就必须再行分析,是初值设的不好(重设),有些量没有纳入拟合(同时考虑),有些量不适合同时拟合(牺牲一些量另拟合另一些量),模型不对(重建模型)等。
如果觉得拟合到最佳程度,就点Done.
附注:
该实例拟合还没完善(Chi^2诊断值偏大)。
谢谢大家
倏忽于2008年12月12日
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- origin6 磁性 拟合 过程