MATLAB及其应用 第三讲 数据处理 授课人:鲍文 在此幻灯片插入公司的徽标 从“插入”菜单 选择图片 找到徽标文件 单击“确定”

Slides:



Advertisements
Similar presentations
简单迭代法的概念与结论 简单迭代法又称逐次迭代法,基本思想是构造不动点 方程,以求得近似根。即由方程 f(x)=0 变换为 x=  (x), 然后建立迭代格式, 返回下一页 则称迭代格式 收敛, 否则称为发散 上一页.
Advertisements

1. 卸下标签 身心松静 关注健康! 2. 坦诚开放 互信互赖 社会支持! 3. 排除干扰 倾心体悟 创造协作! 4. 连接自己 享受成长 和谐社会! 恳请与提醒.
2.5 微分及其应用. 三、可微的条件 一、问题的提出 二、微分的定义 六、微分的形式不变性 四、微分的几何意义 五、微分的求法 八、小结 七、微分在近似计算中的应用.
1/14 练习题 Ex1. 计算球体 V 允许其相对误差限为 1%, 问测量球 半径 R 的相对误差限最大为多少 ? 试分析高度误差对面积计算的影响。 Ex2. 将地球模型取为半径为 R (km) 的球体,赤道上 方高度为 d (km) 的地球同步卫星发射的信号对地球 的覆盖面积计算公式为 Ex3 在计算机上对调和级数逐项求和.
《数值计算》课件 第五章 数值积分与数值微分 第三章 数值积分与数值微分 3.1 引例及 Newton-Cotes 公式 3.2 复合求积公式 3.3 龙贝格求积方法 3.4 数值微分 3.5 引例的 MATLAB 求解.
人的性别遗传 合肥市第四十九中学 丁 艳. 男女成对染色体排序图 1 、男性和女性各 23 对染色体有何异同 ? 哪 一对被称为性染色体 ? 2 、这两幅图中,哪幅 图显示的是男性的染色 体?哪幅图显示的是女 性染色体? 3 、图中哪条染色体是 Y 染色体?它与 X 染色体 在形态上的主要区别是.
《公路纵断面设计》 —— 纵断面设计的要求 道桥系 二○○七年五月. 纵断面设计的一般要求 1 .纵坡设计必须满足《公路工程技术标准》中的各项规定。 2 .为保证汽车能以一定的车速安全舒顺地行驶,纵坡应具有 — 定 的平顺性,起伏不宜过大及过于频繁。尽量避免采用极限纵坡 值.缓和坡段应自然地配合地形设置,在连续采用极限长度的.
第 5 章 基因突变及其他变异 第 3 节 人类遗传病 【思考】 感冒是不是遗传病? 先天性疾病、地方性疾病和遗传 病有什么关系?
安徽7班全真模拟 主讲: 杨洁 时间:6月12日晚.
XX啤酒营销及广告策略.
普陀区税务局 营业税改征增值税试点 最新政策 货物和劳务税科 2013年7月.
第四章:长期股权投资 长期股权投资效果 1、控制:50%以上 有权决定对方财务和经营.
专利技术交底书的撰写方法 ——公司知识产权讲座
小学科学中的化学 武威十九中 刘玉香.
1、一般地说,在生物的体细胞中, 和 都是成对存在的。
辨性别 A B. 辨性别 A B 第三节人类染色体与性别决定 昌邑市龙池初中 杨伟红 学习目标 1.理解人的染色体组成和传递规律。 2.解释人类性别决定的原理。 3.通过探究活动,解读数据了解生男生女的比例。
神州五号、六号的发射和回收都取得了成功 ,圆了几代中国人的航天梦,让全中国人为之骄傲和自豪 神州五号、六号的发射和回收都取得了成功 ,圆了几代中国人的航天梦,让全中国人为之骄傲和自豪!但是你们知道我们的科学家是怎样迅速地找到返回舱着陆的位置的吗? 这全依赖于GPS——卫星全球定位系统”。大家一定觉得很神奇吧!学习了今天的内容,你就会明白其中的奥妙。
这是一个数字的 乐园 这里埋藏着丰富的 宝藏 请跟我一起走进数学的 殿堂.
—— matlab 具有出色的数值计算能力,占据世界上数值计算软件的主导地位
二次函數 高士欽 林國源.
西南科技大学网络教育系列课程 5. 优 化 设 计 5.2 优化方法的数学基础.
欢迎大家来到生命科学课堂.
初中语文总复习 说明文 阅读专题 西安市第六十七中学 潘敏.
班级安全文化建设的思考与实践 夯实安全基础 规范安全行为 培养安全习惯 训练安全能力 尤 学 文 管 理 学 博 士
行政作用法 行政命令.
2016中重卡网络规划 中重卡营销部 2016年6月.
清仓处理 跳楼价 满200返160 5折酬宾.
四种命题 班级:C274 指导教师:钟志勤 任课教师:颜小娟.
1.1.2 四 种 命 题.
高一数学 充分条件与必要条件 教育科学学院03级教育技术2班 刘文平.
色 弱 與 色 盲.
第五章 定积分及其应用.
第4章 数值积分与数值微分 4.1 数值积分概论 4.2 牛顿-柯特斯公式 4.3 复合求积公式 4.4 龙贝格求积公式
微积分基本公式 在上一节我们已经看到,直接用定义计算定积分是十分繁难的,因此我们期望寻求一种计算定积分的简便而又一般的方法。我们将会发现定积分与不定积分之间有着十分密切的联系,从而可以利用不定积分来计算定积分。
宠物之家 我的宠物性别? 雌(♀) or 雄(♂) 第一阶段:我的宠物我做主 第二阶段:宠物“相亲记” 第三阶段:家族诞生
课标教材下教研工作的 实践与思考 山东临沂市教育科学研究中心 郭允远.
北师大版七年级数学 5.5 应用一元一次方程 ——“希望工程”义演 枣庄市第三十四中学 曹馨.
第八章二元一次方程组 8.3实际问题与二元一次方程组.
第八章二元一次方程组 8.3实际问题与二元一次方程组 (第3课时).
2010年桂林理工大学数学建模竞赛暑期培训班 MATLAB编程入门培训课程
運輸與空間的交互作用 運輸發展的階段 一、分散的港口 二、侵入路線 三、發展支線 四、初步相互連結 五、完全相互連結 六、高度優越的幹線
海洋存亡 匹夫有责 ——让我们都来做环保小卫士 XX小学三(3)班.
第七讲 二维连续分布独立性与二维函数分布 本次课讲授:第二章的 ; 下次课讲第三章的 。
MATLAB数学实验 第四章 函数和方程.
Introduction to MATLAB
授課教授:張寶基 助教:梁凱雯 郭千豪 音視訊處理實驗室 2014 / 9 / 30
網路遊戲版 幸福農場168號.
Matlab及其应用 讲座之五 动态系统仿真——Simulink 主讲人:鲍文 在此幻灯片插入公司的徽标 从“插入”菜单 选择图片
第四章 内压薄壁圆筒与封头的强度设计 教学重点: 内压薄壁圆筒的厚度计算 教学难点: 厚度的概念和设计参数的确定.
第二章 函 数 插 值 — Matlab插值函数.
插值与拟合 一、插值的基本原理 二、拟合的基本原理 三、插值与拟合的关系 四、插值的MATLAB实现 五、拟合的Matlab实现.
导数的应用 ——函数的单调性与极值.
数学建模 江西财经大学 数学与管理决策系 制作:华长生 华长生制作.
第二节 极限 一、数列极限 定义:.
医学信号处理的原理和方法 曹银祥 Dept. of Physiology & Pathophysiology
第七章  事业单位支出的核算      §第一节  支出概述     §第二节  拨出款项     §第三节  各项支出     §第四节  成本费用.
Homework 1(1/2) 本頁表格為派密(Peavy)在07年球季各場次ERA的表現,(1)請依此數據完成下頁表格之統計值並說明之;(2)並與其他三名投手之統計值比較之。(請詳述計算過程) 場次 各場次ERA
第8章 MATLAB数值积分与微分 8.1 数值积分 8.2 数值微分
Matlab及其应用 讲座之五 动态系统仿真——Simulink 主讲人:鲍文 在此幻灯片插入公司的徽标 从“插入”菜单 选择图片
西式點心 派的種類 單皮派 雙皮派 油炸派 派的製作 派的烤焙.
第六章 影像幾何 6.1 數據內插法 假設有4 個數值要放大成8 個數值,該怎麼做? 解出線性係數a、b如下:
業務員 傷害險通報作業 新光人壽內網-產險傷害險通報P2~P4 【個人】傷害險通報作業P5~P10 【團體】傷害險通報作業P11~P16
MATLAB 程式設計進階篇 多項式的處理與分析
导数的几何意义及其应用 滨海中学  张乐.
第4章 MATLAB数值运算 编者.
数学题解答 第二章 一元一次方程 2.1从算式到方程 (第1课时) 数学题解答
第7章 MATLAB工程计算.
教学大纲(甲型,54学时 ) 教学大纲(乙型, 36学时 )
实验二 定积分的近似计算.
第二节 偏 导 数 一、 偏导数概念及其计算 二 、高阶偏导数.
Presentation transcript:

