用matlab语言描述下列系统,东南大学信号与系统MATLAB实践第二次作业.doc

. . . .

专业.专注

练习二

实验六

一．用MATLAB语言描述下列系统，并求出极零点、

1.

>> Ns=[1];

Ds=[1,1];

sys1=tf(Ns,Ds)

实验结果：

sys1 =

1

-----

s + 1

>> [z,p,k]=tf2zp([1],[1,1])

z =

Empty matrix: 0-by-1

p =

-1

k =

1

2.

>>Ns=[10]

Ds=[1,-5,0]

sys2=tf(Ns,Ds)

实验结果：

Ns =

10

Ds =

1 -5 0

sys2 =

10

---------

s^2 - 5 s

>>[z,p,k]=tf2zp([10],[1,-5,0])

z =

Empty matrix: 0-by-1

p =

0

5

k =

10

二．已知系统的系统函数如下，用MATLAB描述下列系统。

1．

>> z=[0];

p=[-1,-4];

k=1;

sys1=zpk(z,p,k)

实验结果：

sys1 =

s

-----------

(s+1) (s+4)

Continuous-time zero/pole/gain model.

2.

>> Ns=[1,1]

Ds=[1,0,-1]

sys2=tf(Ns,Ds)

实验结果：

Ns =

1 1

Ds =

1 0 -1

sys2 =

s + 1

-------

s^2 - 1

Continuous-time transfer function.

3．

>> Ns=[1,6,6,0];

Ds=[1,6,8];

sys3=tf(Ns,Ds)

实验结果：

Ns =

1 6 6 0

Ds =

1 6 8

sys3 =

s^3 + 6 s^2 + 6 s

-----------------

s^2 + 6 s + 8

Continuous-time transfer function.

六．已知下列H(s)或H(z)，请分别画出其直角坐标系下的频率特性曲线。

1.

>> clear;

for n = 1:400

w(n) = (n-1)*0.05;

H(n) = (1j*w(n))/(1j*w(n)+1);

end

mag = abs(H);

phase = angle(H);

subplot(2,1,1)

plot(w,mag);title('幅频特性')

subplot(2,1,2)

plot(w,phase);title('相频特性')

实验结果：

2.

>> clear;

for n = 1:400

w(n) = (n-1)*0.05;

H(n) = (2*j*w(n))/((1j*w(n))^2+sqrt(2)*j*w(n)+1);

end

mag = abs(H);

phase = angle(H);

subplot(2,1,1)

plot(w,mag);title('幅频特性')

subplot(2,1,2)

plot(w,phase);title('相频特性')

实验结果：

3.

>>clear;

for n = 1:400

w(n) = (n-1)*0.05;

H(n) = (1j*w(n)+1)^2/((1j*w(n))^2+0.61);

end

mag = abs(H);

phase = angle(H);

subplot(2,1,1)

plot(w,mag);title('幅频特性')

subplot(2,1,2)

plot(w,phase);title('相频特性')

实验结果：

4.

>>clear;

for n = 1:400

w(n) = (n-1)*0.05;

H(n) =3*(1j*w(n)-1)*(1j*w(n)-2)/(1j*w(n)+1)*(1j*w(n)+2);

end

mag = abs(H);

phase = angle(H);

subplot(2,1,1)

plot(w,mag);title('幅频特性')

subplot(2,1,2)

plot(w,phase);title('相频特性')

实验结果：

实验七

三．已知下列传递函数H(s)或H(z)，求其极零点，并画出