- 1、本文档共23页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
数字信号处理实验报告-信号采集与重建
实验二 信号的采样与重建
一.实验目的
通过观察采样信号的混叠现象,进一步理解奈奎斯特采样频率的意义。
通过实验,了解数字信号采样转换过程中的频率特征。
对实际的 音频文件作内插和抽取操作,体会低通滤波器在内插和抽取中的作用。
二.实验内容
采样混叠,对一个模拟信号Va(t)进行等间采样,采样频率为200HZ,得到离散时间信号V(n).Va(t)由频率为30Hz,150Hz,170Hz,250Hz,330Hz的5个正弦信号的加权和构成。
Va(t)=6cos(60pi*t)+3sin(300pi*t)+2cos(340pi*t)+4cos(500pi*t)+10sin(660pi*t)观察采样后信号的混叠效应。
程序:clear,
close all,
t=0:0.1:20;
Ts=1/2;
n=0:Ts:20;
V=8*cos(0.3*pi*t)+5*cos(0.5*pi*t+0.6435)-10*sin(0.7*pi*t);
Vn=8*cos(0.3*pi*n)+5*cos(0.5*pi*n+0.6435)-10*sin(0.7*pi*n);
subplot(221)
plot(t,V),
grid on,
subplot(222)
stem(n,Vn,.),
grid on,
输入信号X(n)为归一化频率f1=0.043,f2=0.31的两个正弦信号相加而成,N=100,按因子M=2作抽取:(1)不适用低通滤波器;(2)使用低通滤波器。分别显示输入输出序列在时域和频域中的特性。
程序:clear;
N=100;
M=2;
f1=0.043;
f2=0.31;
n=0:N-1;
x=sin(2*pi*f1*n)+sin(2*pi*f2*n);
y1=x(1:2:100);
y2=decimate(x,M,fir);
figure(1);
stem(n,x(1:N));
title(input sequence);
xlabel(n);ylabel(fudu);
figure(2);
n=0:N/2-1;
stem(n,y1);
title(output sequence without LP);
xlabel(n);ylabel(fudu);
figure(3);
m=0:N/M-1;
stem(m,y2(1:N/M));
title(output sequence with LP);
xlabel(n);ylabel(fudu);
figure(4);
[h,w]=freqz(x);
plot(w(1:512),abs(h(1:512)));
title(frequency spectrum of the input sequence);
xlabel(w);ylabel(fudu);
figure(5);
[h,w]=freqz(y1);
plot(w(1:512),abs(h(1:512)));
title(frequency spectrum of the output sequence without LP);
xlabel(w);ylabel(fudu);
figure(6);
[h,w]=freqz(y2);
plot(w(1:512),abs(h(1:512)));
title(frequency spectrum of the output sequence without LP);
xlabel(w);ylabel(fudu);
输入信号X(n)为归一化频率f1=0.043,f2=0.31的两个正弦信号相加而成,长度N=50,内插因子为2.(1)不适用低通滤波器;(2)使用低通滤波器。分别显示输入输出序列在时域和频域中的特性。
程序:clear,
close all,
N=50;
L=2;
f1=0.043;
f2=0.31;
n=0:N-1;
x=sin(2*pi*f1*n)+sin(2*pi*f2*n);
figure(1);
stem(n,x(1:N));
title(input sequence);
xlabel(n);ylabel(fudu);
y1=zeros(1,N*2);
y1(1:2:N*2)=x;
figure(2);
m=0:N*L-1;
stem(m,y1(1:N*L));
title(output sequence );
xlabel(n);ylabel(fudu);
y2=interp(x,L);
figure(3);
m=0:N*L-1;
stem(m,y2(1:N*L));
title(output sequence);
xlabel(n);ylabel(fudu);
figure(4);
[h,w]=fr
您可能关注的文档
最近下载
- 最新2023版知识产权贯标GBT29490 18知识产权合规管理控制程序(含表单)[知识产权合规管理体系文件].docx
- 新源(财)Y026-2022国网新源公司工程财务管理办法.pdf VIP
- 二年级上册阅读计划(教材“快乐读书吧”) .pdf VIP
- 《四川省建筑工程施工扬尘防治标准》DBJ51T 231-2023.docx
- 爱心公益基金会章程(精选3篇).docx VIP
- 学校餐厨垃圾处理合同.docx VIP
- 二年级上册语文上册-abb、abcc、aabc、abab、aabb等词语 人教部编版.pdf
- 银行数据治理建设方案.pdf VIP
- 2022类风湿关节炎诊疗规范(全文).pdf
- 内部控制制度固定资产循环.docx VIP
文档评论(0)