MATLAB程序设计.doc

MATLAB程序设计专题实验报告学号：0901010122姓名：史永峰班级：工程力学0901指导老师：周彬实验地点：数学建模实验室实验日期：2010-12-222010-12-28联系电话2、一位长寿老人生于19世纪90年代，有一年他发现自己年龄的平方刚好等于当年的年份，问这位老人出生于哪一年。分析：由年龄平方等于岁数以及现在的年代知此人年龄不超过50程序如下：for i=1:50 n=i2; if n1890 break endend 结果如下：4、求1!+2!+3!+4!+5!程序如下：s=0.0;sum=1;for i=1:5 sum=sum*i; s=s+sum; ends 结果如下：6、设已有数学模型t=a*n2+b*n,数据（n,t）的十个观测值分别为：t为0 20 40 60 80 100 120 140 160 184n为0 1141 2019 2760 3413 4004 4545 5051 5525 6061 ，使用最小二乘法原理估计a,b的值，写出主要的数学推导公式，程序和最后的结果。程序如下：n0=0,1141,2019,2760,3413,4004,4545,5051,5525,6061;t0=0,20,40,60,80,100,120,140,160,184;m=length(n0);p=2; %多项式阶数N=zeros(m,p+1); %构成数据阵Nfor k=1:p N(:,p-k+1)=(n0.k);endN(:,p+1)=ones(m,1);aT=(Nt0); %计算多项式系数disp(a=) %输出a、b的值disp(aT(1)disp(b=)disp(aT(2) 结果如下：8、用筛选法求100以内素数。分析：素数是大于1，且除1和它本身以外不能被其他数整除的整数。要找出2到m的全部素数，首先划去2的倍数（不包括2）首先划去3的倍数（不包括3），由于4被划去，再找出5的倍数（不包括5）,直至划去的倍数，剩下的就是素数。程序如下：m=input(m=); /*从键盘输入m*/p=1:m;p(1)=0;for i=2:sqrt(m) for j=2*i:i:m p(j)=0; endendn=find(p=0);p(n) 结果如下：9710、从键盘上输入一个正整数n(n20)，求n的阶乘并输出结果。（要求结果精确显示，不能表示成科学计数法）程序如下：n=input(n=);sum=1;for i=1:n sum=sum*i;enddisp(sum) 结果如下：二实验任务：MATLAB程序设计专题实验报告三本次实验内容涉及到了以下知识点：1. 循环的运用2. 阶乘3. 最小二乘法原理4. 筛选法求自然数的素数四实验总结：通过本次上机室验收获颇丰，知道matlab语言一些格式、数据类型和基本程序，大体了解了matlab软件，知道如何用运matlab输出输入，编写一些简单的问题和基本程序以及计算，不过实验过程中也存在一些问题和不足，有时候甚至是一个分号导致程序运行错误，以及文件不再工作区等等，知道以后该如何学习这门软件，向什么方向发展，matlab不仅仅是高效率高级语言，而且是图像、科学计数、程序语言共同来表达的数学高级语言。相比c/