《MATLAB之矩阵与多项式运算》由会员分享，可在线阅读，更多相关《MATLAB之矩阵与多项式运算（34页珍藏版）》请在金锄头文库上搜索。
1、第二讲 数值计算功能一、矩阵的生成数组与矩阵的运算1、逐个元素输入法 直接输入法 x=1 2 3 4 5 6 7 8 ;4 5 6 7 8 9 10 11; 1 2 3 4 5 6 7 8; x= ; y=2,4, 5 3 6 8y =2 4 53 6 8 a=1; b=2; c=3; x=5 b c; a*b a+c c/bx =5.0000 2.0000 3.00002.0000 4.0000 1.50002. 快速矩阵生成法（1）数组的冒泡生成法： x=a:inc:b y=1:1:8 y =1 2 3 4 5 6 7 8 y=0:0.2:1 y =0 0.2000 0.4000 0.60
2、00 0.8000 1.0000 x=1:8;4:11x =1 2 3 4 5 6 7 84 5 6 7 8 9 10 11(2) 采用Matlab的矩阵生成函数生成函数法举例说明矩阵生成函数：x=linspace(n1,n2,n) % 在区间上生成线性分度的向量 x=linspace(1,10,10)x =1 2 3 4 5 6 7 8 9 10 a= %生成一个空距阵a = a=logspace(n1,n2,n) %在区间上生成对数分度的向量 a=logspace(1,3,3)a =10 100 1000eyeeye(m,n(m,n); ); eyeeye(m(m) ) %生成单位矩阵生成
3、单位矩阵 eyeeye(2,3)(2,3) ansans= =1 0 0 1 0 0 0 1 0 0 1 0 eye(2) ans=1 0 0 1 zeroszeros(m,n(m,n); ); zeroszeros(m(m) ) % % 生成零矩阵生成零矩阵 zeroszeros(2,3)(2,3) ansans= =0 0 0 0 0 00 0 0 0 0 0 zeros(2) ans=0 0 0 0ones(m,n); ones(m) %生成一矩阵 onesones(2,3(2,3) )ansans= =1 1 1 1 1 11 1 1 1 1 1ones(2) ans=1 11 1V=
4、a1,a2,an A=V=a1,a2,an A=diagdiag(V(V) ) % %生成对角矩阵：对角元素向量生成对角矩阵：对角元素向量 V=5 7 2; V=5 7 2; A=A=diagdiag(V(V) ) A=A=5 0 0 5 0 00 7 0 0 7 00 0 2 0 0 2 A=1 2 3;2 3 4;3 4 5;V=diag(A) V =135rand(m,n) %随机矩阵：产生一个mn的均匀分布随机矩阵 rand(2,4) ans = 0.9501 0.6068 0.8913 0.4565 0.2311 0.4860 0.7621 0.0185 % 2x4随机矩阵 a=1:
5、1:10; b=0.1:0.1:1; c=b a; %组成一个新的数组ans =Columns 1 through 5 1.0000 + 0.1000i 2.0000 + 0.2000i 3.0000 + 0.3000i 4.0000 + 0.4000i 5.0000 + 0.5000iColumns 6 through 10 6.0000 + 0.6000i 7.0000 + 0.7000i 8.0000 + 0.8000i 9.0000 + 0.9000i 10.0000 + 1.0000ic =Columns 1 through 8 0.1000 0.2000 0.3000 0.4000
6、 0.5000 0.6000 0.7000 0.8000Columns 9 through 16 0.9000 1.0000 1.0000 2.0000 3.0000 4.0000 5.0000 6.0000Columns 17 through 20 7.0000 8.0000 9.0000 10.0000 a+b*i %复数数组的生成X,Y=meshprid(-1:0.4:1,-.5:0.2:0.5); %用于3D绘图 X X =-1.0000 -0.6000 -0.2000 0.2000 0.6000 1.0000-1.0000 -0.6000 -0.2000 0.2000 0.6000
7、1.0000-1.0000 -0.6000 -0.2000 0.2000 0.6000 1.0000-1.0000 -0.6000 -0.2000 0.2000 0.6000 1.0000-1.0000 -0.6000 -0.2000 0.2000 0.6000 1.0000-1.0000 -0.6000 -0.2000 0.2000 0.6000 1.0000 Y Y =-0.5000 -0.5000 -0.5000 -0.5000 -0.5000 -0.5000-0.3000 -0.3000 -0.3000 -0.3000 -0.3000 -0.3000-0.1000 -0.1000 -0
8、.1000 -0.1000 -0.1000 -0.10000.1000 0.1000 0.1000 0.1000 0.1000 0.10000.3000 0.3000 0.3000 0.3000 0.3000 0.30000.5000 0.5000 0.5000 0.5000 0.5000 0.5000 X,Y = meshgrid(-2:.2:2, -2:.2:2);Z = X .* exp(-X.2 - Y.2);mesh(Z) %画出3D图二. 矩阵的子矩阵寻访与赋值1.子数组寻址x=1 4 12 3 6 4 7 5 8 6 9 7 10 8 11 x(3) % x的第三个元素ans =
9、12x(1 2 5) % x的第1，2 ，5个元素ans =1 4 6x(1:5) %x的前5个元素ans =1 4 12 3 6x(10:end) %x的第10个元素后的元素ans =6 9 7 10 8 11x(10:-1:2) %x的第10个元素和第二元素的倒排ans =6 8 5 7 4 6 3 12 4x(find(x5) %查找x中大于5的非零元素ans =12 6 7 8 6 9 7 10 8 11A(:) %由A的各列按从左到右的次序，首尾相接而生成“一维长列”数组矩阵寻址A=1 2 3 4;3 4 5 6A(2,3) % A的第二行第三列的元素 A =1 2 3 43 4 5
10、 6ans =5A(2,:) %A的第二行全部列的元素ans =3 4 5 6A(:,3) %A的全部行第3列的元素ans =35 ans =132435462. 数组赋值x=1 4 12 3 6 4 7 5 8 6 9 7 10 8 11x(4)=100 %给x的第4个元素重新赋值为100 x =Columns 1 through 13 1 4 12 100 6 4 7 5 8 6 9 7 10Columns 14 through 15 8 11 x(3)= %删掉掉3个元素 x =Columns 1 through 13 1 4 100 6 4 7 5 8 6 9 7 10 8Column
11、 14 11x(16)=1 %加入第16个元素x =Columns 1 through 13 1 4 100 6 4 7 5 8 6 9 7 10 8Columns 14 through 16 11 0 1矩阵赋值A=1:8;4:11 %直接赋值 A =1 2 3 4 5 6 7 84 5 6 7 8 9 10 11 A(:)=1:4;2:5;3:6;4:7 %全元素赋值，保持A的“行宽，列长”。A和D 两个数组的总元素相等，但“行宽，列长”不一定相同。A =1 3 2 4 3 5 4 62 4 3 5 4 6 5 7数组赋值x=1 4 12 3 6 4 7 5 8 6 9 7 10 8 11x(1 4)=1 1 %把当前 x 数组的第一、四个元素都赋值为 1 x =Columns 1 through 13 1 4 12 1 6 4 7 5 8 6 9 7 10Columns 14 through 15 8 11 D=1:4;2:5;3:6;4:7D =1 2 3 42 3 4 53 4 5 64 5 6 7 D(:) a
