读取图像：用imread函数读取图像文件，文件格式可以是TIFF、JPEG、GIF、BMP、PNG 等。比如
>> f = imread('chestxray.jpg');
读进来的图像数据被保存在变量f中。尾部的分号用来抑制输出。如果图片是彩色的，可以用rgb2gray转换成灰度图：
>> f = rgb2gray(f);
然后可以用size函数看图像的大小
>> size(f)
如果f是灰度图像，则可以用下面的命令把这个图像的大小赋给变量M和N
>> [M, N] = size(f);
用whos命令查看变量的属性
>> whos f
显示图像：用imshow显示图像
imshow(f, G)
其中f是图像矩阵，G是像素的灰度级，G可以省略。比如
>> imshow(f, [100 200])
图像上所有小于等于100的数值都会显示成黑色，所有大于等于200的数值都会显示成白色。pixval命令可以用来查看图像上光标所指位置的像素值。
pixval
例如
>> f = imread('rose_512.tif');
>> whos f
>> imshow(f)
如果要同时显示两幅图像，可以用figure命令，比如
>> figure, imshow(g)
用逗号可以分割一行中的多个命令。imshow的第二个参数用一个空的中括号：
>> imshow(h, [])
可以使动态范围比较窄的图像显示更清楚。
写图像。用imwrite写图像
imwrite(f, 'filename')
文件名必须包括指明格式的扩展名。也可以增加第三个参数，显式指明文件的格式。比如
>> imwrite(f, 'patient10_run1.tif', 'tif')
也可以写成
>> imwrite(f, 'patient10_run1.tif')
还可以有其他参数，比如jepg图像还有质量参数：
>> imwrite(f, 'filename.jpg', 'quality', q)
q是0到100之间的一个整数。对比不同质量的图像效果。用imfinfo命令可以查看一个图像的格式信息，比如
>> imfinfo bubbles25.jpg
可以把图像信息保存到变量中
>> K = imfinfo('bubbles25.jpg');
>> image_bytes = K.Width * K.Height * K.BitDepth / 8;
>> compressed_btyes = K.FileSize;
>> compression_ratio = image_bytes / compressed_bytes
数据类型。MA TLAB的数据类型包括：
double 双精度浮点
uint8 无符号8位整数
uint16 无符号16位整数
uint32 无符号32位整数
int8 有符号8位整数
int16 有符号16位整数
int32 有符号32位整数
single 单精度
char 字符
logical 逻辑型（二值）
数据类型转换
B = data_class_name(A)
比如
>> C = [1.4 1.5]
>> D = uint8(C)
图像类型分为：
Intensity image 灰度图
Binary image 二值图
Indexed image 索引图
RGB image 彩色图
在灰度图中每个像素可以是整型、浮点型或者逻辑型。图像类型的像素类型可以转换
function to from
im2uint8 uint8 logical,uint8,uint16,double
im2uint16 uint16 logical,uint8,uint16,double
mat2gray double double
im2double double logical,uint8,uint16,double
im2bw logical uint8,uint16,double
比如
g = mat2gray(A, [Amin, Amax]);
g = mat2gray(A);
g = im2double(h);
g = im2bw(f, T)
其中A是浮点型的图像，Amin和Amax是浮点数的范围，h和f是任意类型的图像，T 是分割的阈值。
数组（向量）索引：创建向量（数组）：
>> v = [1 3 5 7 9 11 13]
用小括号对向量进行索引（取数组中的某个元素）：
>> v(2)
转置（将行向量通过转置变成列向量）：
>> w = v.'
取向量其中的一部分：
>> v(1:3) 第1个到第3个
>> v(2:4)
>> v(3:end) 第3个到最后一个
>> v(1:end)
>> v(:) 全部
>> v(1:2:end) 第1个到最后一个，每次增加2
>> v(end:-2:1) 最后一个到第1个，每次减2
其中end总是表示最后一个。
>> x = linspace(1, 5, 10)
>> v([1 4 5])
linspace函数产生一个范围内的平均分布。
矩阵索引：创建矩阵
>> A = [1 2 3; 4 5 6; 7 8 9]
取矩阵中的一个元素
>> A(2, 3)
取矩阵中的一行或者一列
>> C3 = A(:, 3)
>> R2 = A(2, :)
取矩阵中某些行某些列
>> T2 = A(1:2, 1:3)
对矩阵中某些元素进行赋值：
>> B = A;
>> B(:, 3) = 0
用end表示最后一行或者最后一列：
>> A(end, end) 最后一行最后一列
>> A(end, end-2) 最后一行倒数第三列
>> A(2:end, end:-2:1) 第2行到最后一行，最后一列到第一列，每次减2 >> E = A([1 3], [2 3]) 第1、3行，第2、3列
>> D = logical([1 0 0; 0 0 1; 0 0 0])
>> A(D) 取A中由D指定的位置上的元素
>> v = T2(:) 把矩阵变成一个向量
>> s = sum(A(:)) 求和
>> sum(sum(A))
可以把矩阵操作用在图像上
>> f = imread('filename');
>> fp = 矩阵列倒转
>> imshow(fp)
>> fc =
>> imshow(fc)
>> fs =
>> imshow(fs)
>> plot(f(512,:))
矩阵可以是多维的，用size看矩阵大小，用ndims命令常看矩阵的维数
>> size(A, 1)
>> ndims(A)
一些常用的矩阵
zeroes(M, N)
ones(M, N)
true(M, N)
false(M, N)
magic(M)
rand(M, N)
randn(M, N)
其中M、N表示矩阵的行数和列数。比如
>> A = 5 * ones(3, 3)
>> magic(3)
>> B = rand(2, 4)
函数：可以把一系列的MA TLAB语句或者一个带参数的函数放在扩展名叫做m的文件中。一个带函数的m文件有一下部分组成
函数定义行
H1行
帮助部分
函数体
注释
函数定义行的格式是
function [outputs] = name(inputs)
比如要写一个函数计算两个图像的和以及乘积
function [s, p] = sumprod(f, g)
其中f和g是输入的图像，而s是和，p是乘积。返回值用中括号括起来，如果返回值只有一个，可以省略中括号。如果函数没有输出，则中括号和等号都可以省略。函数名字的命名规则和C语言是相同的。定义好的函数可以在命令行调用：
>> [s, p] = sumprod(f, g);
也可以被其它函数调用。如果只有一个返回值，调用时中括号也是可以省略的，比如
>> y = sum(x);
H1行是文本的第一行，是一个单行的注释，紧跟在函数定义行后面，之间不能有空行。比如
% SUMPROD Computes the sum and product of two images.
百分号开始的文字表示注释。当使用帮助命令
>> help function_name
时，这个H1行会被首先显示出来。如果使用lookfor命令，则会在所有H1行中查找指定的关键字。这一行应该提供这个函数功能的一个概述。帮助部分是紧跟在H1后的文本块，中间没有空行，用来提供对这个函数更详细的帮助说明。在使用help命令时会显示所有这部分内容。这部分内容由注释语句构成，全部由%开始。接下来第一个非注释语句表示函数体的开始。函数体包含进行计算的语句和给返回值赋值的语句。函数题中的所有注释（百分号开始的行）被认为是普通的注释，不是H1或者帮助部分。m文件可以用任何文本编辑器创建和编辑，只要用.m扩展名保存在MA TLAB可以搜索到的路径里面。另一个创建和编辑函数的方法是在命令行输入edit命令，比如
>> edit sumprod
这命令会编辑已经存在的sumprod.m文件，如果没有则自动在当前目录中创建一个sumprod.m并开始编辑。
运算符。运算符可以分为算术运算符，关系运算符和逻辑运算符。算术运算符分为矩阵算术运算符和数组算术运算符。
+ 矩阵和数组加法plus(A, B) a+b, A+B
- 矩阵和数组减法minus(A,B) a-b, A-B
.* 数组乘法times(A,B) C=A.*B, 意味着C(I,J) = A(I,J)*B(I,J)
* 矩阵乘法mtimes(A,B) A*B, 表示线性代数中的矩阵运算,或者a*A
./ 数组右除rdivide(A,B) C=A./B, 意味着C(I,J)=A(I,J)/B(I,J)
.\ 数组左除ldivide(A,B) C=A.\B, 意味着C(I,J)=B(I,J)/A(I,J)
/ 矩阵右除mrdivide(A,B) A/B 意味着A*inv(B), inv是矩阵求逆
\ 矩阵左除mldivide(A,B) A\B 意味着inv(A)*B
.^ 数组指数power(A, B) C=A.^B,意味着C(I,J)=A(I,J)^B(I,J)
^ 矩阵指数mpower(A,B) 请查看帮助
.' 向量和矩阵转置transpose(A) A.'
' 复数的共轭ctranspose(A) A'
+ 单目加号uplus(A) +A 与0+A相同
- 单目负号uminus(A) -A 与0-A相同
图像处理工具包还提供其他一些算术运算
imadd 两个图像相加，或者一个图像加上一个常量
imsubstract 两个图像相减，或者一个图像减掉一个常量
immultiply 两个图像相乘，或者一个图像乘上一个常量
imdivide 两个图像相除，或者一个图像除以一个常量
imabsdiff 两个图像的差的绝对值
imcomplement 求一个图像的反色图
inlincomb 求一组图像的线性组合
关系运算符包括
<
<=
>
>=
==
~=
关系运算符的结果是逻辑型的矩阵，比如
>> A = [1 2 3; 4 5 6; 7 8 9]
>> B = [0 2 4; 3 5 6; 3 4 9]
>> A == B
>> A >= B
如果关系运算符两边都是矩阵，则要求两边的矩阵是同样大小的。或者一边是矩阵一边
是常数，或者两边都是常数，那也是可以。
>> A > 3
>> 3 ~= 4
逻辑运算符包括与、或、非三个运算
& AND
| OR
~ NOT
在matlab中非0被认为是真，0被认为是假。比如
A = [1 2 0; 0 4 5]
B = [1 -2 3; 0 1 1]
>> A & B
MATLAB还有其它一些逻辑运算函数：
xor 异或
all 如果一整列都是真，则结果是真
any 如果一整列只要有一个是真，则结果是真
比如
>> xor(A, B)
>> all(A)
>> any(A)
>> all(B)
>> any(B)
一些重要的常量
ans
eps
i(或者j)
NaN或者nan
pi
realmax
realmin
computer
version
MATLAB中一般常量的写法
3
-99
0.00001
9.6397238
1.60210e-20
6.02252e23
1i
-3.14159j
3e5j
控制流。包括
if if和else, elseif组合，条件执行一组语句
for 指定次数重复执行一组语句
while 按条件反复执行一组语句
break 终止for或者while循环
continue 马上开始下一次for或者while循环
switch switch和case,otherwise结合，按照条件值的不同执行不同的语句块return 终止当前函数，返回到调用它的地方
try...catch 捕获异常状况
if语句
if expression_r1
statements1
elseif expression_r2
statements2
else
statements3
end
for循环
for index = start:increment:end
statements
end
比如
count = 0;
for k = 0:0.1:1
count = count +1
end
while循环
while expression_r
statements
end
比如
a = 10;
b = 5;
while a
a = a - 1;
while b
b = b - 1;
end
end
switch语句
switch switch_expression_r
case case_expression_r
statements
case {case_expression_r1, case_expression_r2}
statements
otherwise
statements
end
比如
switch newclass
case 'uint8'
g = im2uint8(f);
case 'uint16'
g = im2uint16(f);
case 'double'
g = im2double(f);
otherwise
error('Unknown or improper image class.')
end
例子：写一个函数计算一幅灰度图像所有像素的平均值
function av = average(A)
%A VERAGE Computes the average value of an array
% A V = A VERAGE(A) computes the average value of input array, A, % which must be a 1D or 2D array.
% Check the validity of input.
if ndims(A) > 2
error('The dimensions of the input cannot exceed 2.')
end
% Computes the average
av =
例子：比较各种不同的JPEG质量下的图像质量
for q = 0:5:100
filename = sprintf('series_%3d.jpg', q);
imwrite(f, filename, 'quality', q);
end
其中sprintf语句和c语言的fprintf语句用法类似。例子：写一个函数从一个图像中取出一个矩形的子图。
function s = subdim(f, m, n, rx, cy)
%SUBDIM Extracts a subimage, s, from a given image, f.
% The subimage is of size m-by-n, and the coordinates of its top, left
% corner are (rx, cy).
s = zeros(m, n);
rowhigh = rx + m - 1;
colhigh = cy + n - 1;
xcount = 0;
for r = rx:rowhigh
xcount = xcount + 1;
ycount = 0;
for c = cy:colhigh
ycount = ycount + 1;
s(xcount, ycount) = f(r, c);
end
end
实际上这个功能可以用一个matlab语句就可以实现了。
循环优化。一些循环可以被转换成同样向量计算来代替，比如f(x)=Asin(x/2pi)，生成一个向量包含一组函数的值：
for x = 1:M
f(x) = A*sin((x-1)/(2*pi));
end
可以用下面两个语句来代替
x =
f = A*sin(x/2(*pi))
向量运算要比循环快得多。对于二维的情况，MATLAB提供了一个meshgrid函数
[C, R] = meshgrid(c, r)
比如
>> r = [0 1 2];
>> c = [0 1];
>> [C, R] = meshgrid(c, r)
>> h = R.^2 + C.^2
例子：f(x, y) = Asin(ux + vy)，生成一个矩阵，包含这个函数的值：
function [rt, f, g] = twodsin(A, u0, v0, M, N)
%TWODSIN compare for loops vs vectorization
% The comparision is based on implementing the function
% f(x, y) = Asin(u0x+v0y) for x = 0, 1, 2,..., M-1, and
% y = 0, 1, 2, ..., N-1. The inputs to the function are
% M and N and the constants in the function.
tic
for r = 1:M
u0x = u0*(r-1);
for c = 1:N
v0y = v0*(c-1)
f(r, c) = A*sin(u0x + v0y);
end
end
t1 = toc;
tic
r =
c =
[C,R] = meshgrid(c, r);
g = A*sin(u0*R + v0*C)
t2 = toc;
rt = t1/t2;
运行这个例子可以看出，向量计算至少比循环要快30倍。如果只是针对图像的一部分进行操作，可以简单的提取出来，比如
rowhigh = rx + m - 1;
colhigh = cy + n - 1;
s =
同样功能的程序，这个写法比前面例子看到的写法要快1000倍。
用户交互。函数disp用来提示用户一些信息。比如
>> A = [1 2; 3 4];
>> disp(A)
>> sc = 'Digital Image Processing.';
>> disp(sc)
>> disp('This is another way to display text.')
input函数用来提示用户输入某些值：
t = input('message')
这个函数可以显示message，并等待用户输入一个值，并存到t中。输入的值可以是MATLAB允许的任何类型的值。而如下的格式只接受字符串输入
t = input('message', 's')
如果字符串中包含的都是数字，则可以用函数str2num进行转换
n = str2num(t)
比如
>> t = input('Enter your data: ', 's')
>> class(t)
>> size(t)
>> n = str2num(t)
>> size(n)
>> class(n)
如果混合输入字符串和数值，可以全部按照字符串读入，再使用字符串处理函数strread，比如：
>> t = '12.6, x2y, z';
>> [a, b, c] = strread(t, '%f%q%q', 'delimiter', ',')
>> d = char(b)
其中%f表示浮点数，%q表示字符串。delimiter参数表示分割符。函数strcmp用来比较字符串，如果两个字符串相等返回真，否则返回假。lower函数可以把字符串中的大写字母全部变成小写，upper函数可以把字符串中所有小写字母全部变成大写。
Cell数组和结构体。Cell数组是指数组的元素本身还是一个数组，比如
>> c = {'gause', [1 0; 0 1], 3}
>> c{1}
>> c{2}
>> c{3}
结构体和Cell数组类似，但是其中的元素是用一个名字去访问的，比如
>> S.char_string = 'gause';
>> S.matrix = [1 0; 0 1];
>> S.scalar = 3;
>> S.matrix
实验任务
MA TLAB 没有提供一个函数来确定数组的每个元素是不是整数（即...,-2,-1,0,1,2,...）。写一个函数来实现这个功能，并满足下面的格式：
function I = isinteger(A)
%ISINTEGER Determines which elements of an array are integers.
% I = ISINTEGER(A) returns a logical array, I, of the same size
% as A, with 1s (TRUE) in the locations corresponding to integers
% (i.e., . . . -2 -1 0 1 2 . . . )in A, and 0s (FALSE) elsewhere.
% A must be a numeric array.
要求不能使用while或者for循环。提示：使用floor函数。如果你要做到能够处理复数，那么可以使用real和imag函数。
MATLAB没有提供一个函数用来确定数组的每个元素是不是偶数（即...,-4,-2,0,2,4,...）。写一个函数来实现这个功能，并满足下面的格式：
function E = iseven(A)
%ISEVEN Determines which elements of an array are even numbers.
% E = ISEVEN(A) returns a logical array, E, of the same size as A,
% with 1s (TRUE) in the locations corresponding to even numbers
% (i.e., . . . -3, -1, 0, 2, 4, . . . )in A, and 0s (FALSE) elsewhere.
% A must be a numeric array.
要求不能使用while或者for循环。提示：使用floor函数。
MATLAB没有提供一个函数用来确定数组的每个元素是不是奇数（即...,?3,?1,1,3,...）。写一个函数实现这个功能，并满足下面的格式：
function D = isodd(A)
%ISODD Determines which elements of an array are odd numbers.
% E = ISODD(A) returns a logical array, D, of the same size as A,
% with 1s (TRUE) in the locations corresponding to odd numbers
% (i.e., . . . -3, -1, 1, 3, . . . )in A, and 0s (FALSE) elsewhere.
% A must be a numeric array.
要求不能使用while或者for循环。提示：使用floor函数。
写一个函数满足下面的规定：
function H = imcircle(R, M, N)
%IMCIRCLE Generates a circle inside a rectangle.
% H = IMCIRCLE(R, M, N) generates a circle of radius R centered
% on a rectangle of height M and width N. H is a binary image with
% 1s on the circle and 0s elsewhere. R must be an integer >= 1.
你的程序应该检查R的合理性，使它能够在矩形的内部。要求不能使用for或者while 循环。提示：使用函数meshgrid和floor函数。
这个实验主要目的是学习如何显示和改变目录，并利用目录信息来读取图像。写一个函数满足下面的规定：
function [I, map] = imagein(path)
%IMAGEIN Read image in from current-working or specified directory.
% I = IMAGEIN displays a window containing all the files in the
% current directory, and saves in I the image selected from the
% current directory.
% [I, MAP] = IMAGEIN variable MAP is required to be an output
% argument when the image being read is an indexed image.
% [ . . .] = IMAGEIN('PATH') is used when the image to be read
% resides in a specified directory. For example, the input
% argument 'C:\MY_WORK\MY_IMAGES' opens a window showing
% the contents of directory MY_IMAGES. An image selected from
% that directory is read in as image I.
提示：使用在线帮助来熟悉cd，pwd，uigetfile函数的使用。或者也可以采用fullfile函数来代替cd函数
一、MATLAB常用的基本数学函数 abs(x)：纯量的绝对值或向量的长度 angle(z)：复数z的相角(Phase angle) sqrt(x)：开平方 real(z)：复数z的实部 imag(z)：复数z的虚部 conj(z)：复数z的共轭复数 round(x)：四舍五入至最近整数 fix(x)：无论正负，舍去小数至最近整数 floor(x)：地板函数，即舍去正小数至最近整数ceil(x)：天花板函数，即加入正小数至最近整数rat(x)：将实数x化为分数表示 rats(x)：将实数x化为多项分数展开 sign(x)：符号函数(Signum function)。 当x<0时，sign(x)=-1； 当x=0时，sign(x)=0; 当x>0时，sign(x)=1。 rem(x,y)：求x除以y的馀数 gcd(x,y)：整数x和y的最大公因数 lcm(x,y)：整数x和y的最小公倍数 exp(x)：自然指数 pow2(x)：2的指数 log(x)：以e为底的对数，即自然对数或 log2(x)：以2为底的对数 log10(x)：以10为底的对数 二、MATLAB常用的三角函数 sin(x)：正弦函数 cos(x)：余弦函数
tan(x)：正切函数 asin(x)：反正弦函数 acos(x)：反馀弦函数 atan(x)：反正切函数 atan2(x,y)：四象限的反正切函数 sinh(x)：超越正弦函数 cosh(x)：超越馀弦函数 tanh(x)：超越正切函数 asinh(x)：反超越正弦函数 acosh(x)：反超越馀弦函数 atanh(x)：反超越正切函数 三、适用於向量的常用函数有： min(x): 向量x的元素的最小值 max(x): 向量x的元素的最大值 mean(x): 向量x的元素的平均值 median(x): 向量x的元素的中位数 std(x): 向量x的元素的标准差 diff(x): 向量x的相邻元素的差 sort(x): 对向量x的元素进行排序（Sorting）length(x): 向量x的元素个数 norm(x): 向量x的欧氏（Euclidean）长度sum(x): 向量x的元素总和 prod(x): 向量x的元素总乘积 cumsum(x): 向量x的累计元素总和cumprod(x): 向量x的累计元素总乘积 dot(x, y): 向量x和y的内积 cross(x, y): 向量x和y的外积 四、MATLAB的永久常数
Matlab自定义函数 1、函数文件+调用命令文件：需单独定义一个自定义函数的M文件; 2、函数文件+子函数：定义一个具有多个自定义函数的M文件； 3、Inline:无需M文件，直接定义； 4、Syms+subs:无需M文件,直接定义； 5、字符串+subs：无需M文件,直接定义. 6、匿名函数 7、直接通过@符号定义. 1、函数文件+调用函数文件：定义多个M文件： %调用函数文件:myfile.m clear clc for t=1:10 y=mylfg(t);%调用函数时要注意实参与形参的匹配！ fprintf(‘%4d^(1/3)=%6.4f\n’,t,y); end %自定义函数文件:mylfg.m function y=mylfg(x)%注意：函数名（mylfg）必须与文件名（mylfg.m）一致 Y=x^(1/3); 注：这种方法要求自定义函数必须单独写一个M文件，不能与调用的命令文件写在同一个M文件中。 2、函数文件+子函数：定义一个具有多个子函数的M文件 %函数文件：funtry2.m function[]=funtry2()%可以无自变量()或无因变量[] for t=1:10 y=lfg2(t); fprintf('%4d^(1/3)=%6.4f\n',t,y); end function y=lfg2(x)%%子函数 y=x^(1/3);
%注：自定义函数文件funtry2.m中可以定义多个子函数function。子函数lfg2只能被主函数和主函数中的其他子函数调用。 3、Inline:无需M文件，直接定义； %inline命令用来定义一个内联函数：f=inline(‘函数表达式’,‘变量1’,’变量2’,……)。 调用方式：y=f(数值列表)%注意：代入的数值列表顺序应与inline()定义的变量名顺序一致。 例如： f=inline(‘x^2+y’,’x’,’y’); z=f(2,3) Ans=7 注：这种函数定义方式是将它作为一个内部函数调用。特点是，它是基于Matlab的数值运算内核的，所以它的运算速度较快，程序效率更高。缺点是，该方法只能对数值进行代入，不支持符号代入，且对定义后的函数不能进行求导等符号运算。 内联函数定义方式是将f作为一个内部函数调用。其特点是：调用方式最接近于我们平时对函数的定义，使程序更具可读性。同时由于它是基于Matlab的数值计算内核的，所以它的运算速度较快，程序更有效率。 这种定义方式的缺点： 定义一个内联函数用去的内存空间比相同条件下其他的方法要大得多。 该方法只能对数值进行代入，不支持符号代入，并且对于定义后的函数不能进行求导等符号运算。 例：通过命令clear清除工作空间的所有变量后，执行如下指令 Clear Clc f=’x^2’; Syms x g; g=x^2; h=inline(‘x^2’,’x’); whos 4、Syms+subs:无需M文件,直接定义； 用syms定义一个符号表达式，用subs调用： Syms f x%定义符号 f=1/(1+x^2);%定义符号表达式也是符号
信源函数 randerr 产生比特误差样本 randint 产生均匀分布的随机整数矩阵 randsrc 根据给定的数字表产生随机矩阵 wgn 产生高斯白噪声 信号分析函数 biterr 计算比特误差数和比特误差率 eyediagram 绘制眼图 scatterplot 绘制分布图 symerr 计算符号误差数和符号误差率 信源编码 compand mu律/A律压缩/扩张 dpcmdeco DPCM（差分脉冲编码调制）解码dpcmenco DPCM编码 dpcmopt 优化DPCM参数 lloyds Lloyd法则优化量化器参数 quantiz 给出量化后的级和输出值 误差控制编码 bchpoly 给出二进制BCH码的性能参数和产生多项式convenc 产生卷积码 cyclgen 产生循环码的奇偶校验阵和生成矩阵cyclpoly 产生循环码的生成多项式 decode 分组码解码器 encode 分组码编码器 gen2par 将奇偶校验阵和生成矩阵互相转换gfweight 计算线性分组码的最小距离 hammgen 产生汉明码的奇偶校验阵和生成矩阵rsdecof 对Reed-Solomon编码的ASCII文件解码rsencof 用Reed-Solomon码对ASCII文件编码rspoly 给出Reed-Solomon码的生成多项式syndtable 产生伴随解码表 vitdec 用Viterbi法则解卷积码 （误差控制编码的低级函数） bchdeco BCH解码器 bchenco BCH编码器 rsdeco Reed-Solomon解码器 rsdecode 用指数形式进行Reed-Solomon解码 rsenco Reed-Solomon编码器 rsencode 用指数形式进行Reed-Solomon编码 调制与解调
Matlab工具箱中的BP与RBF函数 Matlab神经网络工具箱中的函数非常丰富，给网络设置合适的属性，可以加快网络的学习速度，缩短网络的学习进程。限于篇幅，仅对本章所用到的函数进行介绍，其它的函数及其用法请读者参考联机文档和帮助。 1 BP与RBF网络创建函数 在Matlab工具箱中有如表1所示的创建网络的函数，作为示例，这里只介绍函数newff、newcf、newrb和newrbe。 表 1 神经网络创建函数 (1) newff函数 功能：创建一个前馈BP神经网络。 调用格式：net = newff(PR,[S1 S2...S Nl],{TF1 TF2...TF Nl},BTF,BLF,PF) 参数说明： ?PR - R个输入的最小、最大值构成的R×2矩阵； ?S i–S NI层网络第i层的神经元个数； ?TF i - 第i层的传递函数，可以是任意可导函数，默认为'tansig'，
可设置为logsig，purelin等； ?BTF -反向传播网络训练函数，默认为'trainlm'，可设置为trainbfg，trainrp，traingd等； ?BLF -反向传播权值、阈值学习函数，默认为'learngdm'； ?PF -功能函数，默认为'mse'； (2) newcf函数 功能：创建一个N层的层叠(cascade)BP网络 调用格式：net = newcf(Pr,[S1 S2...SNl],{TF1 TF2...TFNl},BTF,BLF,PF) 参数同函数newff。 (3) newrb函数 功能：创建一个径向基神经网络。径向基网络可以用来对一个函数进行逼近。newrb函数用来创建一个径向基网络，它可以是两参数网络，也可以是四参数网络。在网络的隐层添加神经元，直到网络满足指定的均方误差要求。 调用格式：net = newrb(P,T,GOAL,SPREAD) 参数说明： ?P：Q个输入向量构成的R×Q矩阵； ?T：Q个期望输出向量构成的S×Q矩阵； ?GOAL：均方误差要求，默认为0。 ?SPREAD：分散度参数，默认值为1。SPREAD越大，网络逼近的函数越平滑，但SPREAD取值过大将导致在逼近变化比较剧烈的函
MatLab函数大全 一、常用对象操作：除了一般windows窗口的常用功能键外。 1、!dir可以查看当前工作目录的文件。 !dir& 可以在dos状态下查看。 2、who可以查看当前工作空间变量名，whos可以查看变量名细节。 3、功能键： 功能键快捷键说明 方向上键Ctrl+P返回前一行输入 方向下键Ctrl+N返回下一行输入 方向左键Ctrl+B光标向后移一个字符 方向右键Ctrl+F光标向前移一个字符 Ctrl+方向右键Ctrl+R光标向右移一个字符 Ctrl+方向左键Ctrl+L光标向左移一个字符 home Ctrl+A光标移到行首 End Ctrl+E光标移到行尾 Esc Ctrl+U清除一行 Del Ctrl+D清除光标所在的字符 Backspace Ctrl+H删除光标前一个字 符Ctrl+K删除到行尾 Ctrl+C中断正在执行的命令 4、clc可以命令窗口显示的内容，但并不清除工作空间。 二、函数及运算 1、运算符： ＋：加，－：减， *：乘， /：除， \：左除 ^：幂，‘：复数的共轭转置，（）：制定运算顺序。 2、常用函数表： sin( ) 正弦（变量为弧度） Cot( ) 余切（变量为弧度）
sind( ) 正弦（变量为度数） Cotd( ) 余切（变量为度数） asin( ) 反正弦（返回弧度） acot( ) 反余切（返回弧度） Asind( ) 反正弦（返回度数） acotd( ) 反余切（返回度数） cos( ) 余弦（变量为弧度） exp( ) 指数 cosd( ) 余弦（变量为度数） log( ) 对数 acos( ) 余正弦（返回弧度） log10( ) 以10为底对数 acosd( ) 余正弦（返回度数） sqrt( ) 开方 tan( ) 正切（变量为弧度） realsqrt( ) 返回非负根 tand( ) 正切（变量为度数） abs( ) 取绝对值 atan( ) 反正切（返回弧度） angle( ) 返回复数的相位角 atand( ) 反正切（返回度数） mod(x,y) 返回x/y的余数 sum( ) 向量元素求和 3、其余函数可以用help elfun和help specfun命令获得。 4、常用常数的值： pi3.1415926……. realmin最小浮点数，2^-1022 i虚数单位 realmax最大浮点数，（2－eps）2^1022
MATLAB主要命令汇总 MATLAB函数参考 附录1.1 管理用命令 函数名功能描述函数名功能描述 addpath 增加一条搜索路径 rmpath 删除一条搜索路径 demo 运行Matlab演示程序 type 列出.M文件 doc 装入超文本文档 version 显示Matlab的版本号 help 启动联机帮助 what 列出当前目录下的有关文件 lasterr 显示最后一条信息 whatsnew 显示Matlab的新特性 lookfor 搜索关键词的帮助 which 造出函数与文件所在的目录 path 设置或查询Matlab路径 附录1.2管理变量与工作空间用命令 函数名功能描述函数名功能描述 clear 删除内存中的变量与函数 pack 整理工作空间内存 disp 显示矩阵与文本 save 将工作空间中的变量存盘 length 查询向量的维数 size 查询矩阵的维数 load 从文件中装入数据 who,whos 列出工作空间中的变量名 附录1.3文件与操作系统处理命令 函数名功能描述函数名功能描述 cd 改变当前工作目录 edit 编辑.M文件 delete 删除文件 matlabroot 获得Matlab的安装根目录 diary 将Matlab运行命令存盘 tempdir 获得系统的缓存目录 dir 列出当前目录的内容 tempname 获得一个缓存(temp)文件 ! 执行操作系统命令 附录1.4窗口控制命令 函数名功能描述函数名功能描述 echo 显示文件中的Matlab中的命令 more 控制命令窗口的输出页面format 设置输出格式 附录1.5启动与退出命令 函数名功能描述函数名功能描述 matlabrc 启动主程序 quit 退出Matlab环境 startup Matlab自启动程序 附录2 运算符号与特殊字符附录 2.1运算符号与特殊字符 函数名功能描述函数名功能描述
matlab遗传算法工具箱函数及实例讲解 最近研究了一下遗传算法，因为要用遗传算法来求解多元非线性模型。还好用遗传算法的工箱予以实现了，期间也遇到了许多问题。借此与大家分享一下。 首先，我们要熟悉遗传算法的基本原理与运算流程。 基本原理：遗传算法是一种典型的启发式算法，属于非数值算法范畴。它是模拟达尔文的自然选择学说和自然界的生物进化过程的一种计算模型。它是采用简单的编码技术来表示各种复杂的结构，并通过对一组编码表示进行简单的遗传操作和优胜劣汰的自然选择来指导学习和确定搜索的方向。遗传算法的操作对象是一群二进制串（称为染色体、个体），即种群，每一个染色体都对应问题的一个解。从初始种群出发，采用基于适应度函数的选择策略在当前种群中选择个体，使用杂交和变异来产生下一代种群。如此模仿生命的进化进行不断演化，直到满足期望的终止条件。 运算流程： Step 1：对遗传算法的运行参数进行赋值。参数包括种群规模、变量个数、交叉概率、变异概率以及遗传运算的终止进化代数。 Step 2：建立区域描述器。根据轨道交通与常规公交运营协调模型的求解变量的约束条件，设置变量的取值范围。 Step 3：在Step 2的变量取值范围内，随机产生初始群体，代入适应度函数计算其适应度值。 Step 4：执行比例选择算子进行选择操作。 Step 5：按交叉概率对交叉算子执行交叉操作。 Step 6：按变异概率执行离散变异操作。 Step 7：计算Step 6得到局部最优解中每个个体的适应值，并执行最优个体保存策略。 Step 8：判断是否满足遗传运算的终止进化代数，不满足则返回Step 4，满足则输出运算结果。 其次，运用遗传算法工具箱。 运用基于Matlab的遗传算法工具箱非常方便，遗传算法工具箱里包括了我们需要的各种函数库。目前，基于Matlab的遗传算法工具箱也很多，比较流行的有英国设菲尔德大学开发的遗传算法工具箱GATBX、GAOT以及Math Works公司推出的GADS。实际上，GADS 就是大家所看到的Matlab中自带的工具箱。我在网上看到有问为什么遗传算法函数不能调用的问题，其实，主要就是因为用的工具箱不同。因为，有些人用的是GATBX带有的函数，但MATLAB自带的遗传算法工具箱是GADS，GADS当然没有GATBX里的函数，因此运行程序时会报错，当你用MATLAB来编写遗传算法代码时，要根据你所安装的工具箱来编写代码。
matlab 常用函数汇总 编程2008-07-10 21:45:20 阅读46 评论0 字号：大中小订阅matlab常用函数 图形注释 Title 图形标题 Xlabel X轴标记 Ylabel Y轴标记 Text 文本注释 Gtext 用鼠标放置文本 Grid 网格线 MATLAB编程语言 Function 增加新的函数 Eval 执行由MA TLAB表达式构成的字串 Feval 执行由字串指定的函数 Global 定义全局变量 程序控制流 If 条件执行语句 Else 与if命令配合使用 Elseif 与if命令配合使用 End For,while和if语句的结束 For 重复执行指定次数（循环） While 重复执行不定次数（循环） Break 终止循环的执行 Return 返回引用的函数 Error 显示信息并终止函数的执行 交互输入 Input 提示用户输入 Keyboard 像底稿文件一样使用键盘输入 Menu 产生由用户输入选择的菜单 Pause 等待用户响应 Uimenu 建立用户界面菜单 Uicontrol 建立用户界面控制 一般字符串函数 Strings MATLAB中有关字符串函数的说明 Abs 变字符串为数值 Setstr 变数值为字符串 Isstr 当变量为字符串时其值为真 Blanks 空串 Deblank 删除尾部的空串 Str2mat 从各个字符串中形成文本矩阵 Eval 执行由MA TLAB表达式组成的串 字符串比较 Strcmp , , , 比较字符串 Findstr 在一字符串中查找另一个子串
Upper 变字符串为大写 Lower 变字符串为小写 Isletter 当变量为字母时，其值为真 Isspace 当变量为空白字符时，其值为真 字符串与数值之间变换 Num2str 变数值为字符串 Int2str 变整数为字符串 Str2num 变字符串为数值 Sprintf 变数值为格式控制下的字符串 Sscanf 变字符串为格式控制下的数值 十进制与十六进制数之间变换 Hex2num 变十六进制为IEEE标准下的浮点数Hex2dec 变十六制数为十进制数 Dec2hex 变十进制数为十六进制数 建模 Append 追加系统动态特性 Augstate 变量状态作为输出 Blkbuild 从方框图中构造状态空间系统Cloop 系统的闭环 Connect 方框图建模 Conv 两个多项式的卷积 Destim 从增益矩阵中形成离散状态估计器Dreg 从增益矩阵中形成离散控制器和估计器Drmodel 产生随机离散模型 Estim 从增益矩阵中形成连续状态估计器Feedback 反馈系统连接 Ord2 产生二阶系统的A、B、C、D Pade 时延的Pade近似 Parallel 并行系统连接 Reg 从增益矩阵中形成连续控制器和估计器Rmodel 产生随机连续模型 Series 串行系统连接 Ssdelete 从模型中删除输入、输出或状态ssselect 从大系统中选择子系统 模型变换 C2d 变连续系统为离散系统 C2dm 利用指定方法变连续为离散系统 C2dt 带一延时变连续为离散系统 D2c 变离散为连续系统 D2cm 利用指定方法变离散为连续系统 Poly 变根值表示为多项式表示 Residue 部分分式展开 Ss2tf 变状态空间表示为传递函数表示 Ss2zp 变状态空间表示为零极点表示
Matlab绘图 强大的绘图功能是Matlab的特点之一，Matlab提供了一系列的绘图函数，用户不需要过多的考虑绘图的细节，只需要给出一些基本参数就能得到所需图形，这类函数称为高层绘图函数。此外，Matlab 还提供了直接对图形句柄进行操作的低层绘图操作。这类操作将图形的每个图形元素（如坐标轴、曲线、文字等）看做一个独立的对象，系统给每个对象分配一个句柄，可以通过句柄对该图形元素进行操作，而不影响其他部分。 本章介绍绘制二维和三维图形的高层绘图函数以及其他图形控制函数的使用方法，在此基础上，再介绍可以操作和控制各种图形对象的低层绘图操作。 一、二维绘图 二维图形是将平面坐标上的数据点连接起来的平面图形。可以采用不同的坐标系，如直角坐标、对数坐标、极坐标等。二维图形的绘制是其他绘图操作的基础。 （一）绘制二维曲线的基本函数 在Matlab中，最基本而且应用最为广泛的绘图函数为plot，利用它可以在二维平面上绘制出不同的曲线。 1．plot函数的基本用法
plot函数用于绘制二维平面上的线性坐标曲线图，要提供一组x 坐标和对应的y坐标，可以绘制分别以x和y为横、纵坐标的二维曲线。plot函数的应用格式 plot(x,y) 其中x,y为长度相同的向量，存储x坐标和y坐标。 例51 在[0 , 2pi]区间，绘制曲线 程序如下：在命令窗口中输入以下命令 >> x=0:pi/100:2*pi; >> y=2*exp(-0.5*x).*sin(2*pi*x); >> plot(x,y) 程序执行后，打开一个图形窗口，在其中绘制出如下曲线 注意：指数函数和正弦函数之间要用点乘运算，因为二者是向量。 例52 绘制曲线 这是以参数形式给出的曲线方程，只要给定参数向量，再分别求出x,y向量即可输出曲线： >> t=-pi:pi/100:pi; >> x=t.*cos(3*t); >> y=t.*sin(t).*sin(t); >> plot(x,y) 程序执行后，打开一个图形窗口，在其中绘制出如下曲线 以上提到plot函数的自变量x,y为长度相同的向量，这是最常见、最基本的用法。实际应用中还有一些变化。
M文件与M函数 Matlab输入命令的常用方式有两种：一种是直接在Matlab的命令窗门中逐条输入Matlab 命令；二是m文件工作方式。当命令行很简单时，使用逐条输入方式还是比较方便的。但当命令行很多时(比如说几十行乃至全成百上千行命令)，显然再使用这种方式输入MATLAB 命令，就会显得杂乱无章，不易于把握程序的具体走向，并且给程序的修改和维护带来了很大的麻烦。这时，建议采用Matlab命令的第二种输入形式m文件工作方式。 m文件工作方式，指的是将要执行的命令全部写在一个文本文件中，这样既能使程序显得简洁明了，又便于对程序的修改与维护。m文件直接采用Matlab命令编写，就像在Matlab 的命令窗口直接输入命令一样，因此调试起来也十分方便，并且增强了程序的交互性。 m文件与其他文本文件一样，可以在任何文本编辑器中进打编辑、存储、修改和读取。利用m文件还可以根据白己的需要编写一些函数，这些函数也可以橡Matlab提供的函数一样进行调用。从某种意义上说，这也是对MATLAB的二次开发。 m文件有两种形式：一种是命令方式或称脚本方式；另一种就是函数文件形式。两种形式的文件扩展名均是.m。 1、M文件 当遇到输入命令较多以及要重复输入命令的情况时，利用命令文件就显得很方便了。将所有要执行的命令按顺序放到一个扩展名为.m的文本文件中，每次运行时只需在MATLAB 的命令窗口输入m文件的文件名就可以了。需要注意的是，m文件最好直接放在Matlab 的默认搜索路径下(一般是Matlab安装目录的子目录work中)，这样就不用设置m文件的路径了，否则应当用路径操作指令path重新设置路径。另外，m文件名不应该与Matlab的内置函数名以及工具箱中的函数重名，以免发生执行错误命令的现象。Matlab对命令文件的执行等价于从命令窗口中顺序执行文件中的所有指令。命令文件可以访问Matlab工作空间里的任何变量及数据。命令文件运行过程中产生的所有变量都等价于从Matlab工作空间中创建这些变量。因此，任何其他命令文件和函数都可以自由地访问这些变量。这些变量一旦产生就一直保存在内存中，只有对它们重新赋值，它们的原有值才会变化。关机后，这里变量也就全部消失了。另外，在命令窗口中运行clear命令，也可以把这些变量从工作空间中删去。当然，在Matlab的工作空间窗口中也可以用鼠标选择想要删除的变量，从而将这些变量从工作空间中删除。 接下来，编写一个名为test.m的命令文件，用来计算矩阵1到100的和，并把它放到变量s中。 第一步创建新的M-文件。在Matlab主菜单上选择菜单命令File→New→M-File
附录1 常用命令 附录1.1 管理用命令 函数名功能描述函数名功能描述 addpath 增加一条搜索路径rmpath 删除一条搜索路径 demo 运行Matlab演示程序type 列出.M文件 doc 装入超文本文档version 显示Matlab的版本号 help 启动联机帮助what 列出当前目录下的有关文件 lasterr 显示最后一条信息whatsnew 显示Matlab的新特性 lookfor 搜索关键词的帮助which 造出函数与文件所在的目录 path 设置或查询Matlab路径 附录1.2管理变量与工作空间用命令 函数名功能描述函数名功能描述 clear 删除内存中的变量与函数p ack 整理工作空间内存 disp 显示矩阵与文本save 将工作空间中的变量存盘 length 查询向量的维数size 查询矩阵的维数 load 从文件中装入数据who,whos 列出工作空间中的变量名 附录1.3文件与操作系统处理命令 函数名功能描述函数名功能描述 cd 改变当前工作目录edit 编辑.M文件 delete 删除文件matlabroot 获得Matlab的安装根目录 diary 将Matlab运行命令存盘t empdir 获得系统的缓存目录 dir 列出当前目录的内容tempname 获得一个缓存(temp)文件 ! 执行操作系统命令 附录1.4窗口控制命令 函数名功能描述函数名功能描述 echo 显示文件中的Matlab中的命令m ore 控制命令窗口的输出页面 format 设置输出格式 附录1.5启动与退出命令 函数名功能描述函数名功能描述 matlabrc 启动主程序quit 退出Matlab环境 startup Matlab自启动程序 附录2 运算符号与特殊字符附录 2.1运算符号与特殊字符 函数名功能描述函数名功能描述 + 加... 续行标志 - 减, 分行符(该行结果不显示) * 矩阵乘; 分行符(该行结果显示) .* 向量乘% 注释标志 ^ 矩阵乘方! 操作系统命令提示符
图像增强 1. 直方图均衡化的Matlab 实现 1.1 imhist 函数 功能：计算和显示图像的色彩直方图 格式：imhist(I,n) imhist(X,map) 说明：imhist(I,n) 其中，n 为指定的灰度级数目，缺省值为256；imhist(X, map) 就算和显示索引色图像X 的直方图，map 为调色板。用stem(x,coun ts) 同样可以显示直方图。 1.2 imcontour 函数 功能：显示图像的等灰度值图 格式：imcontour(I,n),imcontour(I,v) 说明：n 为灰度级的个数，v 是有用户指定所选的等灰度级向量。 1.3 imadjust 函数 功能：通过直方图变换调整对比度 格式：J=imadjust(I,[low high],[bottom top],gamma) newmap=imadjust(map,[low high],[bottom top],gamma) 说明：J=imadjust(I,[low high],[bottom top],gamma) 其中，gamma 为校正量r，[low high] 为原图像中要变换的灰度范围，[bottom top] 指定了变换后的灰度范围；newmap=imadjust(map,[low high],[bottom top],gamm a) 调整索引色图像的调色板map 。此时若[low high] 和[bottom top] 都为2×3的矩阵，则分别调整R、G、B 3个分量。
1.4 histeq 函数 功能：直方图均衡化 格式：J=histeq(I,hgram) J=histeq(I,n) [J,T]=histeq(I,...) newmap=histeq(X,map,hgram) newmap=histeq(X,map) [new,T]=histeq(X,...) 说明：J=histeq(I,hgram) 实现了所谓“直方图规定化”，即将原是图象I 的直方图变换成用户指定的向量hgram 。hgram 中的每一个元素都在[0,1] 中；J=histeq(I,n) 指定均衡化后的灰度级数n ，缺省值为64；[J,T]=histeq(I,...)返回从能将图像I 的灰度直方图变换成图像J 的直方图的变换T ；newma p=histeq(X,map) 和[new,T]=histeq(X,...) 是针对索引色图像调色板的直方图均衡。 2. 噪声及其噪声的Matlab 实现 imnoise 函数 格式：J=imnoise(I,type) J=imnoise(I,type,parameter) 说明：J=imnoise(I,type) 返回对图像I 添加典型噪声后的有噪图像J ，参数type 和parameter 用于确定噪声的类型和相应的参数。 3. 图像滤波的Matlab 实现 3.1 conv2 函数 功能：计算二维卷积
1-1、基本运算与函数 在MATLAB下进行基本数学运算，只需将运算式直接打入提示号（>>）之後，并按入Enter 键即可。例如： >> (5*2+1.3-0.8)*10/25 ans =4.2000 MATLAB会将运算结果直接存入一变数ans，代表MATLAB运算後的答案（Answer）并显示其数值於萤幕上。 小提示： ">>"是MATLAB的提示符号（Prompt），但在PC中文视窗系统下，由於编码方式不同，此提示符号常会消失不见，但这并不会影响到MATLAB的运算结果。 我们也可将上述运算式的结果设定给另一个变数x： x = (5*2+1.3-0.8)*10^2/25 x = 42 此时MATLAB会直接显示x的值。由上例可知，MATLAB认识所有一般常用到的加（+）、减（-）、乘（*）、除（/）的数学运算符号，以及幂次运算（^）。 小提示： MATLAB将所有变数均存成double的形式，所以不需经过变数宣告（Variable declaration）。MATLAB同时也会自动进行记忆体的使用和回收，而不必像C语言,必须由使用者一一指定.这些功能使的MATLAB易学易用，使用者可专心致力於撰写程式，而不必被软体枝节问题所干扰。 若不想让MATLAB每次都显示运算结果，只需在运算式最後加上分号（；）即可，如下例：y = sin(10)*exp(-0.3*4^2); 若要显示变数y的值，直接键入y即可： >>y y =-0.0045 在上例中，sin是正弦函数，exp是指数函数，这些都是MATLAB常用到的数学函数。 下表即为MATLAB常用的基本数学函数及三角函数： 小整理：MATLAB常用的基本数学函数 abs(x)：纯量的绝对值或向量的长度 angle(z)：复数z的相角(Phase angle) sqrt(x)：开平方 real(z)：复数z的实部 imag(z)：复数z的虚部 conj(z)：复数z的共轭复数 round(x)：四舍五入至最近整数 fix(x)：无论正负，舍去小数至最近整数 floor(x)：地板函数，即舍去正小数至最近整数 ceil(x)：天花板函数，即加入正小数至最近整数 rat(x)：将实数x化为分数表示 rats(x)：将实数x化为多项分数展开 sign(x)：符号函数 (Signum function)。 当x<0时，sign(x)=-1； 当x=0时，sign(x)=0; 当x>0时，sign(x)=1。 > 小整理：MATLAB常用的三角函数 sin(x)：正弦函数
A a abs绝对值、模、字符的ASCII码值acos反余弦 acosh反双曲余弦 acot反余切 acoth反双曲余切 acsc反余割 acsch反双曲余割 align启动图形对象几何位置排列工具all所有元素非零为真 angle相角 ans表达式计算结果的缺省变量名any所有元素非全零为真 area面域图 argnames函数M文件宗量名 asec反正割 asech反双曲正割 asin反正弦 asinh反双曲正弦 assignin向变量赋值 atan反正切 atan2四象限反正切 atanh反双曲正切 autumn红黄调秋色图阵 axes创建轴对象的低层指令 axis控制轴刻度和风格的高层指令 B b bar二维直方图 bar3三维直方图 bar3h三维水平直方图 barh二维水平直方图 base2dec X进制转换为十进制 bin2dec二进制转换为十进制blanks创建空格串 bone蓝色调黑白色图阵 box框状坐标轴 break while或for环中断指令brighten亮度控制 C c
capture（3版以前）捕获当前图形 cart2pol直角坐标变为极或柱坐标 cart2sph直角坐标变为球坐标 cat串接成高维数组 caxis色标尺刻度 cd指定当前目录 cdedit启动用户菜单、控件回调函数设计工具cdf2rdf复数特征值对角阵转为实数块对角阵ceil向正无穷取整 cell创建元胞数组 cell2struct元胞数组转换为构架数组 celldisp显示元胞数组内容 cellplot元胞数组内部结构图示 char把数值、符号、内联类转换为字符对象chi2cdf分布累计概率函数 chi2inv分布逆累计概率函数 chi2pdf分布概率密度函数 chi2rnd分布随机数发生器 chol Cholesky分解 clabel等位线标识 cla清除当前轴 class获知对象类别或创建对象 clc清除指令窗 clear清除内存变量和函数 clf清除图对象 clock时钟 colorcube三浓淡多彩交叉色图矩阵 colordef设置色彩缺省值 colormap色图 colspace列空间的基 close关闭指定窗口 colperm列排序置换向量 comet彗星状轨迹图 comet3三维彗星轨迹图 compass射线图 compose求复合函数 cond（逆）条件数 condeig计算特征值、特征向量同时给出条件数condest范-1条件数估计 conj复数共轭 contour等位线 contourf填色等位线 contour3三维等位线
MATLAB函数大全 代充全国移动、联通、电信话费、腾讯QQ业务、网游点卡 淘宝店址：https://www.doczj.com/doc/6e4411809.html,/ 信誉至上，服务第一 A a abs 绝对值、模、字符的ASCII码值 acos 反余弦 acosh 反双曲余弦 acot 反余切 acoth 反双曲余切 acsc 反余割 acsch 反双曲余割 align 启动图形对象几何位置排列工具 all 所有元素非零为真 angle 相角 ans 表达式计算结果的缺省变量名 any 所有元素非全零为真 area 面域图 argnames 函数M文件宗量名 asec 反正割 asech 反双曲正割 asin 反正弦 asinh 反双曲正弦 assignin 向变量赋值 atan 反正切 atan2 四象限反正切 atanh 反双曲正切 autumn 红黄调秋色图阵 axes 创建轴对象的低层指令 axis 控制轴刻度和风格的高层指令 B b bar 二维直方图 bar3 三维直方图 bar3h 三维水平直方图 barh 二维水平直方图 base2dec X进制转换为十进制
bin2dec 二进制转换为十进制 blanks 创建空格串 代充全国移动、联通、电信话费、腾讯QQ业务、网游点卡 淘宝店址：https://www.doczj.com/doc/6e4411809.html,/ 信誉至上，服务第一 bone 蓝色调黑白色图阵 box 框状坐标轴 break while 或for 环中断指令 brighten 亮度控制 C c capture （3版以前）捕获当前图形 cart2pol 直角坐标变为极或柱坐标 cart2sph 直角坐标变为球坐标 cat 串接成高维数组 caxis 色标尺刻度 cd 指定当前目录 cdedit 启动用户菜单、控件回调函数设计工具 cdf2rdf 复数特征值对角阵转为实数块对角阵 ceil 向正无穷取整 cell 创建元胞数组 cell2struct 元胞数组转换为构架数组 celldisp 显示元胞数组内容 cellplot 元胞数组内部结构图示 char 把数值、符号、内联类转换为字符对象 chi2cdf 分布累计概率函数 chi2inv 分布逆累计概率函数 chi2pdf 分布概率密度函数 chi2rnd 分布随机数发生器 chol Cholesky分解 clabel 等位线标识 cla 清除当前轴 class 获知对象类别或创建对象 clc 清除指令窗 clear 清除内存变量和函数 clf 清除图对象 clock 时钟 colorcube 三浓淡多彩交叉色图矩阵
Matlab自定义函数的五种方法 [转] n 1、函数文件+调用命令文件：需单独定义一个自定义函数的M文件; n 2、函数文件+子函数：定义一个具有多个自定义函数的M文件； n 3、Inline:无需M文件，直接定义； n 4、Syms+subs: 无需M文件,直接定义； n 5、字符串+subs：无需M文件,直接定义. 1、函数文件+调用函数文件：定义多个M文件： % 调用函数文件:myfile.m clear clc for t=1:10 y=mylfg(t); fprintf(‘%4d^(1/3)=%6.4f\n’,t,y); end %自定义函数文件: mylfg.m function y=mylfg(x) %注意：函数名（mylfg）必须与文件名（mylfg.m）一致 Y=x^(1/3); 注：这种方法要求自定义函数必须单独写一个M文件，不能与调用的命令文件写在同一个M文件中。 2、函数文件+子函数：定义一个具有多个子函数的M 文件
%命令文件：funtry2.m function []=funtry2() for t=1:10 y=lfg2(t) fprintf(‘%4d^(1/3)=%6.4f\n’); End function y=lfg2(x) Y= x^(1/3); %注：自定义函数文件funtry2.m中可以定义多个子函数function。子函数lfg2只能被主函数和主函数中的其他子函数调用。 3、Inline:无需M文件，直接定义； %inline命令用来定义一个内联函数：f=inline(‘函数表达式’, ‘变量1’,’变量2’,……)。 调用方式：y=f(数值列表) %注意：代入的数值列表顺序应与inline()定义的变量名顺序一致。 例如： f=inline(‘x^2+y’,’x’,’y’); z=f(2,3) Ans=7 注：这种函数定义方式是将它作为一个内部函数调用。特点是，它是基于Matlab 的数值运算内核的，所以它的运算速度较快，程序效率更高。缺点是，该方法只能对数值进行代入，不支持符号代入，且对定义后的函数不能进行求导等符号运算。 例： Clear Clc
MATLAB函数大全 Matlab有没有求矩阵行数/列数/维数的函数？ ndims(A)返回A的维数 size(A)返回A各个维的最大元素个数 length(A)返回max(size(A)) [m,n]=size(A)如果A是二维数组，返回行数和列数nnz(A)返回A中非0元素的个数 MATLAB的取整函数:fix(x), floor(x) :,ceil(x) , round(x) (1)fix(x) : 截尾取整. >> fix( [3.12 -3.12]) ans = 3 -3 (2)floor(x):不超过x 的最大整数.(高斯取整) >> floor( [3.12 -3.12]) ans =
3 -4 (3)ceil(x) : 大于x 的最小整数>> ceil( [3.12 -3.12]) ans = 4 -3 (4)四舍五入取整 >> round(3.12 -3.12) ans = >> round([3.12 -3.12]) ans =
3 -3 >> 如何用matlab生成随机数函数 rand(1) rand(n):生成0到1之间的n阶随机数方阵rand(m,n):生成0到1之间的m×n的随机数矩阵(现成的函数) 另外: Matlab随机数生成函数 betarnd 贝塔分布的随机数生成器 binornd 二项分布的随机数生成器 chi2rnd 卡方分布的随机数生成器 exprnd 指数分布的随机数生成器 frnd f分布的随机数生成器 gamrnd 伽玛分布的随机数生成器 geornd 几何分布的随机数生成器 hygernd 超几何分布的随机数生成器
告: Matlab中插值函数汇总和使用说明收藏 命令1 interp1 功能一维数据插值(表格查找。该命令对数据点之间计算内插值。它找出一元函数f(x在中间点的数值。其中函数f(x由所给数据决定。x:原始数据点 Y:原始数据点 xi:插值点 Yi:插值点 格式 (1yi = interp1(x,Y,xi 返回插值向量yi,每一元素对应于参量xi,同时由向量x 与Y 的内插值决定。参量x 指定数据Y 的点。 若Y 为一矩阵,则按Y 的每列计算。yi 是阶数为length(xi*size(Y,2的输出矩阵。 (2yi = interp1(Y,xi 假定x=1:N,其中N 为向量Y 的长度,或者为矩阵Y 的行数。 (3yi = interp1(x,Y,xi,method 用指定的算法计算插值: ’nearest’:最近邻点插值,直接完成计算; ’linear’:线性插值(缺省方式,直接完成计算;
’spline’:三次样条函数插值。对于该方法,命令interp1 调用函数spline、ppval、mkpp、umkpp。这些命令生成一系列用于分段多项式操作的函 数。命令spline 用它们执行三次样条函数插值; ’pchip’:分段三次Hermite 插值。对于该方法,命令interp1 调用函数p chip,用于对向量x 与y 执行分段三次内插值。该方法保留单调性与数据的外形; ’cubic’:与’pchip’操作相同; ’v5cubic’:在MATLAB 5.0 中的三次插值。 对于超出x 范围的xi 的分量,使用方法’nearest’、’linear’、’v5cubic’的插值算法,相应地将返回NaN。对其他的方法,interp1 将对超出的分量执行外插值算法。 (4yi = interp1(x,Y,xi,method,'extrap' 对于超出x 范围的xi 中的分量将执行特殊的外插值法extrap。 (5yi = interp1(x,Y,xi,method,extrapval 确定超出x 范围的xi 中的分量的外插值extrapval,其值通常取NaN 或0。 例1 1.>>x = 0:10; y = x.*sin(x; 2.>>xx = 0:.25:10; yy = interp1(x,y,xx; 3.>>plot(x,y,'kd',xx,yy 复制代码 例2 1.>> year = 1900:10:2010;