信号时域采样 谱分析(matlab).doc

基于matlab的时域信号采样及频谱分析 一：主要设计方法与步骤： 画出连续时间信号的时域波形及其幅频特性曲线，其中，幅度因子，衰减因子，模拟角频率； 对信号进行采样，得到采样序列，，其中，为采样间隔，通过改变采样频率可改变，画出采样频率分别为，，时的采样序列波形； 对不同采样频率下的采样序列进行频谱分析，绘制其幅频和相频曲线，对各频率下采样序列和的幅频曲线有无差别，如有差别说明原因； 设系统单位抽样响应为，求解当输入为时的系统响应，画出，，的时域波形及幅频特性曲线，并利用结果验证卷积定理的正确性（此内容将参数设置为，，，）； 用对信号，，进行频谱分析，观察与4中结果有无差别； 由采样序列恢复出连续时间信号，画出其时域波形，对比与原来的连续时间信号的时域波形，计算并记录两者最大误差。 二：详细程序及仿真波形分析 1.连续时间信号及其频率抽样信号函数 % 绘制信号x(n)的幅度谱和相位谱 clc clear all close all n=0:50 % 定义序列的长度是50 A=input(请入A的值A:) % 设置信号的有关参数 a=input(请入a的值a:) w0=input(请入w0的值w0:) 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) close all subplot(2,1,1) stem(n,x) % 绘制x(n)的图形 grid on title(离散时间信号) subplot(2,1,2) plot(n,x) grid on title(连续时间信号) figure(2) 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连续时间信号) k=-25:25 W=(pi/12.5)*k w=W/pi Y1=y1*exp(-j*pi/12.5).^(n*k) figure (3) 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/12.5)).^(n*k) figure (4) 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(500Hz理想采样信号序列的相位谱) Y3=y3*(exp(-j*pi/12.5)).^(n*k) figure (5) subplot(2,1,1) plot(w,abs(Y3)) grid xlabel(w) ylabel(幅度) title(1000Hz理想采样信号序列的幅度谱) axis([-2 2 0 1000]) subplot(2,1,2) plot(w,angle(Y3)) grid xlabel(w) ylabel(幅角) title(1000Hz理想采样信号序列的相位谱) 分析：采样频率为时没有失真，时有横线，产生失真，时横线加长，失真加大。说明采样频率越大，失真越小。 2.设系统单位抽样响应，求解当输入为时的系统响应，画出，，的时域波形及幅频特性曲线，并利用结果验证卷积定理的正确性（此内容将参数设置为，，，）。 clc clear all close all n=1:50 % 定义序列的长度是50 hb=zeros(1,50) % 注意：matlab中数组下标从1开始 hb(1)=1 hb(2)=1 hb(3)=1 hb(4)=1 hb(5)=1 close all subplot(3,1,1) stem(hb) title(系统hb[n]) m=1:50 % 设定序列和