实验6 matlab程序设计new

《实验6 matlab程序设计new》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、黔南民族师范学院数学系数学实验实验6Matlab程序设计实验目的：1．掌握脚本式M文件和函数式M文件的编写、调试和运行；2．掌握从分析问题到算法设计再到程序实现的全过程。实验内容：1．编写求[200，300]之间第一个能被21整除的整数的脚本式M文件，并给出运行的方法和结果。2．编写用5次多项式拟合函数y=sin(x)，x∈[0，2π]的脚本式M文件，要求绘图观察拟合的效果。3．编写求多项式P(x)=anxn+an-1xn-1+…+a1x+a0的值的函数式M文件(不能调用函数polyval)，要求通过实

2、例调试通过。4．编写计算第n个Fibonnaci数的函数式M文件，并应用它求第15个Fibonnaci数。5．按“题目分析→程序设计→实现算法的程序→实例应用”的步骤解决下面的问题。救援问题：假设洪水把一个村庄淹没了，只有那些比较高的房屋的屋顶露出水面，人们已经爬上屋顶等待求援。如右图，在图的原点是大本营，救生船每次从大本营出发，救了人之后将人送回大本营。图中的点代表屋顶，每个屋顶由其位置坐标和其上的人数表示。救生船每次从大本营出发，以速度50米/分钟驶向下一个屋顶，达到一个屋顶后，救下其上的所有人，每

3、人上船1分钟，船原路返回，达到大本营，每人下船0.5分钟。假设原点与任意一个屋顶的连线不穿过其它屋顶。假设屋顶数为50，给定屋顶数以及每个屋顶的坐标和人数，求出所有人都到达大本营并登陆所用的时间。并输出所有屋顶的坐标和人数。班级：10级应用2班姓名：刘浩学号：2010051318-9-黔南民族师范学院数学系数学实验6．按“题目分析→程序设计→算法描述→实现算法的程序→实例应用”的步骤解决下面的问题。奖金问题：过年了，村里要庆祝一下。村长对村里的128个村民说：做一个游戏，让每个人把出生年+月+日得到一个

4、数。例如：1968年10月28日=1968+10+28=2006。然后把这个数报上来。村里有一笔钱要作为游戏的奖金，数额为M元。如果有人报上来的数字与M相同，就把这笔钱发给这些人。如果只有一个人得奖，奖金都归这个人。如果有多于一个人得奖，则他们平分这笔钱。请写一个程序算算都有哪些人得到了奖金？得到多少？实验仪器与软件：1．CPU主频在2GHz以上，内存在512Mb以上的PC；2．Matlab2010a及以上版本。实验讲评：实验成绩：评阅教师：2011年月日班级：10级应用2班姓名：刘浩学号：201005

5、1318-9-黔南民族师范学院数学系数学实验实验6Matlab程序设计一、求[200，300]之间第一个能被21整除的整数的脚本式M文件clearforn=200:300ifrem(n,21)~=0continueendbreakendn二、用5次多项式拟合函数y=sin(x)，x∈[0，2π]的脚本式M文件clearx=0:pi/20:2*pi;y=sin(x);p=polyfit(x,y,5)x1=0:pi/30:pi*2;y1=sin(x1);y2=polyval(p,x1);三、求多项式P(x)

6、=anxn+an-1xn-1+…+a1x+a0的值的函数式M文件①程序如下functiony=polynomial(n,x)clearn=input('请依次输入多项式的次数n:')y=0;coe=zeros(1,n+1)fori=1:n+1coe(i)=input('请依次输入多项式的系数：');enddisp('你输入的多项式如下：')poly2sym(coe)x=input('请输入多项式中x的值:')disp('多项式的值为：')fori=1:n+1y=y+coe(i)*x^(n+1-i);en

7、d四、计算第n个Fibonnaci数的函数式M文件clearfunctiony=Fib(n)班级：10级应用2班姓名：刘浩学号：2010051318-9-黔南民族师范学院数学系数学实验ifn<3y=1;elsey=Fib(n-1)+Fib(n-2);end五、救援问题的解决+1．题目分析:这个问题的实际上就是一个简单的求和问题。单独计算出每个屋顶上的人到达大本营的时间，然后将他们相加即可。计算公式如下：totalTime=其中，totalTime是要求的总时间，N是屋顶数（本题中为50），speed是船

8、行驶的速度（本题中为50（米/分钟）），,是第i个屋顶的坐标，是第i个屋顶上的人数。因为要在计算后输出所有屋顶位置和人数，所以需要保存所有有屋顶和人数信息。2．程序设计⑴定义结构类型的数组一存放屋顶信息。⑵依次读入屋顶及其上的人数信息，并保存起来。⑶计算所需营救时间。⑷输出营救时间和所有屋顶坐标及人数。3．实现算法的程序clearnum=50;speen=50;up=1;down=0.5;totalTime=0;disp('请依次输入屋顶横