两个序列的卷积和运算的MATLAB实现（2）

mytomorrow

已知下面两个序列：
1 H; c' C6 c; O$ y  `( g  e3 l
4 p( ]7 ~$ {/ c- ?
( N, H3 H! E1 y; M3 D. g' Q

$ t0 M" O# j! v, K6 x
求这两个序列的卷积。
% E( R) h. d% y6 ~+ O
求卷积的函数是conv，但是使用这个函数有个问题，就是下标问题，也就是求卷积之后的元素值的位置。因此，我们必须要定一个起始点和一个结束点。

方法：
  I1 T' n- R+ E; h, v9 m
$ I8 E/ _+ @0 U

. `/ b% t6 r( K# N4 n9 g

是两个有限长序列，二者卷积的起始点和结束点定义为：


& T! [3 G/ A  G+ O
: e% C# S( K; E9 p

测试脚本：
9 K5 J3 [) M3 v$ ~  Y
4 u& ]5 @& U8 ^/ Fclc
clear
close all

' m5 q; @  h+ Q' G, Y; j
nx = -3:3;
x = [3,11,7,0,-1,4,2];
/ a  @0 o5 ?: ?5 e2 U
nh = -1:4;
h = [2,3,0,-5,2,1];
5 \! I5 a: D3 B/ c# l5 C! _
nyb = nx(1) + nh(1);
4 C5 [9 f  k, N, }4 c( B4 {nye = nx(length(x)) + nh(length(h));
ny = nyb:nye;
y = conv(x,h);
0 N3 \1 p9 X* k* }* j# c
subplot(3,1,1);
stem(nx,x);
/ J0 L5 V2 G. E* N# j9 w2 s3 a, atitle('x(n)');

9 f1 f' g: \7 psubplot(3,1,2);
stem(nh,h);
title('h(n)');


3 C. M: `& V; m# i( |3 a  j4 a* psubplot(3,1,3);
+ V/ F- K5 u" l& c! Wstem(ny, y);
title('y(n)');

# z4 w) V' D1 J& d, K2 ~1 l' K
1 d3 e. Z- A; I  A* |# ~/ d7 K
7 Z* c/ w. Z2 ]- c$ z# L- ^* I0 E
! Y$ \4 b  S) A$ L" S. U7 L" ?昨天，这篇博文就到此结束了，可是呢？你不觉得每次卷积时候都要进行求卷积之后得到的卷积值的位置麻烦吗？

9 j( Z+ V7 v& }8 c; j+ _包括上篇：两个序列的卷积和运算的MATLAB实现（1）

- F  N+ A$ L) K3 n2 Y那我们考虑下把两个信号的卷积简单扩展为一个函数conv_m。
; i! U, Y$ @. A6 P0 k- U
- j3 t; F) l) \- B7 |0 a0 @/ {' C5 ]如下：
/ |7 E% n, Z0 r7 X( {
9 ~% I3 A0 o2 s7 v2 ffunction [y,ny] = conv_m(x,nx,h,nh)
' ^) K+ `) z9 p$ B+ C% Modified convolution routine for signal processing
%___________________________________________________
2 y5 s/ A5 b4 ^& h, H! [# K% [y,ny] = conv_m(x,nx,h,nh)
% [y,ny] = convolution result
  k) l! {! V7 v; `% [x,nx] = first signal
7 W" l2 k8 Q* r! l% [h,nh] = second signal
%
" h) G/ Y4 T' _* H2 B! bnyb = nx(1) + nh(1);
# v2 N4 n3 [# Z6 T5 j2 \. @nye = nx(length(x)) + nh(length(h));
# m, [0 N7 E+ _+ g* Eny = nyb:nye;
+ T2 p) i% v! Cy = conv(x,h);
, G! b1 K3 P9 l# [# ~9 P我们在验证下：

clc
: i: N. S$ x! A! p: O9 jclear
close all

" M& y5 o6 m* P
, N5 v6 r7 H3 s" J1 N0 }nx = -3:3;
! h7 r1 _7 \9 p, kx = [3,11,7,0,-1,4,2];
, E% F" T% t1 h
nh = -1:4;
  `8 w2 v4 `" D/ a6 T4 Eh = [2,3,0,-5,2,1];
; ~& e5 Q8 @1 b8 [$ F  J3 o6 O% ~
[y,ny]=conv_m(x,nx,h,nh);
5 n: s* i/ t/ o$ I# Y
/ D" a6 R! T, k" H1 w$ m+ Y7 R  {8 Qsubplot(3,1,1);
( \* |  n* y4 y. f  \! U# Z3 gstem(nx,x);
title('x(n)');

' `, v* C) N6 b3 c) x+ X2 }subplot(3,1,2);
stem(nh,h);
8 a5 ]. s" E- l* _title('h(n)');
8 s8 I' R/ @2 e

subplot(3,1,3);
- n& A* d8 ?8 b6 H" q) @" Y: F% Jstem(ny, y);
title('y(n)');
- o! h# V. d" J" d


5 C6 P# Y7 R; c: l! K
! ^: z6 z5 ]' L$ Y% u  ?
+ A+ m& u. D! W/ S) Q6 l- |

NNNei256

两个序列的卷积和运算的MATLAB实现