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

mytomorrow

: V  `1 b6 {9 R6 V" j5 M% Z
常用离散时间序列的Matlab产生:https://www.eda365.com/thread-271002-1-1.html

5 A' F" a" Z. |上篇文章介绍了例如单位样值序列、单位阶跃序列产生的脚本以及函数，这篇博文利用到的函数，将直接贴过来：
9 F$ E/ u* Z9 i6 g/ t
单位样值序列：
, s! |' e; s, q0 |3 Z. f( i. |, Q8 p


8 D3 B, L, }9 Y. X  W5 _. R+ R

• 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];
• 
  8 j6 V6 j0 u, |' w+ G  p& Y! k

7 w! d6 B& E/ k2 F* Q0 f单位阶跃序列：

6 a- w0 l. L8 K4 ?2 s
+ \1 c- G3 r8 T$ Z# c6 _5 R6 K

• 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];
  

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

: j4 R* }4 Y8 ~, h例1：
3 \' ^" y# C. F) l
- z( W- u- V6 l
& Y4 x* L9 _- q) E! w$ r
, W/ s+ Z' H) R5 Q0 {. b' z7 K

• 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]);
  

运行即可得到结果：

6 p8 y) X# [- w+ [
* T5 R& a# d% ?' `, e
6 O5 w5 a% G( M+ u0 X& Y% U例2：
  B  h! b5 Z7 @% {, k! M' s( s
, x; Q2 B# d4 `* N' R
( }4 `! n5 n# P& C

• 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]);
  4 v* f- `+ T2 S1 D: @1 H/ q2 p

' }  H& z' \6 X) }- R


! N0 K6 Q% F  T: ^. \3 V例3：

) V! G0 i) \. L1 A9 G$ I% V5 m0 R
( W% O; p9 ?4 U6 N
其中，w(n)是均值为0，方差为1的高斯随机序列。

• 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)');
    }2 H4 h+ }# A' W

0 I( i0 G; w. T' q0 ~
* b4 Q. s  Z; T! |% B, c4 i
' K5 L  A  ]; Z5 W
例4：



5 a4 F7 }! i  a; v

• 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)');
  

) c, `+ i: m2 N0 h) @4 q3 i. p

yin123

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