1-Matlab基本介绍1课件.ppt

什么是数学实验？;先看两个例子; 例2 导弹打击过程仿真; MATLAB程序：---z2 ? k=0; A=[0,0]; %导弹初始位置 B=[0,100]; %飞机初始位置 v=1; dt=1; %离散时间改变量 d=100; %相距距离 while d0.5 plot(A(1),A(2),’ro’); %画导弹位置 hold on plot(B(1),B(2),’b*’); %画飞机位置 pause(0.2); k=k+1; B=B+[v*dt,0]; %飞机移动位置 e=B-A; %导弹指向飞机向量 d=norm(e); e0=e/d; %取向量方向（单位化） A=A+2.0*v*dt*e0; %导弹追击位置 fprintf(k=%.0f A(%.2f,%.2f) B(%.2f,100) d=%.2f\n,k ,A(1),A(2),B(1),d); end ;课程内容 ;《数学实验》课程的考核;什么是MATLAB ？; 20世纪70年代，美国新墨西哥大学计算机科学系主任Cleve Moler为了减轻学生编程的负担，用FORTRAN编写了最早的MATLAB。1984年由Little、Moler、Steve Bangert合作成立了的MathWorks公司正式把MATLAB推向市场。到20世纪90年代，MATLAB已成为国际控制界的标准计算软件。 ;特点;MATLAB命令简单介绍 ;直接点击windows桌面上的matlab图标，或者运行MATLAB，就可创建一个或多个窗口 ;b)?编辑区(Editor\Debugger Window)编制各种M-文件，存盘(Save)、运行（Run）等，在编辑窗口完成并存盘的命令程序默认生成m-文件。程序运行的结果和相关信息显示在命令窗口中。若所写程序不符合要求，则会出现错误提示信息。; 工具箱及命令查询 help topics;MATLAB中基本代数运算符;变量及数组输入 ;2：若数组元素之间有一定的规律，则可用特殊命令格式输入。 例 x=a:p:b ，a,b,p分别表示初值、终值、步长，p为1时可以省略。 还可以利用linspace()函数生成向量。 linspace(a,b,n) 将区间[a,b]等分成n-1段，返回由端点及分段点坐标所产生的n个元素的行向量。 ;矩阵的输入和数组类似。 eye(n,n) 生成n维单位矩阵 ones(m,n) 生成全是1的m行n列矩阵 zeros(m,n) 生成全是0的m行n列矩阵 rand(m,n) 产生均匀分布的m行n列随机矩阵 randn(m,n) 产生正态分布的m行n列随机矩阵 reshape(a,m,n) 在总元素不变的情况下，将矩 阵a改变成m行n列矩阵;数组、矩阵元素的访问;MATLAB中数组、矩阵基本运算符;例1、求方程组Ax=B的解---------q1.m;例2、构造矩阵的方法------q3;例3、 矩阵运算------q2; size(A) %矩阵A的规格,即行数和列数 rank(A) %矩阵A的秩 rref(A) %矩阵A的最简行阶梯型 [Q,R]=qr(A) %矩阵A的QR正交分解 [L,U]= lu(A) %矩阵A的LU分解 trace(A) %矩阵A的迹 ;MATLAB中的数值函数运算 ;1 求函数值：f(2)、u(2,3) 画数值函数的图形： fplot(f,[0,2]); %画一元函数在[0,2]上的图象 3 求函数零点 c=fzero(f,1); %求一元函数在x=1附近的零点c. c=fzero(f,[a,b]); %求一元函数在区间[a,b]内的零点c. c=fzero(u,1,[ ],2) %求二元函数y=2时在x=1附近的零点c. ;4 求数值函数的最大最小值 [x,y]=fminbnd(f,a,b) %求一元函数f在区间[a,b