spss进行籽粒灌浆拟合的详细步骤_第1页
spss进行籽粒灌浆拟合的详细步骤_第2页
spss进行籽粒灌浆拟合的详细步骤_第3页
spss进行籽粒灌浆拟合的详细步骤_第4页
spss进行籽粒灌浆拟合的详细步骤_第5页
已阅读5页,还剩5页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

SPSS进行灌浆拟合的详细过程一、取样:进行灌浆曲线拟合,取样非常重要,务必准确标记抽穗开花时间,同时持续较长时间的取样,一般花后间隔3-5天取样,并持续至充分的黄熟收获取样时间。二、灌浆拟合例如:水稻花后每5天取样,直至成熟,样品80℃烘箱处理48h,质量衡定后脱去颖壳称质量,并换算成千粒质量,数据如下取样时间(天)干物重(g/1000粒)50.921107.5811514.3422016.3922516.0263016.292(一)参数初始值的设定籽粒灌浆拟合实质是曲线回归。曲线回归拟合是否成功,拟合的好坏与参数初始值的定义密切相关。那么灌浆的参数初始值应该定义成多少呢?这个难道了很多人。灌浆拟合采用的是Richards方程,为W=A(1+Be-kt)-1\N。当N=1时,Richards方程就等于Logistics方程:Y=A/(1+Be-kt)。因此在进行Richards方程拟合之前可进行Logistics方程的拟合,以获得其他参数B、K、A的初始值。1、打开spss软件并输入数据,如图点击【分析】、【回归】、【曲线估计】。2、打开曲线估计对话框,将x选入自变量,y选入因变量,同时在模型处选中logistic,在【下限】框内,输入比最大y值略大一些的数值,本处可输入16.5。然后勾选上显示anova表格,如下图所示3、点击确定输出以下结果从系数表中可知b0=1.366,b1=0.756从观测值数据中选取3对观察值,x为3个等间距值,如5、10、15。3对观测值可为(x1,y1)=(5,0.921),(x2,y2)=(10,7.581),(x3,y3)=(15,14.342).按照以下公式计算出A值:A=(y22*(y1+y3)-2*y1*y2*y3)/(y22-y1*y3)=15.293然后B=b0*A=20.89,K=-lnb1=-ln0.756=0.279。因此得到初始值:A=15.293,B=20.89,K=0.279,N=1.(二)Richards方程拟合1、如图所示选择【非线性】2打开非线性对话框,将y选入因变量,同时将Richard方程按照spss的语法规范输入模型。如下图所示然后在非参数对话框内,点击【参数】,打开参数对话框,并输入之前算出的参数初始值3、点击确定,输出结果A=16.255,B=33,K=0.407,N=0.583最后得到Richards方程为:y=16.255(1+33*e-0.407*x)-1/0.583对上述Richards方程求一阶导数,可得下式灌浆速率(GR,即单位时间千粒重的增长量,gd-1)依时间x的灌浆速率方程GR=dy/dx=AkBe-kx/[N(1+Be-kx)

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论