找回密码
 注册
关于网站域名变更的通知
查看: 498|回复: 2
打印 上一主题 下一主题

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

[复制链接]

该用户从未签到

跳转到指定楼层
1#
发表于 2020-1-2 10:49 | 只看该作者 |只看大图 回帖奖励 |倒序浏览 |阅读模式

EDA365欢迎您登录!

您需要 登录 才可以下载或查看,没有帐号?注册

x

6 }! _" g; A$ Y* R) T% i单位样本序列
" @# Y: C7 r5 s- @) H# {+ s0 p7 N& U! \

/ w( Y3 D5 Y( v1 t' }% T; e& r
+ _" C* `5 j: H% P# l+ @5 l- T
  • clc
  • clear
  • close all
  • n1 = 0;
  • n2 = 5;
  • n0 = 3;
  • n = [n1:n2];
  • x = [(n - n0) == 0];
  • stem(n,x,'filled');
  • ylim([-1,2]);. G) R: \. L7 B5 H4 d$ ]" y0 b; I
  2 `& v" M* {2 `& K, g/ J3 z: z' f
1 N0 i/ J. m6 o2 ]

' `7 r9 E0 s5 ^改成一个函数:0 @0 p+ I1 i. k- x2 Z
7 M; c* w5 }, y9 u$ k! ?- O
  • 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];3 }6 E5 t; ^/ t- N6 r7 L0 `7 k
+ h, \4 V: T* @7 D" q6 S( R* I) W
命名为delta.m
2 n  a/ S. w, S% I3 k. H) E% c* h2 B1 w3 T. r8 [7 @
新建一个脚本:
. H6 F! |* [; B1 x1 B% k$ j2 a7 M4 y! @! k+ d
  • n0 = 3;
  • n1 = -3;
  • n2 = 6;
  • [x,n] = delta(n0,n1,n2);
  • stem(n,x,'filled');
  • ylim([-1,2]);
    0 x, v% ~' ]( f9 A6 A6 [0 \5 ~
  
6 E7 F; l% }6 ?( u+ a运行得到:. m2 [6 [& t+ j! t

& @1 O2 @4 w; P0 ?' K* x2 v
, b; v- g0 c5 d- i6 \4 m* [9 y3 A2 v+ _
( \0 ~. x- [; w+ W+ `! @: h1 N6 D
单位阶跃序列
; t+ z: j; P* ?. X# `6 M, s4 n+ a3 Y! x9 ]2 |

( R; k# X3 M8 n- L
0 k8 y8 Y* S2 H6 D0 T; x  n4 y  w直接写成函数形式吧:
7 I- i# ?0 d1 M8 v- b/ M4 `0 C# y8 |9 A; I$ d
  • 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];
    / k: A4 {8 ^/ {
9 B4 i: U( g1 w. A2 [* K- a
新建一个脚本,测试:
* C$ S( l9 N- k$ z3 P# ~  O) |, o% Q4 H+ z8 ?0 j8 s/ o
  • clc
  • clear
  • close all
  • n0 = 3;
  • n1 = -3;
  • n2 = 6;
  • [x,n] = stepseq(n0,n1,n2);
  • stem(n,x,'filled');
  • ylim([-1,2]);* t5 F% d8 d% E1 n! e
   + a* F3 _! K6 L
结果:
4 x5 y0 j1 e. m* V0 d! s( ]5 ]& [
, M5 w/ W! |8 B1 L8 S
" m( |; e' K9 e" ?, D9 M: d% C" l
2 j5 P' |& D  Q& b' j9 Z; x( b1 V8 G1 b6 c  ~! e2 \# R* E
实值指数序列: L! G% j0 {$ }7 u0 ^& E

8 e8 \  G/ J  `% U# y* h 为一般形式。2 E: U7 z6 G0 j$ V

- ]# V: d  W4 H  e, N下面产生一个序列:6 ]# h* z/ c* s! a+ e- M. Q

+ ~2 {" w+ t) r, `; }$ h! F 2 e5 n% l$ G9 H+ X
: M; N! j' l+ ^( R  L' [
脚本代码:
7 Y& B- @* _, K3 ~9 l0 [
7 U) K7 M8 W1 @6 R% s. W) D9 @
  • 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]);& z! v8 l* R: P# C2 a
  
( n/ U% g+ N# F5 M- ]* S$ @ " r3 m9 z5 }( i$ A1 R4 k' k/ C* h: z+ B

$ j$ p, F5 d# g8 L! F+ g% G- [* ^3 D& n3 w. c$ Y) h8 b
复指数序列
6 ?; ^* O- ~4 r5 i
0 u7 {7 Q5 y. X 2 i3 v7 U1 I; N" }

2 b& j. m( I4 O
  • clc
  • clear
  • close all
  • n = 0:10;
  • x = exp( (2+3j)*n );
  • stem(n,x);
  • xlabel('n');
  •   W/ G; L& \& {% _9 E7 M3 l; m
   
) B0 o: C+ H, j+ ~; t+ M) h  v
: f1 @4 [; ~! Z: Y9 c& y7 S! i1 S! \+ T$ i. s  h+ [* [( r4 S

# S& b- D  |1 {* b( A- \0 m周期序列* f3 Y5 K+ O# E9 T# G6 u* W
, V6 ?. n7 N2 y
  • 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');
    3 R4 D1 h: c. o, ?- e
     
3 y/ r: z0 q/ N4 }' r . `( l, X+ E$ C
4 }. u1 Y( d' Q% z: V) I
8 w) F) @1 z, A7 v7 U, e2 Z

该用户从未签到

2#
发表于 2020-1-2 18:43 | 只看该作者
不知道可以用到哪里
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

推荐内容上一条 /1 下一条

EDA365公众号

关于我们|手机版|EDA365电子论坛网 ( 粤ICP备18020198号-1 )

GMT+8, 2025-11-1 20:27 , Processed in 0.156250 second(s), 27 queries , Gzip On.

深圳市墨知创新科技有限公司

地址:深圳市南山区科技生态园2栋A座805 电话:19926409050

快速回复 返回顶部 返回列表