MATLAB及其应用 第三讲 数据处理 授课人:鲍文 在此幻灯片插入公司的徽标 从“插入”菜单 选择图片 找到徽标文件 单击“确定” 2019/2/17 在此幻灯片插入公司的徽标 从“插入”菜单 选择图片 找到徽标文件 单击“确定” 重新设置徽标大小 单击徽标内任意位置。徽标外部出现的方框是“调整控点” 使用这些重新设置对象大小 如果在使用尺寸调整控点前按下 shift 键,则对象改变大小但维持原比例。 MATLAB及其应用 第三讲 数据处理 授课人:鲍文 2019/2/17 哈尔滨工业大学动力工程控制与仿真研究所 > 退出

目录 1 矩阵分析 2 数据分析函数 3 多项式处理 4 曲线拟和与插值 5 数据分析 6 微分方程数值解 < > 主菜单 退出 2019/2/17 哈尔滨工业大学动力工程控制与仿真研究所 < > 主菜单 退出

1 矩阵分析 一、特征值分解 对于方阵a特征值问题:ax=rx,求取a阵的特征值和特征向量使用下面的方法: [v,d]=eig(a) 使用 [v,d]=eig(a,’nobalance’) “平衡” 的作用减少计算误差,不平衡用于A阵大小悬殊的时候。 广义特征值问题:ax=rbx,求解的方式为: [v,d]=eig(a,b) 2019/2/17

