MATLAB与科学计算
实验报告
指导老师：易昆南
实验时间：大二年级上学期
学院：数学科学与计算技术学院
专业班级：统计0902
姓名：丁克明
学号：1304090112
目录
1.学生成绩管理 (3)
2.函数极限问题 (5)
3.蛛网模型 (7)
4.水塔模型 (10)
5.混沌问题 (12)
6.koch曲线 (14)
7.拉格朗日与三次样条插值以及曲线拟合 (17)
8.解线性方程组 (20)
9.矩阵对角化 (23)
10.矩阵标准化 (24)
11.摸球实验 (26)
12.钉板问题 (26)
13.火车问题 (28)
一、学生成绩管理
学号1304090112班级统计0902姓名丁克明指导教师易昆南实验题目学生成绩管理程序评分
1、设计（实习）目的：
1结合实际问题展现MATLAB在生活和学习方面的广泛应用
2学会利用MATLAB编程并求解实际问题
3学会并运用for循环和switch结构，以及MATLAB中已有函数如sum
4了解单元数组cell和结构数组struct的作用，学会创建这些数组.
2、实验内容：
编写一个学生成绩管理的函数程序：将学生成绩按五个等级划分：优秀（90以上），良好（80-90），中等（70-79），及格（60-69），不及格（60以下）对输入的学生成绩按五个等级划分后，打印学生姓名，得分，等级，并计算全班最高分、最低分、平均分。
3、详细设计
clear
for i=1:10
a{i}=89+i;
b{i}=79+i;
c{i}=69+i;
d{i}=59+i;
end
c=[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:n
switch S(i).Score
case100
S(i).Rank='满分';
case a
S(i).Rank='优秀';
case b
S(i).Rank='良好';
case c
S(i).Rank='及格';
case d
S(i).Rank='不及格';
end
end
disp(['学生姓名','得分','等级'])；
for i=1:n
disp([S(i).Name,blanks(6),num2str(S(i).Score),blanks(6),S(i).Rank]);
end
s=0;
for i=1:n
s=S(i).Score+s;
end
averscore=s/n;
t=S(1).Score;
for i=1:(n-1)
if(S(i).Scoret=S(i+1).Score;
end
end
m=S(1).Score;
for i=1:(n-1)
if(S(i).Score>S(i+1).Score)
m=S(i+1).Score;
end
end
disp(['平均成绩']);
disp([averscore]);
disp(['最高分']);
disp(t);
disp(['最低分']);
disp(m);
4、实验结果：
5、实验总结：
这是我第一次接触matlab并运用matlab解决实际问题，在以前学习C++的时候，解决一个问题的代码会很繁琐。但相同的问题，在matlab中却变得很简单。这就是matlab的强大之处。运用它可以方便的解决许多实际问题。
知道了这一点，我决心以后会认真的学习这样一门课程。同时，它也给我们枯燥的数学公式的学习带来了动手解决实际问题的机会。是将理论与实际相结合的方法。让我们更加体会到知识的力量是强大的，我们应该更好的掌握科学技术和相关理论知识，并能够将其运用于实际生活当中。
在解决这个问题的时候也确实遇到了一些问题。比如如何输入符号数组和数值数组的问题，当时一直不理解。当通过查询资料，和同学讨论，最后明白了。输入符号时，符号要用单引号，而数据就不需要了。
通过这样一次实践，我更加明白了动手的必要性。只有理论知识是远远不够的。所以以后一定要加强自己的动手能力，勤动手。
二、函数极限问题
学号1304090112班级统计0902姓名丁克明指导教师易昆南实验题目函数极限问题评分
1、设计（实习）目的：
1结合实际问题展现MATLAB在生活和学习方面的广泛应用
2学会利用MATLAB编程并求解实际问题
3学会并运用movie和moviein函数，以及plot结构。
4理解getframe以及学会运用axis调节坐标抽
2、实验内容：
自选函数，运用plot进行绘图。并使用movie以及moviein，制作函数极限动画。
3、详细设计
clear
m=moviein(16);
hold on
for i=1:50
for j=1:i
t1=(j-1)*5;
t2=j*5;
t=t1:0.01:t2;
axis([0,40,-1,1])%限制动画的坐标显示大小
plot(t,0,'.r')
plot(t,1./t,'.')
end
m(:,i)=getframe;
end
movie(m,50)
clear
m=moviein(16)
