一.順序結構稱程式

1.edit命令建立指令碼

edit filename;

2.順序結構

資料輸入：A = input('提示資訊');
資料輸出：disp(A);
程式暫停：pause(delay time);  注：時間為0，程式停止
         ctrl+c 強行終止程式

例1  有一線段AB，A的座標為(1，1)，B的座標為(4.5，4.5)，求AB的長度，
以及黃金分割點C的座標。
a=input('a=');
b=input('b=');
c=a+0.618*(b-a);
s=abs(a-b);
disp(s)
disp(c)

二.if語句

1.當條件語句為矩陣，矩陣非空且不包含非零元素時為真。

例1  輸入一個整數，若為奇數則輸出其平方根，否則輸出其立方根。（雙分支）  
x=input('請輸入x的值:');
if rem(x,2)==1    %等價if rem(x,2)
    y=sqrt(x);
else
    y=x^(1/3);    %等價y=power(x,1/3);
end
y

例2  輸入一個字元，若為大寫字母，則輸出其對應的小寫字母；若為小寫字母，則輸出
其對應的大寫字母；若為數字字元則輸出其對應數的平方，若為其他字元則原樣輸出。
c=input('請輸入一個字元：','s');
if c>='A' && c<='Z'
 
         %多分支
    disp(lower(c))
elseif c>='a' && c<='z'
    disp(upper(c))
elseif c>='0' && c<='9'
    disp(str2double(c)^2)
else
    disp(c)
end

三.switch語句

例1  輸入一個英文單詞，判斷它是否以母音字母開頭。
%%使用switch語句
c=input('請輸入一個單詞：','s');
switch c(1) 
    case {'A','E','I','O','U'
 
,'a','e','i','o','u'} 
        disp([c,'以母音字母開頭']);
    otherwise 
        disp([c,'以子音字母開頭']);
end

%%使用if語句findstr函式 
c=input('請輸入一個單詞：','s');
if findstr(c(1),'AEIOUaeiou')>0
    disp([c,'以母音字母開頭']);
else
    disp([c,'以子音字母開頭']);
end

