基于matlba的时域信号采样及频谱分析.docVIP

基于 MATLAB ?的时域信号采样及频谱?分析 默认分类 ??2009-05-26? 21:46 ??阅读900???评论0 ? 字号： 大大? 中中? 小小 基于 ?MATLAB 的时域?信号采样及频谱分析 ?? 二、基本要求： ?① 掌握数字信号处理?的基本概念、基本理论?和基本方法； ② 学?会 MATLAB 的?使用，掌握 MATL?AB 的程序设计方法?； ③ 学会用 MA?TLAB 对信号进行?分析和处理； ④ 信?号的各参数需由键盘输?入，输入不同参数即可?得不同的x(t) 和?x(n)； ⑤ 撰写?课程设计论文，用数字?信号处理基本理论分析?结果。 ? 三、设计?方法与步骤： ① 画?出连续时间信号 的时?域波形及其幅频特性曲?线，其中幅度因子A＝?444.128，衰减?因子a＝222.14?4，模拟角频率 ＝2?22.144； ② ?对信号 进行采样，得?到采样序列 ，其中T?＝ 为采样间隔，通过?改变采样频率可改变T?，画出采样频率分别为?200Hz，500 ?Hz，1000 Hz?时的采样序列波形； ?③ 对不同采样频率下?的采样序列进行频谱分?析，绘制其幅频和相频?曲线，对比各频率下采?样序列 和 的幅频曲?线有无差别，如有差别?说明原因。 ④ 设系?统单位抽样响应为 ，?求解当输入为 时的系?统响应 ，画出 , ?, 的时域波形及幅频?特性曲线，并利用结果?验证卷积定理的正确性?（此内容将参数设置为?A＝1，a＝0.4，? ＝2.0734，T?＝1）。 ⑤ 用FF?T 对信号 , , ?进行谱分析，观察与④?中结果有无差别。 ⑥? 由采样序列 恢复出?连续时间信号 ，画出?其时域波形，对比 与?原连续时间信号 的时?域波形，计算并记录两?者最大误差。 ? ?? 四、详细程序及仿真?波形分析： 1、连续?时间信号x(t)及其?200Hz/500H?z/1000Hz频率?抽样信号函数x(n)? %绘制信号x(n)?的幅度谱和相位谱 n?=0:50; ? ? ? %定义序列的长度是?50 A=input?(请输入A的值 A?:); %设?置信号的有关参数 a?=input(请输?入a的值 a:);? w0=input(?请输入w0的值 w?0:); T1=0?.005; T2=0?.002; T3=0?.001; T0=0?.001; x=A*?exp(-a*n*T?0).*sin(w0?*n*T0); %?pi 是MATLAB? 定义的π，信号乘可?采用“.*” y1=?A*exp(-a*n?*T1).*sin(?w0*n*T1); ? %pi 是MATL?AB 定义的π，信号?乘可采用“.*” y?2=A*exp(-a?*n*T2).*si?n(w0*n*T2)?; %pi ?是MATLAB 定义?的π，信号乘可采用“?.*” y3=A*e?xp(-a*n*T3?).*sin(w0*?n*T3); ? %pi 是MATL?AB 定义的π，信号?乘可采用“.*” c?lose all ? ? ? %清除已经绘?制的x(n)图形 s?ubplot(2,1?,1);stem(n?,x),grid o?n %绘制x(n)?的图形 title(?离散时间信号) ?subplot(2,?1,2);plot(?n,x),grid ?on title(?连续时间信号) ?figure(2) ?subplot(3,?1,1);stem(?n,y1),grid? on title(?200Hz理想采样?信号序列); ? %设置结?果图形的标题 sub?plot(3,1,2?);stem(n,y?2),grid on? title(50?0Hz连续时间信号?) subplot(?3,1,3);ste?m(n,y3),gr?id on titl?e(1000Hz连?续时间信号) k=?-25:25; W=?(pi/12.5)*?k; w=W/pi;? Y1=y1*(e?xp(-j*pi/1?2.5)).^(n?*k); figur?e(3) subpl?ot(2,1,1);?plot(w,abs?(Y1));grid?,xlabel(w?),ylabel(?幅度); tit?le(200Hz理?想采样信号序列的幅度?谱); axis(?[-2 2 0 10?00]); subp?lot(2,1,2)?;plot(w,an?gle(Y1));g?rid,xlabel?(w),ylab?el(幅角); ?title (20?0Hz理想采样信号序?列的相位谱) Y2?=y2*(exp(-?j*pi