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

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

[复制链接]

该用户从未签到

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

EDA365欢迎您登录!

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

x

9 i* v/ z( P# b- m1 k) f2 L" S8 @单位样本序列
+ S8 ]3 l  H7 V) s( Z3 u/ x1 h3 z4 K. F/ h0 P( ?
1 E2 P. V. T4 v0 ?. ~
0 G5 r2 Q" |( P3 z  I6 c
  • clc
  • clear
  • close all
  • n1 = 0;
  • n2 = 5;
  • n0 = 3;
  • n = [n1:n2];
  • x = [(n - n0) == 0];
  • stem(n,x,'filled');
  • ylim([-1,2]);% X+ G7 E" H$ L. g+ `& s$ J
  
. m$ A1 R- C2 X. J+ v 6 u2 k& O0 A! v4 h
6 r# E3 z/ ~$ g$ L( X% h& W: _
改成一个函数:
8 I7 ~5 d! W8 R; O* U4 z2 Z6 H! R* h4 \
  • 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];% X2 b6 Q# j& b; f% q
3 p0 i# y3 p$ T
命名为delta.m" w" a7 w2 r( E6 c' x  ?
3 z! ?4 X4 ]- g( f6 t
新建一个脚本:5 c% [8 Y9 N. u* M- l

5 q9 D, g  ^2 [9 ~" Q: d8 V
  • n0 = 3;
  • n1 = -3;
  • n2 = 6;
  • [x,n] = delta(n0,n1,n2);
  • stem(n,x,'filled');
  • ylim([-1,2]);2 r) O# K, x2 h* h
  0 M7 o8 N& d. q6 m- E
运行得到:9 T. D7 A! I# y& `% p1 Q' d
/ Z/ ]3 t! U* S2 d1 ?
" {: U& H7 z, `" g9 |) u6 ^5 B
- T/ Y7 D; Z7 H3 v6 L4 |/ g+ K1 z
0 u/ t( S! I1 n) H
单位阶跃序列- Q$ u( m7 g' U6 u7 ]* \& F
1 ~  k3 A3 [  k2 y6 e% b# Y1 K

' m& u# \+ m# q% A$ }" }, v. ]/ L5 \/ _. q
直接写成函数形式吧:7 ]+ }2 K! q4 ?. t2 ~1 f( Z' F
! N/ {* v1 Q" s1 ~" R
  • 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];
    7 f% G% ]0 U8 [$ Q- b

9 P0 S! I0 Y, K( }新建一个脚本,测试:/ v, ~$ |! m" C' L; i
0 @3 U! i! K: `
  • clc
  • clear
  • close all
  • n0 = 3;
  • n1 = -3;
  • n2 = 6;
  • [x,n] = stepseq(n0,n1,n2);
  • stem(n,x,'filled');
  • ylim([-1,2]);
    & w' N2 ~4 Y. e: z) Y+ {  I
   7 C& Q# Q1 N8 F9 w: ]6 E1 y
结果:) X) O/ h. m; Y

: h* V. W, F- i1 ~# c, G, K 9 Y- |4 X3 ?. N! E
. ?' @9 j/ N+ p
# @4 Z+ P8 U( X
实值指数序列
& O! K+ R" ^1 W, ]6 ]- M* f( z  l2 ^
为一般形式。
5 E+ b) u; b0 l  |2 S7 r, x* a  r2 u( l; D+ `0 u) g0 j, ~& I
下面产生一个序列:" E- o# L( @. ^+ _% S/ E3 q; J
8 y; X. d1 y' r1 C/ }

# k/ G# l- K7 J3 Z' Y0 q( S* i; D+ f+ G
脚本代码:8 D2 ~& c! O# a

) g& @  ?) x% `, k% }, x
  • 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]);8 K& K8 Q+ a7 f# ?( N+ I7 B
  : F9 w6 {% X! [" T( W, y

8 T0 B6 }8 o  p9 f  H. h, y0 r
: }; v# }& R8 Z/ `+ V- y
2 ~# ?/ v8 s4 D+ `/ Q* \' x- A+ ]6 o复指数序列
* y- d, `+ K; `8 t$ O/ S. m
0 m3 {6 l$ n7 L) [$ o# }& \
! F2 E  X  v3 }' S' ]
; a7 e; E2 m9 n
  • clc
  • clear
  • close all
  • n = 0:10;
  • x = exp( (2+3j)*n );
  • stem(n,x);
  • xlabel('n');
  • / v8 k2 c$ {, x! z7 |2 @7 i. D
    , x9 p2 D& J& O& Q  y
& S3 d- B7 B, w3 ]; q; l" y/ A

' R7 \+ {" D5 ?
7 B+ _" D) z3 K* o' l! ^周期序列
' n; C1 K) O. m( ~5 n4 ]* a7 u8 t; X1 e2 [1 _  ~
  • 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');
    / I2 t0 e4 z: B# N! M! \8 r$ \
     & Q& Z7 q  N! d# C0 \3 Q( l

0 S) w) X% q  Y2 F2 t; a2 I$ C& k3 G' e/ T4 d0 P' K
  O* Z6 E& r6 M; ]; r) D  A

该用户从未签到

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

本版积分规则

关闭

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

EDA365公众号

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

GMT+8, 2025-11-23 20:38 , Processed in 0.171875 second(s), 26 queries , Gzip On.

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

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

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