このページの翻訳は最新ではありません。ここをクリックして,英語の最新版を参照してください。
polyfit
多項式の曲線近似
説明
例
三角関数の多項式近似
[0, 4 *π)の区間内で正弦曲線に沿って等間隔に分布する10個の点を生成します。
x = linspace(0, 4 *π,10);y = sin (x);
polyfitを使用して7次の多項式を点にあてはめます。
p = polyfit (x, y, 7);
より細かいグリッド上で多項式を計算して結果をプロットします。
x1 = linspace(0, 4 *π);日元= polyval (p, x1);图绘制(x, y,“o”)举行在情节(x1, y1)从
点の集合の多項式近似
[0, 1]の区間内で等間隔に配置された5つの点のベクトルを作成し,それらの点で
を評価します。
x = linspace (0, 1, 5);y = 1. / (1 + x);
4次の多項式をこの5つの点にあてはめます。一般にn個の点に対して,それらの点を正確に通過するn - 1次の多項式をあてはめられます。
p = polyfit (x, y, 4);
0と2の間のより細かい点のグリッド上で,元の関数と多項式近似を計算します。
x1 = linspace (0, 2);日元= 1. / (1 + x1);f1 = polyval (p, x1);
元の関数の値と多項式近似をより広い区間(0, 2)でプロットし,多項式近似を求めるために使用した点を円で強調表示します。多項式近似は元の[0, 1]の区間では良好ですが,この区間の外側では近似させた関数から急速に離れています。
图绘制(x, y,“o”)举行在情节(x1, y1)情节(x1, f1,“r——”)传说(“y”，“日元”，“f1”）
誤差関数の多項式近似
最初に,区間[0, 2.5]で等間隔に分布するx個の点からなるベクトルを作成し,これらの点で小块土地(x)を評価します。
x =(0:0.1:2.5)”;y =小块土地(x);
6次の近似多項式の係数を求めます。
p = polyfit (x, y, 6)
p =1×70.0084 -0.0983 0.4217 -0.7435 0.1471 1.1064 0.0004
近似の程度を見るために,データ点で多項式を計算し,データ,近似値,誤差を表示する表を生成します。
f = polyval (p (x);T =表(x, y, f, yf,“VariableNames”, {“X”，“Y”，“健康”，“FitError”}）
T =26日×4表X Y Fit fiiterror __________ __________ ___________ 00 0.00044117 -0.00044117 0.1 0.11246 0.11185 0.00060836 0.2 0.2227 0.22231 0.00039189 0.3 0.32863 0.32872 -9.7429e-05 0.4 0.42839 0.4288 -0.00040661 0.5 0.5205 0.52093 -0.00042568 0.6 0.60386 0.60408 -0.00022824 0.7 0.6778 0.67775 4.6383e-05 0.8 0.7421 0.74183 0.00026992 0.9 0.79691 0.796541.0.93401 0.93422 -0.000211 1.4 0.95229 0.95258 -0.00029933 1.5 0.96611 0.96639 -0.00028097⋮
この区間内では,内挿値と実際の値は,かなり一致しています。この区間の外側で外挿値が実際のデータからどの程度急速に離れるのかを見るためにプロットを作成します。
x1 = (0:0.1:5) ';日元=小块土地(x1);f1 = polyval (p, x1);图绘制(x, y,“o”)举行在情节(x1, y1,“- - -”)情节(x1, f1,“r——”)轴([0 5 0 2])保持从
センタリングとスケーリングを使用した数値特性の改善
1750 ~ 2000年の人口のデータのテーブルを作成し,データ点をプロットします。
年= (1750:25:2000)';Pop = 1e6*[791 856 978 1050 1262 1544 1650 2532 6122 8170 11560]';T = table(year, pop)
T =11×2表年流行____ _________ 1750 7.91e+08 1775 8.56e+08 1800 9.78e+08 1825 1.05e+09 1850 1.262 2e+09 1875 1.544e+09 1900 1.65e+09 1925 2.532 2e+09 1950 6.122e+09 1975 8.17e+09 2000 1.156e+10
情节(年,流行,“o”）
3出力のpolyfitを使用し,センタリングとスケーリングを使用して5次の多項式をあてはめます。これにより問題の数値特性が向上します。polyfitは一年のデータを0にセンタリングし,標準偏差が1になるようにスケーリングします。これにより近似計算において悪条件のヴァンデルモンド行列を避けることができます。
(p ~μ)= polyfit (T。年T.pop 5);
polyvalを4入力で使用し,pをスケーリングされた年(year-mu(1)) /μ(2)に対して評価します。結果を元の年に対してプロットします。
f = polyval (p,年,[],μ);持有在情节(一年,f)从
単純な線形回帰
単純な線形回帰モデルを,一連の離散2次元データ点に近似します。
サンプルデータ点のベクトル(x, y)をいくつか作成します。1次多項式をデータに近似します。
x = 1:50;Y = -0.3*x + 2*randn(1,50); / /p = polyfit (x, y, 1);
近似多項式pをxの点で評価します。その結果得られる線形回帰モデルをデータでプロットします。
f = polyval (p (x);情节(x, y,“o”f, x,,“- - -”)传说(“数据”，“线性适应”）
誤差推定を含む線形回帰
線形モデルでデータ点セットを近似し,95%の予測区間の推定を含めて結果をプロットします。
サンプルデータ点のベクトル(x, y)をいくつか作成します。polyfitを使用して1次多項式でデータを近似します。線形近似の係数と誤差推定値の構造体を返す2つの出力を指定します。
x = 1:10 0;Y = -0.3*x + 2*randn(1100);[p, S] = polyfit (x, y, 1);
xの各点でpの1次多項式近似を評価します。polyvalの3番目の入力として誤差推定の構造体を指定し,標準誤差の推定値を計算します。標準誤差の推定値はδに返されます。
[y_fitδ]= polyval (p, x, S);
元のデータ,線形近,似および95%の予測区間 をプロットします。
情节(x, y,“波”)举行在情节(x, y_fit,的r -)情节(x, y_fit + 2 *δ,“m——”, x, y_fit-2 *δ,“m——”)标题(“95%预测区间的数据线性拟合”)传说(“数据”，“线性适应”，“95%的预测区间”）
入力引数
出力引数
制限
多くの点をもつ問題に対して
polyfitを使用して多項式近似の次元数を増やしても,より優れた近似が得られるとは限りません。高次の多項式はデータ点の間で振動する場合があり,データの“より劣った”近似になる場合があります。このような場合は,低次の多項式近似(点の間がより滑らかになる傾向がある)や別の手法を問題に応じて使用します。多項式は本質的に非有界の振動関数です。このため,有界のデータや,単調増加または単調減少するデータの外挿にはあまり適していません。
アルゴリズム
polyfitはxを使用してn + 1列およびm =长度(x)行のヴァンデルモンド行列Vを形成し,次の線形方程式を導きます。
polyfitはこれをp = V \ yによって解きます。ヴァンデルモンド行列の列はベクトルxのべき乗であるため,Vの条件数は高次の近似の場合にはしばしば大きくなり,結果として特異係数行列になります。このような場合,センタリングとスケーリングによって方程式の数値特性が向上し,信頼性の高い近似が得られる場合があります。
拡張機能
你也可以从以下列表中选择一个网站:
