四川农业大学数学实验Matlab

1、数学实验Matlab一、 线性代数1. 函数ones语法：ones(n)生成n行n列的全为1的方阵Ones(m,n) 生成m行n列的全为1的矩阵>> ones(3)ans = 1 1 1 1 1 1 1 1 1>> ones(3,2)ans = 1 1 1 1 1 12函数zeros()语法：zeros(n)生成n行n列的全为零的方阵 Zeros(m,n)生成m行n列的全为零的矩阵>> zeros(4)ans = 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0>> zeros(4,5)ans = 0 0 0 0 0 0 0 0 0

2、 0 0 0 0 0 0 0 0 0 0 03.函数magic语法：magic(n)生成n行n列的魔方阵>> magic(4)ans = 16 2 3 13 5 11 10 8 9 7 6 12 4 14 15 14.函数:eye语法：eye(n)生成n行n列的单位阵>> eye(5)ans = 1 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 15.函数：pascal语法：生成n行n列的PASCAL矩阵>> pascal(5)ans = 1 1 1 1 1 1 2 3 4 5 1 3 6 10 15 1 4 1

3、0 20 35 1 5 15 35 706.函数rand语法：rand(n)返回一个N*N的方阵，每一个值都是0到1之间的随机数 Rand(m,n)返回一个M*N的矩阵，每一个值都是0到1之间的随机数>> rand(2)ans = 0.9501 0.6068 0.2311 0.4860>> rand(3,2)ans = 0.8913 0.0185 0.7621 0.82140.4565 0.44477.函数diag语法：diag(A)求矩阵A的对角线 Diag(V)创建一个对角线为向量V的矩阵>> A=1 2 3;4 5 6;7 8 9;>> d

4、iag(A)ans = 1 5 9>> V=4 3 8;>> diag(V)ans = 4 0 0 0 3 0 0 0 88.函数randn语法：randn(n)返回一个N*N的方阵，每一个数是一个高斯（或正态）随机数，均值为0，方差为1Randn(m,n)返回一个M*N的矩阵，每一个数是一个高斯（或正态）随机数，均值为0，方差为1>> randn(2)ans = -0.4326 0.1253 -1.6656 0.2877>> randn(3,2)ans = -1.1465 -0.0376 1.1909 0.32731.1892 0.17469.

5、函数linspace语法：linspace(a,b,n)创建一个从a到b（包含a,b）的有N个等间距的元素的向量Linspace(a,b)创建一个从a到b（包含a,b）的有100个等间距的元素的向量>> linspace(1,2,5)ans = 1.0000 1.2500 1.5000 1.7500 2.0000>> linspace(1,2)ans = Columns 1 through 12 1.0000 1.0101 1.0202 1.0303 1.0404 1.0505 1.0606 1.0707 1.0808 1.0909 1.1010 1.1111 Colu

6、mns 13 through 24 1.1212 1.1313 1.1414 1.1515 1.1616 1.1717 1.1818 1.1919 1.2020 1.2121 1.2222 1.2323 Columns 25 through 36 1.2424 1.2525 1.2626 1.2727 1.2828 1.2929 1.3030 1.3131 1.3232 1.3333 1.3434 1.3535 Columns 37 through 48 1.3636 1.3737 1.3838 1.3939 1.4040 1.4141 1.4242 1.4343 1.4444 1.4545

7、1.4646 1.4747 Columns 49 through 60 1.4848 1.4949 1.5051 1.5152 1.5253 1.5354 1.5455 1.5556 1.5657 1.5758 1.5859 1.5960 Columns 61 through 72 1.6061 1.6162 1.6263 1.6364 1.6465 1.6566 1.6667 1.6768 1.6869 1.6970 1.7071 1.7172 Columns 73 through 84 1.7273 1.7374 1.7475 1.7576 1.7677 1.7778 1.7879 1.7

8、980 1.8081 1.8182 1.8283 1.8384 Columns 85 through 96 1.8485 1.8586 1.8687 1.8788 1.8889 1.8990 1.9091 1.9192 1.9293 1.9394 1.9495 1.9596 Columns 97 through 100 1.9697 1.9798 1.9899 2.0000二、高等数学10.函数：round语法：round(x)将X四舍五入取整>> round(9.18)ans = 911函数：rem语法：rem(x,y)计算x/y的余数>> rem(13,5)ans

