|
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
: b- [1 [6 n1 S) l0 B单位样本序列
1 C( \+ k# r% s7 C4 z1 n4 W3 d0 r% y- y: L( o( O# L1 Z& {
, d) A [9 y9 @4 v# H7 e6 ~
# R! B. Y1 K. q- O9 @( K8 I
- clc
- clear
- close all
- n1 = 0;
- n2 = 5;
- n0 = 3;
- n = [n1:n2];
- x = [(n - n0) == 0];
- stem(n,x,'filled');
- ylim([-1,2]);
: S5 X8 _2 h, O% X 5 N2 {9 v4 p9 b2 O& T3 H" T7 m2 s
' g# S& s- @" \+ d8 J' m, O6 I" P7 j* W h: D1 H f7 G8 d8 D
改成一个函数:
% E7 @2 n. k, ]5 S4 {
3 N! t+ ~: C: K" Q6 F1 ^( ]- 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];
; J2 D( r+ l1 O5 H/ P" r- L
. p6 M$ g! l ]4 n( F命名为delta.m" x* @' U; `9 M6 W4 E1 _' p
% Z. Y# \+ l5 r) R3 w0 {# r新建一个脚本:
1 i" y8 E/ r& _* g( ^4 n
% C/ w: M8 c! p9 F- n0 = 3;
- n1 = -3;
- n2 = 6;
- [x,n] = delta(n0,n1,n2);
- stem(n,x,'filled');
- ylim([-1,2]);. t6 D2 m# s) X( w$ R3 G6 H, S) V9 w
& @$ s" R& v* ?0 A7 P$ [$ E
运行得到:9 L H5 A# R/ X; s9 ~
5 q6 d q0 U8 F+ R
; U9 c. _# J% N
3 X6 d8 h4 W3 x- i4 x% Q( t
1 _" z7 v: O: i% Z0 W1 Y单位阶跃序列
; Z( y9 z2 v6 p3 t% [! ~8 z! Y$ L. }9 I9 h* h- J4 W* G+ W
G& g1 `, s7 D3 h9 m
5 d6 @% |3 O. Q- q% G3 i
直接写成函数形式吧:
' k& c- g- `; J5 U
0 b! a4 M# H+ f, P! A- 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];: p% `3 q! g0 U( l
+ L' |$ K1 p' P5 m6 M新建一个脚本,测试:7 K$ v% Y- K( Q: i4 H5 a) ]% H
% U- K( u8 V \- clc
- clear
- close all
- n0 = 3;
- n1 = -3;
- n2 = 6;
- [x,n] = stepseq(n0,n1,n2);
- stem(n,x,'filled');
- ylim([-1,2]);9 c7 z/ ^3 D2 b3 X! I
" v/ ]& {9 s) k! g4 t结果:& B) o( H0 ~, J( j: B
6 Q) s( B) `4 \) c
+ j3 | d; W% v9 Y8 K# A7 J! i
- N7 h# r w/ B- M; @. X+ Y2 A) q: E* k2 }
实值指数序列8 ?$ O1 L; f6 o4 e- g) S
V. Q1 Y* R1 Q5 f, }$ }
为一般形式。
) y$ L, I% \# x% u# `& w' f m! Q0 b* Y/ ]* q- \
下面产生一个序列:
% G) ?$ R+ p* g8 _" j, y
6 K+ j% |( z K; ]: {
- E0 I1 c5 {; [4 S/ A4 F
; p* i1 }) S/ y* J/ F" x2 H脚本代码:
, @6 K3 T' {5 X6 k4 j+ L
3 j; \: }% E+ s$ ?2 N, p* N- 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]);
, h1 i0 E$ l' n3 F0 P & r* _$ p0 v C0 }2 P/ }6 H
5 x" C9 E/ I: r+ f2 D0 E4 v
& b+ ?; e9 n) K
8 i' V# z0 ^1 ~$ i W7 F3 ]复指数序列# l& e/ z# u+ T8 V0 P! o7 b
% T, O: j4 ]$ {9 R. @% e
; z6 ]! w% |& h- t; d0 T$ M
% a: I& K5 L( f- V! }/ G/ e- clc
- clear
- close all
- n = 0:10;
- x = exp( (2+3j)*n );
- stem(n,x);
- xlabel('n');
- 3 D. `+ [( N8 q6 @9 c, c9 _. L
+ F% z; v5 t% z# @% C, l
* t# ~' P( U! P( L9 K0 d. M0 X9 x9 B8 |
# J3 B# [) G8 f- U! D# T周期序列
9 j0 W! u1 C/ y6 |6 e1 N/ i& e1 L
5 A& y3 [+ |7 F+ i- 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');
9 \; @5 f* S0 b& k& l
* ]# h: n, M' t7 `
- J& O5 s8 V8 S& y" D+ V$ Q0 F* ?9 D: w$ J: F& Y2 d7 B C" a4 e8 y- S
* M/ h% T! t1 L# h |
|