第5講matlab數(shù)據(jù)分析與多項式計算

1、1第第5講講 MATLAB數(shù)據(jù)分析與多項式計算數(shù)據(jù)分析與多項式計算5.1 數(shù)據(jù)統(tǒng)計處理數(shù)據(jù)統(tǒng)計處理5.2 數(shù)據(jù)插值數(shù)據(jù)插值5.3 曲線擬合曲線擬合5.4 離散傅立葉變換離散傅立葉變換5.5 多項式計算多項式計算25.1 數(shù)據(jù)統(tǒng)計處理數(shù)據(jù)統(tǒng)計處理5.1.1 最大值和最小值最大值和最小值MATLAB提供的求數(shù)據(jù)序列的最大值和最小值的函提供的求數(shù)據(jù)序列的最大值和最小值的函數(shù)分別為數(shù)分別為max和和min，兩個函數(shù)的調(diào)用格式和操，兩個函數(shù)的調(diào)用格式和操作過程類似。作過程類似。1求向量的最大值和最小值求向量的最大值和最小值求一個向量求一個向量X的最大值的函數(shù)有兩種調(diào)用格式，分的最大值的函數(shù)有兩種調(diào)用格式

2、，分別是：別是：(1) y=max(X)：返回向量：返回向量X的最大值存入的最大值存入y，如果，如果X中包含復數(shù)元素，則按模取最大值。中包含復數(shù)元素，則按模取最大值。3(2) y,I=max(X)：返回向量：返回向量X的最大值存入的最大值存入y，最大，最大值的序號存入值的序號存入I，如果，如果X中包含復數(shù)元素，則按模中包含復數(shù)元素，則按模取最大值。取最大值。求向量求向量X的最小值的函數(shù)是的最小值的函數(shù)是min(X)，用法和，用法和max(X)完全相同。完全相同。例例5-1 求向量求向量x的最大值。的最大值。命令如下：命令如下：x=-43,72,9,16,23,47;y=max(x) %求向量求

3、向量x中的最大值中的最大值y,l=max(x) %求向量求向量x中的最大值及其該元素中的最大值及其該元素的位置的位置42求矩陣的最大值和最小值求矩陣的最大值和最小值求矩陣求矩陣A的最大值的函數(shù)有的最大值的函數(shù)有3種調(diào)用格式，分種調(diào)用格式，分別是：別是：(1) max(A)：返回一個行向量，向量的第：返回一個行向量，向量的第i個個元素是矩陣元素是矩陣A的第的第i列上的最大值。列上的最大值。(2) Y,U=max(A)：返回行向量：返回行向量Y和和U，Y向向量記錄量記錄A的每列的最大值，的每列的最大值，U向量記錄每列向量記錄每列最大值的行號。最大值的行號。5(3) max(A,dim)：dim取取

4、1或或2。dim取取1時，時，該函數(shù)和該函數(shù)和max(A)完全相同；完全相同；dim取取2時，該時，該函數(shù)返回一個列向量，其第函數(shù)返回一個列向量，其第i個元素是個元素是A矩矩陣的第陣的第i行上的最大值。行上的最大值。求最小值的函數(shù)是求最小值的函數(shù)是min，其用法和，其用法和max完全相完全相同。同。例例5-2 分別求分別求34矩陣矩陣x中各列和各行元素中中各列和各行元素中的最大值，并求整個矩陣的最大值和最小的最大值，并求整個矩陣的最大值和最小值。值。63兩個向量或矩陣對應(yīng)元素的比較兩個向量或矩陣對應(yīng)元素的比較函數(shù)函數(shù)max和和min還能對兩個同型的向量或矩陣進行還能對兩個同型的向量或矩陣進行比

