1.符號運算

1.1創(chuàng)建符號對象

• sym()
• syms()
  其中syms()可以創(chuàng)建多個對象

a=1/4;
b='1/4';
c=sym(1/4);
d=sym('1/4');
classa=class(a)
classb=class(b)
classc=class(c)
classd=class(d)

classa =

double


classb =

char


classc =

sym


classd =

sym

1.2 例子

1. 計算積分

syms A t tao w
 yf=int(A*exp(-i*w*t),t,-tao/2,tao/2)
>>
(2*A*sin((tao*w)/2))/w

2. 計算矩陣行列式值，特征值，特征根

• det() ：計算行列式
• inv() : 計算特征值
• eig() : 計算特征根

 syms a11 a12 a21 a22
 A=[a11,a12;a21,a22]
 DA=det(A)
 IA=inv(A)
 EA=eig(A)
>>
DA =
 
a11*a22 - a12*a21
 
 
IA =
 
[  a22/(a11*a22 - a12*a21), -a12/(a11*a22 - a12*a21)]
[ -a21/(a11*a22 - a12*a21),  a11/(a11*a22 - a12*a21)]
 
 
EA =
 
 a11/2 + a22/2 - (a11^2 - 2*a11*a22 + a22^2 + 4*a12*a21)^(1/2)/2
 a11/2 + a22/2 + (a11^2 - 2*a11*a22 + a22^2 + 4*a12*a21)^(1/2)/2
 

3. 化簡與展開

• simpilify()

syms fai1 fai2
y=simplify(sin(fai1)*cos(fai2)-cos(fai1)*sin(fai2))

>>
y =
 
sin(fai1 - fai2)

y=sym('(x+1)^3')
 y=expand(y)

>>
y =
 
x^3 + 3*x^2 + 3*x + 1

4. 因式分解

• factor()

syms a x
f1=x^4-5*x^3+5*x^2+5*x-6;factor(f1)
f2=x^2-a^2;factor(f2)
factor(1025)

>>
ans =
 
[ x - 1, x - 2, x - 3, x + 1]
 
 
ans =
 
[ -1, a - x, a + x]
 

ans =

     5     5    41

2. 分支結(jié)構(gòu)

2.1if分支

clear
a=100;
b=20; 
if a<b 
   fprintf ('b>a')          % 在Word中輸入’b>a’單引號不可用，要在Editor中輸入
else 
   fprintf ('a>b')          % 在Word中輸入’b>a’單引號不可用，要在Editor中輸入
end

2.2 switch分支

 clear
n=input('輸入n=');         % 輸入n值
switch mod(n,5)            % mod是求余函數(shù)，余數(shù)為0，得0，余數(shù)不為0，得1
case 0
   fprintf ('%d是5的倍數(shù)',n) 
otherwise
   fprintf('%d不是5的倍數(shù)',n)
end

3.循環(huán)語句

3.1 while循環(huán)

 clear
x=0;                  % 初始化變量x
sum=0;                % 初始化sum變量
   while x<101         % 當x<101執(zhí)行循環(huán)體語句
     sum=sum+x;       % 進行累加
     x=x+2; 
   end                %  while結(jié)構(gòu)的終點
sum                  % 顯示sum

>>
2550

3.2 for循環(huán)

語法規(guī)則如下：

for ii = 初值：增值：終值
語句一
語句二
end

clear
sum=0;                         % 設(shè)置初值(必須要有)
for ii=1:100;                      %  for循環(huán),增量為1
    sum=sum+ii;
end
sum
%end                           % 程序結(jié)束

還有其他的continue，break語句與c語言類似，而input()的用法與python類似。
此外c語言中的”！“在matlab中用”~“

4.函數(shù)

在matlab中對于自帶的函數(shù)可以直接調(diào)用，若是自寫函數(shù)，如冒泡排序算法
新建腳本后，保存為BubbleSort.m，注意文件名要與函數(shù)名一致

function s=BubbleSort(x) 
% 冒泡排序,x為待排序數(shù)組 
n=length(x); 
for i=1:n-1             % 最多做n-1趟排序 
    flag=0;                 % flag為交換標志,本趟排序開始前,交換標志應(yīng)為假 
    for j=1:n-i             % 每次從前向后掃描,j從1到 n-i 
        if x(j)>x(j+1)      % 如果前項大于后項則進行交換 
            t=x(j+1);
            x(j+1)=x(j); 
            x(j)=t; 
            flag=1;         % 當發(fā)生了交換,將交換標志置為真 
        end 
    end 

在調(diào)用的時候，只用新建腳本（與函數(shù)文件在同一路徑下），即可直接調(diào)用

x = rand(1,10)
y = BubbleSort(x)
>>
 1 至 9 列

    0.0292    0.0855    0.2373    0.2625    0.4588    0.4886    0.5785    0.7303    0.8010

  10 列

    0.9289