信号处理作业.docxVIP

信号处理作业 权双璐 3113001126 硕3003班 装备所动态室 问题描述：构造一个多频正弦信号组成的信号，并用FFT分解和小波分解研究不同频带频谱的特点。 构造信号：x(t)=sin(2π*40t)+ sin(2π*170t)+ sin(2π*280t)+ sin(2π*370t)。 抽样频率：fs≧2fmax，这里取800Hz； 抽样点数：n=800； 信号的傅里叶变换分解 用傅里叶变换对信号进行分解，将信号从时域中分解到频域中。其源程序如下： clear all clf %构造多频正弦信号 f1=40;f2=170;f3=280;f4=370;%数据个数采样间隔及格信号频率 fs=800;N=800;dt=1/fs;n=0:N-1;t=n*dt;f=n/(N*dt);%定义时间离散值 x1=sin(2*pi*f1*t);%第一个正弦信号 x2=sin(2*pi*f2*t);%第二个正弦信号 x3=sin(2*pi*f3*t);%第三个正弦信号 x4=sin(2*pi*f4*t);%第四个正弦信号 x=x1+x2+x3+x4;%以上四个正弦信号合成 %Fourier分析 xn=fft(x,N);%对原始信号进行傅里叶变换 magxn=abs(xn);phaxn=angle(xn);%求出Fourier变换的振幅和相位 subplot(2,1,1),plot(t,x);xlabel(时间/s)%绘出原始信号 xlabel(时间/s) ylabel(振幅) title(原始信号(N=800)); k=0:length(magxn)-1; subplot(2,1,2),plot(f(1:N/2),magxn(1:N/2)*2/N); %绘制出Fourier变换的振幅谱 xlabel(频率/hz) ylabel(振幅) title(傅里叶变换结果); 分解后输出结果如下图所示： 图1.多频正弦信号傅里叶变换结果 二、对信号进行二层小波变换，其实现程序为： clear all clf %构造多频正弦信号 f1=40;f2=170;f3=280;f4=370;%数据个数采样间隔及格信号频率 n=800;fs=800;dt=1/fs;N=0:n-1;t=N*dt;f=N/(n*dt);%定义时间离散值 x1=sin(2*pi*f1*t);%第一个正弦信号 x2=sin(2*pi*f2*t);%第二个正弦信号 x3=sin(2*pi*f3*t);%第三个正弦信号 x4=sin(2*pi*f4*t);%第四个正弦信号 x=x1+x2+x3+x4;%以上四个正弦信号合成 subplot(5,2,1),plot(t,x);xlabel(时间/s)%绘出原始信号 title(原始信号); T=wpdec(x,2,db5); %信号y使用db5做两层小波包分析的小波包树 %小波包树各节点处的小波包系数 plot(T); s1=wprcoef(T,[2,0]); s2=wprcoef(T,[2,1]); s3=wprcoef(T,[2,2]); s4=wprcoef(T,[2,3]); subplot(5,2,3),plot(x1); title(小波包[2,0]的重构系数); subplot(5,2,5),plot(x2); title(小波包[2,1]的重构系数); subplot(5,2,7),plot(x3); title(小波包[2,2]的重构系数); subplot(5,2,9),plot(x4); title(小波包[2,3]的重构系数); Y=2*abs(fft(x))/n; Y1=2*abs(fft(s1))/n; Y2=2*abs(fft(s2))/n; Y3=2*abs(fft(s3))/n; Y4=2*abs(fft(s4))/n; subplot(5,2,2),plot(f(1:n/2),Y(1:n/2)); title(原信号幅值谱);ylim([0,1]); xlabel(频率/hz) ylabel(振幅) subplot(5,2,4),plot(f(1:n/8),Y1(1:n/8)); title(小波包分解[2,0]幅值谱); xlabel(频率/hz) ylabel(振幅) subplot(5,2,6),plot(f(n/8:n/4),Y2(n/8:n/4)); title(小波包分解[2,1]幅值谱); xlabel(频率/hz) ylabel(振幅) subplot(5,2,8),plot(f(n/4:3*n/8),Y3(n/4:3*n/8)); title(小波包分解[2,2]幅值谱); xlabel(频率/hz) ylabel(振幅) subplot(5,2,10),plot(f(3*