9、= 312.函数：loglog10语法：log(x)计算In(x)的值，即求以e为底的自然对数Log10(x)计算log10(x)的值，即求以10为底的常用对数>> log(10)ans = 2.3026>> log10(100)ans = 213.函数：floorceilfix语法：floor(x)返回沿负无穷大方向与X最接近的整数Ceil(x)返回沿正无穷大方向与X最接近的整数Fix(x)返回沿零方向与X最接近的整数>> floor(-93.5)ans = -94>> ceil(93.5)ans =94>> fix(-8.3)an

10、s = -814.函数：sign语法：sign(x)若X小于0，返回-1;若X等于0，返回0;若X大于0，返回1>> sign(456)ans = 1>> sign(0)ans = 0>> sign(-456)ans = -115函数：maxmin语法：max(x)返回向量最大值Min(x)返回向量最小值>> X=1 5 98 5 6 8 3 5 4 68 554;>> max(X)ans = 554>> min(X)ans = 116.函数：mean语法：mean(x)计算向量X的平均值，若X为矩阵，则返回一个包含矩阵X

11、的每一列的平均值的向量>> x=1 5 6 8 5;>> mean(x)ans = 5>> X=1 2 3;6 5 4;7 8 9;>> mean(X)ans =4.6667 5.0000 5.333317.函数：median语法：median(x)求向量X中的元素的中值，若X为矩阵，则返回一个包含矩阵X的每一列的中值的行向量>> x=1 5 6 8 4 5.1;>> median(x)ans = 5.0500>> X=1 2 3;6 5 4;7 8 9;>> median(X)ans = 6 5

12、418.函数：sum prod语法：sum(x)计算向量X的元素的和（若X为一个矩阵，则返回一个包含矩阵X的每一列的元素和的行向量）Prod(x)计算向量X的元素的积（若X为一个矩阵，则返回一个包含矩阵X的每一列的元素积的行向量）>> x=1,3,5;>> sum(x)ans = 9>> prod(x)ans = 15>> X=1,3,5;2,4,6;>> sum(X)ans = 3 7 11>> prod(X)ans = 2 12 3019.函数：sort语法：sort(x)对向量X中的元素升序排序（若X是一个矩阵，则对

13、矩阵的每一列按升序排序）>> x=1 5 3 8 9;>> sort(x)ans = 1 3 5 8 9>> X=1 5 6 8 3 45;9 5 85 6 3 4;9 8 56 12 5 54;>> sort(X)ans = 1 5 6 6 3 4 9 5 56 8 3 45 9 8 85 12 5 5420.函数：sphere语法：sphere(n)创建一个半径为n的球体（系统默认为20）>> sphere(50)21.函数：peaks语法：peaks(n)创建一个样本矩阵，用于演示绘图功能（MATLAB提供的三维曲面）>&

14、gt; peaks(50) z = 3*(1-x).2.*exp(-(x.2) - (y+1).2) . - 10*(x/5 - x.3 - y.5).*exp(-x.2-y.2) . - 1/3*exp(-(x+1).2 - y.2)（图形见下图）三、概率论与数理统计22.函数：nchoosek语法：nchoosek(m,n)(m>=n)求组合数>> nchoosek(10,8)ans =4523.函数：factorial语法：factorial(n)求N的阶乘>> factorial(10)ans = 362880024.函数：binopdf语法：binopd

15、f(X,N,P)求该事件在N次中恰好发生X次，且每次的概率均为P（0<=P<=1）的概率>> binopdf(6,10,0.3)ans = 0.036825.函数：binornd语法：binornd(N,P)产生一个参数为N, 概率为P的二项分布的随机数Binornd(N,P,m,n)产生一个矩阵，各元素为参数为N, 概率为P的二项分布的随机数>> binornd(20,0.25)ans = 3>> binornd(20,0.25,3,7)ans = 5 6 5 7 2 5 4 3 2 5 7 5 2 5 2 5 4 6 0 3 526.函数：unidrnd语法：unidrn