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

mytomorrow

常用离散时间序列的Matlab产生:https://www.eda365.com/thread-271002-1-1.html

5 B: f9 L* O5 U7 Y4 P# g' A! @2 V上篇文章介绍了例如单位样值序列、单位阶跃序列产生的脚本以及函数，这篇博文利用到的函数，将直接贴过来：

单位样值序列：
% G. @% u  ^& b
5 r; v4 j6 }! t+ @. k; {

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

6 W* _- v* C5 o3 {( i# D) b单位阶跃序列：



" b. X1 H9 P$ V

• 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];
  $ p2 }7 X" Q# O# y% O

这篇文章给出几个序列，请使用脚本实现：
& h- W; c9 V5 |: J! t4 J# \
5 W) y" l% [4 T% }9 b+ |例1：
5 q6 B- I. {! d+ H7 G* `! r, B7 ^
7 y: z3 F: L% v

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

: F$ G! E# m* |运行即可得到结果：
+ ^1 t, z4 G7 G3 @# x2 V
) H) B8 f4 a& y3 ~- x/ h

) g' N+ G# j& r% ~/ E. l例2：

& A6 m1 g8 M1 V" |. ^9 u. x9 l7 R
7 _" X: j0 `4 @7 l+ e+ z

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

  O4 @& C2 t7 |6 G
. L  c6 _. ~7 O( y! v
% I; }% b* f1 h' z$ R1 B例3：

4 R) Z" i# k- ?+ V9 B! r9 w
4 U' k4 P' D/ W, j" Z* \
0 P) L  h( U8 m, F; v其中，w(n)是均值为0，方差为1的高斯随机序列。

6 m2 ^" `% j5 \4 i. |+ ~

• 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)');
  % ~+ y3 ~) k& A: A! K2 S6 U

5 V9 K! p; P+ Q$ J3 I! D) t" K
% _' `3 O( v# O7 K% a- B" b( `# C

3 A0 Y; w/ o1 Q$ r' x- x4 e; M例4：


3 t- C9 r* W- B2 Z

• 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)');
  " d3 t, `9 s& R/ y# n; r# U* s

1 a8 r7 C2 i) r' \  y

- D2 j# h6 J6 c. U

yin123

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