gaimatlab數(shù)據(jù)分析與多項(xiàng)式計(jì)算-在線瀏覽
2025-07-13 14:03本頁面
【正文】 3125 9765625 ? S=sum(y) ? S =12207031 標(biāo)準(zhǔn)方差與相關(guān)系數(shù) 1．求標(biāo)準(zhǔn)方差 在 MATLAB中，提供了計(jì)算數(shù)據(jù)序列的 標(biāo)準(zhǔn)方差的函數(shù) std。對于矩陣 A， std(A)返回一個(gè)行向量，它的各個(gè)元素便是矩陣 A各列或各行的標(biāo)準(zhǔn)方差。當(dāng) dim=1時(shí)，求各列元素的標(biāo)準(zhǔn)方差；當(dāng) dim=2時(shí)，則求各行元素的標(biāo)準(zhǔn)方差。 缺省 flag=0， dim=1。 corrcoef函數(shù)的調(diào)用格式為： corrcoef(X)： 返回從矩陣 X形成的一個(gè)相關(guān)系數(shù)矩陣。 它把矩陣 X的每列作為一個(gè)變量，然后求它們的相關(guān)系數(shù)。 (P137) 例 68 生成滿足正態(tài)分布的 10 5隨機(jī)矩陣，然后求各列元素的均值和標(biāo)準(zhǔn)方差，再求這 5列隨機(jī)數(shù)據(jù)的相關(guān)系數(shù)矩陣。 M=mean(X) D=std(X) R=corrcoef(X) 排序 MATLAB中對向量 X是 排序函數(shù)是 sort(X)，函數(shù)返回一個(gè)對 X中的元素 按升序排列 的新向量。若dim=1，則按列排；若 dim=2，則按行排。 ? A=[3 1 9。6 2 7] ? B=sort(A) ? C=sort(A,2) ? [X,I]=sort(A) ? A =A=[3 1 9。6 2 7] ? 3 1 9 ? 8 3 1 ? 6 2 7 ? B=sort(A) ? 3 1 1 ? 6 2 7 ? 8 3 9 ? C=sort(A,2) ? 9 3 1 ? 8 3 1 ? 7 6 2 ? [X,I]=sort(A) ? X = ? 3 1 1 ? 6 2 7 ? 8 3 9 ? I = ? 1 1 2 ? 3 3 3 ? 2 2 1 數(shù)據(jù)插值 一維數(shù)據(jù)插值 在 MATLAB中，實(shí)現(xiàn)這些 插值的函數(shù)是interp1，其調(diào)用格式為： Y1=interp1(X,Y,X1,39。) 函數(shù)根據(jù) X,Y的值，計(jì)算函數(shù)在 X1處的值。 method是插值方法，允許的取值有‘ linear’、‘ nearest’、‘ cubic’、‘ spline’。 例 610 用不同的插值方法計(jì)算在 π/2點(diǎn)的值。 例 611 某觀測站測得某日 6:00時(shí)至 18:00時(shí)之間每隔 2小時(shí)的室內(nèi)外溫度 (℃ )，用 3次樣條插值分別求得該日室內(nèi)外 6:30至 17:30時(shí)之間每隔 2小時(shí)各點(diǎn)的近似溫度 (℃ )。命令如下： h =6:2:18。15,19,24,28,34,32,30]39。method39。 Z1是根據(jù)相應(yīng)的插值方法得到的插值結(jié)果。 X,Y,Z也可以是矩陣形式。 例 613 某實(shí)驗(yàn)對一根長 10米的鋼軌進(jìn)行熱源的溫度傳播測試。試用線性插值求出在一分鐘內(nèi)每隔 20秒、鋼軌每隔 1米處的溫度 TI。 h=[0:30:60]39。88,48,32,12,6。 xi=[0:10]。 TI=interp2(x,h,T,xi,hi) 曲線擬合 在 MATLAB中，用 polyfit函數(shù)來求得最小二乘擬合多項(xiàng)式的系數(shù)，再用 polyval函數(shù)按所得的多項(xiàng)式計(jì)算所給出的點(diǎn)上的函數(shù)近似值。其中X,Y是兩個(gè)等長的向量， P是一個(gè)長度為 m+1的向量， P的元素為多項(xiàng)式系數(shù)。 例 614 已知數(shù)據(jù)表 [t,y]，試求 2次擬合多項(xiàng)式 p(t)，然后求ti=1,2,…,10 各點(diǎn)的函數(shù)近似值。設(shè) X的長度 (即元素個(gè)數(shù) )為 N，若 N為 2的冪次，則為以 2為基數(shù)的快速傅立葉變換 ，否則為運(yùn)算 速度很慢的非 2冪次的算法 。 (2) fft(X,N)： 計(jì)算 N點(diǎn)離散傅立葉變換。對于矩陣 X，它同樣應(yīng)用于矩陣的每一列，只是限定了向量的長度為 N。只是當(dāng)參數(shù) dim=1時(shí)，該函數(shù)作用于 X的每一列；當(dāng) dim=2時(shí)，則作用于 X的每一行。這樣，計(jì)算速度將大大加快。ifft(F)返回 F的一維離散傅立葉逆變換；ifft(F,N)為 N點(diǎn)逆變換； ifft(F,[],dim)或ifft(F,N,dim)則由 N或 dim確定逆變換的點(diǎn)數(shù)或操作方向。 在 0~1秒時(shí)間范圍內(nèi)采樣 128點(diǎn)，從而可以確定采樣周期和采樣頻率。又考慮到對離散傅立葉變換來說，其振幅 | F(k)|是關(guān)于 N/2對稱的，
公司管理相關(guān)推薦
鄂ICP備17016276號-1