- 1、本文档共14页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
数字信号处理实验报告
班 级:
姓 名: 霍王迎
学 号:31
实验一 频谱分析与采样定理
一、实验目的
观察模拟信号经理想采样后的频谱变化关系。
验证采样定理,观察欠采样时产生的频谱混叠现象
加深对DFT算法原理和基本性质的理解
熟悉FFT算法原理和FFT的应用
二、实验原理
根据采样定理,对给定信号确定采样频率,观察信号的频谱
三、实验内容和步骤
实验内容
在给定信号为:
1.x(t)=cos(100*π*at)
2.x(t)=exp(-at)
3.x(t)=exp(-at)cos(100*π*at)
其中a为实验者的学号,记录上述各信号的频谱,表明采样条件,分析比较上述信号频谱的区别。
实验步骤
1.复习采样理论、DFT的定义、性质和用DFT作谱分析的有关内容。
2.复习FFT算法原理和基本思想。
3.确定实验给定信号的采样频率,编制对采样后信号进行频谱分析的程序
四、实验设备
计算机、Matlab软件
五、实验程序
T=0.00001;
F=1/T;
L=0.001;
N=L/T;
t=0:T:L;
a=25;
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(两信号相乘频谱);
七.实验结果
实验二 卷积定理
一、实验目的
通过本实验,验证卷积定理,掌握利用DFT和FFT计算线性卷积的方法。
二、 实验原理
时域圆周卷积在频域上相当于两序列DFT的相乘,因而可以采用FFT的算法来计算圆周卷积,当满足时,线性卷积等于圆周卷积,因此可利用FFT计算线性卷积。
三、实验内容和步骤
给定离散信号和,用图解法求出两者的线性卷积和圆周卷积;
编写程序计算线性卷积和圆周卷积;
比较不同列长时的圆周卷积与线性卷积的结果,分析原因。
四、实验设备
计算机、Matlab软件
五、实验程序
x=[1 1 0 1 3];
y=[3 0 0 1 3];
z=conv(x,y);
figure(1);
subplot(3,1,1);stem(x);
axis([1 9 0 4]);
subplot(3,1,2);stem(y);
axis([1 9 0 4]);
subplot(3,1,3);stem(z);
axis([1 9 0 30]);
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(3,2,1);stem(x1);
subplot(3,2,2);stem(real(X1));
subplot(3,2,3);stem(y1);
subplot(3,2,4);stem(real(Y1));
subplot(3,2,5);stem(z1);
subplot(3,2,6);stem(real(Z1));
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);
fig
文档评论(0)