例6.1.2利用MATLAB绘制上述三种信号的波形图 ■clear; ■ n0=0;nf=10;ns- Figure 1 ■ nl=n0:nf: File Edit View Insert Tools Desktop Window Help ■ x1=[zeros(1,ns-r 凸启日是哈®⑨9回早公,同 ■ subplot(2,2,1),st 单位脉冲序列 单位阶跃序列 %用逻辑式产生 ■ nl=n0:nf;x1=[( 0.5 0.5 ■n2=n0:nf, 10 5 10 ■ x2=[zeros(1,ns-r ■ %用逻辑式产生单 ■ nl=n0:nf;x1=[( subplot(2,2,2),st
例6.1.2 利用MATLAB绘制上述三种信号的波形图 n clear; n n0=0; nf=10; ns=3; n n1=n0:nf; n x1=[zeros(1,ns-n0),1,zeros(1,nf-ns)]; %单位脉冲序列的产生 n subplot(2,2,1),stem(n1,x1);title('单位脉冲序列'); n % 用逻辑式产生单位脉冲序列更为简洁 n % n1=n0:nf; x1=[(n1-ns)==0] n n2=n0:nf; n x2=[zeros(1,ns-n0),ones(1,nf-ns+1)]; %单位阶跃序列的产生 n % 用逻辑式产生单位阶跃序列语句 n % n1=n0:nf; x1=[(n1-ns)>=0] n subplot(2,2,2),stem(n2,x2);title('单位阶跃序列');
复指数序列: n3=n0:nf; ■ x3=exp(-0.2+0.5j)*n3·%复指数序列 ☑Figure1 ▣回 subplot(2,2,3),stem(n3 Eile Edit View Insert Tools Desktop Window Help 凸启日®同沙回果公,司口国口四 title(实部指数序列),】 单位脉冲序列 单位阶跃序列 0.5 0.5 subplot(2,2,4),stem(n3 0e-e 10 0 5 10 title(虚部指数序列), 实部指数序列 0.5 橡 0 -0.5 5 10
复指数序列: n n3=n0:nf; n x3=exp((-0.2+0.5j)*n3); %复指数序列 n subplot(2,2,3),stem(n3,real(x3));line([0,10],[0,0]); n title('实部指数序列'),ylabel('实部'); n subplot(2,2,4),stem(n3,imag(x3)); n title('虚部指数序列'), ylabel('虚部');
程序运行结果为: ☑Figure1 File Edit View Insert Iools Desktop Window Help 凸3日曼 ®⊙回¥ 里公同 0图■回 单位脉冲序列 单位阶跃序列 命命 0.5 0.5 0哈ee 5 10 实部指数序列 虚部指数序列 0.6 分 伞 P 0.4 0.5 p 希 希 0.2 0 0 平 -0. -0.2 0 10 0 5 10
程序运行结果为:
(4)正弦序列 定义:x(n)=sin(no+) ⊙。是正弦序列的数字域频率 中为初相位。与连续的正弦 信号不同,正弦序列的自变量必须为整数。可以证明, 2元 只有当 为有理数时,正弦序列具有周期性。 (⑤)单边指数序列 定义: x(n)=a”u(n)
(4)正弦序列 定义: 0 x(n n ) = + sin( ) w f 是正弦序列的数字域频率; 为初相位。与连续的正弦 信号不同,正弦序列的自变量n必须为整数。可以证明, 只有当 为有理数时,正弦序列具有周期性。 w0 f 0 2p w (5)单边指数序列 定义: ( ) ( ) n x n = aun ( ) 1.2 ( ) n x n = u n
例6.1.3 给出MATLAB:绘制正弦序列x(n)=sin( ■ 和单边指数序列x(n)=1.2”u(n)的程序。 ■n=0:39;a=1.2: x1=sin(pi/6*n) x2-a.n; subplot(2,1,1),stem(n,x1,'fill'); xlabel(n'),grid on; title(正弦序列),axis([0,40,-1.5,1.5]); ■ subplot(2,1,2),stem(n,x2,'fill'),grid on; xlabel('n'),title('x(n)=1.2n*u(n));
例6.1.3 n 给出MATLAB绘制正弦序列 n 和单边指数序列 的程序。 n n=0:39; a=1.2; n x1=sin(pi/6*n); n x2=a.^n; n subplot(2,1,1), stem(n,x1, 'fill'); n xlabel('n'), grid on; n title('正弦序列'),axis([0,40,-1.5,1.5]); n subplot(2,1,2),stem(n,x2,'fill'),grid on; n xlabel('n'),title('x(n)=1.2^n*u(n)'); ) 6 ( ) sin( n p x n = ( ) 1.2 ( ) n x n = u n