matlab离散信号与系统的时域分析
离散信号与系统的时域分析
摘要:通过对数字信号处理的深刻理解,在有MATLAB对信号进行仿真前提下,使得离散信号输出更准确更形象。利用MATLAB工具箱中filter函数对差分方程求解,防真出系统响应和系统单位脉冲响应,用conv函数完成线形卷积。在判断系统稳定性时,巧妙的将输入信号改为单位阶跃序列,观察稳态输出是否趋于常数(包括零)来判断,如果输出趋于数,则系统稳定,否则系统不稳定。然而,完整的模拟也会存在细小的瑕疵,所以本文还需要改进。尤其是在卷积编程时,冗余度较大。
报告要求:
(1)简述离散系统时域分析方法;
答:离散系统在时域中,是通过差分方程,单位脉冲响应或系统函数对系统特性进行描写。所以在输入信号已知的情况下,知道其中之一就可以利用MATLAB就可以对其进行求解。倘若一直差分方程就适合用工具箱函数filter求差分方程的解。也可以用conv函数计算输入信号核系统的单位脉冲响应的线形卷积,求出系统响应。
(2)简述通过实验判断系统稳定性的方法;
答:判断系统稳定性,从时域看可以通过检查系统单位脉冲响应是否绝对可和来判断。但实际中,只要用单位阶跃序列作为输入信号,如果稳态输出趋于常数或零,则系统稳定,否则系统不稳定。
(3)完成设计实验,并对结果进行分析和解释;
(4)打印程序清单和要求画出的信号波形;
(5)写出本次课程设计的收获和体会。
程序清单和波形图
定义基础函数:
1 uDT为单位阶跃序列
function y=uDT(n)
y=n>=0; %当参数大于等于0时阶跃为1,否则为0
- impDT为单位采样序列
function y=impDT(n)
y=(n==0); %当参数为0时冲激为1,否则为0
第一题:
1 给定一个低通滤波器的差分方程为
程序一:
%系统响应%
N=8;
a=[1 -0.9];
b=[0.05 0.05];
n=0:30;
subplot(211)
x1=uDT(n)-uDT(n-N);
y1=filter(b,a,x1);
stem(n,y1,'fill','g'),grid on
xlabel('n'),title('系统响应y1(n)')
subplot(212)
x2=uDT(n);
y2=filter(b,a,x2);
stem(n,y2,'fill','r'),grid on
xlabel('n'),title('系统响应y2(n)')
2 系统的单位脉冲响应波形图如下:
程序二:
%系统单位脉冲响应%
a=[1 -0.9];
b=[0.05 0.05];
n=0:30;
impz(b,a,30),grid on
title('系统单位脉冲响应h(n)')
第二题:
用线性卷积法求x1(n)=R8(n)分别对系统h1(n)和h2(n)的输出响应,画出波形如下:
程序三:
%线性卷积%
nx=-1:5; %x(n)向量显示范围(添加了附加的零值)
nh=-2:10; %h(n)向量显示范围(添加了附加的零值)
x=uDT(nx)-uDT(nx-8);
h1=uDT(nh)-uDT(nh-10);
h2=impDT(nh)+2.5.*impDT(nh-1)+2.5.*impDT(nh-2)+impDT(nh-3);
y1=conv(x,h1);
y2=conv(x,h2);
ny1=nx(1)+nh(1); %卷积结果起始点
%卷积结果长度为两序列长度之和减1,即0到(length(nx)+length(nh)-2);
%因此卷积结果的时间范围是将上述长度加上起始点的偏移值
ny=ny1+(0:(length(nx)+length(nh)-2));
subplot(321)
stem(nx,x,'fill'),grid on
xlabel('n'),title('x(n)')
axis([-4 16 0 3])
subplot(322)
stem(nh,h1','fill'),grid on
xlabel('n'),title('h1(n)')
axis([-4 16 0 3])
subplot(323)
stem(nh,h2','fill'),grid on
xlabel('n'),title('h2(n)')
axis([-4 16 0 3])
subplot(324)
stem(ny,y1,'fill'),grid on
xlabel('n'),title('y(n)=x(n)*h1(n)')
axis([-4 16 0 3])
subplot(325)
stem(ny,y2,'fill'),grid on
xlabel('n'),title('y(n)=x(n)*h2(n)')
axis([-4 16 0 3])
第三题:
用实验方法检查系统是否稳定。输入信号为u(n)时,画出系统输出波形如下:
从频域看:
程序四:
B=[1,1.8237,-0.9801];
A=[1/100.49,0,-1/100.49];
zplane(B,A),grid on
legend('零点','极点')
title('零极点分布图')
(由于系统所有极点在单位圆上,则系统稳定)
从时域看:
程序五:
%稳定性判断
a=[1 -1.8237 0.9801];
b=[1/100.49 0 -1/100.49];
n=0:300;
x1=uDT(n);
y1=filter(b,a,x1);
stem(n,y1,'fill''b'),grid on
xlabel('n'),title('系统响应y1(n)')
(由于系统输出趋于零,则系统稳定)
(2) 给定输入信号,求出系统的输出响应,画出其波形图如下:
程序六:
%差分方程对应的系统函数
a=[1 -1.8237 0.9801];
b=[1/100.49 0 1/100.49];
n=0:30;
subplot(211)
x1=uDT(n);
y1=filter(b,a,x1);
stem(n,y1,'fill','g'),grid on
xlabel('n'),title('系统响应y1(n)')
subplot(212)
x2=sin(0.014.*n)+sin(0.4.*n);
y2=filter(b,a,x2);
stem(n,y2,'fill','r'),grid on
xlabel('n'),title('系统响应y2(n)')
参考文献:
[1] 吴礼斌,李伯年,数学实验与建模 北京,国防工业出版社2007
[2] 张威 MATLAB基础与编程入门(第二版)西安, 西安电子科技大学 2008
[3] 高西全,丁玉美,数字信号处理(第三版)西安,西安电子科技大学 2008
[4] Vinay K Lngle,John G Proakis.数字信号处理及其MATLAB实现 陈怀琛 王朝英,高西全, 译.北京:电子工业出版社,1988
matlab离散信号与系统的时域分析相关推荐
- 离散信号内插matlab,离散信号和系统实验报告.doc
离散信号和系统实验报告 三.实验效果分析(包括仪器设备等使用效果) 实验中注意中英文的切换是应注意,特别是标点符号,括号等需要注意.需要仔细键入程序. 教 师 评 语 指导老师 年 月 日 江西师范大 ...
- matlab 复频域,MATLAB信号与系统分析(四)——离散信号与系统的复频域分析及MATLAB实现...
一.系统的z变换和反变换 1.利用部分分式展开求解逆Z变换: 2.例子 3.Z变换的MATLAB函数 clear all f=sym('cos(a*k)'); F=ztrans(f) F=sym('z ...
- 信号与系统 徐亚宁 matlab程序,信号与系统(第4版)
第1章绪论 11信号与系统 12信号的描述与分类 13系统的描述与分类 131连续时间LTI系统及其描述 132离散时间LTI系统及其描述 14系统的基本特性 15信号与系统分析方法 ...
- 海南大学matlab试卷,信号与系统
spContent=你想知道现代通信背后的原理是什么吗? 你想知道现代控制背后的原理是什么吗? 你想知道信息处理背后的原理是什么吗? 请跟随北京交通大学国家级教学名师陈后金教授走进"信号与系 ...
- 实验(2)信号与系统的时域分析
本文的pdf文件:link 一.实验目的 建立线性时不变离散系统的差分方程和系统输入序列的数学模型,产生输入序列:利用matlab信号处理工具箱的差分方程求解库函数设计程序,求解系统的单位脉冲响应,给 ...
- matlab已知ft求f2t,基于Matlab的信号与系统实验指导2
基于Matlab 的信号与系统实验指导 实验一 连续时间信号在Matlab 中的表示 一.实验目的 1.学会运用Matlab 表示常用连续时间信号的方法 2.观察并熟悉这些信号的波形和特性 二.实验原 ...
- 离散时间信号和系统的时域分析
实验目的 加深对离散时间信号和系统的理解: 掌握运用MATLAB求解序列卷积和相关的方法: 掌握运用MATLAB对离散时间信号和系统进行时域分析的方法. 实验内容 1.已知以抽样频率对上述3个信号进行 ...
- 手把手MATLAB 离散信号表示 指数 正弦 单位阶跃序列的表示
MATLAB 离散信号表示 指数 正弦 单位阶跃序列的表示 序号 内容 1 指数序列 2 正弦序列 3 单位阶跃序列 一,指数序列 1,指数序列的数学表达式: f ( t ) = A e α t f( ...
- 基于MATLAB的信号与系统实验指导,《信号与系统及MATLAB实现》实验指导书
内容简介: <信号与系统及MATLAB实现>实验指导书 前 言 长期以来,<信号与系统>课程一直采用单一理论教学方式,同学们依靠做习题来巩固和理解教学内容,虽然手工演算训练了计 ...
- 信号与系统matlab pdf,信号与系统及MATLAB实现.pdf
<信号与系统>实验指导书 <信号与系统及MATLAB 实现> 实验指导书 信息工程学院实验中心 2009 年3 月 1 <信号与系统>实验指导书 前 言 长期以来, ...
最新文章
- Detail-Preserving Pooling in Deep Networks
- 初学Flink,对Watermarks的一些理解和感悟(透彻2)
- 求封闭曲线面积vc代码_圆锥曲线综合5个类型,逐一突破
- 面试官问我:如何设计 QQ、微信等第三方账号登陆 ?
- mysql优化 坑_mysql之我们终将踩过的坑(优化)
- mysql 8.0数据备份恢复_RDS for MySQL8.0物理备份恢复到本地自建数据库
- Java实现Redis的订阅发布功能,亲测可以
- 微信小程序后端Java接口开发
- 禅道项目管理软件配置及使用
- cad高程点怎么抽稀_地形图缩编1:500变成2000,植被符号抽稀功能,
- 单片机c语言编程乘法运算符,单片机C语言编程常用的运算符介绍
- AutoCAD DWG 读模块参数:处理块
- VS间接寻址级别不同
- cad立面索引符号 规范_cad立面索引符号怎么画
- 如何计算机画函数图形,Excel中怎么绘制出好看的函数图像
- python语音通话_如何优雅的用Python玩转语音聊天机器人
- 11:Hive从0到1系列学习:函数之常用的内置函数
- 行波iq调制器_行波型LiNbO3电光调制器的电极优化设计
- 【高等数学如何学,做题方法,期末突击】
- Win10 或者Win7 64位系统安装32位的CAD2010教程