数字信号处理matlab实验3太原理工大学.docx
《数字信号处理matlab实验3太原理工大学.docx》由会员分享，可在线阅读，更多相关《数字信号处理matlab实验3太原理工大学.docx（10页珍藏版）》请在冰豆网上搜索。
数字信号处理matlab实验3太原理工大学
课程名称：
数字信号处理
实验项目：
实验一
实验地点：
化工楼二层
专业班级：
学号：
学生姓名：
指导教师：
2014年月日
实验一
1、单位抽样序列
n=input（'请输入单位抽样序列的n值"x1=[zeros（1,N）1zeros（1,N）]"：
'）
a=input（'请输入单位抽样序列的a值"偏移量a"：
'）
N=n
x=-（n）:
（n）
y=x+a
x1=[zeros（1,N）1zeros（1,N）];
stem（y,x1）;
xlabel（'时间序列n'）;ylabel（'振幅'）;title（'单位抽样序列x1'）;
请输入单位抽样序列的n值"x1=[zeros（1,N）1zeros（1,N）]"：
5
请输入单位抽样序列的a值"偏移量a"：
3
2、单位阶越序列
n=input（'请输入单位抽样序列的n值"x1=[zeros（1,N）1zeros（1,N）]"：
'）
a=input（'请输入单位抽样序列的a值"偏移量a"：
'）
N=n
x=-（n）:
（n）
y=x+a
x1=[zeros（1,n）1ones（1,n）];
stem（y,x1）;
xlabel（'时间序列n'）;ylabel（'振幅'）;title（'单位阶跃序列x1'）;
请输入单位抽样序列的n值"x1=[zeros（1,N）1zeros（1,N）]"：
6
请输入单位抽样序列的a值"偏移量a"：
2
3、正弦序列
N=input（'N='）
A=input（'A='）
n=0:
N-1
f=100
Fs=2000
fai=0
x=A*sin（2*pi*f*n/Fs+fai）
plot（x）
gridon
4、复指数序列
n=input（'请输入复指数序列的n值"x1=R*exp（j*w*N）"：
'）
R=input（'请输入复指数序列的R值"R"：
'）
w=input（'请输入复指数序列的w值"w"：
'）
N=1:
n
x1=R*exp（j*w*N）;
stem（N,x1）;
xlabel（'时间序列n'）;ylabel（'振幅'）;title（'复指数序列x1'）;
请输入复指数序列的n值"x1=R*exp（j*w*N）"：
14
n=
14
请输入复指数序列的R值"R"：
2
R=
2
请输入复指数序列的w值"w"：
4
5、指数序列
n=input（'请输入指数序列的n值"x1=A^N"：
'）
A=input（'请输入指数序列的R值"A"：
'）
N=1:
n
x1=A.^N;
stem（N,x1）;
xlabel（'时间序列n'）;ylabel（'振幅'）;title（'指数序列x1'）;
请输入指数序列的n值"x1=A^N"：
13
请输入指数序列的R值"A"：
3
练习二
2、
（1）单位脉冲响应
N=input（'N='）
b=[1-1]
a=[10.750.125]
x=[zeros（1,N-1）,1,zeros（1,N-1）]
t=（1-N）:
（N-1）
y=filter（b,a,x）
stem（t,y）
gridon
N=20
单位阶跃响应
N=input（'N='）
b=[1-1]
a=[10.750.125]
x=[zeros（1,N）,1,ones（1,N）]
t=（-N）:
（N）
y=filter（b,a,x）
stem（t,y）
gridon
N=20
（2）单位脉冲响应
N=input（'N='）
b=[00.250.250.250.25]
a=[1]
x=[zeros（1,N-1）,1,zeros（1,N-1）]
t=（1-N）:
（N-1）
y=filter（b,a,x）
stem（t,y）
gridon
N=20
单位阶跃响应
N=input（'N='）
b=[00.250.250.250.25]
a=[1]
x=[zeros（1,N）,1,ones（1,N）]
t=（-N）:
（N）
y=filter（b,a,x）
stem（t,y）
gridon
升级会员 