第11章MATLAB71符號計算

1、技術凝聚實力技術凝聚實力專業(yè)創(chuàng)新出版專業(yè)創(chuàng)新出版第11章 符號計算l與一般的數值計算不同，符號計算是對字符串符號進行分析和運算，為了便于理解，讀者可以將符號計算看作“由計算機實現的數學公式推導”。進行符號計算時，matlab負責將計算請求提交給其內置的maple組件并返回maple的計算結果。matlab的符號計算歷經多次的改進和完善，其功能已經非常強大，尤其是在大規(guī)模的簡單公式推導、邏輯推導等應用中有重要應用。本章介紹的內容如下。l符合對象的創(chuàng)建 符號函數l符號對象的復合和化簡 符號對象的運算l符號矩陣l符號微積分l符號方程求解技術凝聚實力技術凝聚實力專業(yè)創(chuàng)新出版專業(yè)創(chuàng)新出版11.1 符號類

2、型與數值類型l符號對象（symbolic objet）是matlab符號工具箱中定義的一種新的數據類型。與一般的簡單數值類型（如logical、int、float、double等）不同，符號對象是一種保存了符號字符串表示的復雜數據結構。本節(jié)主要討論符號類型與普通數值類型的差別。技術凝聚實力技術凝聚實力專業(yè)創(chuàng)新出版專業(yè)創(chuàng)新出版11.2 符號對象及其創(chuàng)建lmatlab符號工具箱定義的數據類型符號對象（symbolic object），是一種存儲了符號字符串表示的復雜數據結構。matlab利用符號對象表示符號變量、符號表達式、符號函數、符號矩陣，本節(jié)介紹符號變量、符號表達式的創(chuàng)建，至于符號函數和符號

3、矩陣，將分別在11.3、11.5節(jié)進行詳細介紹。技術凝聚實力技術凝聚實力專業(yè)創(chuàng)新出版專業(yè)創(chuàng)新出版11.2.1 創(chuàng)建符號變量lsym函數可用于創(chuàng)建符號常量和變量，其調用格式為：lsym_variable = sym(sym_expression)l該指令創(chuàng)建符號變量（或常量）sym_variable，其中sym_expression為符號變量的字符串表示。如果sym_expression為常數（例如pi、1、2.01等），則函數返回符號常量；如果sym_expression是符合命名規(guī)則的字符串（即該字符串只能包含數字、英文字母和下畫線，且不能以數字開頭），則函數返回符號變量。技術凝聚實力技術凝

4、聚實力專業(yè)創(chuàng)新出版專業(yè)創(chuàng)新出版11.2.2 創(chuàng)建符號表達式l符號表達式是符號變量或常量的組合，也可以將符號變量和符號常量認為是符號表達式。符號表達式的創(chuàng)建，依然使用函數sym，其調用格式為：lsym_variable = sym(sym_expression)l這里的輸入字符串參數sym_expression與11.2.1中有所不同，它是由多個子字符串通過數值運算符組合而成的，每個子字符串或者是符合命名規(guī)則的字符串或者是數值常量。也可以事先利用sym或syms創(chuàng)建符號變量和常量，然后通過數值運算符、支持符號計算的函數等將這些符號變量或常量組合起來構成符號表達式。技術凝聚實力技術凝聚實力專業(yè)創(chuàng)新

5、出版專業(yè)創(chuàng)新出版11.2.3 符號對象與數值的轉換lmatlab符號計算工具箱提供如下三類運算操作：lnumberic，即matlab浮點運算；lrational，即maple符號運算；lvpa，即maple精度可變運算。l浮點運算是三者中運算速度最快、內存使用最少的一種，但其結果的精度也較低，format可以控制輸出結果的位數，但其內部表示仍然是8字節(jié)浮點表示，也就是說format不能改變結果的精度。符號運算的結果是精確、無近似的，但運算速度慢、消耗內存最多。vpa是浮點運算與符號運算的折中，用戶可以通過函數digits控制計算結果的精度，默認情況下digits = 32。技術凝聚實力技術凝

6、聚實力專業(yè)創(chuàng)新出版專業(yè)創(chuàng)新出版11.3 符號函數l顯而易見，符號函數與通常的數值函數有很大的不同。本節(jié)從matlab中符號函數的創(chuàng)建開始，分別介紹了符號函數的復合、反函數等內容，然后給出了一些在理論和實際中具有重要應用的符號函數。技術凝聚實力技術凝聚實力專業(yè)創(chuàng)新出版專業(yè)創(chuàng)新出版11.3.1 創(chuàng)建符號函數l可以通過下面兩種方法創(chuàng)建符號函數：l符號表達式；lm文件。l利用符號表達式創(chuàng)建符號函數的方法在例11.4中已經介紹，這種方法適合于創(chuàng)建一些比較簡單的函數，對于比較復雜的符號函數的創(chuàng)建，應該利用m文件。技術凝聚實力技術凝聚實力專業(yè)創(chuàng)新出版專業(yè)創(chuàng)新出版11.3.2 復合函數l已知函數、，復合函數定

