中南大学系统仿真实验报告Word格式文档下载.docx
《中南大学系统仿真实验报告Word格式文档下载.docx》由会员分享，可在线阅读，更多相关《中南大学系统仿真实验报告Word格式文档下载.docx（33页珍藏版）》请在冰豆网上搜索。
0.27850.95750.1576
0.54690.96490.9706
0.95720.80030.4218
0.48540.14190.9157
5、>
diag（5）
5
diag（5,5）
000005
000000
diag（2,3）
0002
6、
（inv（A）为求A的逆矩阵）>
B=[531;
238;
111],inv（B）
B=
531
238
111
0.62500.2500-2.6250
-0.7500-0.50004.7500
0.12500.2500-1.1250
A=[23;
44],B=[53;
38],inv（A）,inv（B）;
A\B,A/B,inv（A）*B,B*inv（A）
A=
23
44
53
38
-1.00000.7500
1.0000-0.5000
-2.75003.0000
3.5000-1.0000
0.22580.2903
0.64520.2581
-2.00002.2500
5.0000-1.7500
7、>
p=[1,-6,-72,-27],roots（p）
p=
1-6-72-27
12.1229
-5.7345
-0.3884
p=[2,3,6],roots（p）
236
-0.7500+1.5612i
-0.7500-1.5612i
8、（A为n*n的方阵）>
A=[010;
-440;
-212],poly（A）,B=sym（A）,poly（B）
010
-440
-212
1-612-8
B=
[0,1,0]
[-4,4,0]
[-2,1,2]
ans=
x^3-6*x^2+12*x-8
9,、（conv是多项式相乘，deconv是多项式相除）
u=[1246],v=[500-67],conv（u,v）
u=
1246
v=
500-67
5102024-5-10-842
v=[1246],u=[500-67],deconv（u,v）
5-10
10、（点乘是数组的运算，没有点的乘是矩阵运算）
a=[25;
34],b=[31;
47],a.*b,a*b
a=
25
34
b=
31
47
65
1228
2637
2531
a=[23];
b=[47];
a.*b=[821];
a*b%错误
a*b'
=29;
11、（who可以看到你用过的一些变量，whos是把该变量及所存储的大小等信息都显示出来了）>
who
Yourvariablesare:
ABaansbpuv
whos
NameSizeBytesClassAttributes
A2x232double
B2x232double
a1x216double
ans1x216double
b1x216double
p1x324double
u1x540double
v1x432double
12、>
A=[253;
654],disp（A）,size（A）,length（A）
253
654
3
实验二MATLAB绘图命令
熟悉MATLAB基本绘图命令，掌握如下绘图方法：
1．坐标系的选择、图形的绘制；
2．图形注解（题目、标号、说明、分格线）的加入；
3．图形线型、符号、颜色的选取。
t=[0:
pi/360:
2*pi];
x=cos（t）+cos（t*4）;
y=sin（t）+sin（t*4）;
xlabel（'
x轴'
）;
ylabel（'
y轴'
plot（y,x）,grid;
2、>
t=0:
0.1:
100;
x=3*t;
y=4*t;
z=sin（2*t）;
plot3（x,y,z,'
g:
'
）
3、>
x=linspace（-2*pi,2*pi,40）;
y=sin（x）;
stairs（x,y）
4、>
x=cos（t）+cos（t*4）+sin（t*4）;
plot（y,x,'
r:
th=[0:
pi/1000:
2*pi]'
;
r=cos（2*th）;
polar（th,r）;
title（'
四叶草图'
6、>
pi/20:
x=exp（j*th）;
plot（real（x）,imag（x）,'
r-.'
）;
grid;
text（0,0,'
中心'
x=-2:
0.01:
2;
8、y=-2:
9、[X,Y]=meshgrid（x,y）;
Z=Y.*exp（-X.^2-Y.^2）;
[C,h]=contour（X,Y,Z）;
set（h,'
ShowText'
'
on'
TextStep'
get（h,'
LevelStep'
）*2）
8、>
x=0:
0.2:
10;
y=2*x+3;
subplot（411）;
plot（x,y）;
grid;
y的原函数'
subplot（412）;
semilogy（x,y）;
对y取对数'
subplot（413）;
semilogx（x,y）;
对x取对数'
subplot（414）;
loglog（x,y）;
对xy均取对数'
9、>
x=-3:
0.3:
3;
bar（x,exp（-x.*x）,'
g'
）
实验三MATLAB程序设计
1．熟悉MATLAB程序设计的方法和思路；
2．掌握循环、分支语句的编写，学会使用lookfor、help命令。
程序举例
f=[1,1];
i=1;
whilef（i）+f（i+1）<
1000
f（i+2）=f（i）+f（i+1）;
i=i+1;
end
f,i
f=
Columns1through14
1123581321345589144233377
Columns15through16
610987
i=
15
m=3;
n=4;
fori=1:
m
forj=1:
n
a（i,j）=1/（i+j-1）;
end
formatrat
a
11/21/31/4
1/21/31/41/5
1/31/41/51/6
（分数格式形式。
用有理数逼近显示数据）
m=5;
end
1/41/51/61/7
1/51/61/71/8
3、程序中没有formatrat命令时，如果上次运行结果没有清除，输出的结果就是上次运行的结果！
但是运用clear命令清楚之前的运行结果之后就会正常运行。
x=input（'
请输入x的值:
ifx==10
y=cos（x+1）+sqrt（x*x+1）;
else
y=x*sqrt（x+sqrt（x））;
y
2
y=
2391/647
y=fprintf（'
不在定义域内，请重新输入：
return
y=1/（x-10）;
y
-1/8
p=[0001302009];
length（p）,ifp
（1）==0,p=p（2:
length（p））;
end;
end;
p
Columns1through5
13020
Columns6through7
09
p（p==0）=[];
1329
e2（500）
lookforffibno
e2-ffibno计算斐波那契亚数列的函数文件
helpe2
ffibno计算斐波那契亚数列的函数文件
n可取任意自然数
程序如下
（用法：
lookfor关键词
在所有M文件中找“关键词”，比如：
lookformax（即寻找关键词“max”）
其实就和我们平时用CTRL+F来查找“关键词”是
升级会员 