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

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

[复制链接]

该用户从未签到

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

EDA365欢迎您登录!

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

x

3 R8 J8 G6 d- y7 E& E9 I% S; X. x+ X单位样本序列) |. L6 n$ Z. K( Q, Q. H  L0 c
4 Z* v0 u- o- l( g
+ f. v( J7 ?& x4 e  }

+ \% N; W6 M7 C6 P
  • clc
  • clear
  • close all
  • n1 = 0;
  • n2 = 5;
  • n0 = 3;
  • n = [n1:n2];
  • x = [(n - n0) == 0];
  • stem(n,x,'filled');
  • ylim([-1,2]);
    # v7 ?% J" S. q) D8 {
  3 ~$ l( D/ U" y5 x1 i5 u
7 ~. O: x5 @/ Z0 k; {6 e

' L+ k1 q8 L: `5 y# q改成一个函数:
' M4 x( v; \; _- m) k1 e( h$ d6 h4 |5 p5 R$ |  A. u/ X
  • 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];
    4 U! y, p1 {0 i2 @9 p

' A  r6 d9 ~# Q2 d* Z& Y' T命名为delta.m& B7 H2 ]. Z( l. P) d

5 r6 U7 v- l% ]9 n( `4 W; E新建一个脚本:8 N* G$ x$ L3 X8 }/ h; h, ~
1 E. F! S; J2 h; q# k2 R# F6 g1 O
  • n0 = 3;
  • n1 = -3;
  • n2 = 6;
  • [x,n] = delta(n0,n1,n2);
  • stem(n,x,'filled');
  • ylim([-1,2]);
    2 G2 h& _. T) f$ y, N
  , h$ I) i4 L" }# {0 b. X( i
运行得到:
) W- c) t# d+ o7 [' ^/ G  d  }* a5 j: G( M* U3 ^" Y# Q
; ?$ c+ K! N1 ]4 ]% r1 e
) p; }1 x. e* G- f

, v% u4 P. H0 l5 }$ q单位阶跃序列
) j% t& I  e! w4 |6 D
9 W3 l; @  ~( _6 H   R5 j4 L: i% I7 \1 O5 Y

! j3 b# `% w2 H+ F直接写成函数形式吧:+ E# g* f2 b) m
( }% C" `" T2 S! {# 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];
    2 E/ i1 h  w+ z$ z+ Y

  {- X: O" p, K" J5 [, j5 P新建一个脚本,测试:
: C) L! R) F! p  R6 [/ @* c. v2 }$ C5 o4 S$ e( i4 r
  • clc
  • clear
  • close all
  • n0 = 3;
  • n1 = -3;
  • n2 = 6;
  • [x,n] = stepseq(n0,n1,n2);
  • stem(n,x,'filled');
  • ylim([-1,2]);
    ' T$ ^7 \. ^& L3 M1 ?  q9 z# d
   9 h5 G# R( g$ t- P, u0 m
结果:
4 O/ l2 X' u9 s
: M9 {0 X+ p. C& J1 `; a+ |9 C . Q# x( \$ @, v3 v0 A, `
. }6 y5 u- P3 j% e. W/ R6 E( b

" M& M7 a+ U( i& m; E4 @实值指数序列
+ E- w4 `0 s6 s+ w* E0 D  D% V+ V% J6 o4 P. C4 h6 f
为一般形式。
  D# g) C$ P  Z' s; z
& d- ~+ _" d8 Y8 ?8 p3 B4 H7 H下面产生一个序列:8 j- d, w. w5 v" V' q( r% H' }+ E
# N& w; A) k2 E  j

/ w2 h$ E% i* h" ]
9 L( N; y3 I: ]; b5 A脚本代码:
8 q- G+ P. ^$ X% s9 z9 g1 l# ?
+ D+ D  p; n# a! M
  • 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]);+ e( E+ l' T0 _3 X- O6 L1 \! p
  
+ R( Z  K/ c# w, n# f   O# _: v$ F) d

2 [0 O- C4 Q) V& x) V* S6 U
- w4 ]( V; T  c复指数序列
4 H: L2 \# ?# j  L
9 a" T  `# T! b% R! c2 b 5 d4 O) j  ^' z- }8 w% O: d, D
  R7 f( w, b9 v) i6 M+ E6 }2 u
  • clc
  • clear
  • close all
  • n = 0:10;
  • x = exp( (2+3j)*n );
  • stem(n,x);
  • xlabel('n');
  • ( i: k- V6 ?0 a
    : Q; R; L. B0 p# p- a( B' T  Q) ^' d
, G& ]# a/ W5 V$ s" B3 c

+ a' j" l9 J! \, S: X
# K0 e' g, W; ~" C- _3 ~7 u2 w2 ?周期序列
$ }& p1 i" T! D0 O* F- [; a' E! ?1 Y6 C. G2 R8 X" b  D
  • 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');8 [5 O3 N* f$ E) L# k! h+ Z9 _
     - W( t( |1 a/ H$ l

) M. D3 W4 Z. u
& R/ C' O- P  @9 i5 n0 M9 ]$ i
6 H( \. v  \# B* V$ y

该用户从未签到

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

本版积分规则

关闭

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

EDA365公众号

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

GMT+8, 2025-10-5 21:18 , Processed in 0.156250 second(s), 26 queries , Gzip On.

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

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

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