数字信号实验报告概要.doc

实验一：快速傅立叶变换的谱分析 、中的，进行信号的谱分析、题目： t1=[0:0.001:0.3]; t2=[0.301:0.001:0.6]; t3=[0.601:0.001:0.9]; t4=[0.901:0.001:1.199]; x1=sin(2*pi*100*t1); x2=sin(2*pi*50*t2); x3=sin(2*pi*25*t3); x4=sin(2*pi*10*t4); 信号s1为 t= 0:0.001:1.199;s1= sin(2*pi*100*t)+sin(2*pi*50*t)+sin(2*pi*25*t)+sin(2*pi*10*t); 信号s2为t=[t1,t2,t3,t4];s2= [x1,x2,x3,x4]; 信号s3为t=[t1,t2,t3,t4];s3= [x1,x4,x2,x3]; t1=[0:0.001:0.3]; t2=[0.301:0.001:0.6]; t3=[0.601:0.001:0.9]; t4=[0.901:0.001:1.199]; x1=sin(2*pi*100*t1); x2=sin(2*pi*50*t2); x3=sin(2*pi*25*t3); x4=sin(2*pi*10*t4); %信号s1的时域图和频域图 t5=0:0.001:1.199; x5=sin(2*pi*100*t5)+sin(2*pi*50*t5)+sin(2*pi*25*t5)+sin(2*pi*10*t5); y5=abs(fft(x5)); f=1000*(0:(1/1200):0.5); figure(1) subplot(2,1,1),plot(t5,x5);grid on;axis tight;%波形图 title(时域波形);xlabel(t(s)); subplot(2,1,2),plot(f,y5(1:601));grid on;axis tight;%幅频图 title(幅频图);xlabel(f(Hz)); %信号s2的时域图和频谱图 t6=[t1,t2,t3,t4]; x6=[x1,x2,x3,x4]; y6=abs(fft(x6)); figure(2) subplot(2,1,1),plot(t6,x6);grid on;axis tight; %波形图 title(时域波形);xlabel(t(s)); subplot(2,1,2),plot(f,y6(1:601));grid on;axis tight; %幅频图 title(幅频图);xlabel(f(Hz)); %信号s3的时域图和频域图 x7=[x1,x4,x2,x3]; y7=abs(fft(x7)); figure(3) subplot(2,1,1),plot(t6,x7);grid on;axis tight; %波形图 title(时域波形);xlabel(t(s)); subplot(2,1,2),plot(f,y7(1:601));grid on;axis tight; %幅频图 title(幅频图);xlabel(f(Hz)); a. 编写程序分别画出信号s1,s2,s3的时域波形和幅频图（参考下图）。 b. 观察信号s1,s2,s3的时域波形和频谱图，思考其幅频图的差别及其原因。 因为信号s1、s2、s3这三个信号中都有频率为hz，hz，hz，的频率分量，可以看出s1含有这四种分量，而s2s3的频谱图相同，除了这四种分量有其他分量因为s1这四种分量存在但s2s3中只是特定时间段出现。 横坐标两相邻点的频率间隔为=5/6hz,所以各频谱图的频率分辨率tp=1/NT=fs/N,因为N=1200,所以s1，s2，s3的抽样频率fs=NFHZ. d .若有信号s4的幅频图与s2的幅频图完全相同（如图2所示），问s4的时域波形和s2是相同的吗，为什么? 若有信号s4的幅频图与s2的幅频图完全相同，s4的时域波形和s2可能是不幅频特性只表明了幅度和频率的关系，即对频率分量的，而其相位特性，时域图幅度和相位图完全一致时才相同。 IR数字滤波器设计 一、 1、熟悉滤波器的基本概念； 2、了解滤波器的分类； 3、熟悉应用matlab设计各种滤波器的方法。 二、 1.数字滤波器分类 实现方法 有限长冲激响应----FIR滤波器 无限长冲激响应----IIR滤波器 功能 低通（LP） 高通（HP） 带通（BP） 带阻（BS） 2、数字滤波器设计步骤 一、给出所需要的滤波器的技术指标 二、设计一个H(Z)使其逼近所需要的技术指标