注册 登录  
 加关注
   显示下一条  |  关闭
温馨提示!由于新浪微博认证机制调整,您的新浪微博帐号绑定已过期,请重新绑定!立即重新绑定新浪微博》  |  关闭

我的博客

 
 
 

日志

 
 

数字信号处理实验报告——IIR  

2007-12-05 23:39:59|  分类: 默认分类 |  标签: |举报 |字号 订阅

  下载LOFTER 我的照片书  |

1.实验目的

2.实验环境(软件、硬件及条件)

   Matlab6.5及其以上版本

3.实验方法

4.实验分析

 

5.实验结论

  思考题解答:

         对双线性变换法,不存在频谱混叠现象,所以T的取值对设计结果没有影响,可以任意选择。      

结论:

用双线性变换法设计;(1)确定数字低通滤波器的技术指标;(2)将数字低通滤波器的技术指标转化成模拟低通滤波器的技术指标,按关系式W=2*tan(w/2)/T转换;(3)按照模拟低通滤波器的技术指报设计模拟低通滤波器,(4)将模拟滤波器H(s),从s平面转换到z平面,得到数字低通滤波器。

按公式s=2*(1-z-1)/(T*(1+z-1))转换。

附件

1.

实验第一步代码及结果

function iir()         %ʵÑéÈýµÚÒ»²½

w=0:0.02*pi:pi/2;

z=exp(j*w);

A=0.09036;

B1=1.2686;

B2=1.0106;

B3=0.9044;

C1=-0.7051;

C2=-0.3583;

C3=-0.2155;

H1=A*(1+2*z.^(-1)+z.^(-2))./(1-B1*z.^(-1)-C1*z.^(-2));

H2=A*(1+2*z.^(-1)+z.^(-2))./(1-B2*z.^(-1)-C2*z.^(-2));

H3=A*(1+2*z.^(-1)+z.^(-2))./(1-B3*z.^(-1)-C3*z.^(-2));

H=H1.*H2.*H3;                %Â˲¨Æ÷µÄ·ùƵÇúÏß

figure;

plot(w,abs(H));

grid on;

title('iir·ùƵÏìÓ¦ÇúÏß');

xlabel('w');

ylabel('|H(exp(jw))|');

%ÐĵçͼÐòÁÐ

x1=[-4,-2,0,-4,-6,-4,-2,-4,-6,-6];

x2=[-4,-4,-6,-6,-2,6,12,8,0,-16];

x3=[-38,-60,-84,-90,-66,-32,-4,-2,-4,8];

x4=[12,12,10,6,6,6,4,0,0,0];

x5=[0,0,-2,-4,0,0,0,-2,-2,0,0,-2,-2,-2,-2,0];

x=[x1,x2,x3,x4,x5];

 

a1=[1,-B1,-C1];

a2=[1,-B2,-C2];

a3=[1,-B3,-C3];

b1=[A,2*A,A];

b2=b1;

b3=b1;

 

figure;

subplot(2,2,1);

stem(x,'.');

title('x(n)²¨ÐÎ');

xlabel('n');

ylabel('x(n)');

%µÚÒ»¼¶Â˲¨

y1=filter(b1,a1,x);

subplot(2,2,2);

stem(y1,'.');

title('y1(n)²¨ÐÎ');

xlabel('n');

ylabel('y1(n)');

%µÚ¶þ¼¶Â˲¨

y2=filter(b2,a2,y1);

subplot(2,2,3);

stem(y2,'.');

title('y2(n)²¨ÐÎ');

xlabel('n');

ylabel('y2(n)');

%µÚÈý¼¶Â˲¨

y=filter(b3,a3,y2);

subplot(2,2,4);

stem(y,'.');

title('y(n)²¨ÐÎ');

xlabel('n');

ylabel('y(n)');

%ƵÓò±È½Ï

figure;

subplot(2,1,1);

plot(abs(fft(x,225)));

subplot(2,1,2);

plot(abs(fft(y,225)));

 

  评论这张
 
阅读(1043)| 评论(1)
推荐 转载

历史上的今天

评论

<#--最新日志,群博日志--> <#--推荐日志--> <#--引用记录--> <#--博主推荐--> <#--随机阅读--> <#--首页推荐--> <#--历史上的今天--> <#--被推荐日志--> <#--上一篇,下一篇--> <#-- 热度 --> <#-- 网易新闻广告 --> <#--右边模块结构--> <#--评论模块结构--> <#--引用模块结构--> <#--博主发起的投票-->
 
 
 
 
 
 
 
 
 
 
 
 
 
 

页脚

网易公司版权所有 ©1997-2017