离散信号与系统的时域分析

  • 离散信号的MATLAB表示
    • 例2-1正弦信号序列(注意序列信号的周期性)
    • 例2-2单位序列信号δ(k)
    • 例2-3单位阶跃信号ε(k)
    • 例2-4实指数序列信号
  • 离散信号的卷积和
    • 例2-5
  • 离散系统的单位响应
    • 例2-6
  • 离散系统的零状态响应
    • 例2-7
  • 课后习题
    • 2
    • 3

离散信号的MATLAB表示

表示离散信号f(k)需要两个行向量,一个表示序号k = [ ],一个表示相应的函数值f = [ ], 画图命令是stem。

例2-1正弦信号序列(注意序列信号的周期性)

k = 0 : 40;
subplot( 2, 1, 1 ),
stem( k, cos( k * pi / 8 ), ‘filled’ ),
title( ‘cos(kpi/8)’ );
subplot( 2, 1, 2 ),
stem( k, cos( k * pi * 2 ), ‘filled’ ),
title( 'cos(k
pi*2)’ );

例2-2单位序列信号δ(k)

编写函数δ(k)m文件

function danweixulie( k1, k2, k0 )%k1、k2为画图范围,k0为脉冲位置
k = k1:k2;
n = length( k );
f = zeros( 1, n );%生成零矩阵
f( 1, -k0 - k1 + 1 ) = 1;%对脉冲位置赋值为1
stem( k, f,‘filled’)
axis( [ k1, k2, 0, 1.5 ] )
title( ‘单位序列δ’)

实现函数danweixulie

danweixulie( -5, 5, 0)

例2-3单位阶跃信号ε(k)

编写函数ε(k)m文件

function jieyuexulie( k1, k2, k0 )
k_zero = k1: -k0-1;
k_one = -k0 : k2;
y_zero = zeros( 1, length( k_zero ) );
y_one = ones( 1, length( k_one ) );
stem(k_one,y_one)
hold on
stem(k_zero,y_zero)
hold off
title( ‘ε(k)’ )
axis( [ k1, k2, 0, 1.5 ] )

实现函数jieyuexulie

jieyuexulie( -3, 8, 0 )

例2-4实指数序列信号

实指数序列m文件

function zhishuxulie( k1, k2, c, a )%k1为起始,k2为终止,c幅度,a底数
k = k1 : k2;
x = c * ( a.^k );
stem( k, x,‘filled’)
hold on
plot( [ k1, k2 ], [ 0, 0 ] )
hold off

实现函数zhishuxulie

zhishuxulie( 0,40,1,5/4 )

离散信号的卷积和

两个有限长的序列f1,f2卷积可调用matlab里的conv,格式f=conv(f1,f2),f为卷积结果,但不显示时间序列,可自己编函数dconv将其显示出来

例2-5

编写dconv.m

function [f,k]=dconv(f1,f2,k1,k2)
%The function of compute f=f1f2
% f: 卷积和序列 f(k)对应的非零样值向量
% k: 序列 f(k)的对应序号向量
% f1: 序列 f1(k)非零样值向量
% f2: 序列 f2(k)的非零样值向量
% k1: 序列 f1(k)的对应序号向量
% k2: 序列 f2(k)的对应序号向量
f=conv(f1,f2) %计算序列 f1 与 f2 的卷积和 f
k0=k1(1)+k2(1); %计算序列 f 非零样值的起点位置
k3=length(f1)+length(f2)-2; %计算卷积和 f 的非零样值的宽度
k=k0:k0+k3 %确定卷积和 f 非零样值的序号向量
subplot(2,2,1)
stem(k1,f1) %在子图 1 绘序列 f1(k)时域波形图
title(‘f1(k)’)
xlabel(‘k’)
ylabel(‘f1(k)’)
subplot(2,2,2)
stem(k2,f2) %在图 2 绘序列 f2(k)时波形图
title(‘f1(k)’)
xlabel(‘k’)
ylabel(‘f2(k)’)
subplot(2,2,3)
stem(k,f); %在子图 3 绘序列 f(k)的波形图
title(‘f(k)f1(k)与 f2(k)的卷积和 f(k)’)
xlabel(‘k’)
ylabel(‘f(k)’)
h=get(gca,‘position’);
h(3)=2.5
h(3);
set(gca,‘position’,h) %将第三个子图的横坐标范围扩为原来的 2.5 倍

运行

f1=[1 2 1];
k1=[-1 0 1];
f2=ones(1,5);
k2=-2:2;
[f, k]=dconv(f1,f2,k1,k2)

  • 由运行结果和代码可以看出卷积f(k)的起点为f1与f2起点之和,终点为f1与f2终点之和,因此长度为f1与f2长度之和减1 *

离散系统的单位响应

MATLAB中系统响应函数impz

impz( b ,a )%b,a是离散系统的行向量
impz( b, a, n )%b,a是离散系统的行向量,n表示时间范围0~n
impz( b, a, n1, n2 )%时间范围n1~n2
y = impz( b, a, n1, n2 )%由y给出数值序列

例2-6

a=[1,-1,0.9];
b=[1];
impz(b,a)
%impz(b,a,60)
%impz(b,a,-10:40)

离散系统的零状态响应

例2-7

a = [ 1 -0.25 0.5 ];
b = [ 1 1 ];
t = 0 : 20;
x = ( 1/2 ).^t;
y = filter( b, a, x )%b为激励系数,a为响应系数,x为激励函数
subplot( 2, 1, 1 )
stem( t, x )
title( ‘输入序列’ )
subplot( 2, 1, 2 )
stem( t, y )
title( ‘响应序列’ )

课后习题

2

a=[ 2, -2, 1 ];
b=[ 1, 3, 2 ];
impz(b,a);

