数字信号实验一

1， 冲击信号 （1）

程序： >> L=20; >> nn=1:20;

>> imp=zeros(20,1); >> imp(5)=1;

>> stem(nn,0.9*imp,'r:') >>

（2） 程序： >> L=31; >> nn=-15:15;

>> imp=zeros(31,1); >> imp(16)=1; >>

stem(nn,0.8*imp,'r:')

（3）

程序： >> L=51; >> nn=300:350; >> imp=zeros(L,1);

>> imp(34)=1;

>> stem(nn,1.5*imp,'r:')

（4）

程序： >> L=11; >> nn=-10:0;

>> imp=zeros(L,1); >> imp(4)=1;

>> stem(nn,4.5*imp,'r:')

2、正弦信号

M文件：

Function x=FC(A,w,q,N1,N2) n=N1:N2;

x=A*sin(w*n+q); stem(n,x)

（1） 程序： >> n=0:25;

>> x=sin(n*pi/17); >> stem(n,x)

（2）

程序：

>> n=-15:15;

>> x=sin(n*pi/17); >> stem(n,x) （3）

程序：

>> n=-10:10;

>> x=sin(n*pi*3+pi/2); >> stem(n,x)

（4）

程序： >> n=0:50;

>> x=cos(n*pi/(23)^(1/2));

>> stem(n,x)

3、指数信号 （a）

M文件：

function y=genexp(b,n0,L) if(L<=0)

error('GENEXP:length not positive') end

nn=n0+[1:L]'-1; y=b.^nn; end

程序：

>> x=genexp(0.9,1,20); >> stem(x)

（4）复值信号 程序： >> nn=0:25;

>> xx=exp(sqrt(-1)*nn/3); >> subplot(211) >> stem(nn,real(xx))

>> title('REAL PART'),xlabel('INDEX(n)') >> subplot(212)

>> stem(nn,imag(xx))

>> title('IMAG PART'),xlabel('INDEX(n)') >>