案例分析报告.docx

案例分析报告.docx

《案例分析报告.docx》由会员分享，可在线阅读，更多相关《案例分析报告.docx（18页珍藏版）》请在冰豆网上搜索。

案例分析报告

案例分析报告

2011—2012学年第2学期

合肥学院卓越工程师班

案例分析报告

课程名称：

工程应用数学B

案例名称：

数据拟合

专业班级：

11级自动化卓越班

起止时间：

8-11周

组别：

第五组

指导教师：

王贵霞成绩：

90

一.小组成员（具体分工）

姓名

学号

具体分工

胡德文

1105011031

实验步骤，实验程序

王佐军

1105012020

实验结果，实验总结

朱自强

1105012015

提炼题目，问题分析

二.案例描述

六．方案实施及结果（实验程序等形式）

案例一：

1.使用三次多项式拟合图形程序代码：

t=[0246810121416182022];

h=[18.520.422.825.630.927.422.015.313.512.713.515.4];

Y=polyfit（t,h,3）;

tt=0:

0.5:

24;

hh=polyval（Y,tt）

plot（tt,hh,'-b',t,h,'k:

'）,text（9.5,29,'实测温度'）;

text（15,16.5,'拟合温度'）;

xlabel（'时间（h）'）;

ylabel（'温度（℃）'）;

title（'从0:

00到23:

00之间的温度变化'）

2.拟合图像：

3.使用五次多项式程序代码为：

t=[0246810121416182022];

h=[18.520.422.825.630.927.422.015.313.512.713.515.4];

Y=polyfit（t,h,5）;

tt=0:

0.5:

24;

hh=polyval（Y,tt）

plot（tt,hh,'-b',t,h,'k:

'）,text（9.5,29,'实测温度'）;

text（15,16.5,'拟合温度'）;

xlabel（'时间（h）'）;

ylabel（'温度（℃）'）;

title（'从0:

00到23:

00之间的温度变化'）

4.五次多项式拟合图形：

案例二

1.程序代码：

clear;closeall;

x=[0.250.511.523468];

y=[19.2118.1515.3614.112.899.327.455.243.01];

X=0:

0.25:

8;

Y=interp1（x,y,X,'cubic'）;

plot（x,y,'ko-',X,Y,'r*:

'）

xlabel（'时间（t）'）;

ylabel（'血药浓度c（ug/ml）'）;

title（'注射药后血药浓度随时间变化'）

2.拟合图形如下：

从拟合图像的曲线变化可以看出，血药浓度（单位体积血液中的药物含量）随着在人体内时间的增加而减少，并且开始减少的快，后来减少的慢。

（1）.设计程序：

x=[0.250.511.523468];

y=[19.2118.1515.3614.112.899.327.455.243.01];

X=-2:

0.25:

8;

Y=interp1（x,y,X,'cubic'）;

plot（x,y,'ko-',X,Y,'r*:

'）;

holdon;

xx=-2:

0.01:

8;

yy1=25;

yy2=10;

plot（xx,yy1）,text（3,26,'yy1=25'）;

plot（xx,yy2）,text（0,11,'yy2=10'）

xlabel（'时间（t）'）;

ylabel（'血药浓度c（ug/ml）'）;

title（'注射药后血药浓度随时间变化'）

（2）.拟合图像：

注入3小时内药效都在有效范围内，给定药物的最小有效浓度和最大治疗浓度，设计给药方案：

向负方向延长基本成正比可以求出比例系数约为15，所以每次注射375mg，持续时间约为4小时。

案例三程序代码：

1.用一维插值拟合曲线程序代码：

clear;closeall;

r=rand

（1）-0.5;

x=[234579121417212756];

y=[35+r42+r47+r53+r59+r65+r58+r73+r76+r82+r86+r99+r];

X=0:

0.5:

60;

Y=interp1（x,y,X,'cubic'）;

plot（x,y,'ko-',X,Y,'r.:

'）

xlabel（'养护时间'）;

ylabel（'抗压强度'）;

2.用多项式拟合曲线图形如下：

3.用非线性回归模型程序代码：

（1）建立非线性回归模型，对得到的模型和系数进行检验，这是曲线拟合中非线性方程用模型CustomExponential——在其中输入y=a+k1*exp（m*x）+k2*exp（-m*x）。

（2）在主界面先输入以下数据：

r=rand

（1）-0.5;

x=[234579121417212756];

y=[35+r42+r47+r53+r59+r65+r58+r73+r76+r82+r86+r99+r];X=0:

0.5:

60;

（3）进入数据拟合工具箱，进行拟合数据设置，新建Newfit，选择Typeoffit,选择自定义CustomEquations进行拟合。

拟合结果如下：

Generalmodel:

f（x）=a+k1*exp（m*x）+k2*exp（-m*x）

Coefficients（with95%confidencebounds）:

a=91.64（70.53,112.8）

k1=0.0829（-0.1122,0.278）

k2=-62.8（-79.11,-46.5）

m=0.08254（0.01263,0.1525）

（4）数据拟合工具箱拟合过程截图：

①创建所要求的函数：

②设置好函数形式后，点击新建Apply即可。

（5）数据拟合工具箱中，非线性回归模型拟合图像如下：

结果：

a=91.64k1=0.0829k2=-62.8m=0.08254

七．总结（围绕心得体会、创新之处、改进方案等方面）

又一次MATLAB实训,自学的又一个阶段结束了。

对MATLAB的认识也增加了不少，从生疏到熟悉，这都是自学的成果。

记得在这段期间曾出现过好多的错误与失败，但是我们没有气馁，尽管这次的题目比较难，我们也不会超越上一次，有比较好的成绩，但是我们确实付出了自己的汗水，也学到了好多实实在在的知识。

从初识MATLAB到现在，对于这一门强大的语言也并不是我们课下这一点时间就能学习的非常好的。

因此对于MATLAB的学习我们依然任重而道远。

接下来谈谈这次我们分析报告的完成之路：

首先，对于老师给我们的题目，当时就傻眼了，这么难，我们该从哪里出发啊？

结果我们分析了四五天，才初步有了规划，从三题中选择了两个难的题目，然后自己查找了一题。

大多数时间还是用于老师那两题的分析。

第二步，在选择好题目之后，我们的进展就非常快了，分析题意，选择拟合方案，编写程序代码，最后再完善Word报告。

当然这只是我们整体的步骤，直至有这样清晰的步骤远远不够，因此我详细讲一下我们的实训过程。

第三步，MATLAB功能强大我也能从我所掌握的方法叙述实训过程：

1.对于编写程序代码，我们仍然选择在MATLAB主界面中打开file再点击new接着在new的菜单中选择script，然后再在script中输入程序代码，这样也很方便，而且只需点击菜单区的生成键，程序对应得图形便会出现。

并且这次我们选择了图形标注于坐标标注，再熟悉了图形标注于坐标标注的使用方法之后，便用于本次实训的大部分图形，这样标注更能使图形曲线的表示意思显而易见。

2.本次实训我们还用到了holdon命令控制，其作用是保持原有图形，主要的过程还是将题目模型化，准备将其转化成什么类型的曲线，这样才更有利于拟合。

3.对于老师布置的第三题，在给定的非线性回归模型模拟曲线，我们也用到了功能强大的数据拟合工具箱，在这里面，可以任意形式的曲线。

快捷方式就是在主界面输入>>cftool即可，对于数据拟合工具箱上面已经说明过了。

4.这些基本上就是我们实训过程中所用到的MATLAB功能，当然本次试验我们也用了基本的对比思想，拟合数据解答问题等。

在步骤说明的过程中我们还引用了红蜻蜓抓图软件的使用，截取自己制作过程的图片，这样也有利于报告的完整性与生动性。

第四步，在编程的时候，我们的改进方案有，对图形标注的说明，对使用拟合函数的修改，对颜色的修改，图形大小的选择，最重要还是我们组成员在写报告时一起讨论，在这其中我们不仅相互交流学到了知识，而且增进了我们之间的感情，在这时期尤其有利于对方案的修改于抉择，不仅对我们现在有用，对将来有一定影响。

最后，在本次MATLAB的学习过程中，我们有进一步认识了它，但这对MATLAB强大的功能的影响下是永远不够的。

从资料书与网上查询可知MATLAB的作用也是取之不尽.用之不竭的。

所以在今后，我们依然要多加学习，相信能够获得更多。

而且在本次学习MATLAB的使用过程中，我们不仅学习到了MATLAB的知识，更加认识到认真仔细学习的重要性，我相信这些对我们将受益终身。

以上是我们这次实训的大概经过，也是我们的心得体会。

相信不久的将来我们一定能掌握MATLAB的大部分功能。

八．教师评语