3

a = [ 1, 1, 0.25 ];
b = [ 1 ];
t = 0 : 15;%范围
x = ( t >= 0 );
y = filter( b, a, x )%b为激励系数,a为响应系数,x为激励函数
subplot( 2, 1, 1 )
stem( t, x )
title( ‘输入序列’ );
subplot( 2, 1, 2 )
stem( t, y )
title( ‘响应序列’ );

信号与系统实验二___MATLAB相关推荐

  1. 信号与系统 实验二 连续时间LTI系统的时域分析

    一.实验目的 1. 学会用MATLAB求解连续系统的零状态响应: 2.学会用MATLAB求解冲激响应及阶跃响应: 3. 学会用MATLAB实现连续信号卷积的方法: 二.实验原理 1.连续时间系统零状态 ...

  2. 信号与系统matlab课设报告,MATLAB信号与系统实验报告

    <MATLAB信号与系统实验报告>由会员分享,可在线阅读,更多相关<MATLAB信号与系统实验报告(9页珍藏版)>请在装配图网上搜索. 1.信号与系统实验报告(5)MATLAB ...

  3. matlab用lism求零输入响应,信号与系统实验报告

    信号与系统实验报告 (14页) 本资源提供全文预览,点击全文预览即可全文预览,如果喜欢文档就下载吧,查找使用更方便哦! 19.90 积分 信号与系统实验信号与系统实验一.实验目的一.实验目的1. 通过 ...

  4. matlab已知ft求f2t,基于Matlab的信号与系统实验指导2

    基于Matlab 的信号与系统实验指导 实验一 连续时间信号在Matlab 中的表示 一.实验目的 1.学会运用Matlab 表示常用连续时间信号的方法 2.观察并熟悉这些信号的波形和特性 二.实验原 ...

  5. matlab卷积实验原理,信号与系统实验报告——卷积(含程序)

    信号与系统实验报告--卷积 电 子 科 技 大 学 实 验 报 告 学生姓名:苏晓菁 学 号:2804301026 指导教师:张鹰 一.实验室名称:信号与系统实验室 二.实验项目名称:离散系统的冲激响 ...

  6. matlab对正弦采样信号采样重建,信号与系统实验报告4

    实验报告 课程名称:信号与系统 实验项目名称:实验4 信号抽样及抽样定理 实 验 目 的: 学会运用MATLAB完成信号抽样及对抽样信号的频谱进行分析:学会运用 MATLAB改变抽样间隔,观察抽样后信 ...

  7. 离散信号内插matlab,离散信号和系统实验报告.doc

    离散信号和系统实验报告 三.实验效果分析(包括仪器设备等使用效果) 实验中注意中英文的切换是应注意,特别是标点符号,括号等需要注意.需要仔细键入程序. 教 师 评 语 指导老师 年 月 日 江西师范大 ...

  8. 信号与系统实验感想 matlab,信号与系统MATLAB实验报告.doc

    信号与系统MATLAB实验报告.doc 成绩 信号与系统实验报告 课 程 名 信号与系统 学 部 机械与电子信息学部 专 业 电子信息工程 学 号 姓 名 曹 禹 指导教师 吴国平 日 期 2013/ ...

  9. ATM系统-实验二:Use Case图与Use Case详细描述

    ATM系统-实验二:Use Case图与Use Case详细描述 实验一文档所在博客:https://blog.csdn.net/aptx1255/article/details/80080013 一 ...

  10. 基于MATLAB的信号与系统实验指导,《信号与系统及MATLAB实现》实验指导书

    内容简介: <信号与系统及MATLAB实现>实验指导书 前 言 长期以来,<信号与系统>课程一直采用单一理论教学方式,同学们依靠做习题来巩固和理解教学内容,虽然手工演算训练了计 ...

最新文章

  1. logstash配置文件
  2. union和union all有什么区别_什么是Python Wheels?为什么要关心它?
  3. cnil在python_在Python中使用cumprod()计算权益曲线
  4. Linux 常用命令与设置
  5. ES5-2 语法、规范、错误、运算符、判断分支、注释
  6. 完整的Web应用程序Tomcat JSF Primefaces JPA Hibernate –第3部分
  7. 100万用户服务器_我的应用在一个月内如何增长超过100万用户
  8. ajax清请求过程,JS深入基础之Ajax的请求过程
  9. python 根据名称获取pid_【Python学习笔记】76、常用第三方模块psutil
  10. 数字图像处理--图像ROI
  11. 快速排序思路(Hoare版),代码实现
  12. mysql修改字段结构_MySQL修改表结构及其添加删除修改字段功能
  13. Playframework项目启动后自动停止问题记录
  14. java多线程设计wait、notify、notifyall、synchronized的使用机制
  15. python 实现SOM: 函数更新
  16. 维纳滤波python 函数_Python中离焦图像的维纳滤波
  17. 成都信息工程大学2019级概率论c类试题
  18. 【游戏运营】【笔记】 谈谈对游戏运营的理解
  19. <Halcon> 区域划分算子partition_dynamic和partition_rectangle
  20. 时间序列matlab的实现

热门文章

  1. 2022最新阿里Java面经,转疯了
  2. 小米路由r2d论坛_维修一台亮黄灯的R2D小米路由器2
  3. python图片处理教程_Python图像入门教程:Python图像入门教程推荐
  4. ImageJ使用手册介绍
  5. 前端使用Aliplayer 播放器 播放rtmp直播流
  6. 随机森林算法的扩展-深度森林(deep forest)
  7. 【网络】TCP报文头部
  8. java代码混淆-ProGuard
  9. 神经动力学模型的建立
  10. leetcode算法88.合并两个有序数组