实验一 Matlab使用方法和程序设计
一、实验目的
1、掌握Matlab软件使用的基本方法；
2、熟悉Matlab的数据表示、基本运算和程序控制语句；
3、熟悉Matlab绘图命令及基本绘图控制；
4、熟悉Matlab程序设计的基本方法。
二、实验内容
1、帮助命令
2、矩阵运算
（1）矩阵的乘法和乘方
已知A=[1 2;3 4]:B=[5 5;7 8]:求A^2*B
( 2 )矩阵除法
已知A=[1 2 3:4 5 6:7 8 9]:B=[1 0 0:0 2 0:0 0 3],求矩阵左除A\B,右除A/B。
( 3 )矩阵的转置及共轭转置
已知A=[5+i,2-i,1;6*i,4,9-i],求该复数矩阵的转置A',共轭转置A.'
( 4 )使用冒号选出指定元素
已知：A=[1 2 3:4 5 6:7 8 9];求A中第3列前2个元素;A中第2、3行元素。
( 5 )方括号[]
用magic函数生成一个4阶魔术矩阵，删除该矩阵的第四列
3、多项式
（1）求多项式P(x)=x3-2x-4的根
( 2 )已知A=[1.2 3 5 0.9;5 1.7 5 6;3 9 0 1;1 2 3 4].,构造多项式，并计算多项式值为20的解。
4、基本绘图命令
( 1 ) 绘制余弦曲线
]
2,0[
),
cos(π
∈
=t
t
y
。
( 2 ) 在同一坐标系中绘制曲线
]
2,0[
),
5.0
sin(
),
25
.0
cos(
yπ
∈
-
=
-
=t
t
y
t
5、基本绘图控制
绘制
]
4,0[π
区间上的y=10sint曲线，并要求：
（1）线形为点划线，颜色为红色，数据点标记为加号；
（2）坐标轴控制：显示范围，刻度线，比例，网络线；
（3）标注控制：坐标轴名称，标题，相应文本。
6、基本程序设计
（1）编写命令文件：计算1+2+....+n<2000时的最大n值；
(2)编写函数文件：分别用n和which循环结构编写程序，求2的0到n次幂的和.
n=input('输入正数n:') ji=1: for i=1:n; ji=ji+2^i: end ji input('输入正数n:') ji-1:i-1: While i<=n ji=ji+2^i; i=i+1; end ji
(3)如果想对一个变量x赋值，当从键盘输入y或Y时，x自动赋为1；当从键盘输入n或N时，x自
动赋为0；输入其他字符时终止程序。
k=input('shuruX:'.'s'): if k=='y' k=='Y' x=1 ; else k=='n' k=='N' x=0; else ruturn end >> n=input('输入正数n:')
输入正数n:20
n =20
>> ji=1;
>> for i=1:n
ji=ji+2^i;
end
>> ji
ji =2097151
>> n=input('输入正数n:')
输入正数n:20
n=20
>> ji=1;
>> i=1
i =1
>> while i<=n
ji=ji+2^i;
i=i+1;
end
>> ji
ji =2097151
>> k=input('shuruX:','s'); shuruX:y
>> if k=='y' k=='Y'
x=1
elseif k=='n' k=='N'
x=0
else
return
end
ans =0
x = 1
>> k=input('shuruX:','s');
shuruX:n
>> if k=='y' k=='Y'
x=1
elseif k=='n' k=='N'
x=0
else
return
end
ans =0
x =0
一．实验目的
1.掌握如何使用Matlab 进行系统的时域分析；
2.掌握如何使用Matlab 进行系统的频域分析。
二．实验内容
1.时域分析
（1）根据下面传递函数模型，绘制其单位阶跃响应曲线并从图上读取最大超调量、调节时间、上升时间，绘制系统的单位脉冲响应、零输入响应曲线（设初始状态x0=[1,0,0]）。 G(s)=8106)
65(5232+++++S S S s s
（2）根据下面传递函数模型，绘制其单位阶跃响应曲线并编程序求该系统的上升时间、
调节时间、峰值时间、超调量和终值（∞）。 G(s)=524
2+++s s s
（3）典型二阶系统传递函数为：G(s)=222
2ωξωω++s s ,当ξ=0.7,wn 取2、4、6、8、10、12的单位阶跃响应。
2.频域分析
（1）典型二阶系统传递函数为：G(s)=222
2ωξωω++s s ,在同一幅图上绘制当ξ=0.7,wn 取2、4、6、8、10、12的伯德图。
运行程序及结果如下：
（2）已知系统开环传递函数为：G(s)H(s)=)1()
1(52++Ts s s τ在同一幅图上绘制当τ=3，T=8和
τ=8，T=3的奈氏图。
