|
|
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
|
|