资源描述：
《Matlab数据处理》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、MATLAB高级应用第一讲数据分析与数据建模授课人：鲍文退出>在此幻灯片插入公司的徽标从“插入”菜单选择图片找到徽标文件单击“确定”重新设置徽标大小单击徽标内任意位置。徽标外部出现的方框是“调整控点”使用这些重新设置对象大小如果在使用尺寸调整控点前按下shift键，则对象改变大小但维持原比例。2021/8/3哈尔滨工业大学动力工程控制与仿真研究所目录数据分析方法数据建模数据拟和工具箱退出主菜单><2021/8/3哈尔滨工业大学动力工程控制与仿真研究所一、绘图最简单的数据分析是绘图。二维绘图是观察数据的最简单的模式。频域数据分析可以从另一个侧面分析数据。2021/8/31.1线图plot操作p
2、lot(x)x为实数/复数自动x/y轴plot(x,y)向量/矩阵plot(x1,y1,x2,y2,…)控制曲线形式、标记类型和颜色颜色控制：y(黄)、m(洋红)、c(青)、r(红)、g(绿)、b(蓝)、w(白)、k(黑)标记控制：.(黑点)、o(圆圈)、x(叉)、+(十字)、*(星)线形：-(实线)、:(虚线)、-.(点划)、--(双划)2021/8/31.2其它二维图形plot迪卡尔坐标图loglog双对数坐标semilogxX轴对数坐标semilogyY轴对数坐标polar极坐标polar(theta,rho)，theta相角，rho半径plotyy两侧显示y坐标的绘图bar(x)棒图
3、hist分布图pie饼图2021/8/31.3三维绘图1.线图plot32.三维面图mesh网格图view改变视角2021/8/31.4自动绘图打开变量名表使用右键绘图修改图形的特性：增加网格线修改坐标增加注释使用不同颜色、线型等等2021/8/3二、时域分析函数名含义max最大值min最小值mean均值std标准方差median中值2021/8/3分析函数函数名含义sum元素的总和prod元素的乘积cumrod元素的累积cumsum元素的累加和diff差分函数:少了一个元素2021/8/3三、频域分析FFT即快速傅立叶变换，是数据分析的基本方法，是x由基2的快速变换算法来计算。如x长度不是
4、精确的2次幂则后面使用0填充，ifft(x)是向量x的离散傅立叶变换的逆变换。在频率轴上绘制FFT曲线，要明确FFT结果与实际频率点的关系。设n个数据点，采样频率为fs，则Nyquist频率或n=N/2+1点与实际频率的关系：f=(num-1)*fs/n2021/8/3FFT需要注意的是fft结果为复数矩阵，为了得到幅频特性，可使用abs函数，使用atan2得到相角，由于有的系统的相角可能大于1800，而相角函数值域在-1800~1800之间，需要使用unwrap函数展开折叠的相角，从而得到相频特性。2021/8/3四、曲线拟和4.1引言在分析试验数据中，常常要面临将试验数据作解析描述的任务
5、，这个问题有曲线拟合和插值两种方法。在曲线拟合中，假定已知曲线的规律，作曲线的最佳逼近，但不需要经过所有的数据点；在插值中，认为数据是准确的，求取其中描述点之间的数据。2021/8/34.2曲线拟合1.多项式的最小二乘曲线拟合使用polyfit，它需要曲线的x、y值，以及曲线的阶数。曲线的阶数：如果曲线的阶数选择的过小，拟合效果不好；如果曲线的阶数过高，虽然数据点上看到效果好，数据点之间会出现有数据振荡的问题，阶数不宜过高，小于5阶。灵活使用拟合2021/8/32.直接最小二乘数据规律并不是多项式形式，直接最小二乘来拟合。最小二乘函数为k=nnls(fx,y)nnls已经被lsqnonneg
6、代替，用法相同计算结果将使得
7、fx*k-y
8、2范数下最小在计算中，fx可以为x的函数。例子：拟合流量特性2021/8/3灵活使用最小二乘拟和多参数的方程，比如三个参数的不高于二阶的方程，y=f(x1,x2,x3)。最小二乘的列可以多一些，穷举所有的二阶以下的变量，二次：x1^2,x2^2,x3^2,x1*x2,x1*x3,x2*x3一次：x1,x2,x3零次：1得到最小二乘的系数。2021/8/3五、三次样条1使用的原因高阶多项式插值出现病态问题，三次样条使用分段多项式，各点上的三次导数相等。它光滑、导数连续。2插值yi=spline(x,y,xi);pp=spline(x,y);分段多项式
9、形式例程：ex442021/8/3三次样条pp形式可以和三次多项式形式转化：[break,coef,np,nc]=unmkpp(pp)断点、三次多项式、多项式数量、系数数量pp=mkpp(break,coef);由于转化为了多项式形式，可以方便的进行积分和微分运算。2021/8/3六、滤波和平滑1.插值和拟合的问题：噪声2.滤波:滞后,filtery=filter(b,a,x)a,b:滤波器的分子分母,x输入
