MATLAB 实验7分析和总结.docx

MATLAB实验报告

一、实验目的

熟悉和掌握MATLAB的m文件程序和程序设计。二、实验内容

for循环结构

fori=1:10;

x(i)=i;

end;x

这是通过“for”循环产生从1到10的是个数字。执行结果如下图：

while循环结构

Fibonaci数组的元素满足Fibonacci规则：a ?a ?a ,(k?1,2, );且a ?a ?1。

k?2 k k?1 1 2

现要求该数组中第一个大于10000的元素。输入命令后执行结果如下图：a(1)=1;a(2)=1;i=2;。

whilea(i)=10000

a(i+1)=a(i-1)+a(i);i=i+1;

end;i,a(i)

if-else-end结构。这种结构就是对于多种可能的情况进行不同情况的划分，然后对

不同情况进行编辑执行不同的命令。输入命令执行结构如下：

cost=10;number=12;ifnumber8

sums=number*0.95*cost;end,sums

4.M文件的一般结构

(A)编写一个画任意半径彩色线型的圆.(B)完整函数文件的基本结构.(C)函数文件各基本组成部分的作用.

如下是M函数，只需要调用即可,执行结果如下：

functionsa=functiondemo(r,s)

%CIRCLE

ifnargin2

error(输入变量太多。);

end

ifnargin==1

s=b;

end;clf;

t=0:pi/100:2*pi;x=r*exp(i*t);

ifnargin==2

plot(x,s);

else

sa=pi*r*r;

fill(real(x),imag(x),s)

end

axis(square)

比如输入“functiondemo(4,red)”即可得到右图。

三、思考题

1．编写函数文件，分别用 for 和 while 循环语句编写程序，求出

K??63

i?0

2i?1?22?23?24?25?26? ?263之和，并比较算法的运算时间。

答：（输入一个n值，就可以得到2的0次方一直加的2的n-1次方的累加值，只需要调用函数输入n的值为63即可完成上题）M函数如下：

for循环：

functions=input(n)tic

a=0;

fori=1:n+1x(i)=i-1;a=a+2^x(i);

endt=toca

while循环：

functioninputs(n)tic

b=0;i=0;

whilei=nb=b+2^i;i=i+1;

endt=tocb

得出结果：

?1,x?1

?

用MATLAB语言实现下面分段函数：y?f(x)??x,|x|?1

???1,x??1

?

x由键盘输入。

答：同上面做法一样，编辑M函数实现，然后输入调用函数即可）

M函数如下：

functionfenduan(n)ifn1

y=1;

elseifn=-1n=1y=n;

elseifn-1y=-1;

end;y

用for-end循环求出从100-200之间的所有素数。

答：编辑M函数实现这一命令，如下：

functionsushu

m=1;

forn=100:200

ifrem(n,2)~=0rem(n,3)~=0rem(n,5)~=0rem(n,7)~=0rem(n,11)~=0rem(n,13)~=0

x(m)=n;m=m+1;

endend

x

输出结果如下：

用while-end循环求1-100的所有数值和。

答：while-end循环M函数如下：：（输入不同的n值就可以得到1~n的累加和。）

functionsb(n)p=0;

i=0;

whilei=np=p+i;i=i+1;

end,p

四、遇到的问题

第一个问题是出现在运算时间的计算上，对于etime函数的使用，计算出来的时间总是0，经过反复的尝试和与同学的讨论，终于找到了好的途径，将“t1=clock”、“t1=clock”命令放在正确的地方才可以执行出正确的结果。

在第二个思考题中的命令“-1=n=1”无法正常执行，这是一个小失误，只需要用到“”罗技操作符，改变命令为“n=-1n=1”就可以正常执行了。

五、体会

这次实习的目的就是熟悉和掌握MATLAB的m文件程序和程序设计。在这里我们要先打开MATLAB，然后建立一个空的M文件程序，在空的M文件中编辑程序就可以了。

在这次实习中我们运用编辑了许多函数，比如“i