实验五实验报告分析解析.docx

文档介绍：综合性、设计性实验报告
姓名陶建敏学号201308002107专业通信工程班级通信1301实验课程名称数字信号处理指导教师及职称李玲香讲师开课学期2015至＿2016学年上学期上课时间2015年12月10、17日湖南科技学院教务处编印
设计任务确定数字滤波器性能指标要求，低通滤波器两个边界频率：
Wp和Ws两个衰减：Rp=1,Rs=60,Wp=2*pi*fp,fp=200,Ws=2*pi*fsw,fsw=300带通滤波器上下线通带截止频率分别为：fp1=250,fp2=400,两个阻带截止频率fs1=150,fs2=500;两个衰减：rp=3,rs=40;：将数字滤波器指标转换成模拟滤波器指标；（6）:归一化模拟滤波器；（7）；求出实际模拟滤波器系统函数，去归一化；（8）：利用冲击响应不变法离散化，将模拟滤波器转化为数字滤波器a:模拟低通到数字低通直接离散化；b数字带通滤波及设计，则要将模拟低通变换为模拟带通，再将其离散化为数字带通滤波器；（9）：用频响特性函数检验是否符合。
6．MATLAB编程实现
clf;
fs1=1000;N=128
%
采样频率和数据点数
n=0:N-1;t=n/fs1;
%
时间序列
x=2*sin(300*pi*t)+sin(600*pi*t);%
信号
y=fft(x,N);
%
对信号进行快速Fourier
变换
mag=abs(y);
%
求得Fourier
变换后的振幅
f=n*fs1/N;
%
频率序列
figure(1),plot(f,mag);%
绘出随频率变化的振幅
xlabel('f/Hz');ylabel('
幅度');title('
原信号x（t）幅度谱');
T=;
%
采样周期
Rp=1;
%Rp
为通带最大衰减,dB
Rs=60;
%Rs
为阻带最小衰减，dB
fp=200;
%
通带截止频率
fsw=300
%
阻带截止频率
fs=1/T;%采用频率Wp=2*pi*fp;%Wp为通带边界频率，rad/sWs=2*pi*fsw;%Ws为阻带边界频率，rad/s%设计模拟巴特沃斯低通滤波fp=200;fs=300;Fs=2000;T=1/Fs;Rp=1;Rs=60;Wp=2*Fs*tan(2*fp*pi*Ts/2);Ws=2*Fs*tan(2*fs*pi*Ts/2);[N,Wc]=buttord(Wp,Ws,Rp,Rs,'s');[Z,P,K]=buttap(N);[Q,P]=zp2tf(Z,P,K)[BT,AT]=butter(N,Wc,'s');[Bd,Ad]=bilinear(BT,AT,Fs);%绘制频率响应曲线[H,W]=freqz(Bd,Ad);figure(2)plot(W*Fs/(2*pi),abs(H));
xlabel('频率/Hz');ylabel('幅度');title('
低通滤波器')
n=0:299;t=n*T;x2=2*sin(300*pi*t)+sin(600*pi*t);%绘制频率响应曲线y=filter(Bd,Ad,x2)figure(3)