- 1、本文档共11页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
数字信号处理实验实验一
数字信号处理实验报告
实验名称:离散时间系统的时域特性分析
学生姓名: z
学生学号:
学生班级:
上课时间: 周二上午
指导老师:
实验目的
线性时不变离散时间系统在时域中可以通过常系数线性差分方程来描述,冲激响应序列可以刻画其时域特性。本实验通过使用MATLAB函数研究离散时间系统的时域特性,以加深对离散时间系统的差分方程、冲激响应和系统的线性和时不变特性的理解。
实验原理
1.线性系统?满足叠加原理的系统称为线性系统,即若某一输入是由N个信号的加权和组成的,则输出就是系统对这几个信号中每一个输入的响应的加权和。?即:如果系统在x1(n)和x2(n)输入时对应的输出分别为y1(n)和y2(n),当对任意常数a1和a2,式T[a1x1(n)+a2x2(n)]=a1T[x1(n)]+a2[x2(n)]=a1y1(n)+a2y2(n)成立,则该系统是线性系统。?2.时不变系统?若输入x(n)的输出为y(n),则将输入序列移动任意位后,其输出序列除了跟着位移外,数值应保持不变。即:当T[x(n)]=y(n),满足T[x(n-m)]=y(n-m)?(m为任意整数)时,则该系统就称为时不变系统。?3.常系数线性差分方程?线性时不变离散系统的输入、输出关系可以用以下常系数线性差分描述:??????y(n)=-?∑aky(n-k)+?∑brx(n-r)?当输入x(n)为单位冲激序列时,输出y(n)即为系统的单位冲击响应h(n)。n=0:299;
clf;
n=0:299;
x1=cos(20/256*pi*n);
x2=cos(200/256*pi*n);
x=x1+x2;
num1=[0.5 0.27 0.77];
num2=[0.45 0.5 0.45];
den1=[1 0 0];
den2=[1 -0.53 0.46];
y1=filter(num1,den1,x1);
y2=filter(num1,den1,x2);
y=filter(num1,den1,x);
yt=y1+y2;
y3=filter(num2,den2,x1);
y4=filter(num2,den2,x2);
yn=filter(num2,den2,x);
ytt=y3+y4;
subplot(4,1,1);
stem(n,x);
ylabel(信号幅度);
title(输入信号);
subplot(4,1,2);
stem(n,yt);
ylabel(信号幅度);
title(输出信号);
subplot(4,1,3);
stem(n,ytt);
ylabel(信号幅度);
title(输出信号);
系统的输入与输出波形如下图形所示:
(2)解:MATLAB程序如下:
N=300;
num1=[0.5 0.27 0.77];
num2=[0.45 0.5 0.45];
den1=[1 0 0];
den2=[1 -0.53 0.46];
y1=impz(num1,den1,N);
subplot(2,1,1);
stem(y1);
xlabel(时间序列);
ylabel(信号幅度);
title(输入信号);y1=impz(num1,den1,N);
y2=impz(num2,den2,N);
subplot(2,1,2);
stem(y2);
xlabel(时间序列);
ylabel(信号幅度);
title(输入信号);
系统1与系统2的冲激响应序列的波形如下图形所示:
(3)ⅰ 解:MATLAB程序如下:
n=0:40;D=5;x=cos(20*pi*n/256)+cos(200*pi*n/256);xd=[zeros(1,D)?x];?%生成新序列xd=x(n-D),延时D个单位num=[0.45?0.5?0.45];den=[1?-0.53?0.46];ic=[0?0];y=filter(num,den,x,ic);yd=filter(num,den,xd,ic);N=length(y);d=y-yd(1+D:N+D);%求y与yd的差值subplot(3,1,1)?;stem(n,y)?;ylabel(信号幅度)?;title(输出y[n])?;grid?;subplot(3,1,2)?;stem(yd(1:length(yd)))?;ylabel(信号幅度)?;subplot(3,1,3)?;stem(n,d)?;xlabel(
您可能关注的文档
最近下载
- 教育、科技、人才一体化与新质生产力发展.pptx VIP
- 2024年湖南水利水电职业技术学院单招职业技能测试题库含答案(研优卷).docx VIP
- 教育、科技、人才一体化视域下发展新质生产力的四重逻辑和实践路径.docx VIP
- 教育、科技、人才一体化视域下发展新质生产力的四重逻辑和实践路径.docx VIP
- 《大棚番茄基质槽式栽培技术规程》.docx
- 立体车库 毕业论文外文翻译.docx VIP
- 教育、科技、人才一体推进新质生产力发展.docx VIP
- 2024年湖南水利水电职业技术学院单招职业技能测试题库含答案(考试直接用).docx VIP
- 社交礼仪教程(第7版 微课版)课件 第五章 餐饮礼仪.pptx
- 2025清华:DeepSeek从入门到精通.pdf
文档评论(0)