- 1、本文档共8页,可阅读全部内容。
- 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
信号与系统的频域分析专题研讨
【目的】
(1) 加深对信号与系统频域分析基本原理和方法的理解。
(2) 学会利用信号抽样的基本原理对信号抽样过程中出现的一些现象的进行分析。
(3) 通过实验初步了解频谱近似计算过程中产生误差的原因。
(4)学会用调制解调的基本原理对系统进行频域分析。
【研讨题目】
1.信号的抽样
频率为f0 Hz的正弦信号可表示为
按抽样频率fsam=1/T对x(t)抽样可得离散正弦序列x[k]
在下面的实验中,取抽样频率fsam=8kHz。
(1)对频率为2kHz, 2.2 kHz, 2.4 kHz和 2.6 kHz正弦信号抽样1 秒钟,利用MATLAB函数 sound(x, fsam)播放这四个不同频率的正弦信号。
(2)对频率为5.4 kHz, 5.6kHz, 5.8 kHz和 6.0kHz正弦信号抽样1 秒钟,利用MATLAB函数 sound(x, fsam)播放这四个不同频率的正弦信号。
(3)比较(1)和(2)的实验结果,解释所出现的现象。
【题目分析】
【信号抽样过程中频谱变化的规律】
【比较研究】
连续的播放两段音频信号,比较函数sound和wavplay的异同。
【仿真结果】
【结果的理论分析和解释】
【自主学习内容】
【阅读文献】
【发现问题】 (专题研讨或相关知识点学习中发现的问题):
【问题探究】
【仿真程序】
2. 连续时间信号Fourier变换的数值近似计算
计算连续信号频谱是对信号和系统进行频域分析的基础,由于实际信号大多无简单的解析表达式,所以要用数值方法进行近似计算。本题要求对频谱近似计算中误差的原因进行初步的分析,希望能在计算实际信号频谱的近似计算中起一定的指导作用。
若信号x(t)的非零值在区间,则可用下面提供的函数ctft1或ctft2近似计算其频谱。函数ctft的调用形式为
[X,f]=ctft1(x,fsam,N)
[X,f]=ctft2(x,fsam,N)
其中调用变量x 存放信号x(t)的抽样值,fsam表示对连续信号x(t)的抽样频率(Hz),N表示用DFT进行近似计算时DFT 的点数,为了能高效的进行计算,N最好取2的整数次幂,如512, 1024等。返回变量X是计算出的信号频谱的抽样值,f(单位Hz)表示对应的频率抽样点。返回变量X一般是复数,可用函数abs(X)计算出幅度谱,函数angle(X) 计算出相位谱。
(1)阅读程序ctft2,叙述该程序的基本原理。该程序中有一处需要产生一个大的2维矩阵,指出该行程序,并评价该方法的优缺点。
(2)取抽样频率fsam=100Hz, 信号抽样长度N=1024, 分别用两个子程序近似计算信号的频谱,比较两种方法的计算时间和误差;
(3)若将信号的时域有效宽度定义为
其中表示信号在时域的最大值。试分析时域有效宽度对近似计算的影响。给出一个由信号时域有效宽度估计近似计算中所需信号长度的经验公式。
(4)定义信号频域有效宽度为
其中表示信号在频域的最大值。给出一个由信号频域有效宽度估计近似计算中所需抽样频率的经验公式。
(5)用计算机录分别一段男生和女生的语音信号,计算其频谱并比较其特点。
(6)讨论:计算误差产生的主要原因?如果不知信号的解析表达式,如何分析计算误差?
【题目分析】
由傅里叶变换,对应1/(1+jw),则幅度为1/√(1+W2);arctan(-w)为相位;
信号时域有效宽度估计近似计算中所需信号长度,必须要包含信号的有效宽度,而该信号为单调递减信号,则若X(t1)=Xmax; 0.1 Xmax = X(t2);=t2-t1;
根据幅度谱的单调性可知,必然对应于两个大于0的频率f1,f2;那么f的范围就是f1:f2;则fs≧2*(f2-f1);
(1):
由题目所给的有ctgt1和ctft2这两个程序的形式:
%近似计算连续信号频谱的函数
Ctft1:
function [X,f]=ctft1(x,Fs,N)
X=fftshift(fft(x,N))/Fs;
f=-Fs/2+(0:N-1)*Fs/N;
ctft2:
function [X,f]=ctft2(x,Fs,N)
tk=(0:N-1)/Fs;% 时域抽样点
dF=Fs/N; %频域抽样间隔
fm=(0:N/2)*dF; %频域抽样点
X=x*exp(-j*2*pi*tk*fm)/Fs;% 近似计算信号频谱
f=[-fliplr(fm(2:end)) fm]; % 增添负频率点
X=[-conj(fliplr(X(2:end))) X];% 增添频率点对应的频谱
用该程序得到的幅度频谱比较准确跟书上的图相似,但是得到的相位相应存在很大的误差。不如ctft1得到的准确。
(2)
源程序:
%The first fun
文档评论(0)