计算机仿真技术课程结业报告.docx
《计算机仿真技术课程结业报告.docx》由会员分享，可在线阅读，更多相关《计算机仿真技术课程结业报告.docx（10页珍藏版）》请在冰豆网上搜索。
计算机仿真技术课程结业报告
学号：
201210203033
计算机仿真技术课程结业报告
学生姓名
学院：
信息工程学院
班级：
电子12-2
二〇一四年十一月
实验四离散时间LTI系统的时域分析
编程练习
1.试用MATLAB命令求解以下离散时间系统的单位取样响应。
（1）
（2）
2.已知某系统的单位取样响应为
，试用MATLAB求当激励信号为
时，系统的零状态响应。
答案：
a=[341];
b=[11];
n=0:
30;
impz（b,a,30）,gridon
titl（'系统单位取样响应h（n）'）
a=[2.5610];
b=[1];
n=0:
30;
impz（b,a,30）,gridon
titl（'系统单位取样响应h（n）'）
nx=-1:
5;%x（n）向量显示范围（添加了附加的零值）
nh=-2:
10;%h（n）向量显示范围（添加了附加的零值）
x=uDT（nx）-uDT（nx-5）;
h=（7/8）.^nh.*（uDT（nh）-uDT（nh-10））;
y=conv（x,h）;
ny1=nx
（1）+nh
（1）;%卷积结果起始点
%卷积结果长度为两序列长度之和减1,即0到（length（nx）+length（nh）-2）
%因此卷积结果的时间范围是将上述长度加上起始点的偏移值
ny=ny1+（0:
（length（nx）+length（nh）-2））;
subplot（311）
stem（nx,x,'fill'）,gridon
xlabel（'n'）,title（'x（n）'）
axis（[-41603]）
subplot（312）
stem（nh,h','fill'）,gridon
xlabel（'n'）,title（'h（n）'）
axis（[-41603]）
subplot（313）
stem（ny,y,'fill'）,gridon
xlabel（'n'）,title（'y（n）=x（n）*h（n）'）
axis（[-41605]）
在MATLAB中，冲激序列可以通过编写uDT.m文件来实现，即
functiony=uDT（n）
y=n>=0;%当参数为非负时输出1
调用该函数时n也同样必须为整数或整数向量。
实验七IIR数字滤波器的MATLAB实现
1、基于chebyshev1型模拟滤波器原型使用冲激不变转换方法设计数字滤波器，要求具有下面的参数指标：
通带截止频率：
通带波动值：
阻带截止频率：
阻带波动值：
2、一个椭圆数字滤波器的设计，要求采用双线性变换方法，指标参数如下：
通带截止频率：
通带波动值：
阻带截止频率：
阻带波动值：
答案：
wp=0.2*pi;
ws=0.3*pi;
Rp=1;
As=15;
[cs,ds]=cheb1ord（wp,ws,Rp,As）;
[num,den]=cheby1（cs,Rp,ds）;
freqz（num,den）;
[num1,den1]=impinvar（num,den）;
freqz（num1,den1）;
Wp1=0.2*pi;
Ws1=0.3*pi;
Rp=1.5;
As=20;
fs=1;
[cs,ds];
[num,den]=ellip（cs,Rp,As,ds）;
[num1,den1]=bilinear（num,den,fs）;
Freqz（num1,den1）;
实验八FIR数字滤波器的MATLAB实现
1.利用hamming窗设计一个48阶的FIR带通滤波器，通带为Wn＝[0.350.65]。
2.用矩形窗设计一个线性相位高通滤波器。
其中
答案:
wn=[0.350.65]%设置阻带范围
b=fir1（48,wn）;%调用fir1函数
freqz（b）;%绘制滤波器的频率响应曲线
figure
stem（b,'.'）;%绘制单位冲激响应序列
line（[0,25],[0,0]）;xlabel（'n'）;ylabel（'h（n）'）;
tr_width=ws-wp;
M=ceil（6.6*pi/tr_width）+1;
%生成hamming窗
w_ham=（hamming（M））'
wn=
0.35000.6500
w_ham=
Columns1through10
0.08000.09860.15300.23880.34890.47450.60550.73110.84120.9270
Columns11through20
0.98141.00000.98140.92700.84120.73110.60550.47450.34890.2388
Columns21through23
0.15300.09860.0800
升级会员 