Matlab数学实验报告.docx

1、Matlab数学实验报告MATLAB与科学计算实验报告指导老师：易 昆 南实验时间：大二年级上学期学院：数学科学与计算技术学院专业班级：统计 0902姓名：丁 克 明学号：13040901121. 学生成绩管理 32. 函数极限问题53. 蛛网模型74. 水塔模型105. 混沌问题126.koch 曲线147. 拉格朗日与三次样条插值以及曲线拟合178. 解线性方程组209. 矩阵对角化2310.矩阵标准化 .2411. 摸球实验 2612. 钉板问题 2613. 火车问题 28一、学生成绩管理学号1304090112班级统计0902姓名丁克明指导教师易昆南实验题目学生成绩管理程序评 分1、设

2、计（实习）目的 ：1 结合实际问题展现 MATLAB在生活和学习方面的广泛应用2 学会利用 MATLAB编程并求解实际问题3 学会并运用 for 循环和 switch 结构，以及 MATLAB中已有函数如4 了解单元数组 cell 和结构数组 struct 的作用，学会创建这些数组sum.2、实验内容：编写一个学生成绩管理的函数程序：将学生成绩按五个等级划分：优秀（ 90 以上），良好（80-90 ），中等（ 70-79 ），及格（ 60-69 ），不及格（ 60 以下）对输入的学生成绩按五个等级划分后，打印学生姓名，得分，等级，并计算全班最高分、最低分、平均分。3、详细设计clearfor

3、i=1:10ai=89+i;bi=79+i;ci=69+i;di=59+i;endc=d,c;Name=input(please input name:);Score=input(please input score:);n=length(Score);Rank=cell(1,n);S=struct(Name,Name,Score,Score,Rank,Rank);for i=1:nswitch S(i).Scorecase 100S(i).Rank=满分 ;case aS(i).Rank=优秀 ;case bS(i).Rank=良好 ;case cS(i).Rank=及格 ;case dS(

4、i).Rank=不及格;endenddisp(学生姓名,得分,等级 )；for i=1:ndisp(S(i).Name,blanks(6),num2str(S(i).Score),blanks(6),S(i).Rank);ends=0;for i=1:ns=S(i).Score+s;endaverscore=s/n;t=S(1).Score;for i=1:(n-1)if(S(i).ScoreS(i+1).Score)m=S(i+1).Score;endenddisp( 平均成绩 );disp(averscore);disp( 最高分 );disp(t);disp( 最低分 );disp(m)

5、;4、实验结果：5、实验总结：这是我第一次接触 matlab 并运用 matlab 解决实际问题，在以前学习 C+的时候，解决一个问题的代码会很繁琐。但相同的问题，在 matlab 中却变得很简单。这就是 matlab 的强大之处。运用它可以方便的解决许多实际问题。知道了这一点，我决心以后会认真的学习这样一门课程。同时，它也给我们枯燥的数学公式的学习带来了动手解决实际问题的机会。是将理论与实际相结合的方法。让我们更加体会到知识的力量是强大的，我们应该更好的掌握科学技术和相关理论知识，并能够将其运用于实际生活当中。在解决这个问题的时候也确实遇到了一些问题。比如如何输入符号数组和数值数组的问题，当

6、时一直不理解。当通过查询资料，和同学讨论，最后明白了。输入符号时，符号要用单引号，而数据就不需要了。通过这样一次实践，我更加明白了动手的必要性。只有理论知识是远远不够的。所以以后一定要加强自己的动手能力，勤动手。二、函数极限问题学号1304090112班级统计0902姓名丁克明指导教师易昆南实验题目函数极限问题评 分1、设计（实习）目的 ：1 结合实际问题展现 MATLAB在生活和学习方面的广泛应用2 学会利用 MATLAB编程并求解实际问题3 学会并运用 movie 和 moviein4 理解 getframe 以及学会运用函数，以及 plotaxis 调节坐标抽结构。2、实验内容：自选函数

