常用matlab语句

setdiff：集合差集
setxor：集合异或(不在交集中的元素)
union：两个集合的并
unique：返回向量作为一个集合所有元素(去掉相同元素)
14、代数式因式分解，展开，极限，倒数，级数求和，泰勒级数展开
>>syms x y a b c t
>>E1 = expand((x-2)*(x-4)*(y-t))//因式展开
3、将向量坐标在matlab里作图：
x=xlsread('C:\Documents and Settings\Administrator\桌面\school\GIS\新建Microsoft Excel工作表.xls',1,'A1:a1000');
y=xlsread('C:\Documents and Settings\Administrator\桌面\school\GIS\新建Microsoft Excel工作表.xls',1,'b1:b1000');
2 3 4 6 7 8
2 3 5 6 7 8
2 4 5 6 7 8
3 4 5 6 7 8
19、设ｎ个人中每个人的生日在一年365天中任一天是等可能的。 求当n为23，40，64时，这ｎ个人中至少有两人生日相同的概率各为多少？
>> n=23;
>> p=1-nchoosek(365,n)*factorial(n)/365^n
Title(‘plot of sin(x)/x’);
Grid on;
fplot与ezplot相类似，但更加精确。前两个参数与函数ezplot中的相同，但是数fplot还有其他优点
1.函数fplot是适应性的，它意味着在自变量范围内函数突然变化显示更多的点。
2.函数fplot支持TEX命令，用来指定坐标图的标题和坐标轴标签，而函数ezplot不能。
c 2x3x296 double//输出两个2x3的数组
>> c
c(:,:,1) =
1 2 3
4 5 6
c(:,:,2) =
12 5 8
9 9 2
8、多维数组的存储方式以列为主导，先存a(1,1)，a(2,1)，……a(1,2)，a(2,2)，……
对于多维数组，数组的第一个下标增长最快，第二个仅次之，一次类推，最后一个变化最慢。例如一个2x2x2的数组中，他的元素内存分配顺序是（1，1，1），（2，1，1），（1，2，1），（2，2，1），（1，1，2），（2，1，2），（1，2，2），
在一般情况下，在画函数图象时，你应当使用函数fplot。
34、pcacov
功能：运用协方差矩阵进行主成分分析
格式：PC=pcacov(X)
[PC,latent,explained]=pcacov(X)
说明：[PC,latent,explained]=pcacov(X)通过协方差矩阵X进行主成分分析，返回主成分(PC)、协方差矩阵X的特征值(latent)和每个特征向量表征在观测量总方差中所占的百分数(explained)。
18、取出1到8,8个数中，6个数的所有组合
combntns(1:8,6)
ans =
1 2 3 4 5 6
1 2 3 4 5 7
1 2 3 4 5 8
1 2 3 4 6 7
1 2 3 4 6 8
1 2 3 4 7 8
1 2 3 5 6 7
1 2 3 5 6 8
1 2 3 5 7 8
1 2 3 6 7 8
end
以上代码效果一样，但是矩阵数组运算比for几乎快15倍还多，而需花更多的内存供中间创建数组使用。
30、tic函数复位内建计时器，而toc函数则从最后一次调用tic以秒开始计时。
>>tic;
>> t=toc
t =
2.8483e-004
31、逻辑数组
>> a=[1,2,3;4,5,6;7,8,9];
>> image(a);
29、for ii = 1:100 ii=1:100
Square(ii)=ii^2; square=ii.^2;
Square_root(ii)=ii^(1/2); square_root =ii.^(1/2);
Cube_root（ii）=ii^(1/3); cube_root =ii.^(1/3);
22、程序运行中用ctrl+c强制结束。
23、查找正在执行文件版本和路径：
>>which cross
D:\MATLAB7\toolbox\MATLAB\specfun\corss.m
24、用zeros（）创建0矩阵：
Zeros(2),zeros(2,3),
C=[1 2;3 4],d=zeros(size(c));
运行结果：
p =
0.5073
>> n=40;
>> p=1-nchoosek(365,n)*factorial(n)/365^n
运行结果：
p =
0.8912
>> n=64;
>>p=1-nchoosek(365,n)*factorial(n)/365^n
运行结果：
p =
0.99ห้องสมุดไป่ตู้2
20、
21、如果程序语句过长，需用…连接，跳行继续编写。
>>E2 = expand(cos(x+y))
>>E3 = expand(exp((a+b)^3))
>>E4 = expand(log(a*b/sqrt(c)))
>>E5 = expand([sin(2*t), cos(2*t)])
>>factor(x^2*y-x^2*t-6*x*y+6*x*t+8*y-8*t)//因式分解
plot(x,y);
//若要在同一图中作多个图，需成对用‘，’隔开即可plot(x,y1,x,y2)
title('函数图像');%设置图像标题
xlabel('X坐标');%设置图像X坐标名称
ylabel('Y坐标');%设置图像Y坐标名称
4、在matlab中寻找自己想找的函数：help(确定的情况下)和lookfor（在不确定函数名的情况下，如下图）
the value of pi is 3.14159>>//%6.5f表示以6个字符，小数点后5位小数的浮点型
但是，fprintf不能显示虚数部分，忽略虚数部分。
13、matlab里有关于集合运算的函数有如下的一些:
intersect：集合交集
ismember：是否集合中元素
issorted：判断集合是否排序(相关函数，sort：排序)
格式B = repmat(A,m,n) %将矩阵A复制m×n块，即B由m×n块A平铺而成。
B = repmat(A,[m n]) %与上面一致
B = repmat(A,[m n p…]) %B由m×n×p×…个A块平铺而成
repmat(A,m,n) %当A是一个数a时，该命令产生一个全由a组成的m×n矩阵
同理ones( ),eye( );1矩阵和单位矩阵
Length（arr）用于返回向量的长度和二维数组最长的那一维数；
Size(arr)返回指定数组的行数和列数。
25、最大值函数，index为第几位
[max index]=max([4 2 4 5])
max =
5
index =
4
26、try
Statement 1
enter am inputdate:1,2,3,5,97//略，如果要输入矩阵，则需要用[ ]括
x =1,2,3,5,97//起来
7、>> c(:,:,1)=[1,2,3;4,5,6];
>> c(:,:,2)=[12,5,8;9,9,2];
>> whos
Name Size Bytes Class Attributes
Statement 2
...
catch
Statement 1
Statement 2
...
end
当程序运行到try/catch语句块，在try语句块中的一些语句将会被执行。如果没有错误出现，catch语句块将会被跳过。另一方面，如果错误发生在一个try语句块，那么程序将中止执行try语句块，并立即执行catch语句块。
>> b=a>5
b =
0 0 0
0 0 1
1 1 1
>>a(b)=sqrt(a(b))或者a(~b)类似
a =
1.0000 2.0000 3.0000
4.0000 5.0000 2.4495
2.6458 2.8284 3.0000
32、计算括号里面的值
a=eval('4^2')
a =
16；
33、ezplot(‘sin(x)/x’,[-4*pi*pi]);
（2，2，2）。
另外：matlaB还允许以一维下标的形式访问多维数组：
如：a=[1 2 3;4 5 6;7 7 8];中，a(5)=5。
End函数在数据提取中的应用：a=[1 2 3 4;5 6 7 8;9 10 11 12],则a(2：end,2：end)将会返回[6 7 8;10 11 12]；而其中第一个end为行数，第二个end为列数。
（注：help命令比lookfor要快得多）
5、搜索函数文件：>> which cross