二、三角分解 三角分解把矩阵分解为上三角矩阵和下三角矩阵,又称为LU分解或者。计算中使用高斯变量消去法。这一分解使用 [l,u]=lu(a)实现。 2019/2/17

三、奇异值分解 [u,s,v]=svd(a) 实现奇异值分解。 分解得到的三个因数有如下关系 a=u*s*v 其中u矩阵和v矩阵是正交矩阵,s矩阵是对角矩阵,它的对角元素是a矩阵的奇异值。 奇异值分解的稳定性很好。 2019/2/17

2 数据分析函数 函数名 含义 max 最大值 min 最小值 mean 均值 std 标准方差 median 中值 2019/2/17

分析函数 函数名 含义 sum 元素的总和 prod 元素的乘积 cumrod 元素的累积 cumsum 元素的累加和 函数名 含义 sum 元素的总和 prod 元素的乘积 cumrod 元素的累积 cumsum 元素的累加和 diff 差分函数:少了一个元素 2019/2/17

例题 求出y=x*sin(x) 在0<x<100的每个峰值 思路: 1、数学上峰值就是导数为零的点 2、导数在matlab中可以使用差分代替 3、差分后怎么求过零点呢? 2019/2/17

3 多项式处理 一、多项式表示 多项式在MATLAB中使用降幂系数的行向量表示。表示中需要包含零系数的项。poly2str:control toolbox中的函数 使用函数roots可找出多项式等于零的根。 规定:多项式用行向量,根用列向量。 给出多项式的根,使用poly函数也可以构造出相应的多项式。 2019/2/17

二、多项式运算 函数conv进行乘法运算,deconv进行除法运算。MATLAB没有提供特别的多项式加减法运算。 多项式除法并不一定能够除尽,很多时候需要有余数多项式。 多项式微分使用polyder(p)函数,估计值使用polyval(p,at)函数。 2019/2/17

4 曲线拟和与插值 在分析试验数据中,常常要面临将试验数据作解析描述的任务,这个问题有曲线拟合和插值两种方法。 在曲线拟合中,假定已知曲线的规律,作曲线的最佳逼近,但不需要经过所有的数据点;在插值中,认为数据是准确的,求取其中描述点之间的数据。 2019/2/17

一、曲线拟合 1、多项式的最小二乘曲线拟合 使用polyfit,它需要曲线的x、y值,以及曲线的阶数。 曲线的阶数:如果曲线的阶数选择的过小,拟合效果不好;如果曲线的阶数过高,虽然数据点上看到效果好,数据点之间会出现有数据振荡的问题,阶数不宜过高,小于5阶。 灵活使用拟合 2019/2/17

2、直接最小二乘 最小二乘函数为k=nnls(fx,y) 数据规律并不是多项式形式,直接最小二乘来拟合。 matlab 数据规律并不是多项式形式,直接最小二乘来拟合。 最小二乘函数为k=nnls(fx,y) 计算结果将使得|fx*k-y|2范数下最小 在计算中,fx可以为x的函数。 例子:拟合 2019/2/17

二、插值函数 1、曲线插值函数interp1 方法 t=interp1(x,y,x0,’method’) matlab 1、曲线插值函数interp1 方法 t=interp1(x,y,x0,’method’) x、y:原始数据点,x0为进行插值的数组,method为插值算法:线性插值('linear'),三次样条插值('spline'),三次多项式插值(‘cubic’). 如果x0出界,则对应值为NaN 例程:ex42.m 2019/2/17