5、較，調(diào)用格式為：比較，調(diào)用格式為：(1) U=max(A,B)：A,B是兩個同型的向量或矩陣，是兩個同型的向量或矩陣，結(jié)果結(jié)果U是與是與A,B同型的向量或矩陣，同型的向量或矩陣，U的每個元素的每個元素等于等于A,B對應(yīng)元素的較大者。對應(yīng)元素的較大者。(2) U=max(A,n)：n是一個標量，結(jié)果是一個標量，結(jié)果U是與是與A同型同型的向量或矩陣，的向量或矩陣，U的每個元素等于的每個元素等于A對應(yīng)元素和對應(yīng)元素和n中的較大者。中的較大者。min函數(shù)的用法和函數(shù)的用法和max完全相同。完全相同。例例5-3 求兩個求兩個23矩陣矩陣x, y所有同一位置上的較大所有同一位置上的較大元素構(gòu)成的新矩陣元素

6、構(gòu)成的新矩陣p。75.1.2 求和與求積求和與求積數(shù)據(jù)序列求和與求積的函數(shù)是數(shù)據(jù)序列求和與求積的函數(shù)是sum和和prod，其使用方法類似。設(shè)其使用方法類似。設(shè)X是一個向量，是一個向量，A是一是一個矩陣，函數(shù)的調(diào)用格式為：個矩陣，函數(shù)的調(diào)用格式為：sum(X)：返回向量：返回向量X各元素的和。各元素的和。prod(X)：返回向量：返回向量X各元素的乘積。各元素的乘積。sum(A)：返回一個行向量，其第：返回一個行向量，其第i個元素是個元素是A的第的第i列的元素和。列的元素和。8prod(A)：返回一個行向量，其第：返回一個行向量，其第i個元素是個元素是A的第的第i列的元素乘積。列的元素乘積。su

7、m(A,dim)：當：當dim為為1時，該函數(shù)等同于時，該函數(shù)等同于sum(A)；當；當dim為為2時，返回一個列向量，時，返回一個列向量，其第其第i個元素是個元素是A的第的第i行的各元素之和。行的各元素之和。prod(A,dim)：當：當dim為為1時，該函數(shù)等同于時，該函數(shù)等同于prod(A)；當；當dim為為2時，返回一個列向量，時，返回一個列向量，其第其第i個元素是個元素是A的第的第i行的各元素乘積。行的各元素乘積。例例5-4 求矩陣求矩陣A的每行元素的乘積和全部元的每行元素的乘積和全部元素的乘積。素的乘積。95.1.3 平均值和中值平均值和中值求數(shù)據(jù)序列平均值的函數(shù)是求數(shù)據(jù)序列平均值

8、的函數(shù)是mean，求數(shù)據(jù)序列中值的函數(shù)，求數(shù)據(jù)序列中值的函數(shù)是是median。兩個函數(shù)的調(diào)用格式為：。兩個函數(shù)的調(diào)用格式為：mean(X)：返回向量：返回向量X的算術(shù)平均值。的算術(shù)平均值。median(X)：返回向量：返回向量X的中值。的中值。mean(A)：返回一個行向量，其第：返回一個行向量，其第i個元素是個元素是A的第的第i列的算術(shù)列的算術(shù)平均值。平均值。median(A)：返回一個行向量，其第：返回一個行向量，其第i個元素是個元素是A的第的第i列的中列的中值。值。mean(A,dim)：當：當dim為為1時，該函數(shù)等同于時，該函數(shù)等同于mean(A)；當；當dim為為2時，返回一個列向

9、量，其第時，返回一個列向量，其第i個元素是個元素是A的第的第i行的算術(shù)行的算術(shù)平均值。平均值。median(A,dim)：當：當dim為為1時，該函數(shù)等同于時，該函數(shù)等同于median(A)；當；當dim為為2時，返回一個列向量，其第時，返回一個列向量，其第i個元素是個元素是A的第的第i行的行的中值。中值。例例5-5 分別求向量分別求向量x與與y的平均值和中值。的平均值和中值。105.1.4 累加和與累乘積累加和與累乘積在在MATLAB中，使用中，使用cumsum和和cumprod函數(shù)能方便地求得函數(shù)能方便地求得向量和矩陣元素的累加和與累乘積向量，函數(shù)的調(diào)用格式向量和矩陣元素的累加和與累乘積向

