matlab華夏colleg第五章-數(shù)值計算(編輯修改稿)
2025-08-20 13:36 本頁面
【文章內容簡介】 YI=interp1(h,t,XI,‘spline’) %用 3次樣條插值計算 二維數(shù)據(jù)插值 在 MATLAB中，提供了解決二維插值問題的函數(shù)interp2，其調用格式為： Z1=interp2(X,Y,Z,X1,Y1,39。method39。) 其中 X,Y是兩個向量，分別描述兩個參數(shù)的采樣點， Z是與參數(shù)采樣點對應的函數(shù)值， X1,Y1是兩個向量或標量，描述欲插值的點。 Z1是根據(jù)相應的插值方法得到的插值結果。 method的取值與一維插值函數(shù)相同。X,Y,Z也可以是矩陣形式。 同樣， X1,Y1的取值范圍不能超出 X,Y的給定范圍，否則，會給出 “ NaN”錯誤。 例 某實驗對一根長 10米的鋼軌進行熱源的溫度傳播測試。用 x表示測量點 0::10(米 )，用 h表示測量時間0:30:60(秒 )，用 T表示測試所得各點的溫度 (℃) 。 T=[95,14,0,0,0。88,48,32,12,6。67,64,54,48,41]。 試用線性插值求出在一分鐘內每隔 20秒、鋼軌每隔 1米處的溫度 TI。 命令如下： x=0::10。 h=[0:30:60]39。 T=[95,14,0,0,0。88,48,32,12,6。67,64,54,48,41]。 xi=[0:10]。 hi=[0:20:60]39。 TI=interp2(x,h,T,xi,hi) 曲線擬合 在 MATLAB中，給定一組數(shù)據(jù)進行曲線擬合，需要考慮兩個問題： ? 用什么曲線來擬合 ? 什么是最佳擬合方案 多項式： x4+8x310 2x2x+3 f(x)=a0xn+a1xn1+a2xn2+...+an1xn+an 多項式是一種基本的數(shù)據(jù)分析工具 ,很多復雜的函數(shù)都可 用多項式逼近。 P=[a0 a1 a2...an1 an] 數(shù)據(jù)插值及曲線擬合 曲線擬合 在 MATLAB中，用 polyfit函數(shù)來求得最小二乘擬合多項式的系數(shù)。 polyfit函數(shù)的調用格式為： [P， S]=polyfit(X,Y,m) 函數(shù)根據(jù)采樣點 X和采樣點函數(shù)值 Y，產生一個 m次多項式 P及其在采樣點的誤差向量 S。其中 X,Y是兩個等長的向量， P是一個長度為 m+1的向量， P的元素為多項式系數(shù)。 曲線擬合 再用 polyval函數(shù)按所得的多項式計算所給出的點上的函數(shù)近似值。 其調用格式是： y=polyval(P,X) 例 :用一個 8次多項式在區(qū)間 [0,3]內逼近函數(shù) sin(x) 例 :比較不同階數(shù)的擬合效果。 例 :比較區(qū)間內外擬合效果。 教學內容： ? 數(shù)據(jù)處理和分析 ? 數(shù)據(jù)插值與曲線擬合 ? 多項式計算 ? 解方程與函數(shù)極值 第五章 數(shù)值計算 多項式的建立 M采用向量表示多項式，將多項式的系數(shù)按降冪次序存放在向量中。 f(x)=a0xn+a1xn1+a2xn2+...+an1xn+an P=[a0 a1 a2...an1 an] 例如： f(x)=x32x5 P=[1 0 2 5] 多項式求根 n次多項式具有 n個根，當然這些根可能是實根，也可能含有若干對共軛復根。 MATLAB提供的 roots函數(shù)用于求多項式的全部根，其調用格式為： x=roots(P) 其中 P為多項式的系數(shù)向量，求得的根賦給向量 x，即 x(1),x(2),… ,x(n)分別代表多項式的 n個根。 P=poly(x) 若 x為具有 n個元素的向量，則 poly(x)建立以 x為其根的多項式，且將該多項式的系數(shù)賦給向量 P。 例：已知 f(x) (1) 計算 f(x)=0 的全部根。 (2) 由方程 f(x)=0的根構造一個多項式 g(x)，并與f(x)進行對比。 命令如下： P=[3,0,4,5,5]。 X=roots(P) %求方程 f(x)=0的根 G=poly(X) %求多項式 g(x) 多項式的求值 polyval函數(shù)用來求代數(shù)多項式的值，其調用格式為： Y=polyval(P,x) 若 x為一數(shù)值，則求多項式在該點的值；若 x為向量或矩陣，則對向量或矩陣中的每個元素求其多項式的值。 例 已知多項式 x4+8x310 ，分別取 x=2 3矩陣為自變量計算該多項式的值。 例：計算多項式的值 多項式的四則運算 1．多項
畢業(yè)設計相關推薦
鄂ICP備17016276號-1