信号与系统实验二___MATLAB
离散信号与系统的时域分析
- 离散信号的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(kpi*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.5h(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相关推荐
- 信号与系统 实验二 连续时间LTI系统的时域分析
一.实验目的 1. 学会用MATLAB求解连续系统的零状态响应: 2.学会用MATLAB求解冲激响应及阶跃响应: 3. 学会用MATLAB实现连续信号卷积的方法: 二.实验原理 1.连续时间系统零状态 ...
- 信号与系统matlab课设报告,MATLAB信号与系统实验报告
<MATLAB信号与系统实验报告>由会员分享,可在线阅读,更多相关<MATLAB信号与系统实验报告(9页珍藏版)>请在装配图网上搜索. 1.信号与系统实验报告(5)MATLAB ...
- matlab用lism求零输入响应,信号与系统实验报告
信号与系统实验报告 (14页) 本资源提供全文预览,点击全文预览即可全文预览,如果喜欢文档就下载吧,查找使用更方便哦! 19.90 积分 信号与系统实验信号与系统实验一.实验目的一.实验目的1. 通过 ...
- matlab已知ft求f2t,基于Matlab的信号与系统实验指导2
基于Matlab 的信号与系统实验指导 实验一 连续时间信号在Matlab 中的表示 一.实验目的 1.学会运用Matlab 表示常用连续时间信号的方法 2.观察并熟悉这些信号的波形和特性 二.实验原 ...
- matlab卷积实验原理,信号与系统实验报告——卷积(含程序)
信号与系统实验报告--卷积 电 子 科 技 大 学 实 验 报 告 学生姓名:苏晓菁 学 号:2804301026 指导教师:张鹰 一.实验室名称:信号与系统实验室 二.实验项目名称:离散系统的冲激响 ...
- matlab对正弦采样信号采样重建,信号与系统实验报告4
实验报告 课程名称:信号与系统 实验项目名称:实验4 信号抽样及抽样定理 实 验 目 的: 学会运用MATLAB完成信号抽样及对抽样信号的频谱进行分析:学会运用 MATLAB改变抽样间隔,观察抽样后信 ...
- 离散信号内插matlab,离散信号和系统实验报告.doc
离散信号和系统实验报告 三.实验效果分析(包括仪器设备等使用效果) 实验中注意中英文的切换是应注意,特别是标点符号,括号等需要注意.需要仔细键入程序. 教 师 评 语 指导老师 年 月 日 江西师范大 ...
- 信号与系统实验感想 matlab,信号与系统MATLAB实验报告.doc
信号与系统MATLAB实验报告.doc 成绩 信号与系统实验报告 课 程 名 信号与系统 学 部 机械与电子信息学部 专 业 电子信息工程 学 号 姓 名 曹 禹 指导教师 吴国平 日 期 2013/ ...
- ATM系统-实验二:Use Case图与Use Case详细描述
ATM系统-实验二:Use Case图与Use Case详细描述 实验一文档所在博客:https://blog.csdn.net/aptx1255/article/details/80080013 一 ...
- 基于MATLAB的信号与系统实验指导,《信号与系统及MATLAB实现》实验指导书
内容简介: <信号与系统及MATLAB实现>实验指导书 前 言 长期以来,<信号与系统>课程一直采用单一理论教学方式,同学们依靠做习题来巩固和理解教学内容,虽然手工演算训练了计 ...
最新文章
- logstash配置文件
- union和union all有什么区别_什么是Python Wheels?为什么要关心它?
- cnil在python_在Python中使用cumprod()计算权益曲线
- Linux 常用命令与设置
- ES5-2 语法、规范、错误、运算符、判断分支、注释
- 完整的Web应用程序Tomcat JSF Primefaces JPA Hibernate –第3部分
- 100万用户服务器_我的应用在一个月内如何增长超过100万用户
- ajax清请求过程,JS深入基础之Ajax的请求过程
- python 根据名称获取pid_【Python学习笔记】76、常用第三方模块psutil
- 数字图像处理--图像ROI
- 快速排序思路(Hoare版),代码实现
- mysql修改字段结构_MySQL修改表结构及其添加删除修改字段功能
- Playframework项目启动后自动停止问题记录
- java多线程设计wait、notify、notifyall、synchronized的使用机制
- python 实现SOM: 函数更新
- 维纳滤波python 函数_Python中离焦图像的维纳滤波
- 成都信息工程大学2019级概率论c类试题
- 【游戏运营】【笔记】 谈谈对游戏运营的理解
- <Halcon> 区域划分算子partition_dynamic和partition_rectangle
- 时间序列matlab的实现