10、量，函數(shù)的調(diào)用格式為：為：cumsum(X)：返回向量：返回向量X累加和向量。累加和向量。cumprod(X)：返回向量：返回向量X累乘積向量。累乘積向量。cumsum(A)：返回一個矩陣，其第：返回一個矩陣，其第i列是列是A的第的第i列的累加和向列的累加和向量。量。cumprod(A)：返回一個矩陣，其第：返回一個矩陣，其第i列是列是A的第的第i列的累乘積列的累乘積向量。向量。cumsum(A,dim)：當：當dim為為1時，該函數(shù)等同于時，該函數(shù)等同于cumsum(A)；當當dim為為2時，返回一個矩陣，其第時，返回一個矩陣，其第i行是行是A的第的第i行的累加行的累加和向量。和向量。cum

11、prod(A,dim)：當：當dim為為1時，該函數(shù)等同于時，該函數(shù)等同于cumprod(A)；當當dim為為2時，返回一個向量，其第時，返回一個向量，其第i行是行是A的第的第i行的累乘行的累乘積向量。積向量。例例5-6 求求s的值。的值。115.1.5 標準方差與相關(guān)系數(shù)標準方差與相關(guān)系數(shù)1求標準方差求標準方差在在MATLAB中，提供了計算數(shù)據(jù)序列的標準方差的函數(shù)中，提供了計算數(shù)據(jù)序列的標準方差的函數(shù)std。對于向量對于向量X，std(X)返回一個標準方差。對于矩陣返回一個標準方差。對于矩陣A，std(A)返回一個行向量，它的各個元素便是矩陣返回一個行向量，它的各個元素便是矩陣A各列或各列或

12、各行的標準方差。各行的標準方差。std函數(shù)的一般調(diào)用格式為：函數(shù)的一般調(diào)用格式為：Y=std(A,flag,dim)其中其中dim取取1或或2。當。當dim=1時，求各列元素的標準方差；當時，求各列元素的標準方差；當dim=2時，則求各行元素的標準方差。時，則求各行元素的標準方差。flag取取0或或1，當，當flag=0時，按時，按1所列公式計算標準方差，當所列公式計算標準方差，當flag=1時，按時，按2所列公式計算標準方差。缺省所列公式計算標準方差。缺省flag=0，dim=1。例例5-7 對二維矩陣對二維矩陣x，從不同維方向求出其標準方差。，從不同維方向求出其標準方差。122相關(guān)系數(shù)相關(guān)

13、系數(shù)MATLAB提供了提供了corrcoef函數(shù)，可以求出數(shù)函數(shù)，可以求出數(shù)據(jù)的相關(guān)系數(shù)矩陣。據(jù)的相關(guān)系數(shù)矩陣。corrcoef函數(shù)的調(diào)用格函數(shù)的調(diào)用格式為：式為：corrcoef(X)：返回從矩陣：返回從矩陣X形成的一個相關(guān)形成的一個相關(guān)系數(shù)矩陣。此相關(guān)系數(shù)矩陣的大小與矩陣系數(shù)矩陣。此相關(guān)系數(shù)矩陣的大小與矩陣X一樣。它把矩陣一樣。它把矩陣X的每列作為一個變量，的每列作為一個變量，然后求它們的相關(guān)系數(shù)。然后求它們的相關(guān)系數(shù)。corrcoef(X,Y)：在這里，：在這里，X,Y是向量，它們是向量，它們與與corrcoef(X,Y)的作用一樣。的作用一樣。13例例5-8 生成滿足正態(tài)分布的生成滿足

