- 1、本文档共14页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
数字信号处理实验二的()
实验二 离散时间系统的时域分析
实验室名称: 实验时间:
姓 名: 学号: 专业: 指导教师:
成绩
教师签名: 年 月 日 一、实验目的
1. 通过MATLAB仿真一些简单的离散时间系统,并研究它们的时域特性。
2. 通过MATLAB仿真滤波器,研究它们的性质。
3. 用MATLAB仿真研究离散时间系统的性质。
二、实验内容
Q2.1 M=2运行内容一的程序代码,生成输入x[n]=s1[n]+s2[n]的输出信号。
Q2.2若线性时不变系统由y[n]=0.5(x[n]+x[n-1])变成y[n]=0.5(x[n]-x[n-1])对输入x[n]=s1[n]+s2[n]的影响是什么?
Q2.4 修改程序P2.1,用一个长度为101、最低频率为0、最高频率为0.5的扫频正弦信号作为输入信号(见程序P1.7),计算其输出信号。
Q2.7 运行程序P2.3,对由加权输入得到的y[n]与在相同权系数下输出y1[n]和y2[n]相加得到的yt[n]进行比较,这两个序列是否相等?该系统是线性系统吗?
Q2.12 运行程序P2.4并比较输出序列y[n]和yd[n-10]。
Q2.19 运行程序P2.5,生成式(2.15)所给离散时间系统的冲击响应。
Q2.20 修改程序P2.5产生如下因果线性时不变系统的冲激响应的前45个样本:
y[n]+0.71y[n-1]-0.46y[n-2]-0.62y[n-3]=0.9x[n]-0.45x[n-1]+0.35x[n-2]+0.002x[n-3]。
Q2.23 运行程序P2.6,计算输出序列y[n]和y2[n]以及差值信号d[n]。y[n]和y2[n]相等吗?
Q2.28 运行程序P2.7,对序列h[n]和x[n]求卷积,生成y[n],并用FIR滤波器h[n]对输入x[n]滤波,求得y1[n]。
三、实验器材及软件
1. 微型计算机1台
2. MATLAB 7.0软件
四、实验原理
1.离散时间系统中信号运算方法和基本性质;
2.离散时间系统的MATLAB表示和基本运算;
3. MATLAB处理离散时间系统的方法。
五、实验步骤
根据实验题目要求进行分析;
2. 运用所学知识用MATLAB编程实现题目要求;
3. 对结果进行分析总结。
六、实验记录(数据、图表、波形、程序等)
Q2.1 程序代码:
%产生输入信号
n = 0:100;
s1 = cos(2*pi*0.05*n); %一个低频正弦
s2 = cos(2*pi*0.47*n); %一个高频正弦
x = s1+s2;
%滑动平均滤波器的实现
M = input(Desired length of the filter = );
num = ones(1,M);
y = filter(num,1,x)/M;
clf;
%显示输入和输出信号
subplot(2,2,1);
plot(n, s1);
axis([0, 100, -2, 2]);
xlabel(Time index n); ylabel(Amplitude);
title(Signal #1);
subplot(2,2,2);
plot(n, s2);
axis([0, 100, -2, 2]);
xlabel(Time index n); ylabel(Amplitude);
title(Signal #2);
subplot(2,2,3);
plot(n, x);
axis([0, 100, -2, 2]);
xlabel(Time index n); ylabel(Amplitude);
title(Input Signal);
subplot(2,2,4);
plot(n, y);
axis([0, 100, -2, 2]);
xlabel(Time index n); ylabel(Amplitude);
title(Output Signal); axis;
显示结果:
Q2.2 代码(将输出信号的代码做出修改,其余代码与Q2.1相同):
n=0:100;
s1=cos(2*pi*0.05*n);
s2=cos(2*pi*0.47*n);
x=s1+s2;
M=input( Desired length of the filter = );
num=(-1).^[0:M-1];
y=filter(num,1,x)/M;
clf;
%显示输入和输出信号
subplot(2,2,1);
plot(n, s1);
axis([0, 100, -2, 2]);
xlabel(Time index n); ylabel(Amplitude);
t
您可能关注的文档
- 提高的调频炉炉衬寿命.ppt
- 提高配电自动化建设与运的行经济性.pptx
- 提高的成交率和转介绍率.pptx
- 搜索引擎的的使用技巧.ppt
- 搜索技巧和文件的的下载.ppt
- 搜索算法-B的FS.ppt
- 搜集的华三面试题的及答案整理.pdf
- 搪玻璃反应釜的URS.doc
- 提高院内危重病人转的运安全率.ppt
- 搜索技巧:资的料搜索.ppt
- 计及电动汽车移动储能动态电价的微电网优化调度研究及解决方案.pdf
- 浅谈电动汽车充电桩绝缘智能化自检装置的设计与应用 .pdf
- 浅谈电动汽车公共充电桩布局方案评价方法.pdf
- 浅谈基于弹性响应的电动汽车快充电价定价策略 汽车充电桩有序充电.pdf
- 浅谈光储充一体化社区的有序充电策略及解决方案.pdf
- 晚期肾透明细胞癌系统性治疗中国专家共识(2024版).pptx
- 中国膀胱癌保膀胱治疗多学科诊治协作共识(2022版).pptx
- 成人心血管外科手术体外循环患者血液管理指南.pptx
- 下尿路修复重建移植物应用规范中国专家共识.pptx
- 中国儿童急性非静脉曲张性上消化道出血诊治指南(2024).pptx
文档评论(0)