资源描述：
《信号时域采样 频谱分析(matlab).doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、基于matlab的时域信号采样及频谱分析一：主要设计方法与步骤：1.画出连续时间信号的时域波形及其幅频特性曲线，其中，幅度因子，衰减因子，模拟角频率；2.对信号进行采样，得到采样序列，，其中，为采样间隔，通过改变采样频率可改变，画出采样频率分别为，，时的采样序列波形；3.对不同采样频率下的采样序列进行频谱分析，绘制其幅频和相频曲线，对各频率下采样序列和的幅频曲线有无差别，如有差别说明原因；4.设系统单位抽样响应为，求解当输入为时的系统响应，画出，，的时域波形及幅频特性曲线，并利用结果验证卷积定理的正确性（此内容将参数设置为，，，）；5.用对信号，，进行频谱分析，观察与4中结果有无
2、差别；6.由采样序列恢复出连续时间信号，画出其时域波形，对比与原来的连续时间信号的时域波形，计算并记录两者最大误差。二：详细程序及仿真波形分析1.连续时间信号及其频率抽样信号函数%绘制信号x(n)的幅度谱和相位谱clcclearallclosealln=0:50%定义序列的长度是50A=input('请入A的值A:')%设置信号的有关参数a=input('请入a的值a:')w0=input('请入w0的值w0:')T1=0.005T2=0.002T3=0.001T0=0.001x=A*exp(-a*n*T0).*sin(w0*n*T0)y1=A*exp(-a*n*T1).*sin
3、(w0*n*T1)y2=A*exp(-a*n*T2).*sin(w0*n*T2)y3=A*exp(-a*n*T3).*sin(w0*n*T3)closeallsubplot(2,1,1)stem(n,x)%绘制x(n)的图形gridontitle('离散时间信号')subplot(2,1,2)plot(n,x)gridontitle('连续时间信号')figure(2)subplot(3,1,1)stem(n,y1)gridontitle('200Hz理想采样信号序列')subplot(3,1,2)stem(n,y2)gridontitle('500Hz连续时间信号')subpl
4、ot(3,1,3)stem(n,y3)gridontitle('1000Hz连续时间信号')k=-25:25W=(pi/12.5)*kw=W/piY1=y1*exp(-j*pi/12.5).^(n'*k)figure(3)subplot(2,1,1)plot(w,abs(Y1))gridxlabel('w')ylabel('幅度')title('200Hz理想采样信号序列的幅度谱')axis([-2201000])subplot(2,1,2)plot(w,angle(Y1))gridxlabel('w')ylabel('幅角')title('200Hz理想采样信号序列的相位谱')
5、Y2=y2*(exp(-j*pi/12.5)).^(n'*k)figure(4)subplot(2,1,1)plot(w,abs(Y2))gridxlabel('w')ylabel('幅度')title('500Hz理想采样信号序列的幅度谱')axis([-2201000])subplot(2,1,2)plot(w,angle(Y2))gridxlabel('w')ylabel('幅角')title('500Hz理想采样信号序列的相位谱')Y3=y3*(exp(-j*pi/12.5)).^(n'*k)figure(5)subplot(2,1,1)plot(w,abs(Y3))gr
6、idxlabel('w')ylabel('幅度')title('1000Hz理想采样信号序列的幅度谱')axis([-2201000])subplot(2,1,2)plot(w,angle(Y3))gridxlabel('w')ylabel('幅角')title('1000Hz理想采样信号序列的相位谱')分析：采样频率为时没有失真，时有横线，产生失真，时横线加长，失真加大。说明采样频率越大，失真越小。2.设系统单位抽样响应，求解当输入为时的系统响应，画出，，的时域波形及幅频特性曲线，并利用结果验证卷积定理的正确性（此内容将参数设置为，，，）。clcclearallcloseall
7、n=1:50%定义序列的长度是50hb=zeros(1,50)%注意：matlab中数组下标从1开始hb(1)=1hb(2)=1hb(3)=1hb(4)=1hb(5)=1closeallsubplot(3,1,1)stem(hb)title('系统hb[n]')m=1:50%设定序列和长度值T=1%设定序列的采样率A=1a=0.4T=1w0=2.0734x=A*exp(-a*m*T).*sin(w0*m*T)subplot(3,1,2)stem(x)title('输入序列x[n]