14、正態(tài)分布的100005隨機矩隨機矩陣，然后求各列元素的均值和標準方差，陣，然后求各列元素的均值和標準方差，再求這再求這5列隨機數(shù)據(jù)的相關(guān)系數(shù)矩陣。列隨機數(shù)據(jù)的相關(guān)系數(shù)矩陣。命令如下：命令如下：X=randn(10000,5);M=mean(X)D=std(X)R=corrcoef(X)145.1.6 排序排序MATLAB中對向量中對向量X是排序函數(shù)是是排序函數(shù)是sort(X)，函數(shù)返，函數(shù)返回一個對回一個對X中的元素按升序排列的新向量。中的元素按升序排列的新向量。sort函數(shù)也可以對矩陣函數(shù)也可以對矩陣A的各列或各行重新排序，的各列或各行重新排序，其調(diào)用格式為：其調(diào)用格式為：Y,I=sort(

15、A,dim)其中其中dim指明對指明對A的列還是行進行排序。若的列還是行進行排序。若dim=1，則按列排；若則按列排；若dim=2，則按行排。，則按行排。Y是排序后的是排序后的矩陣，而矩陣，而I記錄記錄Y中的元素在中的元素在A中位置。中位置。例例5-9 對二維矩陣做各種排序。對二維矩陣做各種排序。155.2 數(shù)據(jù)插值數(shù)據(jù)插值5.2.1 一維數(shù)據(jù)插值一維數(shù)據(jù)插值在在MATLAB中，實現(xiàn)這些插值的函數(shù)是中，實現(xiàn)這些插值的函數(shù)是interp1，其調(diào)用格式為：其調(diào)用格式為：Y1=interp1(X,Y,X1,method)函數(shù)根據(jù)函數(shù)根據(jù)X,Y的值，計算函數(shù)在的值，計算函數(shù)在X1處的值。處的值。X,Y

16、是是兩個等長的已知向量，分別描述采樣點和樣本值，兩個等長的已知向量，分別描述采樣點和樣本值，X1是一個向量或標量，描述欲插值的點，是一個向量或標量，描述欲插值的點，Y1是是一個與一個與X1等長的插值結(jié)果。等長的插值結(jié)果。method是插值方法，是插值方法，允許的取值有允許的取值有l(wèi)inear、nearest、cubic、spline。16注意：注意：X1的取值范圍不能超出的取值范圍不能超出X的給定范圍，的給定范圍，否則，會給出否則，會給出“NaN”錯誤。錯誤。例例5-10 用不同的插值方法計算在用不同的插值方法計算在/2點的值。點的值。MATLAB中有一個專門的中有一個專門的3次樣條插值函數(shù)次

17、樣條插值函數(shù)Y1=spline(X,Y,X1)，其功能及使用方法與，其功能及使用方法與函數(shù)函數(shù)Y1=interp1(X,Y,X1,spline)完全相同。完全相同。17例例5-11 某觀測站測得某日某觀測站測得某日6:00時至時至18:00時之間每隔時之間每隔2小時的小時的室內(nèi)外溫度室內(nèi)外溫度()，用，用3次樣條插值分別求得該日室內(nèi)外次樣條插值分別求得該日室內(nèi)外6:30至至17:30時之間每隔時之間每隔2小時各點的近似溫度小時各點的近似溫度()。設(shè)時間變量設(shè)時間變量h為一行向量，溫度變量為一行向量，溫度變量t為一個兩列矩陣，其中為一個兩列矩陣，其中第一列存放室內(nèi)溫度，第二列儲存室外溫度。命令如

18、下：第一列存放室內(nèi)溫度，第二列儲存室外溫度。命令如下：h =6:2:18;t=18,20,22,25,30,28,24;15,19,24,28,34,32,30;XI =6.5:2:17.5YI=interp1(h,t,XI,spline) %用用3次樣條插值計算次樣條插值計算185.2.2 二維數(shù)據(jù)插值二維數(shù)據(jù)插值在在MATLAB中，提供了解決二維插值問題的函數(shù)中，提供了解決二維插值問題的函數(shù)interp2，其調(diào)用格式為：，其調(diào)用格式為：Z1=interp2(X,Y,Z,X1,Y1,method)其中其中X,Y是兩個向量，分別描述兩個參數(shù)的采樣點，是兩個向量，分別描述兩個參數(shù)的采樣點，Z是與

