整理矩阵基础.docx
- 文档编号:29067879
- 上传时间:2023-07-20
- 格式:DOCX
- 页数:16
- 大小:21.26KB
整理矩阵基础.docx
《整理矩阵基础.docx》由会员分享,可在线阅读,更多相关《整理矩阵基础.docx(16页珍藏版)》请在冰豆网上搜索。
整理矩阵基础
(0)教材第9页的表1.4、表1.5、第11页的表1.6需要重点看。
Matlab程序运行结果
a=2;b=2;
whoa;
Yourvariablesare:
a
(1)who
Yourvariablesare:
ab
(2)whosab
NameSizeBytesClassAttributes
a1x18double
b1x18double
(3)whosa,b
NameSizeBytesClassAttributes
a1x18double
b=
2
(4)i=exist(a)
?
?
?
Errorusing==>exist
Thefirstinputtoexistisastring.
(5)i=exist('a')
i=
1
(6)i=exist('c')
i=
0
(7)helpcos
COSCosineofargumentinradians.
COS(X)isthecosineoftheelementsofX.
Seealsoacos,cosd.
Overloadedmethods:
distributed/cos
codistributed/cos
sym/cos
ReferencepageinHelpbrowser
doccos
(8)lookforcos
返回多个调用cos的信息
(9)helplog
LOGNaturallogarithm.
LOG(X)isthenaturallogarithmoftheelementsofX.
ComplexresultsareproducedifXisnotpositive..
Seealso
log1p,log2,log10,exp,logm,reallog
Overloadedmethods
gf/log,distributed/log,codistributed/log,fints/log,designdev/log,sym/log
(10)a=[123;456;789]
a=
123
456
789
(11)savea
在当前工作目录下,出现一个a.mat文件,但是在脱离matlab环境时,只看到a.mat的快捷方式。
(12)clear
>>a
?
?
?
Undefinedfunctionorvariable'a'.
(13)loada
>>a
a=
123
456
789
(14)path(path,'c:
:
\MyDir')
Warning:
Directoryaccessfailure:
c:
:
\MyDir.
>Inpathat110
(15)what
M-filesinthecurrentdirectoryC:
\DocumentsandSettings\tiger\MyDocuments\MATLAB
Untitledabcmovie1tu1
Untitled2classtempmovie2
Untitled3hilberttest1
MAT-filesinthecurrentdirectoryC:
\DocumentsandSettings\tiger\MyDocuments\MATLAB
asample_file
m
(16)dir
.exp.txtsample_file.mat
..exp1.txttest.dat
Untitled.asvfdat.bintest1.dat
Untitled.mhilbert.asvtest1.fig
Untitled2.mhilbert.mtest1.m
Untitled3.mm.mattext2.txt
a.matmkmatlab.txttu1.asv
abc.mmovie1.mtu1.m
classtemp.asvmovie2.m
classtemp.msample_file
(17)typetu1
x=0:
0.1:
20;
y=1-1/sqrt(1-0.3^2)*exp(-0.3*x).*sin(sqrt(1-0.3^2)*x+acos(0.3));
plot(x,y)
h=line(0,0,'color','red','marker','.','markersize',40,'erasemode','xor');%定义红色的小球
fori=1:
length(x)
set(h,'xdata',x(i),'ydata',y(i));pause(0.005)%暂停0.005秒
drawnow%刷新屏幕
end
(18)whichtu1
C:
\DocumentsandSettings\tiger\MyDocuments\MATLAB\tu1.m
(19)matlabroot
ans=
C:
\ProgramFiles\MATLAB\R2008b
(20)%MATLAB命令窗口的字符和数值的不同颜色显示
a=12.5
b='Hello'
ifa>0c=5,end
(21)%MATLAB中不同标点符号的应用
a=12.5,b='Hello'%逗号表示分隔命令,单引号构成字符串,点号为小数点
c=[12;34;56]%[]表示构成矩阵,分号分隔行,空格分隔元素
d=a*...%...表示续行
c
(22)%MATLAB的通用操作界面各个窗口的使用
a=[123;456;789];
b=[111;222;333];
c='计算';
d=a+b*i
(23)%用from:
step:
to生成向量和矩阵
x1=2:
5
x2=2:
0.5:
4
x3=5:
-1:
2
x4=2:
-1:
3
x5=2:
-1:
0.5
x6=[1:
2:
5;1:
3:
7]
(24)%用linspace和logspace函数生成行向量
x1=linspace(0,2*pi,5)%从0到2π等分成5个点(4等分)
x2=logspace(0,2,3)%从1到100对数等分成3个点
(25)%查看eye函数的功能
X1=eye(2,3)
X2=eye(3,2)
X3=eye(size(X1))
(26)%使用全下标和单下标方式对矩阵的元素进行操作
a=[12;34;56]
a(3,3)%提取a(3,3)的值(提取不出来,因为下标越界)
a(3,3)=9%给赋a(3,3)值(自动对a进行扩充)
(27)%用逻辑矩阵来提取矩阵的元素
a=[12;34;56]
a(3,3)=9
l1=logical([101])%给出逻辑向量l1
l2=logical([110])%给出逻辑向量l2
a(l1,l2)
b=a>1%得出逻辑向量b
a(b)%按单下标顺序排成长列
(28)%用全下标、单下标和全元素方式给矩阵元素赋值
a(1:
2,1:
3)=[111;111]%给第一、二行元素赋值为全1
a(8:
9)=[11]%给第8、9元素赋值为全1
a=[12;34;56]
b=[123;456]
a(:
)=b
(29)%对矩阵元素进行删除
a=[12;34;56]
a(3,3)=9
a(:
3)=[]%删除一列元素
a
(1)=[]%删除一个元素,则矩阵变为行向量
a=[]%删除所有元素为空矩阵
(30)%生成大矩阵
a=[12;34;56]
a(3,3)=9
[a;a]%联接成6×3的矩阵
[aa]%联接成3×6的矩阵
[a(1:
2,1:
2)10*a(1:
2,2:
3)]%计算并联接
(31)%对字符串使用字符串函数和运算
str1='Hello'
str2='Ilike''MATLAB'''%可重复单引号的来输入含有单引号的字符串
str3='你好!
'%支持中文
length(str1)%字符串长度
x1=double(str1)%查看字符串的ASCII码
x2=str1+1%字符串的数值运算
char(x1)%将ASCII码转换成字符串形式
char(x2)%将ASCII码转换成字符串形式
class(str1)
class(x1)
ischar(str1)
(32)%多个字符串构成一个新字符串和执行、显示字符串
str1='Hello'
str2='Ilike''MATLAB'''%可重复单引号的来输入含有单引号的字符串
str3='你好!
'%支持中文
str4=[str1,'!
',str2]%多个字符串并排成一个
str5=[str1;str3,'']%将str3添加两个空格
str6=char(str1,str2,str3)
str7=char(str1,str2,str3)
str8=strvcat(str1,str2)
whos
str9='a=2*5'
eval(str9)%执行字符串
disp('请输入2*2的矩阵a')
(33)%矩阵和数组的加、减和乘法运算
x1=[12;34;56];
x2=eye(3,2)
x1+x2%矩阵相加
x1.*x2%数组相乘
x1*x2%矩阵相乘x1列数不等于x2行数
x3=eye(2,3)
x1*x3%矩阵相乘
(34)%用矩阵除法来解线性方程组
A=[2-13;31-5;4-11]
B=[5;5;9]
X=A\B%求方程解
(35)%矩阵和数组的除法和乘方运算
x1=[12;34];
x2=eye
(2)
x1/x2%矩阵右除
x1\x2%矩阵左除
x1./x2%数组右除
inv(x1)%求逆矩阵
x1.\x2%数组左除
x1^2%矩阵乘方
x1^-1%矩阵乘方,指数为负整数
x1^0.2%矩阵乘方,指数为小数
2^x1%标量乘方(计算规则?
)
2.^x1%数组乘方
x1.^x2%数组乘方
(36)%矩阵和数组转置运算
x1=[12;34];
x2=eye
(2);
x3=x1+x2*i
x3'%矩阵转置
x3.'%数组转置
(37)%数组的算术运算函数
x=linspace(0,2*pi,6)
y=sin(x)%计算正弦
y1=abs(y)%计算绝对值,将正弦曲线变成全波整流
1-exp(-x).*y%计算按指数衰减的正弦曲线
(38)%数组的关系和逻辑运算
t=linspace(0,3*pi,100);
y=sin(t);%计算正弦曲线
t1=(t
y1=t1.*y;%得出0~π和2π~3π的半波整流
plot(t,y1)
axis([03*pi,-1,1])
(39)%用全下标元素赋值方式创建多维数组
a(:
:
2)=[12;34]%创建三维数组
b=[11;22]%先创建二维数组
b(:
:
2)=5%扩展数组
(40)%用函数来生成多维数组
a=[12;34];
b=[11;22];
c=cat(2,a,b)%沿着第二维连接生成数组c
cat(3,a,b)%沿着第三维连接
repmat(a,[222])%放置数组a
reshape(c,[222])%重组二维数组为三维数组
(41)%用函数得出矩阵的大小
a=[12;34;56]
ndims(a)%得出维数
size(a)%得出各维的大小
size(a,2)%得出列的大小
length(a)%得出最大维的大小
(42)%产生稀疏矩阵并转变为全元素矩阵,得出的稀疏矩阵分布图
a=eye(3);
a(4,:
)=[-5-2-3]
b=sparse(a)%创建稀疏矩阵
b=sparse([142434],[112233],[1-51-21-3])%创建与b相同的稀疏矩阵
full(b)
D=[329;249;114]
d=[012]
s=spdiags(D,d,4,3)%这个命令的作用?
理解上移的意思
full(s)
loadspr.dat%从教材59页3开始,装载文本文件
b=spconvert(spr)%转换成稀疏矩阵
savesprb%保存spr.mat文件
clear
loadsprb%装载spr.mat文件
b
nnz(b)
nonzeros(b)
nzmax(b)
spy(b)
(43)%计算多项式的值和根,并计算特征根和部分分式展开
p1=[121200];%降幂排列多项式的系数形成行向量
polyval(p1,2)%计算x=2时多项式的值
x=0:
0.5:
3
polyval(p1,x)%计算x为矩阵时多项式的值
roots(p1)%计算多项式的根
poly([0;-20;-1])%用多项式的根生成多项式,输入的根向量应为列向量
s=[12;34]
p2=poly(s)%计算特征多项式
roots(p2)%计算特征根
p3=[100200];
[r,p,k]=residue(p3,p1)
(44)%计算多项式的乘积和微积分
a1=[10];%对应多项式‘s’
a2=[11];%对应多项式‘s+1’
a3=[120];%对应多项式‘s+20’
p1=conv(a1,a2)
p1=conv(p1,a3)%对应s(s+1)(s+20)
[p2,r]=deconv(p1,a3)%计算多项式除法的商和余子式
conv(p2,a3)+r%用商*除式+余子式验算
p3=polyder(p1)%多项式微分
s=length(p3):
-1:
1
p1=[p3./s,0]%多项式积分,常数k=0
(45)%对多项式进行曲线拟合和插值
x1=1:
10
p=[2-1510];
y0=polyval(p,x1)
p1=polyfit(x1,y0,1)%一阶拟合
p2=polyfit(x1,y0,2)%二阶拟合
p3=polyfit(x1,y0,3)%三阶拟合
y1=polyval(p1,x1)
y2=polyval(p2,x1)
y3=polyval(p3,x1)
figure
(1)
plot(x1,y0,'o')
holdon
plot(x1,y1)
plot(x1,y2,'r')
holdon
plot(x1,y3)
figure
(2)
y01=interp1(x1,y0,9.5)%线性插值(默认)
y02=interp1(x1,y0,9.5,'spline')%三次样条插值
plot(x1,y0)
holdon
plot(9.5,1682,'+r')
plot(9.5,1696,'+b')
(46)%创建元胞数组并对元胞数组的内容进行操作
A={'ThisisthefirstCell.',[12;34];eye(3),{'Tom','Jane'}}
whos
B(1,1)={'ThisisthesecondCell.'}
B(1,2)={5+3*i}
B(1,3)={[12;34;56]}
C{1,1}='ThisisthethirdCell.';
C{2,1}=magic(4)
celldisp(A)
celldisp(B)
celldisp(C)
cellplot(A)
x1=A{1,2}%取A(1,2)元胞元素的内容
x2=A{1,2}(2,2)%取A(1,2)元胞元素的矩阵第二行第二列内容
x3=A(1,2)
[x4,x5,x6]=deal(A{[2,3,4]})
(47)%创建结构数组
ps
(1).name='曲线1'
ps
(1).color='red'
ps
(1).position=[0,0,300,300]
ps
(2).name='曲线2';
ps
(2).color='blue';
ps
(2).position=[100,100,300,300]
ps
(1)=struct('name','曲线1','color','red','position',[0,0,300,300]);
ps
(2)=struct('name','曲线2','color','blue','position',[100,100,300,300])
x1=ps
(1)
x2=ps
(1).position
x3=ps
(1).position(1,3)
x4=getfield(ps,{1},'color')
x5=getfield(ps,{1},'color',{1})
ps=setfield(ps,{1},'color','green');
ps
(1)
x6=fieldnames(ps)
all_x=[ps.name]
cat(1,ps.position)%沿第一维排列
cat(2,ps.position)%沿第二维排列
cat(3,ps.position)%沿第三维排列
(48)%矩阵的微分和积分运算
t=0:
0.5:
10;
y=exp(-0.2).*sin(t);
y=exp(-0.2).*sin(t)
d=[0diff(y)]%计算微分
plot(t,y)
holdon
plot(t,d,':
')
s1=0.5*cumsum(y)%用矩形法计算积分
plot(t,s1,'-.')
s2=cumtrapz(t,y)%用梯形法计算积分
plot(t,s2,'--')
环境的两个特点:
(49)%利用傅立叶变换和卷积公式计算离散序列的卷积
A=ones(1,10);
A
(1)=0
安全评价可针对一个特定的对象,也可针对一定的区域范围。
B=ones(1,9);
B([123])=0
(3)专项规划环境影响报告书的内容。
除包括上述内容外,还应当包括环境影响评价结论。
主要包括规划草案的环境合理性和可行性,预防或者减轻不良环境影响的对策和措施的合理性与有效性,以及规划草案的调整建议。
C=conv(A,B)%计算卷积
(3)环境影响分析、预测和评估的可靠性;N=32;%序列长度为32
一、环境影响评价的发展与管理体系、相关法律法规体系和技术导则的应用AF=fft(A,N);%傅立叶变换
BF=fft(B,N);
(5)建设项目对环境影响的经济损益分析。
CF=AF.*BF;
CC=real(ifft(CF))%过滤掉虚部
(5)污染防止措施能否达到要求。
plot(CC)
环境影响经济损益分析一般按以下四个步骤进行:
holdon
plot(C,'r')
(一)规划环境影响评价的适用范围和责任主体(50)%计算向量的叉乘和点乘
a=[123];
一、环境影响评价的发展与管理体系、相关法律法规体系和技术导则的应用b=[456];
c=cross(a,b)
c=dot(a,b)
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 整理 矩阵 基础