预览加载中,请您耐心等待几秒...
1/10
2/10
3/10
4/10
5/10
6/10
7/10
8/10
9/10
10/10

亲,该文档总共100页,到这已经超出免费预览范围,如果喜欢就直接下载吧~

如果您无法下载资料,请参考说明:

1、部分资料下载需要金币,请确保您的账户上有足够的金币

2、已购买过的文档,再次下载不重复扣费

3、资料包下载后请先用软件解压,在使用对应软件打开

§1.MATLAB操作环境一、MATLAB的安装、启动和退出图1MATLAB桌面工作环境二、MATLAB的窗口构成三、MATLAB的语言要素【注】①几个特殊常量:3.变量、函数、表达式4.语句(命令)§2.矩阵及其运算一、矩阵的输入或>>A=[123 456](回车)例2.>>z=zeros(2,3),u=ones(3)例3.>>e=eye(3,4),x=rand(1,3)3.利用M-文件生成矩阵4.载入外部数据文件生成矩阵(load,fread)二、向量(数组)的输入例7.>>A=1:5;%A=12345(公差缺省为1)三、矩阵的裁剪、拼接与修改2.矩阵的拼接3.矩阵的编辑修改A= 1210 5601 0101112四、矩阵的数值运算2.标量与矩阵的运算例11.>>A=[123;312],B=[210;302],C=[12;30]>>X1=A+B,X2=C*B,X3=A\B>>X4=A',X5=C^2>>Z1=A.*B,Z2=A.\B,Z3=A.^B§3.语句、函数和窗口环境一、语句(命令)例1.>>a=[12345];b=[13579];... >>c=a.*b,d=a*b',e=a'*b2.几个特殊常量二、函数例3.设x=[00.250.50.751],计算函数②简单函数可以用inline函数形式输入: y=feval(inline(‘函数表达式'),x)2.向量函数(2)构造矩阵函数: zeros(n,m)0阵ones(n,m)1阵eye(n,n)单位阵 rand(n,m)(0,1)均匀分布随机阵 diag(A)(构造或提取对角阵) triu(A)/tril(A)(提取上/下三角阵) reshape(A,m,n)(重新整理A的维数)例5.熟悉函数max(x),sum(x),length(x),size(X),diag(X)的功能。>>D=diag(A),DD=diag(B)三、在线帮助与数据显示格式2.数据显示格式四、MATLAB工作区§4.图形功能一、2维图形例1.绘制正玄曲线2.画多重线(在同一个坐标系下画多条曲线)例2.在同一个坐标系下绘制正玄曲线及余玄曲线3.网格与标记(在图上加网格、标题、轴标记)例3.绘制正玄,余玄曲线,并加上标题,网格,注释等。4.多幅图形>>x=linspace(0,2*pi,30);y=sin(x); >>subplot(2,2,1),plot(x,y) >>axis([02*pi-11]),title('y=sin(x)')5.其它有关函数二、三维图形2.绘制网格曲面(分4步)例7.绘制“阔边帽”曲面:§5.MATLAB程序设计MATLAB提供了一个完善的程序设计环境,使用户能够方便地编制复杂的程序,完成各种计算。本节先介绍关系运算、逻辑运算,再介绍M-文件(即程序文件)的结构及MATLAB的程序控制流语句。一、关系运算和逻辑运算D=00100②当作用于两个同样大小矩阵时,则分别对两个矩阵的对应元素运算,结果是一个0-1矩阵。(真值表见P27)d=000111111find(x)1个输出参数:返回向量x中非零元下标,若x 是矩阵,则视其为一长列向量。 2个输出参数:分别返回矩阵x中非零元的行 标和列标。例3.对向量或矩阵进行逻辑运算,或调用逻辑函数。y=134二、M-文件的结构例4.给定矩阵C=524 468>>A=[123;456];B=[401;012];%输入参数赋值三、M-文件的建立、编辑与调用其中,实际参数和形式参数可以不同;输入参数要用圆括号将其括起,输出参数多于一个时,要用方括号“[]”将其括起;参数间用逗号隔开。四、MATLAB的控制流语句解:编写函数文件(fenduan.m)如下:switch<表达式> case值1 <语句体1> … case值n <语句体n> otherwise <语句体n+1> end例7.编写命令文件,从键盘输入某学生成绩(5分制), 由计算机判断并输出该生成绩的优良等级.基本功能:循环循环变量依次取矩阵的每一列,然后执行一次语句体;遍历矩阵各列后,然后执行end的后续命令.。>>%调用函数 >>h=hhilb(3,4) h=11/21/31/4 1/21/31/41/5 1/31/41/51/6例8*.编写函数文件,计算并输出矩阵A的各行向量的 元素和.(2)while语句>>eeps%运行程序 n=53 EPS=2.2204e-0163.其它控制流语句X=1013 2831pause/pause(n)暂停程序执行(n秒),按任意键 继续;functionX=anysum(A,B) try X=A+B; catch k=min([size(A);size(B)]); X=A(1:k(1),1:k(2))+B(1:k(1),1:k(2));