7、義為。matlab符號計算工具箱提供函數compose用于符號函數復合，其調用格式為：lcompose(f,g)lcompose(f,g,z)lcompose(f,g,x,z)lcompose(f,g,x,y,z)技術凝聚實力技術凝聚實力專業(yè)創(chuàng)新出版專業(yè)創(chuàng)新出版11.3.3 反函數lmatlab符號工具箱的函數finverse用于計算函數的反函數，其調用格式為：lg = finverse(f)lg = finverse(f,v)l其中f為符號表達式，該函數求f關于v的反函數，如果未指定v，那么v = findsym(f)。技術凝聚實力技術凝聚實力專業(yè)創(chuàng)新出版專業(yè)創(chuàng)新出版11.3.4 特殊符號函

8、數lmatlab符號工具箱對一些重要的數學函數提供了支持，技術凝聚實力技術凝聚實力專業(yè)創(chuàng)新出版專業(yè)創(chuàng)新出版11.4 符號表達式的復合和化簡l符號表達式的顯示、展開、分解、化簡等是符號運算的重要內容，matlab符號工具箱對這些內容提供了強大的支持。本節(jié)將詳細介紹matlab符號工具相中符號表達式的顯示、展開、分解、替換、化簡等操作。技術凝聚實力技術凝聚實力專業(yè)創(chuàng)新出版專業(yè)創(chuàng)新出版11.4.1 顯示符號表達式lmatlab顯示符號表達式時，默認的方式是matlab式的顯示，例如在matlab命令窗口輸入以下代碼。l f = a*x2 + b*x + cl運行結果如下。lf =la*x2+b*x+

9、c技術凝聚實力技術凝聚實力專業(yè)創(chuàng)新出版專業(yè)創(chuàng)新出版11.4.2 展開、分解及化簡符號表達式l對數學表達式的展開、分解、化簡是學習數學和應用數學中經常要做的工作，matlab符號工具箱這些工作提供了強大的支持，下面列出了一些相關的函數及其使用方法。l函數collect用于同類項合并，其調用格式為：lcollect(f, s)lcollect(f)lcollect(f, s)將符號表達式f看作符號變量s的多項式（如果未指定s，則默認s為findsym(f, 1)），對具有相同階數的項進行合并。技術凝聚實力技術凝聚實力專業(yè)創(chuàng)新出版專業(yè)創(chuàng)新出版11.4.3 符號表達式替換l符號表達式的替換有兩種，一是

10、利用符號變量或符號表達式替換目標符號表達式中的符號變量（包括符號常量），利用函數subs實現；二是用符號變量替代目標符號表達式中的某個子符號表達式，利用函數subexpr實現。函數subs的一般調用格式為：lr = subs(s) lr = subs(s, new)lr = subs(s,old,new)技術凝聚實力技術凝聚實力專業(yè)創(chuàng)新出版專業(yè)創(chuàng)新出版11.5 符號矩陣l符號矩陣是數值矩陣的擴展，符號矩陣的每個元素為符號表達式。與數值矩陣一樣，符號矩陣也有代數運算、邏輯運算等；還可以求符號矩陣的奇異值，對符號進行奇異值分解。對于符號方陣，可以求行列式、特征值；符號方陣可逆時，還可以求符號方陣的

11、逆。本節(jié)涉及的許多符號矩陣運算與數值矩陣相同，讀者在學習本節(jié)內容時可以參考第7章。技術凝聚實力技術凝聚實力專業(yè)創(chuàng)新出版專業(yè)創(chuàng)新出版11.5.1 創(chuàng)建符號矩陣lmatlab符號工具箱中沒有提供專門用于創(chuàng)建符號矩陣的函數，用戶可以通過以下三種方式創(chuàng)建符號矩陣：l直接輸入符號矩陣的各個元素符號表達式；l利用sym函數將數值矩陣轉化為符號矩陣，這種方法只能創(chuàng)建符號常量矩陣；l簡單符號矩陣的組合創(chuàng)建復雜符號矩陣。技術凝聚實力技術凝聚實力專業(yè)創(chuàng)新出版專業(yè)創(chuàng)新出版11.5.2 符號矩陣的代數運算l同數值矩陣一樣，符號矩陣也能進行加、減、乘、除四則運算，并且具有與數值矩陣幾乎相同的運算規(guī)則，只是矩陣元素由數值

