数据与函数的可视化
plot指令的基本调用方式
（1）plot（x）
x为向量时，以该元素的下标为横坐标、元素值为纵坐标绘出曲线。
x为实数二维数组时，则按列绘制每列元素值相对其下标的曲线，曲线数等于x数组的列数。
x为复数二维数组时，则按列分别以数组的实部和虚部位横、纵坐标绘制多条曲线。
（2）plot（x，y）
x、y为同维数组时，绘制以x、y元素为横纵坐标的曲线。
x为向量，y为二维数组、且其列数或行数等于x的元素数时，绘制多条不同颜色的曲线。
x为二维数组，y为向量时，情况与上相同，只是y仍为纵坐标。
（3）plot（x1，y1，x2，y2...）
绘制以x1为横坐标、y1位纵左边的曲线1，以x2位横坐标，y2位纵坐标的曲线2，等等。
其中x为横坐标，y为纵坐标，绘制y = f（x）函数曲线。
x = 0:pi/100:2*pi;
y = 2*exp(-0.5*x).*sin(2*pi*x);
plot（x，y）多次叠绘、双纵坐标和多子图
多次调用plot命令在一幅图上绘制多条曲线，需要hold指令的配合
hold on 保持当前坐标轴和图形，并可以接受下一次绘制
hold off 取消当前坐标轴和图形保持，调用plot绘制完全新的图形，不保留以前的坐标格式、曲线
多子图
matlab允许在同一窗口布置几幅独立的子图，具体指令
subplot（m，n，k）
subplot（‘position’，【left，bottom，width，height】）
subplot（m，n，k）的含义，图形窗口包含（mxn）个子图，k为指定的当前子图的编号。
subplot(‘postion’, [left, bottom, width, height])用于手工指定子图位 置，指定位置的四元组采用归一化的标称单位，即认为整个图形窗 口绘图区域的高、宽的取值范围都是[0， 1]，而左下角为（0,0）坐标。
绘制图形的辅助操作
曲线颜色控制符
数据点型控制符
?曲线的线性控制符、颜色控制符、数据点形控制符可以组合使用。
t=(0:15)*2*pi/15;
y=sin(t);
subplot(3,2,1), plot(t, y); title('Lins style is default')
subplot(3,2,2), plot(t, y, 'o'); title('Lins style is o')
subplot(3,2,3), plot(t, y, 'k:'); title('Lins style is k:')
subplot(3,2,4), plot(t, y, 'k:*'); title('Lins style is k:*')
subplot(3,2,5), plot(t, y, 'm-d'); title('Lins style is m-d')
subplot(3,2,6), plot(t, y, 'r-p'); title('Lins style is r-p')
常用的坐标指令
?
?刻度、分格线和坐标框
grid on? ? ? 画出分格线?
grid off? ? ?不画出分格线
box on? ? ? ?控制加边框线
box off? ? ? ? 控制不加边框线
刻度设置：指令格式：set(gca, ‘xtick’, xs, ‘ytick’, ys)
图形标识
图形标识包括：图名（title）坐标轴名（xlabel、ylabel）图形文本注释（text）图例（legend）
简洁使用格式：title（s） xlabel（s） ylabel（） legend（s） text（x，y，s）
字体样式设置：\fontname{arg} \arg \fontsize {arg} string
? ? ? ? ? ? ? ? ? ? ? ? string为要输出的字符串，其前边均为属性控制
上下角标的控制
希腊字母与特殊字符
t=(0:100)/100*2*pi;
y=sin(t);
plot(t, y)
text(3*pi/4,sin(3*pi/4), '\fontsize{16}\leftarrowsin(t) = .707 ')
text(pi, sin(pi), '\fontsize{16}\leftarrowsin(t) = 0 ')
text(5*pi/4, sin(5*pi/4), '\fontsize{16}sin(t)
= -.707\rightarrow','HorizontalAlignment','right')
其中：?‘HorizontalAlignment’, ‘right’设置图形标识 为水平右对齐
特殊图形
直方图（柱形图）---累计式? ?---分组式
水平直方图? ? ? ----累计实? ? ? ---分组式?
饼图pie，饼图用来表示各元素占总和的百分数。该指令的第二输出变量是与第一变量通常的0-1变量，使对应饼块突出。
离散杆图stem
极坐标图polar
三维绘图的基本操作
三维绘图指令plot3，三维绘图指令中，plot3最易于理解，它的使用参数和plot十分相似
图像文件的读写与图像显示
imread指令 -----读取图像文件
imshow指令-----显示图像
imwrite指令------保存图像
matlab命令的执行方式
交互式命令执行方式：逐条输入、逐条执行、操作简单、直观、但是速度慢
M文件的程序执行方式：将命令编成程序储存在一个文件中，依次运行文件中的命令，可以重复进行。
用matlab编写的程序称为M文件，M文件根据调用方式不同分为两类
? ? ? ? 命令文件（script file）? ? ? ?
? ? ? ? 函数文件（function file）
它们的扩展名都是.m
命令文件和函数文件的区别：
????????命令文件没有输入参数，也不返回输出参数，函数文件可以带输入参数，也可以返回输出参数。
????????命令文件对工作空间中的变量进行操作，文件中所有命令的执行结果也返回工作空间中；函数文件中定义的变量为局部变量，当函数文件执行完毕时，这些变量也被清楚。
? ? ? ? 命令文件可以直接运行；函数文件不能直接运行，要以函数调用的方式来调用它。
程序控制结构
任何复杂的程序都可以由这三种基本结构构成：顺序结构、选择结构、循环结构
顺序结构：指按照程序中的语句的排列顺序依次执行，直到最后一个语句
1.数据的输入，从键盘输入数据，可以使用input函数来进行
A = input （提示信息，选项）
2.数据的输出，命令窗口输出函数主要有disp函数，其调用格式为disp（输出项）
A = ‘hello tom’
disp（A）
3.程序的暂停
程序执行过程中暂停，可用pause函数，其调用格式为：pause（延迟描述）
选择结构：根据给定的条件成立或不成立，分别执行不同的语句
1.if语句
（1）单分支if格式
if fix（x） == x?
? ? ? ? disp（x）；
end
（2）双分支格式
if 条件
? ? ? ? 语句组1
else?
? ? ? ? 与剧组2
end
（3）多分支if语句
if 条件1
? ? ? ? 语句组1
elseif 条件2
? ? ? ? 语句组2
elseif 条件m
? ? ? ? 语句组m
else
? ? ? ? 语句组n
end
2.switch语句
switch 表达式
case 表达式1 语句组1
case 表达式2 语句组2 …
case 表达式m 语句组m
otherwise 语句组 n
end
3.try语句
try
语句组1
catch
语句组2
end
循环结构
1.for语句
for 循环变量 = 表达式1：表达式2：表达式3?
? ? ? ? 循环体语句
end
2.while语句
while条件
? ? ? ? 循环体语句
end
3.break和continue语句
一般配合if使用? break 用来终止循环执行，continue用来跳过循环中的某些语句
函数文件
函数文件是另一种形式的M文件，每一个函数文件都定义一个函数
函数文件的基本结构
函数文件由function语句引导，其基本结构为：
? ? ? ? function 输出形参表 = 函数名（输入形参表）
? ? ? ? ? ? ? ? 注释说明部分
? ? ? ? ? ? ? ? 函数体语句
注意：
1. 关于函数文件名
函数文件名通常由函数名再加上扩展名.m组成。 当函数文件名与函数名不同时，Matlab将忽略函数名而确认文件名 因此调用时使用函数文件名。
2. 关于注释说明部分 注释说明包括3部分：
① 紧随引导行之后以%开头的第一注释行。 这一行一般包括大写的函数文件名和函数功能简要描述，供lookfor 关键词查询和help在线帮助时使用。
② 第一注释行及之后连续的注释行。 通常包括函数输入/输出参数的含义及调用格式说明等信息，构成全 部在线帮助文本。
③ 与在线帮助文本相隔一空行的注释行。 包括函数文件编写和修改的信息，如作者和版本等。
3.关于return语句
如果函数文件中插入了return语句，则执行到该语句就结束函数的执行，流程转到调用该函数的位置。通常也不适用return语句。
function [s,p] = fcircle(r)
% FCIRCLE calculate the area and perimeter of a circle of radii r
% r 圆半径
% s 圆面积
% p 圆周长
s = pi*r*r;
p = 2*pi*r;
# 调用
[s,p] = fcircle(10)函数的嵌套使用
function f = factor(n)
if n<=1
f = 1;
else
f = factor(n-1)*n; %递归调用求(n-1)!
end函数参数的可调性
函数所传递参数数目的可调性，即参数的数量可以改变。
在调用函数时，Matlab用两个预定义变量nargin和nargout分别记录 调用该函数时的输入实参和输出实参的个数。
function fout = charray(a,b,c)
if nargin == 1
fout = a;end
if nargin == 2
fout = a+b;end
if nargin == 3
fout = (a*b*c)/2;
end全局变量与局部变量
函数文件中的变量是局部变量
全局变量的作用域是整个matlab工作空间，所有函数都可以对它进行存取和修改
矩阵的分析与处理
常用的产生通用特殊矩阵的函数有
zeros：产生全0矩阵
ones：产生全1矩阵
eye：产生单位矩阵
rand：产生0~1之间均匀分布的随机矩阵
randn：产生均值为0，方差为1的标准正态分布随机矩阵
矩阵结构变换
1.对角阵和三角阵
只有对角线上有非零元素的矩阵称为对角矩阵，在研究矩阵时， 有时候需要将矩阵的对角线上的元素提取出来形成一个列向量，有 时也需要用一个向量构造一个对角阵。
（1）提取矩阵的对角线元素函数：diag
A = [1,2,3;4,5,6];
D = diag(A)
D =
1
5另一种形式：diag（A，k）提取第k条对角线的元素
D1 = diag(A,1)
D =
2
6（2）构造对角矩阵
如果V是一个m个元素的向量，diag(V)将产生一个m×m对角矩阵， 其主对角线元素即为向量V的元素。
diag([1,2,-1,4])
ans =
1 0 0 0
0 2 0 0
0 0 -1 0
0 0 0 4
例如： diag(1:3,-1)
ans =
0 0 0 0
1 0 0 0
0 2 0 0
0 0 3 0矩阵求逆与线性方程组求解
对于一个方阵A，如果存在一个与其同阶的方阵B，使得： (I为单位矩阵） 则称B为A的逆矩阵，当然，A也是B的逆矩阵。求方阵A的逆矩阵可 调用函数inv(A)。
A = [1,-1,1;5,-4,3;2,1,1];
B = inv(A);
A*B
ans =
1.0000 0 0
-0.0000 1.0000 0
-0.0000 0 1.0000将包含n个未知数，由n个方程构成的线性方程组表示为：
A = [1,2,3;1,4,9;1,8,27];
b = [5,-2,6]’;
x = inv(A)*b %x = A\b
x =
23.0000
-14.5000
3.6667矩阵行列式值
把一个方程看做一个行列式，并按行列式的规则求值，称为行列式 的值。在Matlab中，使用函数det(A)得到。
A = rand(5)
A =
0.9501 0.7621 0.6154 0.4057 0.0579
0.2311 0.4565 0.7919 0.9355 0.3529
0.6068 0.0185 0.9218 0.9169 0.8132
0.4860 0.8214 0.7382 0.4103 0.0099
0.8913 0.4447 0.1763 0.8936 0.1389
B = det(A)
B =
-0.0071
多项式
多项式由一个行向量表示
? ? ? ? 该向量元素是该多项式的系数
? ? ? ? 且按照降幂次序排列
如：多项式x 4－12x3＋25x＋116由行向量： p=[1 -12 0 25 116]表示。
>>p=[1 -12 0 25 116]
p =
1 -12 0 25 116
>>r=roots(p)
r =
11.7473
2.7028
-1.2251 + 1.4672i
-1.2251 - 1.4672i注意：matlab规定，多项式是行向量，根是列向量
已知多项式的根，求解多项式？ 使用poly指令
>> pp=poly(r)
pp =
1.0000 -12.0000 -0.0000 25.0000
116.0000多项式乘法（conv指令）
举例：多项式a(x)=x3＋2x2＋3x＋4 和b(x)= x3＋4x2＋9x＋16的乘积。
>> a=[1 2 3 4] ; b=[1 4 9 16];
>> c=conv(a, b)
c =
1 6 20 50 75 84 64多项式除法（deconv）
c(x)=x6＋6x5＋20x4＋50x3＋75x2＋84x＋64 除以b(x)= x3＋4x2＋9x＋16
>> c=[1 6 20 50 75 84 64];
>> b=[1 4 9 16];
>> [q , r]=deconv(c , b)
q =
1 2 3 4
r =
0 0 0 0 0 0 0多项式的导数
举例：求b(x)= x3＋4x2＋9x＋16的导数。
>> b=[1 4 9 16];
>> d=polyder(b)
d =
3 8 9
结果为： 3x2＋8x＋9函数的数值导数
