gaimatlab數(shù)據(jù)分析與多項式計算(文件)
【正文】 oef函數(shù)，可以求出數(shù)據(jù)的相關(guān)系數(shù)矩陣。 corrcoef(X,Y)： 在這里， X,Y是向量，它們與corrcoef([X,Y])的作用一樣。 sort函數(shù)也可以對矩陣 A的各列或各行重新排序，其調(diào)用格式為： [Y,I]=sort(A,dim) 其中 dim指明對 A的列還是行進(jìn)行排序。8 3 1。method39。 注意： X1的取值范圍不能超出 X的給定范圍，否則，會給出“ NaN”錯誤。 設(shè)時間變量 h為一行向量，溫度變量 t為一個兩列矩陣，其中第一列存放室內(nèi)溫度，第二列儲存室外溫度。 XI =:2: YI=interp1(h,t,XI,‘spline’) % 用 3次樣條插值計算 二維數(shù)據(jù)插值 在 MATLAB中，提供了解決 二維插值問題的函數(shù) interp2，其調(diào)用格式為： Z1=interp2(X,Y,Z,X1,Y1,39。 method的取值與一維插值函數(shù)相同。用 x表示測量點 0::10(米 )，用 h表示測量時間 0:30:60(秒 )，用 T表示測試所得各點的溫度 (℃ )。 T=[95,14,0,0,0。 hi=[0:20:60]39。 polyval函數(shù)的功能是按多項式的系數(shù)計算 x點多項式的值，將在 。對于 矩陣 X，fft(X)應(yīng)用于矩陣的每一列。 (3) fft(X,[],dim)或 fft(X,N,dim)： 這是對于矩陣而言的函數(shù)調(diào)用格式，前者的功能與 FFT(X)基本相同，而后者則與 FFT(X,N)基本相同。 相應(yīng)地，一維離散傅立葉 逆變換函數(shù)是 ifft。由于離散傅立葉變換時的下標(biāo)應(yīng)是從 0到 N1，故在實際應(yīng)用時下標(biāo)應(yīng)該前移 1。 % 采樣時間終點 t=linspace(0,T,N)。 % 采樣頻率 (Hz) X=fft(x)。*39。 ylabel(39。 例 616 求多項式 x4+8x310與多項式 2x2x+3的乘積。其中 Q返回多項式 P1除以P2的商式， r返回 P1除以 P2的余式。B=[2 1 3]。 例 618 求有理分式的導(dǎo)數(shù)。兩者的區(qū)別在于前者是 代數(shù)多項式求值，而后者是矩陣多項式求值。 ? X=。 多項式求根 n次多項式具有 n個根，當(dāng)然這些根可能是實根，也可能含有若干對共軛復(fù)根。 x=roots(A) 若已知多項式的全部根，則可以用 poly函數(shù)建立起該多項式，其調(diào)用格式為： P=poly(x) 若 x為具有 n個元素的向量，則 poly(x)建立以 x為其根的多項式，且將該多項式的系數(shù)賦給向量 P。 命令如下： P=[3,0,4,5,5]。 ? x=roots(A) ? x = ? ? ? + ? ? P=poly(x) P = ? 例 622 已知 f(x)= (1) 計算 f(x)=0 的全部根。 例 621 求多項式 x4+8x310的根。polyvalm函數(shù)要求 x為方陣，它以方陣為自變量求多項式的值。 例 619 已知多項式 x4+8x310，分別取 x=個 2 3矩陣為自變量計算該多項式的值。 Q=[1,0,5]。Q 的導(dǎo)函數(shù) [p,q]=polyder(P,Q)：求 P/Q的導(dǎo)函數(shù)，導(dǎo)函數(shù)的分子存入 p，分母存入 q。 deconv是 conv的逆函數(shù)，即有P1=conv(P2,Q)+r。B=[2 1 3]。) 多項式計算 多項式的四則運(yùn)算 1．多項式的加減運(yùn)算 (注意 :同冪次的系數(shù)對齊 ) 2．多項式乘法運(yùn)算，函數(shù) conv(P1,P2)用于求多項式 P1和 P2的乘積。Frequency39。 % F(k)=X(k)(k=1:N/2+1) f=f*(0:N/2)/N。 % 求各采樣點樣本值 x dt=t(2)t(1)。 程序如下： N=128。 例 615 給定數(shù)學(xué)函數(shù) x(t)=12sin(2π 10t+π/4)+5cos(2π 40t) 取 N=128，試對 t從 0~1秒采樣，用 fft作快速傅立葉變換，繪制相應(yīng)的振幅 頻率圖。 注意 :當(dāng)已知 給出的樣本數(shù) N0不是 2的冪次時，可以取一個 N使它大于 N0且是 2的冪次， 然后利用函數(shù)格式 fft(X,N)或 fft(X,N,dim)便可進(jìn)行快速傅立葉變換。它限定 向量的長度為 N，若 X的長度小于 N，則不足部分補(bǔ)上零；若大于 N，則刪去超
公司管理相關(guān)推薦
鄂ICP備17016276號-1