例2  PM2.5是指大氣中直徑小於或等於2.5微米的可入肺顆粒物，是衡量空氣質量的重要標。
假定空氣質量等級以PM2.5數值劃分為6級。 PM2.5數值在[0，35）空氣質量為優，[35，75）
為良，[75，115）為輕度汙染，[115，150）為中度汙染，[150，250）為重度汙染，大於等
於250為嚴重汙染。編寫程式，輸入PM2.5數值，輸出空氣質量等級。
g=input('請輸入PM2.5值：');
switch fix(g) 
    case num2cell(0:34)   % num2cell：將數值矩陣轉換為單元矩陣
        disp('空氣質量優');
    case num2cell(35:74) 
        disp('空氣質量良好');
    case num2cell(75:114) 
        disp('空氣質量輕度汙染');
    case num2cell(115:149) 
        disp('空氣質量中度汙染'); 
    case num2cell(150:249) 
        disp('空氣質量重度汙染');  
    otherwise 
        disp('空氣質量嚴重汙染');
end

四.for語句

1.說明

    1)for語句針對向量的每一個元素執行一次迴圈體
    2)退出迴圈後，迴圈變數的值就是向量中最後元素的值
    3)向量為空迴圈體不執行
    4)對於矩陣，迴圈變數=矩陣的每一列，此時迴圈變數為向量
     （例：5行4列矩陣執行4次）

例：計算圓周率π。
（1）利用無窮級數展開式求π的近似值。
%%用if迴圈語句%
y=0;
g=-1;
n=input('n=?'); 
for i=1:n
    g=-g;
    y=y+g/(2*i-1); 
end
pai=4*y

%%用向量求和的方法實現程式%
n=input('n=?');
x=1:2:(2*n-1);
y=(-1).^(2:n+1)./x;
pai=sum(y)*4


（2）利用定積分的近似值求π的近似值。
a=0;
b=1;
n=input('n=?');
h=(b-a)/n;
x=a:h:b;
f=sqrt(1-x.*x);
s=[];
for k=1:n
    s1=(f(k)+f(k+1))*h/2;
    s=[s,s1];
end
pai=4*sum(s)


（3）利用蒙特卡洛法求π的近似值。
s=0;
n=input('n=?');
for i=1:n
    x=rand(1);
    y=rand(1);
    if x*x+y*y<=1
        s=s+1;
    end
end
pai=s/n*4

五.while語句

例1  從鍵盤輸入若干個數，當輸入0時結束輸入，求這些數的平均值和它們之和。
msum=0;
n=0;
x=input('Enter a number (end in 0):');
while x~=0
    msum=msum+x;
    n=n+1;
    x=input('Enter a number (end in 0):');
end
if n>0
    msum
    mean=msum/n
end

例2  求[100，200]之間第一個能被21整除的整數。
for n=100:200
    if rem(n,21)~=0
        continue      %continue跳出本次迴圈，繼續執行下一次
    end
    n
    break          %break跳出迴圈
end

例3  用篩選法求某自然數範圍內的全部素數。
m=input('m='); 
p=1:m;
p(1)=0;
for i=2:sqrt(m)
    for j=2*i:i:m
        p(j)=0;
    end
end
n=find(p~=0);
p(n)

六.函式檔案定義與呼叫

1.定義、呼叫

    return：表示結束函式的執行，可不使用，程式會自動返回

例1  編寫函式檔案，求半徑為r的圓的面積和周長。
%%functoin 輸出形參表=函式名(輸入形參表)%
function [s,p]=fcircle(r)
s=pi*r*r;
p=2*pi*r;

2.匿名函式

    基本格式：
    函式控制代碼變數[email protected](匿名函式輸入引數)匿名函式表示式   @：函式控制代碼運算名
    函式控制代碼變數[email protected]函式名

例2  已知=(40)/((30)+(20))
①當()=+10 ln⁡(^2+5)時，的值是多少。
②當()=1×2+2×3+3×4+⋯+×(+1)時，的值是多少。
分別用匿名函式和函式檔案定義函式。
%%第二問的函式檔案f2.m。%
function f=f2(n)
f=0;
for k=1:n
    f=f+k*(k+1);
end
%%第一問指令碼檔案mf.m。%
[email protected](n) n+10*log(n*n+5);
y1=f1(40)/(f1(30)+f1(20))
y2=f2(40)/(f2(30)+f2(20))

七.函式的遞迴呼叫

例1  利用函式的遞迴呼叫，求n！。
函式檔案fact.m如下：
function f=fact (n)
if n<=1
    f=1;
else
    f=fact (n-1)*n;    %遞迴呼叫求(n-1)!
end
在指令碼檔案a.m中呼叫函式檔案fact.m，求n!。
n=input('Please input n=');
s=fact (n);
disp(s)

例2  Fibonacci數列定義如下：
                f1=1
                f2=1
                fn=fn-1+fn-2   （n>2） 
編寫遞迴呼叫函式求Fibonacci數列的第n項，然後呼叫該函式驗證Fibonacci數列
的如下性質：
                  f1^2+f2^2+f3^2+……+fn^2=fn×fn+1
首先建立函式檔案ffib.m。
function f=ffib(n)
if n>2
    f=ffib(n-1)+ffib(n-2);
else
    f=1;
end
建立程式檔案test.m。
F=[];
for k=1:20
    F=[F,ffib(k)*ffib(k)];
end
sum(F)
ffib(20)*ffib(21)

八.函式引數與變數的作用域

1.函式引數的可調性

    nargin： 記錄輸入實參的個數
    nargout：記錄輸出實參的個數

例1  利用nargin和nargout建立函式檔案test.m。
function fout=test(a,b,c)
if nargin==1       %輸入引數為一個時
    fout=a;
elseif nargin==2   %輸入引數為兩個時
    fout=a+b;
elseif nargin==3
    fout=(a*b*c)/2;
end

2.全域性變數與區域性變數

    全域性變數定義格式：global 變數名

例2  利用全域性變數建立函式檔案wad.m。
function f=wad(x,y)
global ALPHA BETA
f=ALPHA*x+BETA*y;

%%在工作區使用全域性變數
>>global ALPHA BETA
>>ALPHA=1
>>BETA=2
>>s=wad(2)
s = 
     5