Matlab学习笔记2011-09-16【数据分析】

1. 插值（Interpolation）
   一维插值命令interp1，二维插值命令interp2
    1) 一维插值命令
    yi=interp1(x,y,xi)  x必须是向量，y可以是向量也可以是矩阵。如果y是向量，则必须与变量x具有相同的长度，这时xi可以是标量，也可以是向量或者矩阵。
    yi=interp1(y,xi)  在默认情况下，x变量选择为1：n，其中n是向量y的长度
    yi=interp1(x,y,xi,method)  输入变量method用于制定插值的方法
    yi=interp1(x,y,xi,method,'extrap')  对超出数据范围的插值运算使用外推方法
    yi=interp1(x,y,xi,method,extrapval)  对超出数据范围的插值数据返回extrapval数值，一般为NaN或0
    pp=interp1(x,y,method,'pp') 返回数值pp为数据y的分段多项式形式，method制定产生分段多项式pp的方法
   还提供了interp1q命令作为interp1的快速命令。对于数据及那个不均的数据，interp1q命令运行的速度比interp1速度要快，主要原因在于interp1q命令没有对输入数据进行检测。


   具体的插值方法：
    nearest：最邻近插值方法（Nearest neighbor interpolation） 
    linear:线性插值（Linear interpolation）
    spline:三次样条插值（Cubic spline interpolation）
    pchip:分段三次厄米多项式插值（Piecewise cubic Hermite interpolation）
    cubic:  三次多项式插值，与分段三次厄米多项式插值方法相同
    v5cubic :Matlab5中使用的三次多项式插值
 -----
%一维插值
%已知数据
t=1900:10:1990;
p=[75.995 91.972 105.711 123.203 131.222 150.222 179.222 203.232 226.222 249.633];
%使用不同的方法进行插值运算
x=1900:0.01:1990;
yi_linear=interp1(t,p,x);
yi_spline=interp1(t,p,x,'spline');
yi_cubic=interp1(t,p,x,'cubic');
yi_v5cubic=interp1(t,p,x,'v5cubic');
%绘制对比图形
subplot(2,1,1);
plot(t,p,'ko');hold on;
plot(x,yi_linear,'g','LineWidth',1.5); hold on;
plot(x,yi_spline,'y','LineWidth',1.5);
 grid on;
%创建新的图形窗口
subplot(2,1,2);
 yi_nearest=interp1(t,p,x,'nearest');
 plot(t,p,'ko'); hold on;
 plot(x,yi_nearest);
 grid on;
---
 
各种插值方法的比较
  Nearest   最快的插值方法，但是数据平滑方面最差，数据是不连续的
  Linear        执行速度较快，有足够的精度。最为常用，而且是默认设置
  Cubic       较慢，精度高，平滑度好，当希望得到平滑的曲线时可以使用
  Spline     执行速度最慢，精度高，最平滑
  2） 二维插值命令
   ZI=interp2(X,Y,Z,XI,YI)  原始数据x,y和z决定插值函数z=f(x,y),返回的数值ZI是（XI，YI）在函数z=f(x,y)上的数值
  2. 曲线拟合
    1）多项式拟合
      matlab提供polyfit函数来计算多项式拟合系数，其设定曲线拟合的目标是最小方差（Least Squares）或者被成为最小二乘法。
    [p,S,mu]=ployfit(x,y,n)
        [y,delta]=polyval(p,x,S,mu)
      其中，x和y表示的是已知测量数据，n是多项式拟合的阶数。
----
>> x=(0:0.1:5)';
>> y=erf(x);
>> help erf
 ERF Error function.
    Y = ERF(X) is the error function for each element of X.  X must be
    real. The error function is defined as:
 
      erf(x) = 2/sqrt(pi) * integral from 0 to x of exp(-t^2) dt.
 
    Class support for input X:
       float: double, single
 
    See also erfc, erfcx, erfinv, erfcinv.


    Overloaded methods:
       sym/erf


    Reference page in Help browser
       doc erf


>> %计算多项式拟合的参数
>> [p,s]=polyfit(x,y,6);
>> [yp,delta]=polyval(p,x,s);
>> plot(x,y,'+',x,yp,'g-',x,yp+2*delta,'r:',x,yp-2*delta,'r:'),grid on
>> legend on
>> x=[1 3 5 7 9 2 4 6 8 10];
>> y=[0.0881 0.9290 2.4987 4.9292 7.8888 0.9766 3.4444 3.4444 6.4444 20.44444];
>> plot(x,y,'ro')
>> grid


----