1 Matlab初步.pptVIP

数学建模实验 Matlab Lingo Matlab软件 MATLAB是矩阵实验室（Matrix　Laboratory）之意。除具备卓越的数值计算能力外，它还提供了专业水平的符号计算，文字处理，可视化建模仿真和实时控制等功能。 MATLAB的基本数据单位是矩阵，它的指令表达式与数学,工程中常用的形式十分相似,故用MATLAB来解算问题要比用C,FORTRAN等语言完相同的事情简捷得多. 当前流行的MATLAB 包括拥有数百个内部函数的主包和三十几种工具包(Toolbox).工具包又可以分为功能性工具包和学科工具包.功能工具包用来扩充MATLAB的符号计算,可视化建模仿真,文字处理及实时控制等功能.学科工具包是专业性比较强的工具包,控制工具包,信号处理工具包,通信工具包等都属于此类. 教学计划 符号函数运算 数与形 矩阵运算初步 数值函数差分与积分 数据处理 统计分析初步 回归分析与函数拟合 数据补充与函数插值 优化问题与函数极值 线性规划与Lingo软件 动态仿真 随机现象的蒙特卡洛模拟 概率模型与解析模拟 *存储问题与马尔可夫决策 *物流管理与图论算法 数学建模实验（一） Matlab初步 §1.1 函数及其图像 一、函数值及图像 初等函数数值直接表示： 幂函数 sqrt sin asin tan exp log 1. 求连续函数在一点函数值 sin(2) 2.求函数在一列点上的值 x=1:6;y=x.^2+sin(x); 3. 画函数的图像 plot(x,y),grid 进一步: 二元函数的定义、求值、画图 x=-3:0.1:3;y=x; [X,Y]=meshgrid(x,y); Z=X.^2+cos(4.*Y); surf(X,Y,Z); mesh(X,Y,Z); 尝试meshc，meshz，surfc，surfl ezsurf(’x.^2 +sin(y)’) 二、自定义函数 (1)定义内联函数 fun=inline(‘x^2+y^2’,’x’,’y’); (2)函数句柄 fun=@(x,y) x.^2+sin(y); (3)建立M函数文件 在Matlab命令窗口的File菜单中选择New命令，出现子菜单， 选择m-file的命令，将得到M文件的编辑窗口 编辑M文件：可以是一段程序或一个函数。 编辑完成后，在此窗口的File菜单中选择save as的命令，得到对话 框，输入文件名 ***.m 三、函数运算求函数值 fun(3) x=[-0.5 0 0.8 18]； fun(x); feval(fun,3) 求初等函数在一列点上的函数值 function y=fun(x) y=x^2+sin(x); Save as fun.m x=1:10; eval(‘x.^2+sin(x)’) fvec=vectorize(fun) feval(fvec,x) 四、矩阵运算——Matlab软件的核心 1. 矩阵的生成 a=[1 2 3 -4; 5 6 7 8; -1 3 -5 7] %生成矩阵 b=[a; 2 4 6 8] % 添加一行 c=[b,b(4,:)’] % 将b的第四行转置为一列添在后面 d=ones(2) %2×2的全 1 矩阵 f=ones(3;1) %全1向量 I=eye(3) %3×3的单位矩阵 三维矩阵生成方法：repmat 复制矩阵，cat联合矩阵 v=repmat(0,[3,5,7]) C=cat(3,zeros(3,5),ones(3,5),randn(3,5),rand(3,5)) 2. 矩阵基本运算 h=A’ %矩阵转置 A±B %矩阵（与数）加减 A.*B , A./B %矩阵对应元素相乘，除 X=A*B %矩阵相乘 X=A.^2 %矩阵各元素乘方 X=inv(A) %方阵A求逆 X=A\B %左除，A*X=B, X=inv(A)*B, 解方程 X=A/B %右除，X*B=A, X=A*inv(B) size(A) %矩阵A的维数 length(a) %向量 a 元素的个数 3. 矩阵函数 sin(a)+2.*a b.^2 b^2 det(b) rank(b) a1=sum(a) %矩阵 a 各列元素的和 sum(a1) %向量 a1 各分量的和 a1=max(a) %矩阵 a 各列元素的最大元 max(a1) %向量 a1 各分量的的最大元 4.矩阵的特征值，特征向量 AX=λX eig(A); %返回值就是特征值矩阵和特征向量矩阵,注意该命令要求A为方阵。 5.矩阵的基础解系 null(A) 图像——矩阵 Matlab图像读入： X