2、曲面插值 插值函数: interp2,基本形式: zi=interp2(x,y,z,xi,yi,method) method包括 linear:线性 cubic:三次多项式 nearest:粗略估计数据 例程:ex43 2019/2/17

三、三次样条 1、使用的原因 高阶多项式插值出现病态问题,三次样条使用分段多项式,各点上的三次导数相等。它光滑、导数连续。 2、插值 yi=spline(x,y,xi); pp=spline(x,y); 分段多项式形式 例程:ex44 2019/2/17

三次样条 pp形式可以和三次多项式形式转化: [break,coef,np,nc]=unmkpp(pp) 断点、三次多项式、多项式数量、系数数量 pp=mkpp(break,coef); 由于转化为了多项式形式,可以方便的进行积分和微分运算。 2019/2/17

四、滤波和平滑 1、插值和拟合的问题:噪声 2、滤波: 滞后,filter y=filter(b,a,x) a,b:滤波器的分子分母,x输入 a(1)*y(n) = b(1)*x(n) + b(2)*x(n-1) + ... + b(nb+1)*x(n-nb) - a(2)*y(n-1) - ... - a(na+1)*y(n-na) 例程:ex46 2019/2/17

3、平滑 yi=csaps(x,y,P,xi) yi=csaps(x,y,P) 其中P为平滑因子0~1 0: 最小二乘 1:平滑近似 0: 最小二乘 1:平滑近似 ex46 ex45 2019/2/17

5 数据分析 1、极小化 MATLAB提供了fmin和fmins两个函数来求极值,它们分别寻找一维和n维函数的极值。它使用的单纯性法搜索。函数计算量大,或搜索区内有多极值,搜索的过程较长,也可能找不到极值。如找不到极值,将停止运行并提供解释。 寻找极大值点,重定义函数为-f(x)即可。 2019/2/17

2、求零点 函数fzero可以寻找一维函数的过零点。 应用:使用bode图判断控制系统稳定性,要看幅频特性过零点和相频特性过1800点。 fzero函数也可以寻找函数值等于常值点,只要重新定于函数为f(x)-c即可 2019/2/17

3、积分 有限区域内积分函数:trapz、quad和quad8。 函数trapz通过计算梯形面积的和近似函数的积分,函数的分割是人为地。 quad使用Simpson递归方法,quad8使用Newton-costes递归方法进行数值积分。为了获得更精确的结果,它们在所需的区间都计算被积函数。quad8比quad更精确。 2019/2/17

4、微分 微分描述了函数在一点处的斜率,是函数的微观性质,它对函数的微小变化十分敏感,函数的很小的变化,容易产生相邻点斜率的巨大变化。 尽量避免使用数值微分,尤其是试验数据的微分。如果迫切需要,最好先将试验数据进行最小二乘拟合伙这三次样条拟合,然后对拟合函数进行微分。 2019/2/17

5、FFT变换 FFT即快速傅立叶变换,是数据分析的基本方法,是x由基2的快速变换算法来计算。如x长度不是精确的2次幂则后面使用0填充,ifft(x)是向量x的离散傅立叶变换的逆变换。 在频率轴上绘制FFT曲线,要明确FFT结果与实际频率点的关系。设n个数据点,采样频率为fs,则Nyquist频率或n=N/2+1点与实际频率的关系:f=(num-1)*fs/n 2019/2/17

FFT 需要注意的是fft结果为复数矩阵,为了得到幅频特性,可使用abs函数,使用atan2得到相角,由于有的系统的相角可能大于1800,而相角函数值域在-1800~1800之间,需要使用unwrap函数展开折叠的相角,从而得到相频特性。 2019/2/17

6 微分方程数值解 常微分方程数值解用逐步积分方法实现,Runge-Kutta法是应用最多的微分方程数值解的方法。两种Runge-Kutta法函数: [t,x]=ode23(‘xfun’,t0,tf,x0,tol,trace) [t,x]=ode45(‘xfun’, ,t0,tf,x0,tol,trace) 这两种方法格式相同。其中xfun为定义的常微分方程函数名,该函数必须以为输出,以t、x为输入。 2019/2/17

微分方程 输入变量t0、tf为积分的启始和中止时间,单位是秒。x0为初始的状态向量。tol控制结果的精度,可以缺省。一般来说,ode45比ode23运算速度快一些。 Var der Pol微分方程 重新定义变量,令 x1=x x2=dx/dt 则 dx1/dt=x2 dx2/dt=u(1-x12)x2-x1 2019/2/17

应用举例 一、特性拟合 title(‘string’) ^上标 _下标 二、模型辨识的阶数确定 三、数值积分:已知加速度求速度, 2019/2/17