|
|
EDA365欢迎您登录！
您需要 登录 才可以下载或查看，没有帐号？注册
x
8 U" j$ @! `) K/ n. i3 n9 t单位样本序列
& N3 \: U! u" `; U
5 R! K' A( O0 C' V
& }5 s- y0 d' F
; [) E+ p5 {; s- clc
- clear
- close all
- n1 = 0;
- n2 = 5;
- n0 = 3;
- n = [n1:n2];
- x = [(n - n0) == 0];
- stem(n,x,'filled');
- ylim([-1,2]);
5 d, s/ r4 t* L4 D4 G
9 H" }8 V: @; ?* G9 |2 H5 X2 T
5 i- S: C/ U0 }+ y$ v1 ?6 m0 \) a$ g9 r
改成一个函数：
! f6 b) V: V, p p+ W: S
4 B/ S$ n |3 N- V; J! d V- 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];* ?: [# m7 o* x! x
# Z8 W$ O% o* C* C+ w命名为delta.m
: c0 p; v" k$ \0 d$ z6 Z5 `. E, B4 v6 b z2 X' M1 E* [
新建一个脚本：
5 K2 U: v6 q- j. J" h0 T& Y' x4 M3 F% {& k) u0 g# b! E
- n0 = 3;
- n1 = -3;
- n2 = 6;
- [x,n] = delta(n0,n1,n2);
- stem(n,x,'filled');
- ylim([-1,2]);
: I# ?# z b9 `- }: w: F! P+ d4 d
' i* q. X+ m! a' {6 F4 s运行得到：
" j/ @$ t$ F. r4 C2 i$ ?$ ^4 N! c2 ?# V; y4 n7 Q( `
8 Y# }4 }, t0 E# }' j0 I7 D( d/ H8 N) v, O& o
. j) H/ N( K) ~3 s3 `单位阶跃序列- V' C( [& a, @$ g' K$ _3 H& g
- n, n* L1 G2 D% ]8 ]8 g
& s3 v7 q9 r; q* K$ ?
0 x8 s9 A7 A) t# R直接写成函数形式吧：" F9 F9 |$ \8 T
) R/ d2 k2 S. W# z/ e: `) r! L9 _
- 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];0 l0 A, g8 T% `" A
2 x! g% R0 C" r# s, v" T, r- e
新建一个脚本，测试：" M1 O1 z; F9 G3 E; ^8 p5 O, a1 n/ \
y( ~5 S% u2 P$ a
- clc
- clear
- close all
- n0 = 3;
- n1 = -3;
- n2 = 6;
- [x,n] = stepseq(n0,n1,n2);
- stem(n,x,'filled');
- ylim([-1,2]);! }0 C1 N+ r/ K" Q7 t
2 |* r6 O, o& d
结果：
1 ~! G9 z& ^6 Y8 z% \# K. N: T3 y0 `; d
- _9 x4 b" Y! J$ |# k l
+ p+ q" D* r* B
0 w7 I1 a& d$ y' n8 [5 W实值指数序列1 O. t* d H6 y8 U2 E1 ]
5 d7 }7 ~8 H5 K- Z
为一般形式。
! {2 ?* n! C$ k. Q$ z2 [7 `( J
6 X* `- |4 P) z5 K* L; n) e下面产生一个序列：
6 q0 ^. _7 o; x, w/ H6 Q. R# Q& I0 e3 n1 Z! |! ]! a0 k- G
. g8 ^7 }) h" H5 z% v! d1 d
% W A) Z# b3 s: l/ D. X$ ]脚本代码：
P X& v+ d8 v( B5 Z8 F% t
+ D* Q/ i L8 E# S* I" }! c- v) y" _- 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]);
- Z3 R2 X; o( |4 \* S
l Q; _- f1 W1 l" v
9 c0 e1 l/ ? z; Y( o" g* L' S+ R' V- e1 V6 J- n J9 p
U* Z: C2 c6 U5 V+ I复指数序列
0 H: N0 X2 @6 ]
& e+ R7 V, l+ F' P) e! q
0 o; \8 G0 m% l% G
" ?4 M. }8 f# ]! e5 H- clc
- clear
- close all
- n = 0:10;
- x = exp( (2+3j)*n );
- stem(n,x);
- xlabel('n');
6 t( y+ O4 m7 d$ x! s/ J
4 T0 E7 p; M) S* }& n
6 w& m; E K- U4 j0 D) l
+ C8 U2 X- F* Y8 l" c; }2 c5 t5 L% }! Q! u1 V5 [; r0 i' h! F
周期序列
2 S5 z) B) q0 k0 M0 q5 v% s. Z( O$ t/ S2 K o' ?1 C4 X5 S, D; r& r1 t
- 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');
3 j; E! I0 l1 D! A6 u
4 C) u/ l% f3 C% _% e, a6 t8 a
P7 M" J8 L0 a- `$ H
% |5 n, [: u3 F/ l5 I$ C/ v; d& Y$ u, R) t6 S! Q# j# v) F
|
|
/1 
发表于 2020-1-2 10:49
收藏
淘帖
支持！
反对！