数学建模论文Word格式.docx
- 文档编号:22511553
- 上传时间:2023-02-04
- 格式:DOCX
- 页数:10
- 大小:18.74KB
数学建模论文Word格式.docx
《数学建模论文Word格式.docx》由会员分享,可在线阅读,更多相关《数学建模论文Word格式.docx(10页珍藏版)》请在冰豆网上搜索。
马元凯
201409822姓名:
李炳毅
日期:
2016年12月15日
TOC\o"
1-1"
\h\uHYPERLINK\l_Toc30914摘要PAGEREF_Toc309141
HYPERLINK\l_Toc21875一、问题描述PAGEREF_Toc218751
HYPERLINK\l_Toc26869二、问题分析PAGEREF_Toc268692
HYPERLINK\l_Toc13851三、模型假设PAGEREF_Toc138513
HYPERLINK\l_Toc27219四、模型建立PAGEREF_Toc272193
HYPERLINK\l_Toc8865五、模型求解PAGEREF_Toc88656
HYPERLINK\l_Toc17160六、模型分析PAGEREF_Toc171606
HYPERLINK\l_Toc27977七、模型改进PAGEREF_Toc279779
HYPERLINK\l_Toc21970八、模型评注与推广PAGEREF_Toc2197011
HYPERLINK\l_Toc2372九、参考文献PAGEREF_Toc237212
HYPERLINK\l_Toc32507十、附录PAGEREF_Toc3250712
摘要
某医药公司的新药研究部门为了掌握一种新止痛剂的疗效,进行了相应的研究,根据医药公司给出的相关实验数据,通过MATLAB软件进行数据处理工作,根据拟合曲线和散点图建立基本模型(模型
(1)).用MATLAB软件解出该基本模型,再利用残差图两次剔除数据并且回归,得到最佳模型(模型
(2)),该模型的拟合度为85.14%.拟合度不是很高,针对不同的性别,引入独立变量的交互作用,对模型进行改进,得到两个拟合度较高的模型(4)和(5),针对男性建立的模型,拟合度为90.87%,针对女性建立的模型,拟合度为97.72%.
本文给出的模型较为科学,若实验数据真实可信,则可以进行药物的推广.
关键字:
MATLAB软件,残差图,拟合度,统计回归模型
问题描述
一个医药公司的新药研究部门为了掌握一种新止痛剂的疗效,设计了一个药物试验,给患有同种疾病的病人使用这种新止痛剂的以下4个剂量中的某一个:
2g,5g,7g和10g,并记录每个病人病痛明显减轻的时间(以分钟计).为了解新药的疗效与病人性别和血压有什么关系,试验过程中研究人员把病人按性别及血压的低、中、高三档平均分配来进行测试.通过比较每个病人血压的历史数据,从低到高分成3组,分别记作0.25,0.50和0.75.实验结束后,公司的记录结果见下表1(性别以0表示女,1表示男).
请你为该公司建立一个数学模型,根据病人用药的剂量、性别和血压组别,预测出服药后病痛明显减轻的时间.
表1
问题分析
一般来说,药物的疗效可以直观的用服药后病痛明显减轻的时间来衡量.在新药推广中,医药公司的新药研究部门设计了一种药物给患有同种疾病的病人使用后,根据病人的用药剂量、性别和血压组别,预测病痛减轻时间的多少来预测止痛药的疗效,这是一个统计回归问题,针对这个问题,我们需要给出合理的假设,尤其是变量的选取,进而预测出药物的疗效.
先将因变量分别与变量进行单独分析,得出两者间的大致函数关系,进一步整合这些关系,得出一个因变量与各个变量间的关系函数模型并进行求解,得出最终结论。
模型假设
假设病人只服用了新型止痛药,未服用其它药物.
假设题中给出的实验数据真实可信,误差很小.
假设24名病人都是在服用新型止痛药的人群中随机选取的.
假设病人在实验阶段吃的食物对新型止痛药无影响.
假设模型中出现的符号含义如下表2所示.
表2
模型建立
为了大致地分析y与,,之间的关系,首先利用表1-1的数据分别作出y对,和的散点图(见图4-1,图4-2和图4-3的圆点).
如图1为y对的散点图,图2为y对的散点图,图3为y对的散点图.
图1
图2
图3
由上图可知:
y对可用二次函数拟合,拟合后如图4所示.
图4
根据对散点图图1和图4的分析可得出y对的二次函数模型;
根据图2可得出y对的线型模型;
根据图3可得出y对的线型模型;
结合上述的3个模型可建立如下
(1)的多元线性回归模型:
(1)
上式右端的,,称为回归变量(自变量),给定,,时,病痛减轻时间y的平均值为;
由表1的数据估计,影响y的其他因素作用都包含在随机误差中,如果模型选的合适,应大致服从均值为0的正态分布.
模型求解
直接利用MATLAB统计工具箱中的regress求解,进一步求出回归系数估计值及其置信区间(置信水平=0.05)、检验统计量R2,F,p,s2,详细数据见表3所示.
表3
由表3得出的模型为:
模型分析
表3显示,R2=0.8275指因变量y(病痛减轻时间)的82.75%可由模型确定,F值远远超过F检验的临界值,p远小于,因而模型
(1)从整体来看是可用的.表3的回归系数给出了模型
(1)中参数的估计值,,的置信区间包含零点,因而对这两个系数的解释是不可靠的,所以需要残差分析.
首次回归所得残差图如下图5所示,可以看出,第3个和第24个数据存在异常,剔除,进行第二次回归.
图5
第二次回归所得残差图如下图6所示,可以看出,第5个数据存在异常,剔除,进行第三次回归.
图6
第三次回归所得残差图如下图7所示.
图7
第三次回归的结果如表4所示.
表4
由图7可知数据没有异常项,因此模型基本可用,此时得出的最佳模型应为:
(2)
易知因变量y(病痛减轻时间)的85.14%可由模型确定.
模型改进
从以上的分析可以看出,模型的拟合度最大为85.14%,拟合度不是很高,需要进行模型的改进;
因为服药期间要考虑生理反应,而性别对生理反应有直观的影响,故改进的模型需要对男女分开进行讨论;
模型
(2)中回归变量和对因变量y的影响是相互独立的.根据直觉和经验可以猜想和之间的相互作用会对y有影响,不妨简单地用和的乘积代表它们的交互作用,于是在模型
(1)中增加一项,得到
(3)
针对男性,利用表1的数据估计模型(3)的系数,利用MATLAB得到如下表5的结果.
表5
根据表5,可得出此时的最佳模型为:
(4)
对表5的数据进行残差分析,得到如下图8的结果.
图8
从图8可以看出数据无异常项,则由此得出的模型(4)基本可信,病痛减轻时间的90.87%可由模型确定.
针对女性,利用表1的数据估计模型(3)的系数,利用MATLAB得到如下表6的结果.
表6
对表6的数据进行残差分析,得到如图9所示的结果,
图9
从图9可以看出,第8项数据存在异常,但影响不大,所以不用剔除数据,模型基本可信.根据表5得出此时的最佳模型为:
(5)
病痛减轻时间的97.72%可由模型确定.通过以上分析可以确定,经过对模型的改进,针对不同的性别,得出的两种模型(4)和(5)的拟合程度都达到了百分之九十以上,达到了模型改进的目的.
模型评注与推广
评注从这个实例我们可以看出,建立回归模型要根据已知的数据,从常识和经验进行分析,辅以作图,决定最终的函数模型.用MATLAB软件求解后,作统计分析R2的值是对模型的直观评价,决定模型的拟合程度.本次建立的模型,优点在于,由简到繁,先单独考虑变量对结果的影响,再综合考虑,最后引入交互项对模型进行改进,为病痛减轻时间的预测提供了简洁方便的工具,用简单的形式表示出了事物之间复杂的关系;
缺点在于,该模型不能实际的去观测数据,而是使用实验给出的数据,存在一定的误差,还有改进模型中针对不同性别得出的模型拟合程度差别较大,不是很理想;
其实在模型建立时,还可以引进其他的数据,如心率等因素,使模型更加真实可信.
推广我们建立的这个预测模型还可以用于软件开发人员的薪金预测或其它的薪资预测模型.
参考文献
[1]姜启源,谢金星,叶俊.数学模型.第4版[M].高等教育出版社,2011.
[2]宋来忠,王志明,数学建模与实验,北京:
科学出版社,2005
[3]姜启源,谢金星,叶俊.数学模型(第四版)习题参考解答[M].高等教育出版社,2011.(第十章习题参考答案)
附录
1、图1,图2,图3,图4的MATLAB程序:
y对的散点图
>
x1=[222222555555777777101010101010];
y=[354355474357262728292229191114232022138327265];
scatter(x1,y,'
r'
);
x2=[000111000111000111000111];
scatter(x2,y,'
x3=[0.250.500.750.250.500.750.250.500.750.250.500.750.250.500.750.250.500.750.250.500.750.250.500.75];
scatter(x3,y,'
y对的拟合曲线
p=polyfit(x1,y,2);
x1x1=linspace(min(x1),max(x1));
yy=polyval(p,x1x1);
plot(x1,y,'
o'
x1x1,yy);
2、表3的MATLAB程序:
x=[ones(24,1),x1'
x2'
x3'
(x1.^2)'
];
[b,bint,r,rint,stats]=regress(y'
x)
图5,图6,图7的MATLAB程序:
首次回归分析
x,0.05)
rcoplot(r,rint)
第二次回归分析
x1=[222225555557777771010101010];
x2=[0011100011100011100011];
x3=[0.250.500.250.500.750.250.500.750.250.500.750.250.500.750.250.500.750.250.500.750.250.50];
y=[354347435726272829222919111423202213832726];
x=[ones(22,1),x1'
第三次回归
x1=[22225555557777771010101010];
x2=[001100011100011100011];
x3=[0.250.500.250.500.250.500.750.250.500.750.250.500.750.250.500.750.250.500.750.250.50];
y=[3543474326272829222919111423202213832726];
x=[ones(21,1),x1'
模型改进后,相关的MATLAB程序
表5,图8的MATLAB程序
x1=[222555777101010];
x3=[0.250.500.750.250.500.750.250.500.750.250.500.75];
y=[47435729222923202227265];
x=[ones(12,1),x1'
(x1.*x3)'
表6,图9的MATLAB程序
x1=[222555777101010];
x3=[0.250.500.750.250.500.750.250.500.750.250.500.75];
y=[3543552627281911141383];
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 数学 建模 论文
![提示](https://static.bdocx.com/images/bang_tan.gif)