数字信号处理仿真设计.docVIP

数字信号处理仿真设计 课程设计题目： 基于 MATLAB 的数字信号处理基本概念、基本理论和基本方法学会用MATLAB对信号进行分析和处理的时域波形，其中幅度因子A＝10，衰 减因子a＝8，模拟角频率＝π； 2、对信号进行采样，得到采样序列，其中T＝为采样间隔，通过改变采样频率可改变T，画出采样频率分别为200Hz，500 Hz，1000 Hz时的采样序列波形； 3、分别用傅里叶变换和DFT对不同采样频率下的采样序列进行频谱分析，绘制其幅频和相频曲线，对比各频率下采样序列和的幅频曲线的差别。 4、考察序列：x (n ) = cos(0 .48πn ) + cos(0 .52πn ) 0 ≤ n ≤ 10时，用 DFT 估计 x (n )的频谱；将 x (n )补零加长到长度为 100点序列用 DFT估计 x (n )的频谱。要求画出相应波形。 0 ≤ n ≤ 100时，用 DFT 估计 x (n )的频谱，并画出波形 利用DFT计算线性卷积，并分析线性卷积等于循环卷积的条件。 四、详细程序及仿真波形： 1、连续时间信号x(t)及其200Hz/500Hz/1000Hz频率抽样信号函数x(n)的波形和频谱分析。 %1.绘制信号x(t)的时域波形 t=0:0.001:0.4; %设置连续信号的有关参数 A=input(请输入A的值 A:); a=input(请输入a的值 a:); w0=input(请输入w0的值 w0:); x=A*exp(-a*t).*sin(w0*t); %不是周期函数 figure(1) subplot(2,1,1);plot(t,x),grid on title(x(t)连续时间信号) xlabel(时间 t); ylabel(幅度 x(t)); subplot(2,1,2);plot(t,x),grid on axis([0 0.045 -2 160]); title(调整后的x(t)连续时间信号) xlabel(时间 t); ylabel(幅度 x(t)); %2.绘制各采样序列x(n)的时域波形 n=0:64; %定义序列的长度是64 T1=0.005; T2=0.002; T3=0.001; T0=0.001; x=A*exp(-a*n*T0).*sin(w0*n*T0); y1=A*exp(-a*n*T1).*sin(w0*n*T1); %是周期函数 y2=A*exp(-a*n*T2).*sin(w0*n*T2); y3=A*exp(-a*n*T3).*sin(w0*n*T3); figure(2) subplot(2,1,2);plot(n,x),grid on title(x(n)连续时间信号) subplot(2,1,1);stem(n,x),grid on %绘制序列x(n)的图形 title(x(n)离散时间信号) figure(3) subplot(3,1,1);stem(n,y1),grid on title(200Hz理想采样信号序列); subplot(3,1,2);stem(n,y2),grid on title(500Hz理想采样信号序列) subplot(3,1,3);stem(n,y3),grid on title(1000Hz理想采样信号序列) %3.绘制各采样序列x(n)的幅度谱和相位谱 k=-32:32; W=(pi/16)*k; %周期序列的傅里叶级数系数中旋转因子中的W w=W/pi; %坐标上的w，即此值是冲击函数刚好为1的点 Y1=y1*(exp(-j*pi/16)).^(n*k); %周期序列的傅里叶变换函数,n是n的转置 figure(4) subplot(2,1,1);plot(w,abs(Y1));grid,xlabel(w),ylabel(幅度); title(200Hz理想采样信号序列的幅度谱); axis([-2 2 0 1000]); subplot(2,1,2);plot(w,angle(Y1));grid,xlabel(w),ylabel(幅角); title (200Hz理想采样信号序列的相位谱) Y2=y2*(exp(-j*pi/16)).^(n*k); figure(5) subplot(2,1,1);plot(w,abs(Y2));grid,xlabel(w),ylabel(幅度); title(500Hz理想采样信号序列的幅度谱); axis([-2 2 0 1000]); subplot(2,1,2);plot(w,angle(Y2));grid,xlabel(w),ylabel(幅角); title (500