【 MATLAB 】两个序列的卷积和运算的MATLAB实现（1）

设矩形脉冲 $x(n) = u(n) - u(n -10)$ 是脉冲响应 $h(n) = (0.9)^nu(n)$ 的LTI系统的输入，求输出 y(n).

下面的脚本中用到了一个自定义的函数，也就是两个信号相加的函数：

相关博文：【 MATLAB 】基本序列运算及其MATLAB的等效表示

function [y,n] = sigadd(x1,n1,x2,n2)
% implements y(n) = x1(n) + x2(n)
% [y,n] = sigadd(x1,n1,x2,n2)
%——————————————————————————————
% y = sum sequence over n, which includes n1 and n2
% x1 = first sequence over n1
% x2 = second sequence over n2( n2 can be different from n1)
%
n = min( min(n1), min(n2) ):max( max(n1), max(n2) ); %duration of y(n)
y1 = zeros(1,length(n)); y2 = y1; %initialization
y1( find( ( n >= min(n1) )&( n <= max(n1) ) == 1  )  ) = x1; %x1 with duration of y1
y2( find( ( n >= min(n2) )&( n <= max(n2) ) == 1  )  ) = x1; %x2 with duration of y2
y = y1 + y2;

直接给出MATLAB脚本：

clc
clear
close all% help stepseq
%   generate x(n) = u(n - n0); n1 <= n <= n2
%  _____________________________________________
%  [x,n] = stepseq(n0, n1, n2);
[u1,n1] = stepseq(0,-5,45);
[u2,n2] = stepseq(10,-5,45);% generate signal x(n)
[x,n] = sigadd(u1,n1,-u2,n2);% generate signal h(n)
m = -5:45;
h = ( (0.9).^m ).* u1;% the convolution of x(h) and h(n)
y = conv(x,h);
% ensure the index
nyb = n(1)+ m(1);
nye = n(length(x)) +n(length(h));
ny = nyb:nye;subplot(3,1,1);
stem(n,x);
title('x(n)');
xlabel('n')subplot(3,1,2);
stem(m,h);
title('h(n)');
xlabel('n')subplot(3,1,3);
stem(ny,y);
title('the conv of x(n) and h(n)');
xlabel('n')
xlim([-5,45]);

【 MATLAB 】两个序列的卷积和运算的MATLAB实现（1）相关推荐

【 MATLAB 】两个序列的卷积和运算的MATLAB实现（2）
已知下面两个序列: 求这两个序列的卷积. 求卷积的函数是conv,但是使用这个函数有个问题,就是下标问题,也就是求卷积之后的元素值的位置.因此,我们必须要定一个起始点和一个结束点. 方法: 是两个有限 ...
序列的卷积运算与相关运算——MATLAB
一.实验目的 1.掌握有限长序列线性卷积的编程计算原理,并能够利用Matlab或C语言编写算法程序进行线性卷积运算的程序实现; 2.学会线性卷积函数和线性相关函数的使用方法,并能利用二者进行有限长序列 ...
卷积、相关（matlab）
卷积.相关本次博客主要是图示化卷积过程,能够进一步加深学者在学习过程中对数学卷积的理解.首先,再次回顾一下利用MATLAB产生指数序列 x[k]=Kαku[k], a=input('a='); K= ...
matlab重叠相加法求卷积,通过重叠相加法实现卷积的报告.doc
数字信号处理课程设计题目:通过重叠相加法实现卷积院系:自动化与信息工程学院专业:通信工程班级: 通信091 学号: 3090432028 姓名: 雷帛川指导教师: 李建勋职称: 副教授 ...
基于matlab实现信号的线性卷积与循环卷积
系列文章目录数字信号处理(DSP:Digital Signal Process)是电子通信领域非常重要的研究方向,博主汇总了数字信号处理(DSP)中常用的经典案例分析,主要基于算法分析.MATLAB ...
matlab作卷积的公式,卷积相关公式的matlab代码
取半径=3 用matlab代码实现上式公式: length=3; for Ki = 1:length for Kj = 1:length for Kk = 1:length Ksigma(Ki,Kj, ...
matlab伪随机码,Matlab生成M序列的伪随机码
图3.5 线性反馈移位寄存器 MATLAB生成M序列伪随机码函数如下: % 在MATLAB命令窗口输入以下: % fbconnection=[0 1 1 0 1 0 0 0 0 0 0 0 0 0 0 ...
matlab由频率响应计算差分方程,现代线性系统：使用MATLAB
中译本出版者的话译者的话出版者的话符号一览表全书内容简介前言第1章信号与序列概述基本概念与解说题信号. 序列和系统 IP1. 1 描述连续时间信号 IP1. 2 序列 ...
基于MATLAB完成卷积运算,基于Matlab的离散卷积
基于Matlab 的离散卷积刘国良 (洛阳理工学院河南洛阳 471000) 摘要:卷积运算广泛用于通讯.电子.自动化等领域的线性系统的仿真.分析及数字信号处理等方面.在Matlab 中可以使用线 ...

【 MATLAB 】两个序列的卷积和运算的MATLAB实现（1）

【 MATLAB 】两个序列的卷积和运算的MATLAB实现（1）相关推荐

最新文章

热门文章