matlab上机练习1.doc

实验1 MATLAB7.0基本编程实验目的： 掌握脚本和函数； 掌握matlab中的变量和M文件的流控制语句； 掌握函数的设计和实现。实验要求:给出程序和实验结果。实验内容：一、一个三位整数各位数字的立方和等于该数本身则称该数为水仙花数。输出全部水仙花数。二、从键盘输入若干个数，当输入10时结束输入，求这些数的平均值和它们之和。三、求1000，2000之间第一个能被17整除的整数。四、若一个数等于它的各个真因子之和，则称该数为完数，如6=1+2+3，所以6是完数。求1,10000之间的全部完数。五、Fibonacci数列定义如下：=1=1求Fibonacci数列的第20项。六、设计一个猜数游戏。首先由计算机产生1,100之间的随机整数，然后由用户猜测所产生的随机数。根据用户猜测的情况给出不同提示，如猜测的数大于产生的数，则显示“High”，小于则显示“Low”，等于则显示“You won”，同时退出游戏。用户最多可以猜7次。实验结果： 一、 一个三位整数各位数字的立方和等于该数本身则称该数为水仙花数。输出全部水仙花数。 程序如下：for m=100:999 m1=fix(m/100); %求m的百位数字 m2=rem(fix(m/10),10); %求m的十位数字 m3=rem(m,10); %求m的个位数字 if m=m1*m1*m1+m2*m2*m2+m3*m3*m3 disp(m) endend 153 370 371 407二、 从键盘输入若干个数，当输入0时结束输入，求这些数的平均值和它们之和。程序如下：sum=0;cnt=0;val=input(Enter a number (end in 0):);while (val=0) sum=sum+val; cnt=cnt+1; val=input(Enter a number (end in 0):);endif (cnt 0) sum mean=sum/cntend Enter a number (end in 0):45Enter a number (end in 0):65Enter a number (end in 0):342Enter a number (end in 0):45Enter a number (end in 0):67Enter a number (end in 0):12Enter a number (end in 0):0sum = 576mean =96三、 求100，200之间第一个能被21整除的整数。程序如下：for n=100:200 if rem(n,21)=0 continue end breakendnn = 105四、 若一个数等于它的各个真因子之和，则称该数为完数，如6=1+2+3，所以6是完数。求1,500之间的全部完数。for m=1:500 s=0; for k=1:m/2 if rem(m,k)=0 s=s+k; end end if m=s disp(m); endend 6 28 496五、Fibonacci数列定义如下：=1=1求Fibonacci数列的第20项。编写程序为：for i=3:20; f(1)=1; f(2)=1; f(i)=f(i-1)+f(i-2);end ff = Columns 1 through 7 1 1 2 3 5 8 13 Columns 8 through 14 21 34 55 89 144 233 377 Columns 15 through 20 610 987 1597 2584 4181 6765六、设计一个猜数游戏。首先由计算机产生1,100之间的随机整数，然后由用户猜测所产生的随机数。根据用户猜测的情况给出不同提示，如猜测的数大于产生的数，则显示“High”，小于则显示“Low”，等于则显示“You won”，同时退出游戏。用户最多可以猜7次。命令如下：a=fix(rand(1)*100+1);for k=1：7n=input(输入一个1-100的整数:n=);if na disp(High);elseif n clear;y=input(Input Fahrenheit temperature：);x=f2c(y)Input Fahrenheit temperature：70c = 21.1111x = 21.1111例3-2 输入x,y的值，并将它们的值互换后输出。 程序如下： x=input(Input x please.); y=input(Input y please.); z=x; x=y; y=z; disp(x);disp(y);Input x please.34Input y please.54 54 34例3-3 求一元二次方程ax2 +bx+c=0的根。程序如下：a=input(a=?);b=input(b=?);c=input(c=?);d=b*b-4*a*c;x=(-b+sqrt(d)/(2*a),(-b-sqrt(d)/(2*a);disp(x1=,num2str(x(1),x2=,num2str(x(2);a=?1b=?7c=?9x1=-1.6972,x2=-5.3028例3-4 计算分段函数的值。 程序如下：x=input(请输入x的值:);if x=A & c=a& c=0& c=9 disp(abs(c)-abs(0);else disp(c);end请输入一个字符OO例3-6 某商场对顾客所购买的商品实行打折销售，标准如下(商品价格用price来表示)： price200 没有折扣 200price500 3%折扣 500price1000 5%折扣 1000price2500 8%折扣 2500price eg请输入商品价格23090price = 1.9857e+004例3-7 矩阵乘法运算要求两矩阵的维数相容，否则会出错。先求两矩阵的乘积，若出错，则自动转去求两矩阵的点乘。程序如下：A=1,2,3;4,5,6; B=7,8,9;10,11,12;try C=A*B;catch C=A.*B;endClasterr %显示出错原因C = 7 16 27 40 55 72ans =Error using = mtimesInner matrix dimensions must agree.已知 ，当n=100时，求y的值。 程序如下：y=0;n=100;for i=1:n y=y+1/(2*i-1);endyy =3.2843例3-10 写出下列程序的执行结果。 s=0; a=12,13,14;15,16,17;18,19,20;21,22,23; for k=a s=s+k; enddisp(s); 39 48 57 66例3-14 编写函数文件求半径为r的圆的面积和周长。 编写m文件并命名为fcircle.m如下： function s,p=fcircle(r)%FCIRCLE calculate the area and perimeter of a circle of radii r%r 圆半径%s 圆面积%p 圆周长s=pi*r*r;p=2*pi*r;circle(2)ans = 12.5664例3-15 利用函数文件，实现直角坐标(x,y)与极坐标(,)之间的转换。 函数文件tran.m：function rho,theta=tran(x,y)rho=sqrt(x*x+y*y);theta=atan(y/x); 在命令窗体输入命令：x=input(Please input x=:);y=input(Please input y=:);rho,the=tran(x,y);rhothePlease input x=:5Please input y=:7rho = 8.6023the =0.9505例3-16 利用函数的递归调用，求n！。 n!本身就是以递归的形式定义的：显然，求n!需要求(n-1)!，这时可采用递归调用。递归调用函数文件factor.m如下：function f=factor(n)if n factor(67)ans = 3.6471e+094例3-17 nargin用法示例。 函数文件charray.m：function fout=charray(a,b,c)if nargin=1 fout=a;elseif nargin=2 fout=a+b;elseif nargin=3 fout=(a*b*c)/2;end在命令窗体中输入：x=1:3; y=1;2;3; charray(x)ans = 1 2 3 charray(x,y)ans