常用离散时间序列的Matlab产生

mytomorrow

, C" S  Y1 ]7 t' U单位样本序列

$ P6 k9 c" ]1 i3 \. s+ }
) i( \5 ~5 k  L& T
5 ~1 ?' g5 r  i7 W  M6 v: N. @- g- r. K

• clc
• clear
• close all
• n1 = 0;
• n2 = 5;
• n0 = 3;
• n = [n1:n2];
• x = [(n - n0) == 0];
• stem(n,x,'filled');
• ylim([-1,2]);
  

  

5 Y/ ?4 g& h9 ?0 P2 K) s
改成一个函数：

+ u$ [9 M: w+ n6 l3 a* ~

• 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];
  * G9 x2 E  Z$ z! m3 l7 w# Z1 }

9 s4 U* M. J, L5 z命名为delta.m
7 f! `0 t) _6 y" x& Z; B! a
# p' Y- Q8 ~0 A" Z, L9 z新建一个脚本：

• n0 = 3;
• n1 = -3;
• n2 = 6;
• [x,n] = delta(n0,n1,n2);
• stem(n,x,'filled');
• ylim([-1,2]);
  

  
$ ~4 }& e. Q& f7 ^( I0 r$ H运行得到：

+ H8 E9 h& z1 i9 G9 [
8 M: ]% ?) }3 T/ t8 f. i

单位阶跃序列
0 r8 ]) {5 f' Z. D2 S: K- ~8 ]


6 r  ?; @: H  p3 u: M直接写成函数形式吧：
+ k8 `& V7 W/ B' T% O/ 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];
    D/ O" h2 q, ?! ^+ ~" S

新建一个脚本，测试：

5 k! a9 \( L! M, s$ ~

• clc
• clear
• close all
• n0 = 3;
• n1 = -3;
• n2 = 6;
• [x,n] = stepseq(n0,n1,n2);
• stem(n,x,'filled');
• ylim([-1,2]);
  # o, G/ X  i; c% W

   
% }6 V. e% v4 W; p  ^2 P结果：

3 `5 n7 p5 }% K' e! h4 A& L3 i

$ v: k3 C. a9 B, {& p8 L7 w( I# v# Q
# V% G# }" Z" D, n. ~) [实值指数序列

为一般形式。

3 g* g. G; s: F* ~( D, t下面产生一个序列：
% P1 l3 A- c* A  O7 J" B- p

6 C) Z: k0 Z; v# M! R
脚本代码：

• clc
• clear
• close all
• n = -3:10;
• x = 0.9.^n;
• stem(n,x,'filled');
• xlabel('n');
• ylabel('0.9^n');
• xlim([-5,13]);
• ylim([0,2]);
  % ^, L6 n  x$ |% Q

  



2 C2 M% r: P6 D$ n* O3 ?( H复指数序列
) {3 j6 {$ p9 R% F

; _: k% Y0 R4 |. m' y, S. |9 h
  F& t) f/ T) Q! v6 i& h

• clc
• clear
• close all
• n = 0:10;
• x = exp( (2+3j)*n );
• stem(n,x);
• xlabel('n');
• 
  3 _- a, P" I7 }) n4 R, V

   
8 b9 x% Z' {7 i9 q/ s6 p
' k/ r; S0 Q2 L; L" K* Z/ [5 X
$ E$ Y! }: i( M. L1 k0 O
% x7 z. Q2 j. W  B" b' `5 I' F周期序列

: O9 T* d# C% b$ N6 C2 b  U# }

• clc
• clear
• close all
• n = 0:10;
• x = exp( (2+3j)*n );
• subplot(2,1,1)
• stem(n,x);
• xlabel('n');
• title('Not Period Sequence');
• % the number of period
• P = 5;
• xtilde = x' * ones(1,P);
• xtilde = xtilde(：）';
• subplot(2,1,2);
• stem(xtilde);
• title('Period Sequence');
  

     
# i/ c- b8 s3 P/ K

8 n  w. ~# E3 O5 H# Y9 U$ M: N1 f, C
: {) [; T6 q, D! z% o9 P6 z

ExxNEN

不知道可以用到哪里

xixihahaheihei

收藏一下