原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
数字信号处理实验报告96186.doc
重 庆 交 通 大 学
学 生 实 验 报 告
实验课程名称 数字信号处理 开课实验室 数学实验室 学 院 理学院 年级 11 专业班 信息与计算科学2
学 生 姓 名 王家南 学 号 座机电话号码0203 开 课 时 间 2013 至 2014 学年第 1 学期
程序代码：
Convexam.m：
nx 1:120
x impseq 1,0,119 +impseq 2,40,120 +impseq 2,70,120 +impseq 1,80,120 ;
nh 0:24;
h 0.95.^nh;
[y,ny] my_conv x,h ;
subplot 2,1,1 ,stem nx,x ;title 输入信号’
subplot 2,1,2 ,stem ny,y ;title 卷积滤波
impseq.m：
function[x] impseq n0,n1,n2
if n1 n2 error
end
x n0.*[zeros 1, n1-1 ,1,zeros 1, n2-n1 ];
my_conv.m:
function[y,ny] my_conv x,h
ny [0: length x +length h -2 ];
y conv x,h ;
实验结果：
心得体会： 通过本次实验，了解到了conv函数的用法；同时了解到函数名不得和自己取的m文件不能重名否则会出现错误，同时学会了stem函数。
程序代码：
Casfiltr.m:
function y casfiltr b0,B,A,x
[K,L] size B ;
N length x ;
w zeros K+1,N ;
w 1,: x;
for i 1:1:K w i+1,: filter B i,: ,A i,: ,w i,: ;
end
y b0*w K+1,: ;
cas2dir.m:
function[b,a] cas2dir b0,B,A
[K,L] size B ;
b [1];
a [1];
for i 1:1:K b conv b,B i,: ; a conv a,A i,: ;
end
b b*b0;
casfiltrexam.m:
b0 2;
B [1,-0.5,0;1,2,1];
A [1,-0.75,0;1,0,0.81];
N 64;
n 1:N;
x [ n-1 0];
h_cas casfiltr b0,B,A,x ;
subplot 2,1,1
stem n,h_cas,filled,k ;
axis [0,70,-4,6] ;
title ??áaDíμ￥????3??ìó|
[b,a] cas2dir b0,B,A ;
hdir filter b,a,x ;
subplot 2,1,2
stem n,hdir,filled,k ;
axis [0,70,-4,6] ;
title ?±?óDíμ￥????3??ìó|
实验结果：
心得体会：
通过本次实验了解到了级联型转直接性的转变方法，在函数层面上了解了级联转直连的过程，对这两种脉冲响应加深了印象。
程序代码：
Ideal_lp.m:
function hd ideal_lp wc,M
alpha M-1 /2;
n [0:1: M-1 ];
m n-alpha+eps;%加一个小数以避免零做除数没有意义
hd sin wc*m ./ pi*m ;
[db,mag,pha,grd,w] freqz_m hd,[1]
Freqz_m.m:
function [db,mag,pha,grd,w] freqz_m b,a
[H,w] freqz b,a,1000,whole ;
H H 1:1:501 ;w w 1:1:501 ;
mag abs H ;
db 20*log10 mag+eps /max mag ;
pha angle H ;
grd grpdelay b,a,w ;
chuangexam.m
M 121;
n [0:1:M-1];
wc 0.45*pi;
hd ideal_lp wc,M ;
w_rec boxcar M ;
w_ham hamming M ;
h_rec hd.*w_rec;
h_ham hd.*w_ham;
[db1,mag1,pha1,grd1,w1] freqz_m h_rec,[1] ;
[db2,mag2,pha2,grd2,w2] freqz_m h_ham,[1] ;
figure 1
subplot 2,2,1 ;stem n,hd,k ;grid
title 加窗前脉冲响应
axis [0,M-1,-0.2,0.5] ;ylabel hd n
s
文档评论（0）