信号与系统-实验6-离散时间系统的 Z 域分析

一、实验目的

1、掌握 z 变换及其性质;了解常用序列的 z 变换、逆 z 变换;

2、掌握利用 MATLAB 的符号运算实现 z 变换;

3、掌握利用 MATLAB 绘制离散系统零、极点图的方法;

4、掌握利用 MATLAB 分析离散系统零、极点的方法;

5、掌握利用 MATLAB 分析离散系统频率响应的方法。

二、实验内容

1、分别求cos(kn)和 的 z 变换。

Matlab 复制代码
clear;
syms k n a;
f=cos(k*n);
f1=a^k;
disp("cos(kn)的z变换为");
y=ztrans(f)
disp("a^k的z变换为" );
y1=ztrans(f1)

2、系统函数为

分析并绘制出离散系统的零、极点图。

Matlab 复制代码
clear;
num=[0.2,0.1,0.3,0.1,0.2];
den=[1,-1.1,1.5,-0.7,0.3];
subplot(121)
zplane(num,den)

num1=[3,-5,10,0];
den1=[1,5,2,-10];
subplot(122)
zplane(num1,den1)

3、绘制 的频率响应图形。

Matlab 复制代码
clear;
num=[1,1];
den=[2 2 1];

[H,w]=freqz(num,den);
Hf=abs(H);  %取幅度值实部
Hx=angle(H);  %取相位值对应相位角
clf
subplot(121)
plot(w,20*log10(Hf))  %幅值变换为分贝单位
title('离散系统幅频特性曲线')
subplot(122)
plot(w,Hx)
title('离散系统相频特性曲线')

4、已知某一离散系统的系统函数为 ,试用 MATLAB 绘出该系统的零、极点分布图及幅频特性曲线,并分析该系统 的频率特性。

Matlab 复制代码
clear;
B=[1 5 5 1];
A=[1 0 0 0];
N=60;
[h,t]=impz(B,A,N);

subplot(211)
plot(t,h)
title("h(t)");
xlabel("N");ylabel("h(t)");
subplot(212)
zplane(B, A) 

figure
freqz(B,A);

5、已知某离散系统的系统函数为 ,试用 MATLAB 求出 该系统的单位响应h(k)。

Matlab 复制代码
clear;
num=[1,0,0]
den=[1,5,6]
impz(num,den)

6、已知某序列的 z 变换为 ,试用 MATLAB 求 F(z) 的逆变换。

Matlab 复制代码
clear;
syms z;
F=(z^2+2*z)/((z-0.1)*(z-0.2)*(z-0.3));
X=iztrans(F)

7、已知某离散系统的系统函数为 ,试用 MATLAB 分析 该系统的频率特性,绘制其幅频及相频特性曲线。

Matlab 复制代码
B=[1.25 1.25];
A=[1 -0.25];
freqz(B, A) 
相关推荐
XuX032 小时前
手搓雷达图(MATLAB)
matlab·贴图
freexyn6 小时前
Matlab自学笔记五十一:(推荐)输入参数的数量和可变数量的输入
笔记·算法·matlab
不吃酸的柠檬6 小时前
MATLAB 中的图形绘制
人工智能·机器学习·matlab
studyer_domi6 小时前
Matlab 复合模糊PID
开发语言·matlab
DarrenPig1 天前
【新能源科学与技术】MATALB/Simulink小白教程(一)实验文档【新能源电力转换与控制仿真】
matlab·开源·github·simulink·交流
简简单单做算法1 天前
基于GA遗传优化TCN-BiGRU注意力机制网络模型的时间序列预测算法matlab仿真
matlab·tcn-bigru·时间序列预测·注意力机制·ga遗传优化
Cc小跟班1 天前
MATLAB小技巧记录(特殊符号、图例位置...)
开发语言·算法·matlab
随风飘摇的土木狗2 天前
【MATLAB第116期】基于MATLAB的NBRO-XGBoost的SHAP可解释回归模型(敏感性分析方法)
matlab·gsa·敏感性分析·特征排序·灵敏度分析·shap·可解释
海天一色y2 天前
matlab设置不同颜色的柱状图
matlab·信息可视化
海洋与大气科学2 天前
【matlab|python】矢量棍棒图应用场景和代码
开发语言·python·matlab