19、參數(shù)采樣點對應(yīng)的函數(shù)值，是與參數(shù)采樣點對應(yīng)的函數(shù)值，X1,Y1是兩個是兩個向量或標量，描述欲插值的點。向量或標量，描述欲插值的點。Z1是根據(jù)相應(yīng)的是根據(jù)相應(yīng)的插值方法得到的插值結(jié)果。插值方法得到的插值結(jié)果。 method的取值與一的取值與一維插值函數(shù)相同。維插值函數(shù)相同。X,Y,Z也可以是矩陣形式。也可以是矩陣形式。同樣，同樣，X1,Y1的取值范圍不能超出的取值范圍不能超出X,Y的給定范圍，的給定范圍，否則，會給出否則，會給出“NaN”錯誤。錯誤。19例例5-12 設(shè)設(shè)z=x2+y2，對，對z函數(shù)在函數(shù)在0,10,2區(qū)域內(nèi)進行插值。區(qū)域內(nèi)進行插值。例例5-13 某實驗對一根長某實驗對一根長10

20、米的鋼軌進行熱源的溫度傳播測米的鋼軌進行熱源的溫度傳播測試。用試。用x表示測量點表示測量點0:2.5:10(米米)，用，用h表示測量時間表示測量時間0:30:60(秒秒)，用，用T表示測試所得各點的溫度表示測試所得各點的溫度()。試用線。試用線性插值求出在一分鐘內(nèi)每隔性插值求出在一分鐘內(nèi)每隔20秒、鋼軌每隔秒、鋼軌每隔1米處的溫度米處的溫度TI。命令如下：命令如下：x=0:2.5:10;h=0:30:60;T=95,14,0,0,0;88,48,32,12,6;67,64,54,48,41;xi=0:10;hi=0:20:60;TI=interp2(x,h,T,xi,hi)205.3 曲線擬合

21、曲線擬合在在MATLAB中，用中，用polyfit函數(shù)來求得最小二乘擬合多項式函數(shù)來求得最小二乘擬合多項式的系數(shù)，再用的系數(shù)，再用polyval函數(shù)按所得的多項式計算所給出的函數(shù)按所得的多項式計算所給出的點上的函數(shù)近似值。點上的函數(shù)近似值。polyfit函數(shù)的調(diào)用格式為：函數(shù)的調(diào)用格式為：P,S=polyfit(X,Y,m)函數(shù)根據(jù)采樣點函數(shù)根據(jù)采樣點X和采樣點函數(shù)值和采樣點函數(shù)值Y，產(chǎn)生一個，產(chǎn)生一個m次多項式次多項式P及其在采樣點的誤差向量及其在采樣點的誤差向量S。其中。其中X,Y是兩個等長的向量，是兩個等長的向量，P是一個長度為是一個長度為m+1的向量，的向量，P的元素為多項式系數(shù)。的元

22、素為多項式系數(shù)。polyval函數(shù)的功能是按多項式的系數(shù)計算函數(shù)的功能是按多項式的系數(shù)計算x點多項式的值，點多項式的值，將在將在5.5.3節(jié)中詳細介紹。節(jié)中詳細介紹。21例例5-14 已知數(shù)據(jù)表已知數(shù)據(jù)表t,y，試求，試求2次擬合多項式次擬合多項式p(t)，然后求，然后求ti=1,1.5,2,2.5,9.5,10各點的函數(shù)近似值。各點的函數(shù)近似值。225.4 離散傅立葉變換離散傅立葉變換5.4.1 離散傅立葉變換算法簡要離散傅立葉變換算法簡要5.4.2 離散傅立葉變換的實現(xiàn)離散傅立葉變換的實現(xiàn)一維離散傅立葉變換函數(shù)，其調(diào)用格式與功一維離散傅立葉變換函數(shù)，其調(diào)用格式與功能為：能為：(1) fft

