matlab教学内容.pptVIP

第二章 MATLAB 7的安装和用户界面 第4章 MATLAB 的其他函数库 教学目标:理解数据分析和傅里叶变换函数库， 多项式函数库，了解矩阵的分解与变换函数库，理解函数功能与数值函数库，了解字符串、符合数学与系统仿真函数库 教学重点:数据分析和傅里叶变换函数库，多项 式函数库 max(A)：各列最大值 min(A)：各列最小值 mean(A)：各列平均值 std(A)：各列标准差 median(A)：各列中间元素 sum(A)：各列元素和 sum(A(:))：全部元素和 trapz(A)：梯形法求积分，可看做求和 处理结果不是标量而是列向量的函数 例：说明trapz和cumtrapz的区别 例:求梯度 x=rand(1,1000);hist(x) y=randn(1,1000);hist(y,50) z=randi(10,1,10)，hist(z) filter(b,a,x)：根据输入信号x和线性系统特性b，a求输出信号 X=fft(x,N)：输入信号x的N点离散傅里叶变换 x=ifft(X)：傅里叶反变换 sound(u,s)：在音箱中产生u对应的声音，s为重放速度 (2) 当nm（有效方程数大于未知数数目），超定方程组，A\B仍然合法，得到解为最小二乘解。此时inv(A)不存在，x=inv(A’*A)*(A’*B) A\B= pinv(A)*B (除数A可以不是方阵) (3) nm（有效方程数小于未知数），x=A\B仍然有效，所得解不是全解，是令x中m-n个元素为0的一个特殊解。 最小二乘法（又称最小平方法）是一种数学优化技术。它通过最小化误差的平方和寻找数据的最佳函数匹配。 利用最小二乘法可以简便地求得未知的数据，并使得这些求得的数据与实际数据之间误差的平方和为最小。 [l,u]=lu(a); [q r]=qr(a); [u s v]=svd(a) 3 矩阵的特征值分析 [e, r]=eig(a): e表示特征向量，r表示特征根。 2 多项式求导、求根、求值 a为多项式系数向量 求导：e=polyder(a) 求根：ra=roots(a) 由根求多项式系数：a=poly(ra) 当a是向量时，poly把它看作根来组成多项式；当a是方阵时，poly把它组成方阵的特征多项式 多项式求值：将多项式a中的自变量x赋予值xv时，该多项式的值可用F=polyval(a,xv) 3 多项式拟合 4 多项式插值 5 线性微分方程的解 具体步骤 1 用Lapalace变换表示出常微分方程的解形式 Y(s)=B(s)/A(s)，其中B(s)和A(s)都是s的多项式，分母多项式的次数n通常高于分子多项式的次数m。 2 用residue函数求出Y(s)的极点数组和留数数组。 Y(s)=r(1)/(s-p(1))+ r(2)/(s-p(2))+…… 3 通过求Y(s)的拉普拉斯反变换得到常微分方程的解 y(t)=r(1)*exp(p(1)*t)+ r(2)*exp(p(2)*t)+… 实验任务 实验名称：用MATLAB进行函数库的仿真 实验任务： 1 MATLAB用于数据分析(基本、场论和随机数据分析) 2 MATLAB用于矩阵的特征分析及其分解。 3 MATLAB用于多项式函数库（四则运算、求导、求根、求值、拟合和插值） 4.4 函数功能和数值分析函数库 例：求多项式p(x)=2x4-6x3+3x2+0x+7的微分。 p=[2 -6 3 0 7]; q=polyder(p) p=[2 -6 3 0 7]; x=roots(p) 例：求多项式p(x)=2x4-6x3+3x2+0x+7的根 a=poly(x) 例：求上例中x=1时的值。 p=[2 -6 3 0 7]; polyval(p,1) 注意：polyval 中的xv可以是复数，也可以是矩阵 或数组，此时对输入变量作元素群运算 应用举例: 求多项式p= x4+3x3-5x+6的一阶导数，求x=1，3，5时多项式p的值。 p=[1 3 0 -5 6]; d=polyder(p); poly2str(d,x) x=1:2:5; y=polyval(p,x) ans = 4 x^3 + 9 x^2 - 5 y = 5 153 981 例2：将式(x-4)(x+5)(x2-6x+9)展开为多项式形 式，并求其对应的一元n次方程的根。 p=conv([1,-4],conv([1 5],[1,-6,9])) px=poly2str(p,x) x=roots(p) 例3：已知一元四次方程所对应的四个根为 -5.0000 4.0000 3.0000 3.0000 求这个方程所对应的