工程信号处理作业.doc

1. 对于信号x(t)=Asin(2πft)+0.1*u(t)，其中A=1，f=50 Hz，u(t)是一个均值为零的噪声信号（该信号可以由randn 函数产生）。以两种不同的采样频率(一个不满足采样定理，一个满足采样定理)分别对其进行采样得到两组离散信号，信号采样总时间自定。利用fft 命令对其进行离散傅立叶变换，作出频谱图形，体会频率混淆现象。 Matlab程序如下： T1=0.005; T2=0.01; fs1=1/T1; fs2=1/T2; t1=0:T1:1; t2=0:T2:1; A=1;f=50;u1=randn(size(t1));u2=randn(size(t2)); x1=A*sin(2*pi*f*t1)+0.1*u1; x2=A*sin(2*pi*f*t2)+0.1*u2; X1=fft(x1,256); X2=fft(x2,128); y1=abs(X1); y2=abs(X2); f1=(0:255)*fs1/256; f2=(0:127)*fs2/128; plot(f1(1:128),y1(1:128)*2/256); xlabel(频率/Hz); ylabel(幅值); title(满足采样定理频谱图形); figure; plot(f2(1:64),y2(1:64)*2/128); xlabel(频率/Hz); ylabel(幅值); title(不满足采样定理频谱图形); 运行结果如下图所示： 2. 对于信号x(t)=A1sin(2πf1t)+ A2sin(2πf2t)，其中A1=2，f1=20 Hz；A2=3，f2=22 Hz。以一定的采样频率采集不同时间长度的两组离散信号，分别对其进行离散傅叶变换，作出频谱图形，体会数据长度对频率分辨率的影响。 程序如下： A1=2,f1=20; A2=3,f2=22; T=0.01; fs=1/T; t1=0:T:0.2; t2=0:T:1; x1=A1*sin(2*pi*f1*t1)+A2*sin(2*pi*f2*t1); x2=A1*sin(2*pi*f1*t2)+A2*sin(2*pi*f2*t2); X1=fft(x1,32); X2=fft(x2,128); y1=abs(X1); y2=abs(X2); f1=(0:31)*fs/32; f2=(0:127)*fs/128; plot(f1(1:16),y1(1:16)*2/32); xlabel(频率/Hz); ylabel(幅值); title(长度过短对频率分辨率的影响); figure; plot(f2(1:64),y2(1:64)*2/128); xlabel(频率/Hz); ylabel(幅值); title(长度正常对频率分辨率的影响); 3. 对于信号x(t)= sin(2πft)，其中f=20 Hz。分别以矩形窗和hann 窗对其进行截取一段时间内的信号（采样频率自定），得到两组离散信号，分别对其进行离散傅立叶变换，作出频谱图形，体会窗函数对能量泄露的影响。 程序如下： f=20; N=512; fs=50; n=0:N-1; t=n/fs; x=sin(2*pi*f*t); w_han=(hanning(N)); w_rect=(rectwin(N)); x1=x.*w_han; x2=x.*w_rect; y1=fft(x1,512); z1=abs(y1); y2=fft(x2,512); z2=abs(y2); y3=fft(x,512); z3=abs(y3); f1=(0:511)*fs/512; subplot(3,1,1) plot(f1(1:256),z1(1:256)*2/512); xlabel(频率/Hz); ylabel(幅值); title(hanning窗对能量泄露的影响); subplot(3,1,2) plot(f1(1:256),z2(1:256)*2/512); xlabel(频率/Hz); ylabel(幅值); title(矩形窗对能量泄露的影响); subplot(3,1,3) plot(f1(1:256),z3(1:256)*2/512); xlabel(频率/Hz); ylabel(幅值); title(未加窗); 运行结果如图所示： 4. 对于线性扫频信号x(t)=exp(jαt)，α 是一个常数，可以通过chirp.m 文件产生该信号，如t=0:0.001:1;x=chirp(t,0,1,200); % 产生一个长度为1 s 的频率为0 Hz—200 Hz 扫频信号对该信号进行离散傅立叶变换，作出频谱图形，并验证Parseval 等式： 程序如下： t=0:0.001:1; fs=1000; x=chirp(t,0,1,20