23、(X)：返回向量：返回向量X的離散傅立葉變換。設(shè)的離散傅立葉變換。設(shè)X的長度的長度(即元素個數(shù)即元素個數(shù))為為N，若，若N為為2的冪次，的冪次，則為以則為以2為基數(shù)的快速傅立葉變換，否則為為基數(shù)的快速傅立葉變換，否則為運算速度很慢的非運算速度很慢的非2冪次的算法。對于矩陣冪次的算法。對于矩陣X，fft(X)應(yīng)用于矩陣的每一列。應(yīng)用于矩陣的每一列。23(2) fft(X,N)：計算：計算N點離散傅立葉變換。它限點離散傅立葉變換。它限定向量的長度為定向量的長度為N，若，若X的長度小于的長度小于N，則，則不足部分補上零；若大于不足部分補上零；若大于N，則刪去超出，則刪去超出N的那些元素。對于矩陣的那

24、些元素。對于矩陣X，它同樣應(yīng)用于，它同樣應(yīng)用于矩陣的每一列，只是限定了向量的長度為矩陣的每一列，只是限定了向量的長度為N。(3) fft(X,dim)或或fft(X,N,dim)：這是對于矩：這是對于矩陣而言的函數(shù)調(diào)用格式，前者的功能與陣而言的函數(shù)調(diào)用格式，前者的功能與FFT(X)基本相同，而后者則與基本相同，而后者則與FFT(X,N)基基本相同。只是當參數(shù)本相同。只是當參數(shù)dim=1時，該函數(shù)作時，該函數(shù)作用于用于X的每一列；當?shù)拿恳涣?；當dim=2時，則作用于時，則作用于X的每一行。的每一行。24值得一提的是，當已知給出的樣本數(shù)值得一提的是，當已知給出的樣本數(shù)N0不是不是2的冪次時，可以取

25、一個的冪次時，可以取一個N使它大于使它大于N0且是且是2的冪次，然后利用函數(shù)格式的冪次，然后利用函數(shù)格式fft(X,N)或或fft(X,N,dim)便可進行快速傅立葉變換。這便可進行快速傅立葉變換。這樣，計算速度將大大加快。樣，計算速度將大大加快。相應(yīng)地，一維離散傅立葉逆變換函數(shù)是相應(yīng)地，一維離散傅立葉逆變換函數(shù)是ifft。ifft(F)返回返回F的一維離散傅立葉逆變換；的一維離散傅立葉逆變換；ifft(F,N)為為N點逆變換；點逆變換；ifft(F,dim)或或ifft(F,N,dim)則由則由N或或dim確定逆變換的點確定逆變換的點數(shù)或操作方向。數(shù)或操作方向。25例例5-15 給定數(shù)學函數(shù)

26、給定數(shù)學函數(shù)x(t)=12sin(210t+/4)+5cos(240t)取取N=128，試對，試對t從從01秒采樣，用秒采樣，用fft作快速作快速傅立葉變換，繪制相應(yīng)的振幅傅立葉變換，繪制相應(yīng)的振幅-頻率圖。頻率圖。在在01秒時間范圍內(nèi)采樣秒時間范圍內(nèi)采樣128點，從而可以確點，從而可以確定采樣周期和采樣頻率。由于離散傅立葉定采樣周期和采樣頻率。由于離散傅立葉變換時的下標應(yīng)是從變換時的下標應(yīng)是從0到到N-1，故在實際應(yīng)，故在實際應(yīng)用時下標應(yīng)該前移用時下標應(yīng)該前移1。又考慮到對離散傅立。又考慮到對離散傅立葉變換來說，其振幅葉變換來說，其振幅| F(k)|是關(guān)于是關(guān)于N/2對稱對稱的，故只須使的，

