天天看点

初学Matlab编程

计算机程序就是计算机指令的集合,不同的编程语言指令与功能是不一样的.MATLAB语言是一种面向对象的高级语言,它具有编程效率高、易学易用的优点.

    MATLAB与其它大部分高级语言一样,有它自己的控制流语句.控制流极其重要,因为它使过去的计算影响将来的运算。MATLAB提供如下几种控制流结构:For循环,While循环,If-Else-End结构和switch-case-end结构。由于这些结构经常包含大量的MATLAB命令,故经常出现在M文件中.MATLAB支持的控制流语句和C语言支持的控制流语句在调用格式上非常相似.

1.For 循环

    For循环允许一条语句或一组语句被重复执行预先指定的次数。For循环的一般形式是:

    for  x = array

        语句

    end

在for和end语句之间的语句按数组中的每一列执行一次。在每一次迭代中,x被指定为数组的下一列,即在第n次循环中,x=array(:, n)。例如,

       for n=1:10

           x(n)=sin(n*pi/10);

       end

       x

       x =

    0.3090    0.5878    0.8090    0.9511    1.0000    0.9511    0.8090    0.5878    0.3090    0.0000

    换句话,第一语句是说:对n等于1到10,执行所有语句,直至下一个end语句。第一次通过For循环n=1,第二次,n=2,如此继续,直至n=10。在n=10以后,For循环结束,然后执行end语句后面的任何命令.注意,该循环结束后,n=10.

    For循环的其它重要方面是:

    (1)For循环不能用For循环内重新赋值循环变量n来终止。

       for n=1:10

           x(n)=sin(n*pi/10);

           n=9;

       end

       x

       x =

    0.3090    0.5878    0.8090    0.9511    1.0000    0.9511    0.8090    0.5878    0.3090    0.0000

        n

        n=

           9

执行过程是这样的:

n=1,

x(1)=sin(pi/10),

n=9,

n=2,

x(2)=sin(2*pi/10),

n=9,

n=3,

...,

n=10,

x(10)=sin(10*pi/10),

n=9.

循环结束后,n=9.

    (2)在For循环内接受任何有效的MATLAB数组。

    data=[3  9  45  6;  7  16  -1  5]

    data =

          3     9    45     6

          7    16    -1     5

    for n=data

        x=n(1)-n(2)

    end

       x =

          -4

       x =

          -7

       x =

          46

       x =

          1

    (3)For循环可按需要嵌套。

      for n=1:5

          for m=1:5

              A(n,m)=n^2+m^2;

          end

          disp(n)

      end

      1

      2

      3

      4

      5

      A

      A =

          2     5    10    17    26

          5     8    13    20    29

          10    13    18    25    34

          17    20    25    32    41

          26    29    34    41    50

    (4)当有一个等效的数组方法来解给定的问题时,应避免用For循环。例如,上面的第一个例子可被重写为

       n=1:10;

       x=sin(n*pi/10)

       x =

    0.3090    0.5878    0.8090    0.9511    1.0000    0.9511    0.8090    0.5878    0.3090    0.0000

两种方法得出同样的结果,而后者执行更快,更直观,要求较少的输入。

    (5)为了得到最大的速度,在For循环(While循环)被执行之前,应预先分配数组。例如,前面所考虑的第一种情况,在For循环内每执行一次命令,变量x的大小增加1。迫使MATLAB每通过一次循环要花费时间对x分配更多的内存。为了消去这个步骤,For循环的例子应重写为

    x=zeros(1,10); 

    for n=1:10

        x(n)=sin(n*pi/10);

    end

现在,只有x(n)的值需要改变。

    例1 相传古代印度国王要褒奖他的聪明能干的宰相达依尔(国际象棋发明者),问他要什么?达依尔回答:“陛下只要在国际象棋棋盘的第一个格子上放一粒麦子,第二个格子上放二粒麦子,以后每个格子的麦子数都按前一格的两倍计算。如果陛下按此法给我64格的麦子,就感激不尽,其他什么也不要了。”国王想:“这还不容易!”让人扛了一袋麦子,但很快用光了,再扛出一袋还不够,请你为国王算一下共要给达依尔多少小麦?(1小麦约颗)

    解 麦粒总数为

    程序如下:

a=1;

s=0

for i=1:64

    s=s+a;

    a=2*a;

end

s=s/1.4/10^8

运行后得:

