根据采样序列重建原始模拟信号的仿真Word文件下载.doc
- 文档编号:14181949
- 上传时间:2022-10-19
- 格式:DOC
- 页数:8
- 大小:105.50KB
根据采样序列重建原始模拟信号的仿真Word文件下载.doc
《根据采样序列重建原始模拟信号的仿真Word文件下载.doc》由会员分享,可在线阅读,更多相关《根据采样序列重建原始模拟信号的仿真Word文件下载.doc(8页珍藏版)》请在冰豆网上搜索。
2.设计一个10倍内插算法,将信号x(n)采样率由100HZ变为10000HZ,画出内插后信号的时域波形。
分步实施:
1)熟悉多速率采样;
熟悉内插算法;
熟悉滤波器的设计。
2)编写主要程序;
3)画出画出采样后信号的时域波形x(n),内插后信号的时域波形。
二、实验原理
整数I倍内插是在已知的相邻两个原采样点之间等间隔的插入I-1个新的采样值。
由于这I-1个采样值并非已知值,所以关键问题是如何求出这I-1个采样值。
整数内插是先在已知采样序列x(n)的相邻两个样点之间等间隔插入I-1个0值点,然后进行低通滤波,即可求出I倍内插的结果。
整数倍内插后,相当于再原函数的频谱之间插入了I-1个波形相同的频谱。
只需通过低通滤波器后即可将低频分量去除还原出原波形。
三、设计流程
对20hz信号进行100hz采样
进行10倍插零
低通滤波器设计
还原信号
将信号与理论结果进行对比
结束
开始
四、程序运行结果
图1中最后一幅图是调用系统内插函数后理论波形。
图2中最后一幅图是经过低通滤波器后还原的波形。
将两幅图进行比较
图中.为理论结果,*为滤波后还原波形,可见两波形基本一致。
五、实验源程序
见附录
六、总结与体会
通过这次通信原理课程设计实践,我复习了MATLAB编程语言的基本概念、语法、语义和数据类型的使用特点,加深了对课堂所学理论知识的理解,掌握了运用结构化程序设计的基本思想和方法,更重要的是培养了自己的自学能力。
因为这是我们第二次接触MATLAB编程语言,在编写程序以及调试的过程中遇到了很多困难,但是我通过去图书馆查找资料,请教同学老师,再自己一点点改善程序,最终编写出一个比较完善的程序,实现了所有要求功能,这是最值得我欣慰的一点。
七、参考文献
【1】高西全.丁玉美.数字信号处理(第三版)西安:
西安电子科技大学出版社2008.8
【2】樊昌信.曹丽娜.通信原理(第六版)北京:
国防工业出版社2010.6
【3】冯象初.甘小冰.数值泛函与小波理论西安:
西安电子科技大学出版社2003.5
附录(源程序)
clc
T=2;
f0=20;
f1=100;
d=10000;
M=10;
t=0:
1/d:
T/f0;
t1=0:
1/f1:
t2=0:
1/f1/M:
xt=sin(2*pi*f0*t);
%20hz信号
xn=sin(2*pi*f0*t1);
%100hz采样后信号
figure
(1);
subplot(331)
plot(t,xt);
title('
20Hz信号'
);
subplot(332)
stem(t1,xn);
100Hz采样'
%-------插零----------------%
L=length(xn);
xn2=zeros(1,(L-1)*M+1);
forj=1:
L
xn2(M*(j-1)+1)=xn(j);
end
%----------------------------%
subplot(333)
stem(t2,xn2,'
.'
10倍内插'
fr=fft(xt,1024);
freq=(0:
length(fr)-1)*d/length(fr);
subplot(334)
plot(freq,abs(fr));
%xt频谱
xlim([0200])
frn=fft(xn,1024);
freqn=(0:
length(frn)-1)*f1/length(frn);
subplot(335)
plot(freqn,abs(frn));
%xn频谱
frn2=fft(xn2,1024);
freqn2=(0:
length(frn2)-1)*f1*10/length(frn2);
subplot(336)
plot(freqn2,abs(frn2));
%xn2频谱
%%%%%%%%%%10倍内插%%%%%%%%%%%%%%%%%%%%%%%%
xn3=interp(xn,M);
tt=0:
T/f0+1/f1-1/f1/M;
subplot(313)
stem(tt,xn3,'
)
%%%%%滤波器%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
figure
(2);
wd=50;
N=100;
M=(N-1)/2;
nn=-M:
M;
n=nn+eps;
hd=sin(wd*n)./(pi*n);
w1=hanning(N)'
;
h1=hd.*w1;
H11=(abs(fft(h1,1024)));
HH11=[H11(1:
1024)];
w=(0:
1023);
subplot(311)
plot(w,HH11)
yy=H11.*frn2;
%xn2频谱与滤波器频谱相乘
xx=0:
1023;
subplot(312)
plot(xx,abs(yy))
hy=ifft(yy);
%对滤波后频谱进行傅立叶反变换
rr=0:
rrr=rr/1000;
subplot(313)
stem(rrr,10*hy,'
)
xlim([00.05*T+0.05])
figure(3);
%两图进行对比
hold
k.'
*'
-7-
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 根据 采样 序列 重建 原始 模拟 信号 仿真