数字信号matlab实验报告,matlab数字信号实验报告

matlab数字信号实验报告

数字信号处理实验报告班 级： 1105074143 姓 名： 路晓冬 学 号： 1105074143 实验一 频谱分析与采样定理一、实验目的1、观察模拟信号经理想采样后的频谱变化关系。2、验证采样定理，观察欠采样时产生的频谱混叠现象。3、加深对 DFT 算法原理和基本性质的理解。4、熟悉 FFT 算法原理和 FFT 的应用。二、实验原理根据采样定理，对给定信号确定采样频率，观察信号的频谱。三、实验内容和步骤1)实验内容在给定信号为：1．x(t)=cos(100*π*at)2．x(t)=exp(-at)3．x(t)=exp(-at)cos(100*π*at)其中 a 为实验者的学号，记录上述各信号的频谱，表明采样条件，分析比较上述信号频谱的区别。2)实验步骤1．根据采样理论、DFT 的定义、性质和用 DFT 作谱分析的有关内容。2．根据 FFT 算法原理和基本思想。3．确定实验给定信号的采样频率，编制对采样后信号进行频谱分析的程序四、实验过程%实验一:频谱分析与采样定理T=0.0001; %采样间隔T=0.0001F=1/T; %采样频率为F=1/TL=0.02; %记录长度L=0.02 N=L/T; t=0:T:L; a=36;f1=0:F/N:F;f2=-F/2:F/N:F/2;%%%%%%%%%%%%%%%%%%%%%%%%%x1=cos(100*pi*a*t);y1=T*abs(fft(x1)); % 求复数实部与虚部的平方和的算术平方根 y11=fftshift(y1);figure(1),subplot(3,1,1),plot(t,x1);title( 正弦信号 );subplot(3,1,2),stem(y1);title( 正弦信号频谱 );subplot(3,1,3),plot(f2,y11);title( 正弦信号频谱 );%%%%%%%%%%%%%%%%%%%%%x2=exp(-a*t);y2=T*abs(fft(x2));y21=fftshift(y2);figure(2),subplot(3,1,1),stem(t,x2);title( 指数信号 );subplot(3,1,2),stem(f1,y2);title( 指数信号频谱 );subplot(3,1,3),plot(f2,y21);title( 指数信号频谱 );%%%%%%%%%%%%%%%%%%%%%x3=x1.*x2;y3=T*abs(fft(x3));y31=fftshift(y3);figure(3),subplot(3,1,1),stem(t,x3);title( 两信号相乘 );subplot(3,1,2),stem(f1,y3);title( 两信号相乘频谱 );subplot(3,1,3),plot(f2,y31);title( 两信号相乘频谱 );正弦信号频谱：指数信号频谱：两信号相乘频谱：5、实验结果及分析奈奎斯特抽样定理为抽样频率必须大于或等于信号频谱最高频率的 2 倍，即 。取采样间隔 T=0.0001s，满足奈奎斯特抽样定理。实验中正弦信hs2号的频率分量为 ，指数信号的频率分量为 0，两信号相Hz140280乘后频率分量为 1400Hz，所以三幅频谱图中相应出现 1400Hz、0Hz、1400Hz 的频率分量，与实际计算结果一致。频谱图中显示区域的最高频率即为在该采样频率下的最高可分析频率，等于信号采样频率的一半，即 。Hz5021.0六、实验体会通过实验一我熟悉了 matlab 的使用，加深了对奈奎斯特定理的理解。实验二 卷积定理一、实验目的通过本实验，验证卷积定理，掌握利用 DFT 和 FFT 计算线性卷积的方法。二、 实验原理时域圆周卷积在频域上相当于两序列 DFT 的相乘，因而可以采用 FFT 的算法来计算圆周卷积，当满足 时，线性卷积等于圆周卷积，因此12LN可利用 FFT 计算线性卷积。三、实验内容和步骤1、给定离散信号 和 ，用图解法求出两者的线性卷积和圆周卷积；()xnh2、编写程序计算线性卷积和圆周卷积；3、比较不同列长时的圆周卷积与线性卷积的结果，分析原因。四、实验过程%实验二:卷积定理x=[3 0 2 1 3]; %原始序列y=[3 0 2 1 3];%直接计算线性卷积z=conv(x,y);figure(1),subplot(311),stem(x);axis([1 9 0 4]); subplot(312),stem(y);axis([1 9 0 4]);subplot(313),stem(z);axis([1 9 0 30]);%利用 FFT 计算N=10;%N=10 时x1=[x zeros(1,N-length(x))];y1=[y zeros(1,N-length(y))];X1=fft(x1);Y1=fft(y1);Z1=X1.*Y1;z1=ifft(Z1);figure(2),subplot(321),stem(x1);subplot(322),stem(real(X1));subplot(323),stem(y1);subplot(324),stem(real(X1));subplot(325),stem(z1);subplot(326),stem(real(Z1));N=5;%N=5 时x2=[x zeros(1,N-length(x))];y2=[y zeros(1,N-length(y))];X2=fft(x2);Y2=fft(y2);Z2=X2.*Y2;z2=ifft(Z2);figure(3),subplot(321),stem(x2);subplot(322),stem(real(X2));subplot(323),stem(y2);subplot(324),stem(real(X2));subplot(325),stem(z2);subplot(326),stem(real(Z2));直接计算线性卷积：利用圆周卷积计算，取 N=10：利用圆周卷积计算，取 N=5：