原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
3.9 MATLAB语言在 科学运算中的应用 一、 矩阵分解 1 特征值分解 [V,D]=eig(A) 例： a = 9 8 6 8 [v,d]=eig(a) v = d = 0.7787 -0.7320 15.4462 0 0.6274 0.6813 0 1.5538 2 三角分解 [L,U]=lu(A) 将A做对角线分解，使得A=L*U,其中L为下三角矩阵，U为上三角矩阵。 注意：L实际上是一个置换矩阵P的逆矩阵与一个真正下三角矩阵L1（其对角线元素为1）的乘积。 [L1,U1,P]=lu(A) 二、 多项式处理 1、多项式的表示方法 *在MATLAB中，多项式使用降幂系数的行向量表示 表示为：p=[1 -12 0 25 116] ，使用函数roots可以求出多项式等于0的根，根用列向量表示。若已知多项式等于0的根，函数poly可以求出相应多项式。 2、多项式的运算 求多项式的微分多项式polyder 求多项式函数值polyval(p,n)：将值n代入多项式求解。 三、 曲线拟合与插值 根据样本进行拟合，找出满足样本分布的多项式。这在实验数据分析时非常有用。 命令格式：p=polyfit(x,y,n)，其中x和y为样本向量，n为所求多项式的阶数，p为求出的多项式。 多项式插值是指根据给定的有限个样本，产生另外的估计点以达到数据更为平滑的效果。在信号处理与图像处理上应用广泛。 一维内插interp1 y=interp1(x,y,xi,method) 四、 数据分析 绘制函数图形：fplot() 求极值：fmin,fmins 求零点：寻找一维函数的过零点fzero() 频谱分析：y= fft(x)；unwrap()；abs；angle画出幅频和相频曲线 数据分析函数：max,min,mean,sum 积分运算：trap2,quad,quad8 五、 常微分方程数值解 [t,x]=ode23(‘xfun’,[t0,tf],x0,选项) [t,x]=ode45(‘xfun’,[t0,tf],x0,选项) * 求矩阵A的特征向量V及特征值D，满足A*V=V*D。其中D的对角线元素为特征值，V的列为对应的特征向量。如果D=eig(A)则只返回特征值。 r=roots(p) r = 11.7473 2.7028 -1.2251 + 1.4672i -1.2251 - 1.4672i p=poly(r) p = 1 -12 -0 25 116 conv a=[1 2 3] ; b=[4 5 6] c=conv(a,b)=4 13 28 27 18 conv 指令可以嵌套使用，如conv (conv (a,b),c) deconv [q,r]=deconv(c,b) q=1 2 3 ％商多项式 r=0 0 0 0 0 ％余多项式 *
文档评论（0）