s=

   1.3176e+011

    例2 公元前五世纪我国古代数学家张丘建在《算经》一书中提出了“百鸡问题”:鸡翁一值钱五,鸡母一值钱三,鸡雏三值钱一。百钱买百鸡,问鸡翁、母、雏各几何?

     解 设 x:鸡翁数,则x的范围:0~19

           y:鸡母数,则y的范围:0~33

           z:鸡雏数,则z的范围:0~100

则:

     x+y+z=100

     5x+3y+z/3=100

这是一个不定方程。

for x=0:19

    for y=0:33

        for z=0:100

            if (x+y+z==100)&(5*x+3*y+z/3==100)

                d=[x,y,z]

            end

        end

    end

end

运行后得结果:

d =

     0    25    75

d =

     4    18    78

d =

     8    11    81

d =

    12     4    84

2.While  循环

    与For循环以固定次数求一组命令的值相反,While 循环以不定的次数重复执行一组语句。While循环的一般形式是:

    while 表达式1

          语句1

    end

只要表达式1里的所有元素为真,就执行while和end之间的语句1,否则,就结束循环。通常,表达式的值给出一个标量值,但数组值也同样有效。在数组情况下,当数组的所有元素为真(值不等零)时,就执行语句1,数组中有一个元素为假(值为零),就结束循环。

   例3 按下面的公式计算:

使误差小于给定的.

   解 把作为误差,程序如下:

error=input('请输入误差:');

x=1;

y=0;

n=1;

while x>error

    y=y+1;

    x=x/n;

    n=n+1;

end

e=y

运行如下:

请输入误差:0.001

e =

    2.7181

3.IF-ELSE-END 结构

    很多情况下,命令的序列必须根据关系的检验有条件地执行。在编程语言里,这种逻辑由某种If-Else-End结构来提供。最简单的If-Else-End结构是:

   if  表达式1

       语句1

   end

    如果在表达式1中的所有元素为真(非零),那么就执行if和end语言之间的语句1。

    假如有两个选择,If-Else-End结构是:

    if  表达式1

        语句1

     else

        语句2

     end

在这里,如果表达式1为真,则执行语句1;如果表达式是假,则执行语句2。

    当有三个或更多的选择时,If-Else-End结构采用形式

     if 表达式1

        语句1

     elseif 表达式2

        语句2

     elseif 表达式3

        语句3

     elseif 表达式4

        语句4

     elseif ……

            .

            .

            .

     else

         语句

     end

    如果表达式1为真,则执行语句1,结束循环;如果表达式1为假,则检验表达式2,如果表达式2为真,则执行语句2,结束循环;如果表达式2为假,则检验表达式3,如此下去,如果所有表达式都为假时,则执行最后的语句。即只执行第一个真值表达式相关的语句;接下来的表达式不检验,跳过其余的If-Else-End结构。而且,最后的else命令可有可无。

4.switch-case-end结构

    如果在一个程序中,必须针对某个变量不同取值情况进行相应操作,switch语句比if else语句更方便。switch语句的一般形式为:

    switch 分支条件(数值或字符串)

    case  数值(或字符串)条件1

          语句1

    case  数值(或字符串)条件2

          语句2

    case  数值(或字符串)条件3

          语句3

    case  ...

          ...

    otherwise

          语句

    end

其中分支条件可以是一个函数、变量或者表达式.如果条件1与分支条件匹配就执行语句1,退出循环;否则,检验条件2,如果条件2与分支条件匹配执行语句2,退出循环;否则,检验条件3,...,当所有条件都不与分支条件匹配时就执行最后的语句。注意otherwise是可以省略的。

    例4 在图形界面上放置一个弹出式菜,点击弹出式菜单可以设置曲线颜色。

    解:程序如下:

x=linspace(-6,6,50);

y=sin(x);

h=plot(x,y,'linewidth',4);

h0=uicontrol('position',[600,500,100,30],...

    'style','popupmenu',...

    'string','红色|蓝色|黄色|黑色|青色',...

    'callback',...

    ['v=get(h0,''value'');',...

        'switch v;',...

        'case 1;',...

        'set(h,''color'',''r'');',...

        'case 2;',...

        'set(h,''color'',''b'');',...

        'case 3;',...

        'set(h,''color'',''y'');',...

        'case 4;',...

        'set(h,''color'',''k'');',...

        'case 5;',...

        'set(h,''color'',''c'');',...

        'end;']);

继续阅读