能够用利用三次和五次多项式分别拟合温度观测问题,比较其优劣；对图形标注的说明，对使用拟合函数的修改，对颜色的修改，图形大小的选择，通过拟合图像的曲线变化解决药品临床实验问题；一维插值拟合曲线和多项式拟合曲线及非线性回归模型处理混凝土的抗压强度，达到本次实验的目的，是一篇水平很高的实验报告。

教师签名：

王贵霞

2012年5月2日

2011—2012学年第2学期

合肥学院卓越工程师班

案例分析报告

课程名称：

工程应用数学B

案例名称：

随高度变化的压强

专业班级：

机制卓越班

起止时间：

2012/4/13---2012/5/12

组别：

第四组

指导教师：

张霞

一.小组成员（具体分工）

姓名

学号

具体分工

成绩

张林

1106012036

查阅相关资料调试检测

85

韩贤跃

106012026

查阅相关资料计算核对检测

80

胡利有

1106013003

查阅相关资料编写程序调试检测

90

二.案例描述

随高度变化的大气压强

大气压强/随高度

变化的理论公式为

。

为验证这一公式，测得某地大气压强随高度变化的一组数据如表2-2，试用拟合法进行计算并绘图。

表2-2大气压强随高度变化的一些数据

高度

0

300

600

1000

1500

2000

压强

0.9689

0.9322

0.8969

0.8519

0.7989

0.7491

三.提炼问题（将实际问题数学化）

表2-2大气压强随高度变化的一些数据

高度

0

300

600

1000

1500

2000

压强

0.9689

0.9322

0.8969

0.8519

0.7989

0.7491

利用上述数据验证是否满足公式

?

四.问题分析

结合该题目与所学的内容，我们决定用polyfit,polyval等函数解决本组问题。

考虑我到们所用的函数是解决数据之间的函数关系是线性关系而我们要解决的函数是非线性关系，于是我们将函数

化成

的形式期中

;在令

得到新函数

。

最后由我来编写程序，主要采取两种方案，其他的成员负责调试。

五．解决方案（提供多种备选）

方案一：

利用matlab软件自身带的功能进行拟合。

方案二：

编写程序进行拟合并验证。

六．方案实施及结果（实验程序等形式）

方案一：

x=[0300600100015002000]';

m=[0.9689,0.9322,0.8969,0.8519,0.7989,0.7491]';

y=log（m）;

plot（x,y,'o'）,gridon

xlabel（'x'）

ylabel（'y'）

方案二：

x=[0300600100015002000]';

m=[0.9689,0.9322,0.8969,0.8519,0.7989,0.7491]';

y=log（m）;

p2=polyfit（x,y,1）;

disp（'一阶拟合'）,f1=poly2str（p2,'x'）

x1=[0,300,600,1000,1500,2000];

y1=polyval（p2,x1）;

wch=abs（y-polyval（p2,x））./y

pjwch=mean（wch）

minwch=min（wch）

maxwch=max（wch）

plot（x,y,'rp',x1,y1）

xlabel（'x'）

ylabel（'y'）

legend（'拟合点','一次拟合'）

gridon

七．总结（围绕心得体会、创新之处、改进方案等方面）

从这次的实验中，我们对matlab软件有了更深刻的了解，特别是它在数学方面的应用。

学好一门软件光靠理论知识是不行的，必须与计算机结合起来，多动手。

在这次的实验中我们学会了用matlab绘图和和用它进行高级数值运算，更重要的是除了我们学会了一门课外知识外，我们认识到团队合作的重要性，仅凭个人的能力是很难完成此次任务，所以我们在以后的学习中多于别人交流与沟通，使我们有更高的进步。

我们这次实验创新之处在于将数据之间的非线性关系转换成线性关系和利用软件自带的功能进行曲线拟合免了写程序的烦劳，而且在方案二中有误差分析，使结果更加准确；方案一巧妙利用软件本身的功能进行解题，给我们带来了很大的方面。

我们的实验还有许多改进之处。

我们没有用解决分线性关系的函数解决我们的课题，在用线性关系函数解决我们我们的问题时，有许多复杂的计算。

我们如果能编出解决指函数与对数函数互换，复杂的计算程序或用非线性函数解决我们的课题，我们的实验会更加的完美。

八．教师评语

能正确提炼出数学问题，利用所学的知识分析问题，给出两种方案解决问题，在总结中对两种方案进行比较，并指出了改进的方向。

应用数学解决问题的能力得到了有效的锻炼，汇报流畅，达到了本次自主学习的目的。

教师签名：

张霞

2012年5月11日