MATLAB在动物科学实验数据处理中的应用.pdf

《MATLAB在动物科学实验数据处理中的应用.pdf》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、第16卷第4期湖南文理学院学报（自然科学版）VoI.16No.42004年12月JournaIofHunanUniversityofArtsandScience（NaturaIScienceEdition）Dec.2004文章编号：1672-6146（2004）04-0043-03MATLAB在动物科学实验数据处理中的应用!王京仁，赵东海，李淑红，金霞（湖南文理学院生命科学系，湖南常德415000）摘要：在动物科学实验数据处理中，往往需要对所测数据2k!1=!，!2=!，⋯，!k=!，则上式化为：进行统计分析，建立拟合方程.将MAT

2、LAB语言编写程序，用"y=60+61!1+62!2+63!3+⋯6k!k（2）于实验数据处理，计算结果准确可靠、符合实验要求.通过实此式为一般多元线性回归方程，可采用矩阵方法求例介绍MATLAB在求解多项式回归方程中的应用，具有方法解.即由简捷、直观等特点.关键词：MatIab语言；实验数据处理；多项式回归é1!11!21⋯!k1ù中图分类号：O332文献标示码：Aê1!12!22⋯!k2úX=êú=ê⋯úêú在动物科学专业的实验研究中，数据多，处理麻ë1!1n!2n⋯!knû烦.常用TSP、SAS、SPSS等的软件进行统计分析.

3、这1!2⋯ké11!11!11ùé"1ù里利用MatIab强大的统计分析功能（仅统计工具箱ê2kúêúê1!12!12⋯!12úê"2úStatisticTooIboX中的功能函数就达200多个），对动êú和Y=⋯êú物科学专业《生物统计附试验设计》教材中的两个典êúêúë1!2⋯kûë"û1n!1n!1nn型实例，通过几行简单的命令编程，即可解出多项式直接用MATLAB中提供的专用拟合函数poIyfit，其调回归方程，展示了MatIab在数据处理中简洁、直观与用格式如下：快捷的特点.pIoyfit（X，Y，n）其中X、Y为拟合数据

4、，n为拟1多项式回归模型及回归方程式合多项式的阶数.coef=poIyfit（X，Y，n），则coef（1）=60，在处理动物科学的实验数据时，经常遇到研究资料的一个依变量与一个自变量间的曲线关系很难coef（2）=61，⋯，coef（n+1）=6k.确定时，这时可以使用一个解析函数来描述依变量注意上式对k阶的多项式会有k+1项的系数.与自变量之间的关系.pIoyvaI函数可计算出二次多项式函数的预测值，由2［2］一元k次多项式的一般形式为：自编相关系数R可计算出相关系数."23kyk=60+61!+62!+63!+⋯6k!（1）3

5、实例程序这是一个具有k-1个弯曲（k-1个极值）和k-2个拐点的曲线.式中6，6，6，⋯，6是待定的系数；例1给动物口服某种药物A1000mg，每间隔1012k小时测定血药浓度（g/mI），得表1数据，试建立血药k是多项式的次数.浓度（依变量"）对服药时间（自变量）回归方!2多项式回归分析的方法程［3］.多项式回归问题一般仍采用类似于多元线性回表1服药时间与血药浓度关系归的方法求解多项式回归的统计数.令上式中：服药时间!血药浓度"服药时间!血药浓度"/1-1/1-1/g·mI/g·mI121.89666.36收稿日期：2004-07

6、-21247.13750.34修回日期：2004-09-24361.86825.31!基金项目：湖南文理学院院教研教改课题资助（JG2003-27）.作者简介：王京仁（1963-），男，副教授，硕士，研究方向为动物科470.7893.17学饲养及繁殖.572.8144湖南文理学院学报（自然科学版）2004年（1）输入数据并保存为p225exceSe9-3.mat文件系数x=119，由散点图知，变量间的关系更接近与二次函数，y=［21.89，47.13，61.86，70.78，72.81，66.36，拟合次数以2次为好.故取l=2，程

7、序为：50.34，25.31，3.17］，cIearaII；（2）作散点图，确定初步确定拟合函数的性质及x=1：1：12，多项式的次数.y=［3833.433811.583769.473565.74newx=’服药时间1’，newy=’血药浓度/g·3481.993372.823476.763466.22-1mI’，3395.423807.083817.033884.52］；pIo（tx，y，’0’），［b，s］=poIyfi（tx，y，2），xIabe（I’newx’），yIabe（I’newy’），运行得：（3）执行poIyfi

8、t命令b=［b，s］=poIyfi（tx，y，2），1.0e+003!得结果：0.0158-0.20494.1172b=-3.762434.8269-8.3655s=即回归方程为：y^2=-3.7624!+34.8269!-R：［3x