- 1、本文档共29页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
现代信号处理课程设计报告报告帮助课程设计
现代信号处理课程设计报告
目录
课程设计目的…………………………………2
课程设计内容…………………………………2
课程设计流程及分析…………………………4
3.1熟悉MATLAB基本操作…………………………. 4
3.2 MATLAB数字信号处理图形界面…………………. 5
3.3用MATLAB绘出信号波形………………………. .5
3.4对信号频谱分析………………………. . . . . . . . .7
3.5设计滤波器并对信号滤波……………………….11
设计经验总结…………………………………27
参考资料………………………………………27
一、课程设计目的
1.所学知识
2.提高综合运用所学知识独立分析和解决问题的能力;
(一).信号分析
1、单位抽样序列 2)矩形序列 ???? 三角波序列
反三角波序列
5)Gaussian序列
???
6) 正弦序列 ?? ? 进行采样,可得到测试序列
。令A=50,采样周期T=1ms,即fs=1000Hz,f0=62.5,a=100。
2. 对上述信号完成下列信号分析
1)对三角波序列和反三角波序列,作N=8点的FFT,观察比较它们的幅频特性,说明它们有什么异同?绘出两序列及其它们的幅频特性曲线。
在和的尾部补零,作N=16点的FFT,观察它们的幅频特性发生了什么变化?分析说明原因。
? ?2)、观察高斯序列,固定信号中的参数p=8,令q分别等于2,4,8,观察它们的时域和幅频特性,了解当q取不同值时,对信号序列的时域幅频特性的影响;固定q=8,令p分别等于8,13,14,观察参数p变化对信号序列的时域及幅频特性的影响,观察p等于多少时,会发生明显的泄漏现象,混叠是否也随之出现?记录实验中观察到的现象,绘出相应的时域序列和幅频特性曲线。
3)对于正弦序列
4)、观察衰减正弦序列的时域和幅频特性,绘出幅频特性曲线,改变采样频率fs,使fs=300Hz, 观察此时的频谱的形状和谱峰出现位置?说明产生现象的原因。
3.设有一连续时间信号s(t),其由20Hz、220Hz和750Hz的正弦信号叠加而成,分析确定采样频率及数据分析长度,计算并绘出信号的频谱,指出各个频率份量。
(二)、数字滤波器设计
3.2 MATLAB数字信号处理图形界面
3.3用MATLAB绘出信号波形
global xy
popup_sel_index = get(handles.popupmenu2, Value);
switch popup_sel_index
case 1
n0=0;n1=-10;
n2=10;n=[n1:n2];
x=[(n-n0)==0];
stem(n,x,.);
xlabel(n);ylabel(x(n));
title(单位抽样序列);
case 2
n=-1:10;
n0=2;width=5;
x=rectpuls(n-n0,width);
stem(n,x,.);
xlabel(n);ylabel(x(n));
axis([-1 10 0 2]);
title(矩形序列);
case 3
n1=0:3;
n2=4:7;
x1=n1+1;
x2=8-n2;
n=[(-10:0) n1 n2 (7:17)];
x=[zeros(1,11) x1 x2 zeros(1,11)];
stem(n,x,.);
xlabel(n);ylabel(x(n));
title(三角波序列);
case 4
n1=0:3;
n2=4:7;
x1=4-n1;
x2=n2-3;
n=[(-10:0) n1 n2 (7:17)];
x=[zeros(1,11) x1 x2 zeros(1,11)];
stem(n,x,.);
xlabel(n);ylabel(x(n));
title(反三角波序列);
case 5
n=0:15;p=8;q=2;
x=exp(-(n-p).^2/q);
stem(n,x,.);
xlabel(n);ylabel(x(n));
title(Gaussian序列);
case 6
t=1/64;
n=0:t:16;
phase=0;A=1;
arg=2*pi*n-phase;
x=A*sin(arg);
stem(n,x,.);
xlabel(n);ylabel(x(n));
title(正弦序列);
case 7
t=0.001;
n=0:t:0.2;
x1=sin(
文档评论(0)