matlab p1xp2,MATLAB应用第6章第2版

MATLAB应用第6章第2版

第 6 章 MATLAB 数值计算6．1 数据处理最大值和最小值MATLAB 提供的求数据序列的最大值和最小值的函数分别为 max 和min，两个函数的调用格式和操作过程类似。向量的最大值和最小值求一个向量 X 的最大值的函数有两种调用格式，分别是：(1) y=max(X)：将向量 X 的最大值 max(X)存入 y，如果 X 中包含复数元素，则按模取最大值。例 1 求向量 x=[-43,72,9,16,23,47]的最大值x=[-43,72,9,16,23,47]y=max(x) %求向量 x 中的最大值(2) [y,k]=max(X)：将向量 X 的最大值 max(X)存入 y，最大值的序号存入k，如果 X 中包含复数元素，则按模取最大值。例 2 求向量 x=[-43,72,9,16,23,47]的最大值及其该元素的位置[y,k]=max(x) %求向量 x 中的最大值及其该元素的位置求向量 X 的最小值的函数是 min(X)，用法和 max(X)完全相同。例 3 求向量 x=[-43,72,9,16,23,47]的最小值及其该元素的位置z=min(x) %求向量 x 中的最小值[y,m]= min (x) %求向量 x 中的最小值及其该元素的位置矩阵的最大值和最小值求矩阵 A 的最大值的函数有 3 种调用格式，分别是：(1) max(A)：给出一个行向量，向量的第 i 个元素是矩阵 A 的第 i 列上的最大值 max(a1i,a2i,…ami)。例 4 求矩阵 A=[13, -56, 78; 25, 63, -235; 78, 25, 563; 1, 0, -1]各列的最大值及其整个矩阵的的最大值max(A),max(max(A)), %相当于求向量最大值(2) [Y,U]=max(A)：给出行向量 Y 和 U，Y 向量记录 A 的每列的最大值，U 向量记录每列最大值的行号。例 5 求矩阵 A=[13, -56, 78; 25, 63, -235; 78, 25, 563; 1, 0, -1]各列的最大值及其行号(3) max(A,[],dim)：dim 取 1 或 2。dim 取 1 时，该函数和 max(A)完全相同；dim 取 2 时，该函数返回一个列向量，其第 i 个元素是 A 矩阵的第 i 行上的最大值。例 6 用 max(A,[],dim)求矩阵 A=[13, -56, 78; 25, 63, -235; 78, 25, 563; 1, 0, -1]每行及每列的最大值求最小值的函数是 min，其用法和 max 完全相同。例 7 求矩阵 A=[13, -56, 78; 25, 63, -235; 78, 25, 563; 1, 0, -1]每行及每列的最小值及整个矩阵的最小值两个向量或矩阵对应元素的比较函数 max 和 min 还能对两个同型的向量或矩阵进行比较，调用格式为：(1) U=max(A,B)：A,B 是两个同型的向量或矩阵，结果 U 是与 A,B 同型的向量或矩阵，U 的每个元素等于 A,B 对应元素的较大者。例 8 求矩阵 A=[4, 5, 6; 1, 4, 8]及 B=[1, 7, 5; 4, 5, 7]所有同一位置上较大元素构成的新矩阵 P(2) U=max(A,n)：n 是一个标量，结果 U 是与 A 同型的向量或矩阵，U 的每个元素等于 A 对应元素和 n 中的较大者。例 8 求矩阵 A=[4, 5, 6; 1, 4, 8]的所有元素与 f=4.5 比较后较大者构成的新矩阵P1min 函数的用法和 max 完全相同。元素排序sort(X)，它给出一个对向量 X 中的元素按升序排列的新向量。例 9 对向量 x=[-43,72,9,16,23,47]进行升序排列sort 函数也可以对矩阵 A 的各列或各行重新排序，其调用格式为：[Y,k]=sort(A,dim)其中 dim 指明对 A 的列还是行进行排序。若 dim=1，则按列排，可省略；若 dim=2，则按行排。Y 是排序后的矩阵，而 k 记录 Y 中的元素在 A 中位置。例 10 对矩阵 A=[1, -8, 5; 4, 12, 6; 13, 7, -13]进行各种排序A=[1, -8, 5; 4, 12, 6; 13, 7, -13],sort(A) %对 A 的每列按升序排列sort(A,2) %对 A 的每行按升序排列-sort(-A,2) %对 A 的每行按降序排列[Y,k]=sort(A) %对 A 的按列升序排列，并将每个元素所在的行号送到矩阵 k标准方差设有 N 个数据组成数据序列 x1,x2,…xN，这些数据的标准方差为或 其中，iixS121 iixS12Nix1在 MATLAB 中，提供了计算数据序列的标准方差的函数 std。对于向量 X，std(X)给出一个标准方差。例 11 随机地取 8 个活塞环，测得它们的直径为(以 mm 计)74.001, 74.005, 74,003,74.001, 74.000, 73.998, 74.002，求它的标准方差对于矩阵 A，std(A)给出一个行向量，它的各个元素便是矩阵 A 各列或各行的标准方差。std 函数的一般调用格式为：Y=std(A,flag,dim)其中 dim 取 1 或 2。当 dim=1 时，求各列元素的标准方差；当 dim=2 时，则求各行元素的标准方差。flag 取 0 或 1，当 flag=0 时，按公式 1 所列公式计算标准方差，当 flag=1 时，按公式 2 所列公式计算标准方差。缺省flag=0，dim=1。例 12 对矩阵 A=[4, 5, 6; 1, 4, 8]，从不同维方向求其标准方差A=[4, 5, 6; 1, 4, 8], Y1=std(A,0,1); %按按公式 1 各列元素的标准方差Y2=std(A,1,1); %按按公式 2 各列元素的标准方差Y3=std(A,0,2); %按按公式 1 各行元素的标准方差Y4=std(A,1,2); %按按公式 2 各行元素的标准方差多项式计算n 次多项式表达为： ，是 n+1 项之和nax xa p(x)1-1-n0在 MATLAB 中，n 次多项式可以用长度为 n+1 的行向量表示[a0, a1,……an-1,an]多项式的加减运算设有两个多项式 和nax xa p1(x)1-1-n0。它们的加减运算实际上就是它们的对应mbb xp2() 1--n0系数的加减运算。当它们的次数相同时，可以直接对多项式的系数向量进行加减运算。当它们的次数不同时，应该把次数低的多项式无高次项部分用 0 系数表示。例 2 计算 16352xxa=[1, -2, 5, 3]; b=[0, 0, 6, -1]; c=a+b例 3 设 ， ，求 f(x)+g(x)6572345f 352xxgf=[3, -5, 2, -7, 5, 6];