ans =为3-1次展开taylor级数
76+2*x+(x-1)^2
15、数组乘法A.*B(对应元素相乘),数组右除Ahttps://m.360docs.net/doc/B及数组左除A.B(数组对应相除),数组指数运算A.^B(A(i,j)^B(i,j))。
16、阶乘计算
factorial(10)%10的阶乘
ans =
3628800
17、计算组合数的命令nchoosek和求所有组合的命令combntns
运行结果:
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.9972
20、
21、如果程序语句过长,需用…连接,跳行继续编写。
Statement 2
...
catch
Statement 1
Statement 2
...
end
当程序运行到try/catch语句块,在try语句块中的一些语句将会被执行。如果没有错误出现,catch语句块将会被跳过。另一方面,如果错误发生在一个try语句块,那么程序将中止执行try语句块,并立即执行catch语句块。
(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为列数。
同理ones( ),eye( );1矩阵和单位矩阵
Length(arr)用于返回向量的长度和二维数组最长的那一维数;
Size(arr)返回指定数组的行数和列数。
25、最大值函数,index为第几位
[max index]=max([4 2 4 5])
max =
5
index =
4
26、try
Statement 1
>> 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);
在一般情况下,在画函数图象时,你应当使用函数fplot。
34、pcacov
功能:运用协方差矩阵进行主成分分析
格式:PC=pcacov(X)
[PC,latent,explained]=pcacov(X)
说明:[PC,latent,explained]=pcacov(X)通过协方差矩阵X进行主成分分析,返回主成分(PC)、协方差矩阵X的特征值(latent)和每个特征向量表征在观测量总方差中所占的百分数(explained)。
3、将向量坐标在matlab里作图:
x=xlsread('C:Documents and SettingsAdministrator桌面schoolGIS新建Microsoft Excel工作表.xls',1,'A1:a1000');
y=xlsread('C:Documents and SettingsAdministrator桌面schoolGIS新建Microsoft Excel工作表.xls',1,'b1:b1000');
x=-pi:pi/20:pi;
y=sin(x);
plot(x,y);
title('subplot 1 title');
subplot(2,1,2);
x=-pi:pi/20:pi;
y=cos(x);
plot(x,y);
title('subplot2 title');
28、>> a=imread('F:桌面图片2011-07-04-0.bmp');
1 2 4 5 6 7
1 2 4 5 6 8
1 2 4 5 7 8
1 2 4 6 7 8
1 2 5 6 7 8
1 3 4 5 6 7
1 3 4 5 6 8
1 3 4 5 7 8
1 3 4 6 7 8
1 3 5 6 7 8
1 4 5 6 7 8
2 3 4 5 6 7
2 3 4 5 6 8
2 3 4 5 7 8
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),
Title(‘plot of sin(x)/x’);
Grid on;
fplot与ezplot相类似,但更加精确。前两个参数与函数ezplot中的相同,但是数fplot还有其他优点
1.函数fplot是适应性的,它意味着在自变量范围内函数突然变化显示更多的点。
2.函数fplot支持TEX命令,用来指定坐标图的标题和坐标轴标签,而函数ezplot不能。
(注:help命令比lookfor要快得多)
5、搜索函数文件:>> which cross
D:MATLAB7 oolboxMATLABspecfuncross.m
6、关键字input用来提示直接从键盘输入初始化变量。
x=input('enter am inputdate:','s')//s表示以字符串形式输入,无需是省
>> date
ans =10-Mar-2011
10、
11、disp显示函数
>> str=['the value of pi=' num2str(pi)];//num2str字符到字符串
>> disp(str)
the value of pi=3.1416
12、fprintf函数格式化输出数据
>> fprintf('the value of pi is%6.5f',pi)
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、设n个人中每个人的生日在一年365天中任一天是等可能的。 求当n为23,40,64时,这n个人中至少有两人生日相同的概率各为多少?
>> n=23;
>> p=1-nchoosek(365,n)*factorial(n)/365^n
enter am inputdate:1,2,3,5,97//略,如果要输入矩阵,则需要用[ ]括
x =Biblioteka Baidu,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
格式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矩阵
27、figure函数同时绘多个图像
>> figure(1);
>> x=1:0.5:3;
>> y1=exp(x);
>> plot(x,y1);
>> figure(2);
>> y2=exp(-x);
>> plot(x,y2);
figure(1);
subplot(2,1,1);%2*1个窗口中第一个窗口
plot(x,y);
//若要在同一图中作多个图,需成对用‘,’隔开即可plot(x,y1,x,y2)
title('函数图像');%设置图像标题
xlabel('X坐标');%设置图像X坐标名称
ylabel('Y坐标');%设置图像Y坐标名称
4、在matlab中寻找自己想找的函数:help(确定的情况下)和lookfor(在不确定函数名的情况下,如下图)
>>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)//因式分解
setdiff:集合差集
setxor:集合异或(不在交集中的元素)
union:两个集合的并
unique:返回向量作为一个集合所有元素(去掉相同元素)
14、代数式因式分解,展开,极限,倒数,级数求和,泰勒级数展开
>>syms x y a b c t
>>E1 = expand((x-2)*(x-4)*(y-t))//因式展开
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];
>>collect(a^2+2.3*a^2,'a')//合并同类项
>> limit(x^2,x,0,’right’or’left’)//表示x^2在0点对x求左or右极限,x,0,right可默认为0和左右极限相等
syms k
>> symsum(k^2,1,k-1)//k^2级数1到k-1项求和
>> taylor(x^2+77,3,'x',1)//将x^2+77按x在1处,以最高次数
9、Inf表示无穷大,一般是除以0产生的;
NaN表示没有这个数,一般由数学运算得到,如0除以0
Clock为包含年,月,日,时,分,秒的特殊变量,是一个6元素行向量;
>> clock
ans = 1.0e+003 *
2.0110 0.0030 0.0100 0.0210 0.0520 0.0496
Date输出当前日期,如:
22、程序运行中用ctrl+c强制结束。
23、查找正在执行文件版本和路径:
>>which cross
D:MATLAB7 oolboxMATLABspecfuncorss.m
24、用zeros()创建0矩阵:
Zeros(2),zeros(2,3),
C=[1 2;3 4],d=zeros(size(c));
>> 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]);
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
1、将txt文件数据导入matlab里面:
textread('D:MATLAB7worka.txt')或者用load('D:MATLAB7worka.txt')
2、将excel文件数据导入matlab里面:
x=xlsread('a.xls','sheet1','a1:a10')//excel文件名是a.xls,sheet1是表名,a1:a10是a列的1到10行
the value of pi is 3.14159>>//%6.5f表示以6个字符,小数点后5位小数的浮点型
但是,fprintf不能显示虚数部分,忽略虚数部分。
13、matlab里有关于集合运算的函数有如下的一些:
intersect:集合交集
ismember:是否集合中元素
issorted:判断集合是否排序(相关函数,sort:排序)
版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 举报,一经查实,本站将立刻删除。
如需转载请保留出处:https://bianchenghao.cn/bian-cheng-ri-ji/33880.html