Matlab的其他函数库培训讲义.ppt

第4章 Matlab的其他函数库 第1节 数据分析数据库4.1.1 基本的数据分析 基本数据处理功能是按列进行，行向则表示数据的不同的样本。P63表 例 data=[154 49 83 67;158 99 81 75;155 100 68 86;145 63 75 96;145 63 75 96;... 141 55 65 75;155 56 64 85;147 89 87 77;147 96 54 100;145 60 76 67] std(标准差)指各列中N个元素与该列平均值的平方和开方 trapz（求积分）梯形法求积分近似求元素和 4.1.2 用于场论的数据分析 cradient 求二维场和三维场的近似梯度 del2求二维场和三维场的拉普拉斯算子 cross 两个相量的矢量积（叉乘） cot 两个相量的数量积（点乘） 4.1.3 用于随机数据分析的函数 rand（m，n） 产生在0和1之间的m行n列的随机数矩阵 randn（m,n） 产生正态分布的m行n列随机数矩阵 x=rand(1,1000)；hist(x) y=randn(1,1000)；hist(y,50) 4.1.4 用于相关分析和傅立叶分析的函数 用于信号的时域分析和频域处理。 corrcoef 给出两个同长信号的相关系数 cov（x,y）给出协方差矩阵 conv（x,y）给出x,y的卷积 Filter（b,a,y）根据输入信号x和线性系统求输出信号的函数 X＝fft（x，N）求出时域信号x的离散傅立叶变换X t=0:.001:3;u=sin(300*t)+2*cos(200*t);U=fft(u); plot(abs(U)) 第2节 矩阵的分解与变换（matfun函数库） 4.2.1 线性方程组的系数矩阵 det（a）求方阵a的行列式 rank（a）秩 trace（a）矩阵主对角线上元素的和（迹） inv（a）逆阵 pinv(a) 4.2.2 矩阵的分解p68 正交分解（qr分解）[q,r]=qr(b) 奇异值分解（svd分解）[u,s,v]=svd(b) 常用的矩阵函数 第3节 多项式函数库（polyfun） 一元高次代数多项式 4.3.1 多项式简介 多项式是形如下式的式子： P(x)=a0xn+a1xn-1+…+an-1x+an 在MATLAB中，多项式用行向量表示： P＝[a0 a1 a2…an-1 an] 多项式行向量的构造 直接键入，无该次幂则取零 用poly(A)命令来生成矩阵A的特征多项式。 ＊如果A是形如[a0 a1…an-1 an]的向量，由命令poly(A)可生成(x-a0)(x-a1)…(x-an-1)(x-an)所对应的多项式。 例：已知向量A=[1 –34 –80]，用此向量构造一多项式并显示结果。 A=[1 -34 -80]; PA=poly(A) PA = 1 113 2606 -2720 PAX=poly2sym(PA,X) PAX = X^3+113*X^2+2606*X-2720 4.3.2 多项式运算 多项式的加减运算 运算符　＋　－ 进行加减运算的多项式应该具有相同的阶次，如果阶次不同，低阶的多项式必须用零添补至高阶多项式的阶次。 例：求两个多项式a(x)=5x4+4x3+3x2+2x+1和b(x)=3x2+0x+1的和。 a=[5 4 3 2 1];b=[3 0 1]; c=a+[0 0 b] c = 5 4 6 2 2 多项式乘法 多项式乘法采用conv函数。 据上例求积： d=conv(a,b) d = 15 12 14 10 6 2 1 多项式除法 由deconv函数完成，结果包括商和余数两部分。 [div,rest]=deconv(d,a) div = 3 0 1 rest = 0 0 0 0 0 0 0 微分 提供了微分函数　polyder 例：求多项式p(x)=2x4-6x3+3x2+0x+7的微分。 p=[2 -6 3 0 7]; q=polyder(p) q = 8 -18 6 0 pd=poly2sym(q,x) pd = 8*x^3-18*x^2+6*x 求根 多项式求根可用roots函数。 例：求多项式p(x)=2x4-6x3+3x2+0x+7的根。 p=[2 -6 3 0 7]; x=roots(p) x = 1.9322 + 0.4714i 1.9322 - 0.4714i