27、故只須使k從從0到到N/2即可。即可。26程序如下：程序如下：N=128; % 采樣點數(shù)采樣點數(shù)T=1; % 采樣時間終點采樣時間終點t=linspace(0,T,N); % 給出給出N個采樣時間個采樣時間ti(I=1:N)x=12*sin(2*pi*10*t+pi/4)+5*cos(2*pi*40*t); % 求各采樣點求各采樣點樣本值樣本值xdt=t(2)-t(1); % 采樣周期采樣周期f=1/dt; % 采樣頻率采樣頻率(Hz)X=fft(x); % 計算計算x的快速傅立葉變換的快速傅立葉變換XF=X(1:N/2+1); % F(k)=X(k)(k=1:N/2+1)f=f*(0:N/2

28、)/N; % 使頻率軸使頻率軸f從零開始從零開始plot(f,abs(F),-*) % 繪制振幅繪制振幅-頻率圖頻率圖xlabel(Frequency);ylabel(|F(k)|)275.5 多項式計算多項式計算5.5.1 多項式的四則運算多項式的四則運算1多項式的加減運算多項式的加減運算2多項式乘法運算多項式乘法運算函數(shù)函數(shù)conv(P1,P2)用于求多項式用于求多項式P1和和P2的乘積。的乘積。這里，這里，P1、P2是兩個多項式系數(shù)向量。是兩個多項式系數(shù)向量。例例5-16 求多項式求多項式x4+8x3-10與多項式與多項式2x2-x+3的乘積的乘積。283多項式除法多項式除法函數(shù)函數(shù)Q,

29、r=deconv(P1,P2)用于對多項式用于對多項式P1和和P2作除作除法運算。其中法運算。其中Q返回多項式返回多項式P1除以除以P2的商式，的商式，r返回返回P1除以除以P2的余式。這里，的余式。這里，Q和和r仍是多項式仍是多項式系數(shù)向量。系數(shù)向量。deconv是是conv的逆函數(shù)，即有的逆函數(shù)，即有P1=conv(P2,Q)+r。例例5-17 求多項式求多項式x4+8x3-10除以多項式除以多項式2x2-x+3的結(jié)的結(jié)果。果。295.5.2 多項式的導函數(shù)多項式的導函數(shù)對多項式求導數(shù)的函數(shù)是：對多項式求導數(shù)的函數(shù)是：p=polyder(P)：求多項式：求多項式P的導函數(shù)的導函數(shù)p=pol

30、yder(P,Q)：求：求PQ的導函數(shù)的導函數(shù)p,q=polyder(P,Q)：求：求P/Q的導函數(shù)，導函數(shù)的分的導函數(shù)，導函數(shù)的分子存入子存入p，分母存入，分母存入q。上述函數(shù)中，參數(shù)上述函數(shù)中，參數(shù)P,Q是多項式的向量表示，結(jié)果是多項式的向量表示，結(jié)果p,q也是多項式的向量表示。也是多項式的向量表示。30例例5-18 求有理分式的導數(shù)。求有理分式的導數(shù)。命令如下：命令如下：P=1;Q=1,0,5;p,q=polyder(P,Q)315.5.3 多項式的求值多項式的求值MATLAB提供了兩種求多項式值的函數(shù)：提供了兩種求多項式值的函數(shù)：polyval與與polyvalm，它們的輸入?yún)?shù)均為多

31、項式系數(shù)向，它們的輸入?yún)?shù)均為多項式系數(shù)向量量P和自變量和自變量x。兩者的區(qū)別在于前者是代數(shù)多項。兩者的區(qū)別在于前者是代數(shù)多項式求值，而后者是矩陣多項式求值。式求值，而后者是矩陣多項式求值。321代數(shù)多項式求值代數(shù)多項式求值polyval函數(shù)用來求代數(shù)多項式的值，其調(diào)用函數(shù)用來求代數(shù)多項式的值，其調(diào)用格式為：格式為：Y=polyval(P,x)若若x為一數(shù)值，則求多項式在該點的值；若為一數(shù)值，則求多項式在該點的值；若x為向量或矩陣，則對向量或矩陣中的每個為向量或矩陣，則對向量或矩陣中的每個元素求其多項式的值。元素求其多項式的值。例例5-19 已知多項式已知多項式x4+8x3-10，分別取，分別取x=1.2和一個和一個23矩陣為自變量計算該多