- 1、本文档共6页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
文档全是精心选出来的精品资料,绝对值得下载!
南昌航空大学实验报告
2012 年 04 月 26 日
课程名称: 数字信号处理 实验名称:系统及系统响应
班级: 090423班 学号: 姓名: 张 丽
实验二 系统及系统响应
实验目的
观察离散系统的频率响应;
熟悉连续信号经理想采样前后的频谱变化关系,加深对时域采样定理的理解;
利用序列的FT对连续信号、离散信号及系统响应进行频域分析;
利用卷积方法观察分析系统的时域特性。
实验内容
给定一因果系统H(z)=1+2z-1+z-21-0.67z-1+0.9z-2 ,求出并绘制H(z)的幅频响应与相频响应;
对信号xa(t)=Ae-atsinΩ0πTu(n) 0≤n50 其中A=444.128,a=502π,Ω0=502π, 实现下列实验内容:
取采样频率fs=10KHZ,观察所得采样xa(n)的幅频特性|X(ejω)|和图中的|Xa(jΩ)|在折叠频率附近有无明显差别。
改变采样频率fs=1KHZ,观察|X(ejω)|的变化,并作记录:进一步降低采样频率,fs=300HZ,观察频谱混叠是否明显存在,说明原因,并记录这时的|X(ejω)|曲线。
给定系统的单位抽样响应为h1(n)=R10(n)
利用线性卷积求信号x1(n)=δ(n),通过该系统的响应y1(n)。比较所求响应y1(n)和h1(n)之间有无差别,绘图说明,并用所学理论解释所得结果。
利用线性卷积求信号x2(n)=R10(n),通过该系统的响应y2(n),并判断y2(n)图形及其非0值序列长度是否与理论结果一致,改变x2(n)的长度,取N=5,重复该试验。注意参数变化的影响,说明变化前后的差异,并解释所得结果。
求x(n)=11δ(n+2)+7δ(n+1)- δ(n-1)+4δ(n-2)+2δ(n-3)通过系统h(n)=2δ(n+1)+3δ(n)-5δ(n-2)+2δ(n-3)+ δ(n-4)的响应y(n)。
实验程序及解析
程序
clear; close all;
b=[1,sqrt(2),1];a=[1,-0.67,0.9];
[h,w]=freqz(b,a);
am=20*log10(abs(h));% am=20*log10(abs(h))为幅频响应取dB
subplot(2,1,1);%将窗口划分为2*1的小窗口并选择第一个显示
plot(w,abs(h));xlabel(w);ylabel(幅频响应);title(系统响应)
ph=angle(h);
subplot(2,1,2); %选择第二个窗口显示
plot(w,ph);xlabel(w);ylabel(相频响应);
2、系统响应结果
图1 因果系统的H(z)的系统响应
3、结果分析
分析z域系统的特性主要是由系统的零点和极点的分布得出结论的。首先求出系统的极点和零点,再画出零极图,由理论得出的结论是:
系统的幅频特性是零极???上各零点到单位圆上某点的向量的积与各极点到单位圆上该点的向量的积之比,而幅频特性曲线就是在单位圆上的那一点从0到2π变化的结果。
很显然,从系统的幅频特性响应来看,该系统明显是一个带通滤波器。
(2)1、程序
function X=dtft(x,w)
X=x*exp(-j*[0:(length(x)-1)]*w);%离散傅立叶变换的函数定义
fs=10000;fs1=1000;fs2=300;%设置三种采样频率
t=0:1/fs:0.1; %采集信号长度为0.1s
A=444.128;a=50*sqrt(2)*pi;b=a;
xa=A*exp(-a*t).*sin(b*t);
k=0:511;f=fs*k/112; %序列长度为512,即N=512,还有模拟频率f
Xa=dtft(xa,2*pi*k/512); %离散傅立叶变换
T1=1/fs1;t1=0:T1:0.1; %采集信号长度为0.1s
x1=A*exp(-a.*t1).*sin(b*t1); %1KHz的采样序列x1(n)
X1=dtft(x1,2*pi*k/512); %x1(n)的512点傅立叶变换
T2=1/fs2;t2=0:T2:0.1; %采集信号长度为0.1s
x2=A*exp(-a.*t2).*sin(b*t2);
X2=dtft(x2,2*pi*k/512);
figure(1);
subplot(6,1,1);
plot(t,xa); %画出原始波形
axis([0,max(t),min(xa),max(xa)]);
title(模拟信号);
xlabel(t(s));yl
您可能关注的文档
最近下载
- 四川省红十字应急救护知识竞赛考试题库资料(含答案).pdf
- 新沪教版高中信息技术必修1《解决温标转换问题——认识程序和程序设计语言》说课稿.doc
- 基于功能磁共振成像的脑网络研究.pptx
- 码头水利技术规定条文2010-09.pdf
- 中医院临床药师培训理论考核试题(抗感染专业).doc.docx
- 学校生活-备战高考英语话题通关精准练(词汇+阅读+写作)(解析版).docx VIP
- 2024年初级会计职称《初级会计实务》精讲课件 第1-5章.pdf
- 史陶比尔机器人TX2-90用户手册.pdf
- 广西医科大学附属肿瘤医院临床药师培训基地学员理论考核(抗感染专业).docx
- 2024年-煤矿安全规程(精编).doc
文档评论(0)