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