数字信号处理中的几个常用序列的MATLAB代码产生

mytomorrow

, [( E7 C8 ]) P' w( _
3 f/ l0 w) G  w" }9 p2 ^) r0 s常用离散时间序列的Matlab产生:https://www.eda365.com/thread-271002-1-1.html

) a+ {' j$ U) s* b上篇文章介绍了例如单位样值序列、单位阶跃序列产生的脚本以及函数，这篇博文利用到的函数，将直接贴过来：

+ n9 t2 C0 l& f: p0 E' q6 H单位样值序列：

$ H" |( C% ^  ]* [

• function [x,n]=delta(n0,n1,n2);
• % generate x(n) = delta(n - n0); n1 <= n <= n2
• %_____________________________________________
• %[x,n] = delta(n0, n1, n2);
• %
• n = [n1:n2];
• x = [(n-n0) == 0];
• 
  

" W5 r7 B5 R7 a1 |2 M单位阶跃序列：
! S+ L9 @* A$ f
/ G% R2 v2 U1 \4 o
, u; U( y! n  q. }" h: \8 a$ G# w

• function [x,n]=stepseq(n0,n1,n2);
• % generate x(n) = u(n - n0); n1 <= n <= n2
• %_____________________________________________
• %[x,n] = stepseq(n0, n1, n2);
• %
• n = [n1:n2];
• x = [(n-n0) >= 0];
  

这篇文章给出几个序列，请使用脚本实现：

例1：
3 T4 [& {; G/ Z' ~" K" p. i% c
8 e: P  J: S% o1 ]
& d' h6 E9 @0 \9 }
; k4 ]+ P6 y# w- v& G

• clc
• clear
• close all
• n = [-5:5];
• x = 2*delta(-2,-5,5)-delta(4,-5,5);
• stem(n,x);
• xlabel('n');ylabel('x(n)');
• ylim([-2,3]);
  % j3 h. ^% x3 l8 E, _

9 D% z; _' h8 Z9 j5 j
运行即可得到结果：
- K+ n# ]5 S: O1 }" ]
  K, I2 v3 w( e* G$ z* C, c
5 L( R$ o( x3 ?: Y
例2：
% ^+ d; T! j. @! T3 g8 ?) |% @

• clc
• clear
• close all
• n = [0:20];
• x1 = n.*(stepseq(0,0,20) - stepseq(10,0,20));
• x2 = 10*exp(-0.3*(n-10)).*(stepseq(10,0,20) - stepseq(20,0,20));
• x = x1+x2;
• stem(n,x);
• xlabel('n');ylabel('x(n)');
• ylim([-1,10]);
  

0 t8 ~1 x2 Y' Y% U. m& ~2 R
例3：

( M. u& i5 `  L6 e( w2 w; ?
4 z& I% O4 x5 _- r
其中，w(n)是均值为0，方差为1的高斯随机序列。
- Z8 I  j2 q, c. p5 Z7 E+ o: r
* ]4 _9 p" p6 S" p6 m

• clc
• clear
• close all
• n = [0:50];
• x = cos(0.04 * pi * n) + 0.2*randn(size(n));
• stem(n,x);
• xlabel('n');ylabel('x(n)');
  " c0 f- X( k2 @! I0 C! G

例4：

• clc
• clear
• close all
• n = [-10:9];
• x = [5,4,3,2,1];
• xtilde = x' * ones(1,4);
• xtilde = xtilde(：)';
• stem(n,xtilde);
• xlabel('n');ylabel('xtilde(n)');
  

5 {5 _; [' K( M+ q& M
9 c* j$ H. G$ v' w

yin123

数字信号处理中的几个常用序列的MATLAB代码产生