第七讲 数值数组及向量化运算.ppt

第七讲 数值数组及向量化运算 Xidian University * 第一页，共二十四页，2022年，8月28日 Xidian University * 数组创建 逐个构造 等差数列型 :, linspace 等比数列型 logspace, semilogx/y 工作空间浏览器 load/save命令 第二页，共二十四页，2022年，8月28日 Xidian University * 全1矩阵、全0矩阵、空矩阵、随机矩阵(均匀分布、正态分布) 对角阵、单位阵、魔方数组 第三页，共二十四页，2022年，8月28日 Xidian University * 数组的访问 A, A(:), A(:, :) A(r, c), A(r, :), A(:, c) A(I) A(L) %sizeof(A) == sizeof(L) A(2:end, 3:end) 第四页，共二十四页，2022年，8月28日 Xidian University * reshape, repmat flipud, fliplr, rot90 第五页，共二十四页，2022年，8月28日 Xidian University * 数组运算 A’, A.’, A+B, A-B A*B, A.*B, A/B, A./B, A\B, A.\B A.^a, a.^A 第六页，共二十四页，2022年，8月28日 Xidian University * 非数 NaN, inf, isnan, isinf 第七页，共二十四页，2022年，8月28日 Xidian University * 特殊操作 find ind2sub %单下标换全下标 第八页，共二十四页，2022年，8月28日 Xidian University * 空数组 a=[] b=rand(0) c=rand(1, 0) size(a) size(b) size(c) length(a) length(b) length(c) isempty(a) isempty(c) 第九页，共二十四页，2022年，8月28日 Xidian University * 关系操作 , =, , =, ==, ~= 第十页，共二十四页，2022年，8月28日 Xidian University * 逻辑操作 , |, ~, xor , || all, any false, true, logical isempty, isfinite, isprime, isreal, isletter, isspace, isinf, isnan isa, ishandle, ischar, islogical, isglobal, isnumeric 第十一页，共二十四页，2022年，8月28日 Xidian University * 数值计算 差分 diff 梯度 gradient 求和 sum: sum(A, 1), sum(A, 2) 累计求和 cumsum 近似数值积分 trapz 累计近似数值积分 cumtrapz 第十二页，共二十四页，2022年，8月28日 Xidian University * 数值积分 trapz(x, y, tol) quad(fun, a, b, tol) dblquad(fun, xmin, ymin, xmax, ymax, tol) 第十三页，共二十四页，2022年，8月28日 Xidian University * 求极值 一元函数极值 fminbnd(fun, x1, x2, options)fminbnd(@(x)x^2+2*x-3, -1, 1) 单纯形法求多元函数的极值fminsearch(fun, x0, options) 第十四页，共二十四页，2022年，8月28日 Xidian University * 常微分方程 数值略 第十五页，共二十四页，2022年，8月28日 Xidian University * 矩阵和代数方程 高斯列主元消去法 [R, ci] = rref(A) 特征值、特征向量、行列式 X=null(A); A矩阵零空间的全部正交基 AX=0 Z=orth(A); A矩阵值空间的全部正交基 span(A) = span(Z) 第十六页，共二十四页，2022年，8月28日 Xidian University * 代数方程的求解 fsolve(fun, x0) fzero(fun, x0) 对分法 第十七页，共二十四页，2022年，8月28日