整理矩阵基础.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=（t2*pi）%表示什么？
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）
升级会员 