12、變成了符號表達式。技術凝聚實力技術凝聚實力專業(yè)創(chuàng)新出版專業(yè)創(chuàng)新出版11.5.3 符號矩陣的邏輯運算lmatlab符號工具箱對符號的邏輯運算支持很有限，僅支持=和=。l符號的邏輯運算。 l運行結果如下。l? function gt is not defined for values of class sym.lerror in = gt at 16lbuiltin(gt, varargin:);技術凝聚實力技術凝聚實力專業(yè)創(chuàng)新出版專業(yè)創(chuàng)新出版11.5.4 符號矩陣的行列式l符號矩陣的行列式與數值矩陣的行列式具有相同的定義，只是符號矩陣的行列式為符號函數。同樣，行列式僅對符號方陣才有意義。lmat

13、lab利用同名函數det計算符號矩陣的行列式，該函數在符號工具箱中被重載。l符號矩陣的行列式。l矩陣 被稱為單位旋轉矩陣，試求該矩陣的行列式。 sincoscossinxxxx技術凝聚實力技術凝聚實力專業(yè)創(chuàng)新出版專業(yè)創(chuàng)新出版11.5.5 符號矩陣的逆l對階符號方陣，滿足的符號方陣稱為的逆，記為。其中為單位符號方陣，對角線上元素為符號常量，其余為0。lmatlab符號工具箱提供同名函數inv用于計算符號方陣的逆，這里的inv在符號工具箱中被重載。下面是符號矩陣求逆應用的例子。技術凝聚實力技術凝聚實力專業(yè)創(chuàng)新出版專業(yè)創(chuàng)新出版11.5.6 符號矩陣的秩l符號矩陣秩的含義與數值矩陣相似，對符號矩陣，取

14、其任意行列得到的階子式， 的秩定義為使“至少有一個階子式的行列式不為0”的最大，即所有的階子式行列式為0，且至少有一個階子式行列式不為0。lmatlab提供同名函數rank計算符號矩陣的秩，這里的rank在符號工具箱中被重載。技術凝聚實力技術凝聚實力專業(yè)創(chuàng)新出版專業(yè)創(chuàng)新出版11.5.7 符號矩陣的特征分解l特征值反映了符號矩陣的一些重要的性質，對符號矩陣特征值的分析在信號處理、自動控制中有重要的應用。符號方陣的特征值是滿足的，其中為單位符號矩陣。由于即為a的特征多項式，所以的特征值即為特征多項式的根?；貞浀?0章中利用函數poly求數值矩陣的特征多項式，事實上，在符號工具箱中該函數被重載，利用

15、poly(a)可以得到符號方陣的特征多項式，對特征多項式求根：lsolve(poly(a)技術凝聚實力技術凝聚實力專業(yè)創(chuàng)新出版專業(yè)創(chuàng)新出版11.5.8 符號矩陣的svd分解l關于符號矩陣svd分解的相關概念，可以參考第7章數值矩陣svd分解，這里不再贅述。需要強調的是，matlab符號工具箱僅支持常量符號矩陣（元素為符號常量的符號矩陣）的svd分解，這其中的原因可能是對一般的符號矩陣來說，svd分解過于復雜而變得不可用。lmatlab中利用同名函數svd對符號矩陣進行svd分解，其一般的調用格式為：lu，s，v = svd(a)技術凝聚實力技術凝聚實力專業(yè)創(chuàng)新出版專業(yè)創(chuàng)新出版11.6 符號微積

16、分l第7章已經介紹了數值微積分，數值微積分的適用范圍廣泛，但對問題的描述一般是近似的，而符號微積分可能得到問題的精確描述，因此在條件允許的情況下利用符號微積分解決理論和實際問題是一個很好的選擇。本節(jié)內容包括符號的極限、微分、積分、taylor級數以及級數求和等內容。技術凝聚實力技術凝聚實力專業(yè)創(chuàng)新出版專業(yè)創(chuàng)新出版11.6.1 極限l極限是當變量無限接近特定值時函數的值，例如，一元函數的導數的定義為下面的極限：。matlab符號工具箱利用函數limit計算符號極限，技術凝聚實力技術凝聚實力專業(yè)創(chuàng)新出版專業(yè)創(chuàng)新出版11.6.2 微分技術凝聚實力技術凝聚實力專業(yè)創(chuàng)新出版專業(yè)創(chuàng)新出版11.6.3 積分l函數積分是微分的反運算，對符號函數，的積分是使的，用 表示。同diff相似，如果為多變量符號表達式，用戶應該指定對某個自變量求積分，即f = int(f, v)；如果未指定，則v設定為由函數fingsym(f, 1)決定的默認符號變量。技術凝聚實力技術凝聚實力專業(yè)創(chuàng)新出版專業(yè)創(chuàng)新出版11.6.4 taylor展開技術凝聚實力技術凝聚實力專業(yè)創(chuàng)新出版專業(yè)創(chuàng)新出版11