基于MATLAB的数值分析(1)(精品·公开课件).ppt

基于MATLAB的数值分析 第一章 Matlab入门 1.1 引言 1.2 MATLAB的操作桌面和启动 1.3 指令窗简介 1.4 历史指令窗和实录指令diary 1.5 用M文件开发程序 1.6 当前目录浏览器和文件管理 1.7 工作空间浏览器和数组编辑器 1.8交互界面分类目录窗和帮助系统 1.9 数组的创建和运算 一、一维数组变量 2.冒号创建向量 x1=3.4:6.7, x2=3.4:2:6.7 x3=2.6:-0.8:0 运算结果 x1 = 3.4000 4.4000 5.4000 6.4000 x2 = 3.4000 5.4000 x3 = 2.6000 1.8000 1.0000 0.2000 3.生成线性等分向量 指令x=linspace(a,b,n) 在[a,b]区间产生 n 个等分点(包括端点) x=linspace(0,1,5) 结果 x = 0 0.2500 0.5000 0.7500 1.0000 5.向量的运算 设x=[x1 x2 x3]; y=[y1 y2 y3];为两个三维向量，a,b为标量。 向量的数乘：a*x=[a*x1 a*x2 a*x3] 向量的平移： x+b=[x1+b x2+b x3+b] 向量和： x+y=[x1+y1 x2+y2 x3+y3] 向量差： x-y=[x1-y1 x2-y2 x3-y3] x/a=[x1/a x2/a x3/a] 向量的数组运算(四则运算） x.*y=[x1*y1 x2*y2 x3*y3] x./y=[x1/y1 x2/y2 x3/y3] x.\y=[y1/x1 y2/x2 y3/x3] x.^5=[x1^5 x2^5 x3^5] 2.^x=[2^x1 2^x2 2^x3] x.^y=[x1^y1 x2^y2 x3^y3] 向量的数组运算（函数计算） Matlab有许多内部函数，可直接作用于向量产 生一个同维的函数向量。 x=linspace(0,4*pi,100);（产生100维向量x） y=sin(x); (y也自动为100维向量) y1=sin(x).^2; y2=exp(-x).*sin(x); 观察结果 二、二维数组变量 取出A的1、3行和1、3列的交叉处元素构成新矩阵A1 程序 A=[1 0 1 1 2;0 1 -1 2 3; 3 0 5 1 0;2 3 1 2 1], vr=[1, 3];vc=[1, 3]; A1=A(vr, vc) 观察结果 将A分为四块，并把它们赋值到矩阵B中，观察运行后的结果。 程序 A11=A(1:2,1:2),A12=A(1:2,3:5), A21=A(3:4,1:2),A22=A(3:4,3:5) B=[A11 A12;A21 A22] 结果 4.矩阵的修改和提取 修改矩阵A，将它的第1行变为0。 程序： A=[1 0 1 1 2;0 1 -1 2 3; 3 0 5 1 0;2 3 1 2 1] A(1,:)=[0 0 0 0 0]; A 删除上面矩阵A的第1、3行。 程序： A([1,3],:)=[ ] 结果 5.生成特殊矩阵 全1阵 ones(n), ones(m,n), ones(size(A)) 全零阵： zeros(n) ，zeros(m,n), zeros(size(A)) 常常用于对某个矩阵或向量赋0初值 单位阵： eye(n)，eye(m,n) 随机阵： rand(m,n)， rand(n)=rand(n,n)用于随机模拟，常和rand(seed,k)配合使用。 将rand指令运行多次，观察结果。 程序： y1=rand(1,5), y2=rand(1,5), rand(seed,3), x1=rand(1,5), rand(seed,3), x2=rand(1,5) 结果 7.数据的简单分析 1.当数据为行向量或列向量时，函数对整个向量进行计算. 2.当数据为矩阵时，命令对列进行计算，即把每一列数据当成同一变量的不同观察值。 max(求最大)、min(求最小)、mean(求平均值)、sum(求和)、std(求标准差)、cumsum(求累积和)、median(求中值)、diff(差分)、sort(升序排列)、sortrows(行升序排列)等等。 观察：生成一个3×6的随机数矩阵，并将其各列排序、求各列的最大值与各列元素之和。 程序 rand(seed,1);A=rand(3,6), Asort=sort(A), Amax=m