7、，运用plot进行绘图。并使用movie以及moviein，制作函数极限动画。3、详细设计clearm=moviein(16);hold onfor i=1:50for j=1:it1=(j-1)*5;t2=j*5;t=t1:0.01:t2;axis(0,40,-1,1)%限制动画的坐标显示大小plot(t,0,.r)plot(t,1./t,.)endm(:,i)=getframe;endmovie(m,50)clearm=moviein(16)hold onfor i=1:10for j=1:it1=(j-1)*pi/10;t2=j*pi/10;t=t1:0.01:t2;axis(0,5*p

8、i/2,-1,1)% 限制动画的坐标显示大小plot(pi-0.1,t,.k)plot(pi-0.1,-t,.k)plot(pi+0.1,-t,.k)plot(pi+0.1,t,.k)plot(t,0,.r)plot(t,sin(t),.)t1=(20-j)*pi/10;t2=(21-j)*pi/10;t=t1:0.01:t2;axis(0,5*pi/2,-1,1)% 限制动画的坐标显示大小plot(t,0,.r)plot(t,sin(t),.)endm(:,i)=getframe;endmovie(m,50)4、实验结果：5、实验总结：使用 matlab 进行绘图觉得很有意思， 简单的代码，

9、 却能绘出各种颜色的图形动画。 在学习的过程中还是有很多的问题。还有很多方面的东西不够了解，只能边尝试边查询资料，让自己能够更了解，更明白。以便绘出更好的图形。三、蛛网模型学号1304090112 实验题目班级统计蛛网模型0902姓名丁克明指导教师评 分易昆南1、设计（实习）目的 ：1.学会运用 Matlab 解决实际问题2.进一步掌握 Matlab 的一些基本操作3.通过范例体会 Matlab 的初步建模过程4学会用做动画2、实验内容：已知前两年的猪肉的产量和猪肉的价格分别为： 39 吨， 28 吨， 12 元 / 公斤， 17 元/ 公斤，根据前一年的猪肉价格影响后一年猪肉产量当年猪肉产量

10、影响当年猪肉价格的线性关系，编写程序，利用动画原理画出前十年猪肉的产量价格的动态图形。（参数设置为，c1=39,c2=28,c3=36,r1=12,r2=17,k=16) 。3、详细设计clear%c1为产量%肉价 1, r21, c2 为产量为肉价 2, k2, c3 为产量 3, r1 为 K 年后产量与肉价为%是否稳定hold offc1=39;c2=28;c3=36;r1=12;r2=17;k=16;a1=c1 1;c2 1;% 系数矩阵b1=r1,r2;% 列向量a2=r1 1;r2 1;b2=c2,c3;a=a1b1;b=a2b2;%x0(1)=c1;for n=1:30y0(n)

11、=a(1)*x0(n)+a(2); %a(1)a(2)x0(n+1)=b(1)*y0(n)+b(2);x(n)=x0(n);y(n)=x0(n+1);为矩阵元素的引用，下行类似endplot(x,y0,-g,y,y0,-b)hold onm=moviein(100);for n=1:kfor j=1:30t1=x0(n)+(j-1)*(x0(n+1)-x0(n)/30;t2=x0(n)+j*(x0(n+1)-x0(n)/30;if t2t1t=t1;t1=t2;t2=t;% 这样做方便画图elseendt=t1:0.01:t2;plot(t,y0(n),.r)% 划横线t1=y0(n)+(j-

12、1)*(y0(n+1)-y0(n)/30;t2=y0(n)+j*(y0(n+1)-y0(n)/30;if t2t1t=t1;t1=t2;t2=t;elseendt=t1:0.01:t2;plot(x(n+1),t,.r)% 划竖线endm(:,n)=getframe;endmovie(m,20)4、实验结果：ans =Columns 1 through 639.0000 28.0000 36.0000 30.1818 34.4132 31.3358Columns 7 through 1233.5739 31.9462 33.1300 32.2691 32.8952 32.4398Columns 13 through 1832.7710 32.5302 32.7053 32.5779 32.6706 3