matlab数值计算命令与示例Word下载.docx
《matlab数值计算命令与示例Word下载.docx》由会员分享,可在线阅读,更多相关《matlab数值计算命令与示例Word下载.docx(19页珍藏版)》请在冰点文库上搜索。
2.4.2行列式det8
2.4.3特征值eig8
2.4.4转置'
和.'
9
2.4.5变形reshape10
2.4.6翻转rot90,fliplr,flipud11
2.4.7抽取diag,tril,triu12
2.5数组运算12
2.5.1乘12
N].*[M×
N]12
2.5.2除13
N]./[M×
N]14
N].\[M×
2.5.3乘方14
N].^[M×
N]15
a.^[M×
2.6除法15
2.6.1求解线性方程组15
3多项式16
3.1系数表示法poly16
3.2求根roots16
3.3乘法conv16
3.4除法deconv17
3.5求值polyval17
3.6微分polyder18
1创建矩阵
1.1直接输入
>
x=[123;
456]%实数
x=
1
2
3
4
5
6
y=[3*i
3+4*j;
sqrt(-1)+2
sin(3)]%复数、表达式
y=
0+3.0000i
3.0000+4.0000i
2.0000+1.0000i
0.1411
1.2向量
线性分布
x=linspace(1,10,9)
%初值,终值,元素个数
1.0000
2.1250
3.2500
4.3750
5.5000
6.6250
7.75008.8750
10.0000
y=linspace(1,4);
%初值,终值
100个元素
1.2.2冒号法
x=1:
0.8:
3
%初值:
间隔:
终值
1.0000
1.8000
2.6000
y=1:
4.6
终值
间隔为1
3
444
1.3函数创建
单位矩阵
x=eye(3),y=eye(2,4)
0
0
1
随机矩阵
x=rand(3),y=rand(2,4)
0.9501
0.4860
0.4565
0.2311
0.8913
0.0185
0.6068
0.7621
0.8214
0.4447
0.7919
0.7382
0.4057
0.6154
0.9218
0.1763
0.9355
全0矩阵
x=zeros(3),y=zeros(2,4)
全1矩阵
x=ones(3),y=ones(2,4)
2矩阵运算
2.1加减
N]
A=[1:
3;
4:
6],
B=ones(2,3),C=A+B,D=A-B
A=
B=
C=
4
6
7
D=
2
5
2.2乘
N]*a
B=A*3
12
15
18
M]
B=[11;
22;
33],C=A*B
14
14
32
32
2.3乘方
M]^a
6;
7:
9],
B=A^2
7
8
30
36
42
66
81
96
102
126
150
B=2^A
1.0e+004*
0.7962
0.9782
1.1603
1.8029
2.2154
2.6276
2.8097
3.4523
4.0950
2.4特殊运算
2.4.1求逆inv
429;
701],
B=inv(A)
0.0256
-0.0256
0.1538
0.7564
-0.2564
0.0385
-0.1795
0.1795
-0.0769
2.4.2行列式det
x=det(A)
78
2.4.3特征值eig
C=eig(A)
8.1313
-2.0656+2.3078i
-2.0656-2.3078i
M=A'
N=A.'
M=
9
N=
701],B=2*ones(3),C=A+B*i,
M=C'
N=C.'
1.0000+2.0000i
2.0000+2.0000i
3.0000+2.0000i
4.0000+2.0000i
9.0000+2.0000i
7.0000+2.0000i
0+2.0000i
1.0000+2.0000i
1.0000-2.0000i
4.0000-2.0000i
7.0000-2.0000i
2.0000-2.0000i
0-2.0000i
3.0000-2.0000i
9.0000-2.0000i
1.0000-2.0000i
7.0000+2.0000i
0+2.0000i
3.0000+2.0000i
9.0000+2.0000i
2.4.5变形reshape
A=1:
12,
B=reshape(A,3,4),C=reshape(A,2,6)
10
11
12
10
11
2.4.6翻转rot90,fliplr,flipud
B=rot90(A),C=rot90(A,-1)
%90*k,k=1,2...逆时针,-1,-2...顺时针
8
B=fliplr(A),C=flipud(A)
%fliplr左右翻转,flipup上下翻转
2.4.7抽取diag,tril,triu
B=diag(A),C=triu(A),D=tril(A)
%diag对角线,triu上三角矩阵,tril下三角矩阵
2.5数组运算
2.5.1乘
B=[111;
222],C=A.*B
%同维矩阵,对应元素相乘
2.5.2除
222],C=A./B,D=A.\B
%同维矩阵,对应元素相除
%除号,想象成杠杆撬动,上面是被除数
2.0000
3.0000
2.5000
0.5000
0.3333
0.4000
2.5.3乘方
222],C=A.^B
16
25
36
C=2.^A
64
2.6除法
2.6.1求解线性方程组
x+2y+2z=3
2x-y-z=6
x+y-z=2
·
AX=B
-1
-1
X=A\B
X=
-0.5000
0.5000
XA=B
X=B/A
3.0000
-0.5000
3多项式
3.1系数表示法poly
%由根确定多项式系数,按降幂排列
x=[123],p=poly(x),poly2str(p,'
x'
)
p=
-6
-6
ans=
x^3-6x^2+11x-6
3.2求根roots
%求x^3-6x^2+11x-6=0的根
p=[1-611-6],x=roots(p)
2.0000
1.0000
3.3乘法conv
%求多项式相乘a(x)=x^2+2x+3,b(x)=4x^2+5x+6,c(x)=a(x)*b(x)
a=[123],b=[4
6],c=conv(a,b)
a=
b=
c=
13
28
27
poly2str(c,'
)
4x^4+13x^3+28x^2+27x+18
3.4除法deconv
%求多项式相除a(x)=x^2+2,c(x)=4x^4+13x^3+28x^2+27x+18,b(x)=c(x)/a(x)
c=[4
18],a=[102]
[b,f]=deconv(c,a)
%b为商,f为余式
20
f=
-22
3.5求值polyval
%求f(x)=x^4+4x^3+5x^2+7x-3在x=2.6的值
f=polyval(c,2.6)
165.0016
3.6微分polyder
%求多项式的一阶导数
c=[1457-3],fc=poly2str(c,'
-3
fc=
x^4+4x^3+5x^2+7x-3
d=polyder(c),fd=poly2str(d,'
d=
fd=
4x^3+12x^2+10x+7
a=[123],b=[13],d=polyder(a,b),e=polyder(conv(a,b))
%求多项式的积的微分
e=
c=[1234],a=[13],[p,q]=polyder(c,a)
q=
%(a/b)'
=p/q