matlab讲解

3
While 循环,与 for 循环以固定次数求一组命令相反， while 循环以不定的次 数求一组语句的值. while expression {commands} end 只要在表达式(expression)里的所有元素为真， 就执行 while 和 end 语句之 间的命令串{commands}. 例：设银行年利率为 11.25%。将 10000 元钱存入银行，问多长时间会连本 带利翻一番？ 第三种控制： If-Else-End 结构（有三种格式） 一般形式是： if expression {commands} end 如果在表达式(expression)里的所有元素为真，就执行 if 和 end 语句之间 的命令串{commands}.
x 2 +1 x > 1 , 求f (2), f (−1) 例：设 f ( x) = 2 x x ≤ 1
先建立 M 文件 fun1.m 定义函数 f（x） ，再在 Matlab 命令窗口输入 fun1(2),fun1(-1)即可。 例:研究下述程序段 while 1 n = input(’Enter n, negative quits...’); if n <= 0 break; end while n > 1 if rem(n,2) == 0 n = n/2; else n = 3*n + 1; end end end
5
x=-3:0.1:3; y=1:0.1:5; [X,Y]=meshgrid(x,y); Z=(X+Y).^2; mesh(X,Y,Z) meshz(X,Y,Z) 在网格周围画一个 curtain 图(如,参考平面)
（6）hold hold on 与 hold off 保持当前图形, 以便继续画图到当前图上 2 例：将 y=sin(x)、y=cos(x )分别用点和线画出在同一屏幕上 x=linspace(0,2*pi,30); y=sin(x); z=cos(x.^2) plot(x,z,:) hold on Plot(x,y) （7）subplot 格 式 为 h=subplot(mrows,ncols,thisplot) ， 划 分 整 个 作 图 区 域 为 mrows*ncols 块（逐行对块访问）并激活第 thisplot 块，其后的作图语句将图 形画在该块上。 例 ： 将 屏 幕 分 割 为 四 块 ， 并 分 别 画 出 y=sin(x) ， z=cos(x) ， a=sin(x)*cos(x),b=sin(x)/cos(x)。 程序为： x=linspace(0,2*pi,100); y=sin(x); z=cos(x); a=sin(x).*cos(x);b=sin(x)./(cos(x)+eps) subplot(2,2,1);plot(x,y),title(‘sin(x)’) subplot(2,2,2);plot(x,z),title(‘cos(x)’) subplot(2,2,3);plot(x,a),title(‘sin(x)cos(x)’) subplot(2,2,4);plot(x,b),title(‘sin(x)/cos(x)’) （8）scatter 绘制散点图: scatter （X,Y,S,C） ， 在向量 X 和 Y 的指定位置显示彩色圈． X 和 Y 必须大小相同。 （9）contour（ contour（usful point） point） 绘制平面等值线图： contour (x,y,z,n) 绘制 n 个等值线的二维等值线图。 例：在范围-2<x<2,-2<y<3 内绘 z = xe 程序如下： ［X,Y］=meshgeid(-2:.2:2,-2:.2:3); Z=X.*exp(-X.^2-Y.^2); [C,h]=contour(X,Y,Z);
（3） plot3
pl*pi;
plot3(sin(t),cos(t),t) rotate3d %旋转 （4）surf 格式是 surf(x,y,z), 画出数据点（x，y，z）表示的曲面。 例 画函数 Z=(X+Y).^2 的图形. x=-3:0.1:3; y=1:0.1:5; [X,Y]=meshgrid(x,y); Z=(X+Y).^2; surf(X,Y,Z) shading flat %将当前图形变得平滑 (5)mesh 格式为 Mesh(x,y,z) ，用来画网格曲面 例：画出曲面 Z=(X+Y).^2 在不同视角的网格图.
二、
Matlab 的基本操作
1. 基本的算术功能都可以在 command window 内实现； 而且 command history 记录命令的历史。但是，实验室的 C 盘保护，大家注意保存各类数据，包括你的 word 文件，文件名必须清晰易辩，一种简单的命名法是：功能+名字缩写+日期。 2. Matlab 的两类文件 就像其他软件的文件存取一样，Matlab 在程序设计中也需要对文件操作。 但是它的文件和一般的软件不同，主要分为命令性型文件（角码文件）和函数型 文件两类。 （1） 命令型文件 打开白板，直接输入命令行即可，保存文件名自拟。譬如： a=b; b=c; c=a 保存为 jiaohuan.mat。 调用时，在 command window 窗口先定义 a=1, b=2, c=3,jiaohuan，回车得 结果。 （2） 函数型文件 函数性文件顾名思义，肯定是定以后可以象一个函数一样调用它，哪怕是这 个函数的定义异常复杂，经过特殊的变成而来。譬如： f ( x, y ) = 100( y − x 2 ) 2 + (1− x) 2 在（1,2）的数值。
三、Matlab 的绘图功能 绘图功能
Matlab 的绘图功能显然是一般的程序序言难以具备的，而且图形本本事便 于修饰美化、适合程序操作。但是，就图形的美观程度上来说还是应该和其它的 绘图软件结合使用，譬如 SMARTDRAW、几何画板、origin 等等。 （1） 基本命令 plot
4
Plot(x,y,s),X,y 用来定义点，s 表示 c-s 设计。 颜色：y 黄色;m 洋红;c 蓝绿色; r 红色 线型：. 点 - 连线 o 圈 : 短虚线 -.长短线 + 加号 -- 长虚线 等等 例：x=linspace(0,2*pi,30); y=sin(x); z=cos(x); plot(x,y,'r',x,z,’g0')
MATLAB 软件使用简介
主讲人 赵国喜
面向对象：2012 建模队员 编写原则：简洁、实用 一、
Matlab 的宏观认识
1. 源自 Mathworks 公司，要比 C 语言更直观、功能模块更多、且语法容易 掌握，但 Matlab 的源程序还是 C 语言，二者有接口，可以相互转化。 2. 基本功能包括数值计算、符号推导、图形处理等，当然加上扩展的功能 模块就非常强大，但是我们用的较多的还是前三者。 3. 数学软件的四架马车就是 Matlab、Mathematics、Maple、Mathcad，四 者各有千秋。譬如 Mathematics 更长于符号推导，而 Matlab 的矩阵处理功能无 以伦比, Maple 在绘图上也很有特色。 但就某些专项功能而言，上述四种软件并非都是佼佼者。譬如，统计功能 的实现用 R、SPSS、TSP 都很有特色，处理非线性规划 Lingo 做得也很好。 4. 学习 Matlab 的一些网站主要有： /sxjm/ 5.help 帮助功能和软件自带的 demo 功能要非常熟悉
2
删除 A 的第 j1 到 j2 列，构成新矩阵:A(： ， j1:j2)=[ ] 将矩阵 A 和 B 拼接成新矩阵：[A B]；[A；B],这个功能至关重要，可以用 来记录程序生成的数据。 （3）一些简单的矩阵预算 Inv, eig, det，lu, qr, svd, sum 等等 5.数组 （1）某些意义下矩阵和数组的运算一致 >> x = [1 2 3]; >> y = [4 5 6]; >> x + y ans = 5 7 9 （2）数组的乘法 >> z = x. * y z = 4 10 18 （3）数组的除法 Z1 = x. \ y Z1 = 4.0000 2.5000 2.0000 Z2 = x. / y Z2 = 0.2500 0.4000 0.5000 （4）数组的乘幂 >> z3 = x. ˆ y Z3 = 1 32 729 Z4 = x. ˆ 2 4 = 1 4 9 Z5 = 2. ˆ y Z5 = 16 32 64 6.逻辑语句与控制
1
例 1、计算函数
function ff=fun22(x) f=100*(x(2)-x(1)^2)^2+(1-x(1))^2 回车保存文件，文件名自动生成为 fun22.mat,千万不要自己更改文件名， 保存地址自动置于 C:matlab7\work 下。 调用程序，在 command window 输入 x=[1 2], fun22(x) y1 ' = y2 例2、定义则微分方程变为一阶微分方程组： y2 ' = 1000(1 − y12 ) y2 − y1 y (0) = 2, y (0) = 0 1 2 建立 M 文件 vdp1000．m 如下： function dy=vdp1000(t,y) dy=zeros(2,1); dy(1)=y(2); dy(2)=1000*(1-y(1)^2)*y(2)-y(1); 在数值解分析时用 3.变量 MATLAB 中变量的命名规则是： （1）变量名必须是不含空格的单个词； （2）变量名区分大小写； （3）变量 名最多不超过 19 个字符； （4）变量名必须以字母打头，之后可以是任意字母、数字或下划线，变量 名中不允许使用标点符号； （5） （1）MATLAB 的每条命令后，若为逗号或无标点符号，则显示命令的结 果；若命令后为分号，则禁止显示结果； （6） “%” 后面所有文字为注释； （7） “... ...”表示续行（数据太长了，这时就方便了） 。 ... 4.矩阵 （1）逗号或空格用于分隔某一行的元素，分号用于区分不同的行. 除了分 号，在输入矩阵时，按 Enter 键也表示开始一新行. 输入矩阵时，严格要求所有 行有相同的列. 例 m=[1 2 3 4 ；5 6 7 8；9 10 11 12] p=[1 1 1 1 2 2 2 2 3 3 3 3] （2）特殊矩阵生成（编程中使用） c=ones(m，n) 产生一个 m 行、n 列的元素全为 1 的矩阵 c=eyess(m，n)产生一个 m 行、n 列的单位矩阵 b=zeros(m，n)产生一个 m 行、n 列的零矩阵 （3） 矩阵元素的操作 取出矩阵 A 的第 r 行：A（r， ： ） 取出矩阵 A 的第 r 列：A（： ，r） 删除 A 的第 i1 到 i2 行，构成新矩阵:A(i1:i2， ：)=[ ]