资源描述:
《常见离散信号产生与实现+grid on》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库。
1、实验1常见离散信号产生和实现一、实验目的1、加深对常用离散信号的理解;2、熟悉使用MATLAB在时域中产生一些基本的离散时间信号。二、实验原理1、单位抽样序列在MATLAB中可以利用函数实现。k1=input('k1='),k2=input('k2='),k0=input('k0=')k=k1:k2;n=length(k);x=zeros(1,n);x(1,k0-k1+1)=1;stem(k,x,'filled')axis([k1,k2,0,1.5])title('单位抽样序列')gridon2、单位阶跃序列在MATLAB中可
2、以利用函数实现:k1=input('k1='),k2=input('k2='),k0=input('k0=')k=k1:k0-1;kk=k0:k2;n=length(k);nn=length(kk);u=zeros(1,n);uu=ones(1,nn);stem(kk,uu,'filled')holdonstem(k,u,'filled')holdofftitle('单位阶跃序列')axis([k1,k2,0,1.5])gridon3、正弦序列在MATLAB中实现过程如下:k1=input('k1='),k2=input('k
3、2='),A=input('A='),f=input('f='),fai=input('fai='),Fs=input('Fs=')k=k1:k2;y=A*sin(2*pi*f*k/Fs+fai)stem(k,y,'filled')title('离散时间正弦序列f(n)=Asin(2*pi*f*n/Fs+fai)')xlabel('n')ylabel('f(n)')gridon参考初值:k1=-30;k2=30;A=3;f=1/(8*pi);fai=0;Fs=44、复指数序列在MATLAB中实现过程如下:k1=input('k
4、1='),k2=input('k2='),c=input('c='),w=input('w=')k=k1:k2;x=c*exp(j*w*k);subplot(2,2,1);stem(k,real(x),'filled')title('指数序列实部')gridonxlabel('n');ylabel('f(n)')subplot(2,2,2);stem(k,imag(x),'filled')title('指数序列虚部')gridonsubplot(2,2,3);stem(k,angle(x),'filled')title('指数
5、序列辐角')gridonxlabel('n');ylabel('f(n)')subplot(2,2,4);stem(k,abs(x),'filled')title('指数序列幅值')gridonstem(k,x,'filled')holdonplot([k1,k2],[0,0])holdofftitle('复指数序列')xlabel('n')ylabel('f(n)')gridon5、指数序列在MATLAB中实现过程如下:k1=input('k1='),k2=input('k2='),a=input('a=')k=k1:k2;
6、x=(a.^k);stem(k,x,'filled')holdonplot([k1,k2],[0,0])holdofftitle('指数序列')gridonxlabel('n')ylabel('f(n)')三、预习要求1、预先阅读实验讲义(MATLAB基础介绍);2、讨论正弦序列、复指数序列的性质。A.绘出信号,当、时、、时的信号实部和虚部图;当时呢?此时信号周期为多少?当时,信号实部和虚部图如下z1=-1/12+j*pi/6;k=0:20;x=exp(z1*k);subplot(1,2,1);stem(k,real(x),'
7、filled')title('指数序列实部')xlabel('n');ylabel('f(n)')gridonsubplot(1,2,2);stem(k,real(x),'filled')title('指数序列虚部')gridonxlabel('n');ylabel('f(n)')当时,信号实部和虚部图如下z2=1/12+j*pi/6;k=0:20;x=exp(z2*k);subplot(1,2,1);stem(k,real(x),'filled')title('指数序列实部')gridonxlabel('n');ylabel
8、('f(n)')subplot(1,2,2);stem(k,imag(x),'filled')title('指数序列虚部')gridonxlabel('n');ylabel('f(n)')当时,信号实部和虚部图如下z3=1/12;k=0:20;x=exp(z3*k);subp