欢迎来到天天文库
浏览记录
ID:54016856
大小:209.73 KB
页数:5页
时间:2020-04-28
《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
此文档下载收益归作者所有