重庆大学光电学院数字信号处理实验
五、实验过程原始记录(实验源程序、图标、计算等) 1、用matlab 编制信号产生子程序,产生以下典型信号序列。
()100处有一单位脉冲,1在,)(单位脉冲序列1≤≤=n n
n δ
()()100,,1处及以后为4,在0前为4在,单位阶跃序列2≤≤==n n n
n u
()()n
R
8
矩形序列
3
()()9
3
5
sin
3
正弦序列
4≤
≤
??
?
?
?
?
+
?
=n
n
n
x
π
π
()()90)4(5)3(4)2(3)1(2)(序列5≤≤-+-+-+-+=n n n n n
n n x δδδδδ
()()90)3(2)2()1(2)(序列6≤≤-+-+-+=n n n
n n n h δδδδ
2、利用MATLAB 编程完成序列(5)和(6)的移位、反转、加法、乘法等运算,并绘制运算后序列的波形。
()())4(5)3(4)2(3)1(2)(序列5-+-+-+-+=n n n n n n x δδδδδ ()())3(2)2()1(2)(序列6-+-+-+=n n n
n n h δδδδ
3、利用MATLAB编制一个计算两个序列线性卷积的通用程序,计算上述两序列,并绘制卷积后序列的波形。
六、实验结果分析
1.如何产生方波信号序列和锯齿波信号序列?
答:方波信号由square函数可得,锯齿波信号由循环语句产生。
%方波
T=0:0.001:2*pi; %方波的时间向量
y=square(6*T); %方波w=6,周期为T=2*pi/6
subplot(2,1,1);
plot(T,y); %画方波
title('方波');
axis([T(1)-1 T(end)+1 -2 2]); %坐标轴区域
hold on
一个序列
x (n )的离散时间傅里叶变换就是它的频谱函数。
()()
()1
20
N j nk N k
n N
X k X e
x n W ω
πω-====∑
四、使用仪器、器材
计算机,MATLAB
五、实验过程原始记录(实验源程序、图标、计算等)
(1)用matlab 编程产生并画出信号x1(n)、x2(n)、x3(n)、x4(n)、x5(n)。
??
?
??≤≤-≤≤+==n
n n n n n x n R n x 其它,074,
830,1)()
()(241???
??≤≤-≤≤-=n
n n n n n x 其它,
074,
330,
4)(3 n n x 4
cos
)(4π
=n n x 8
sin
)(5π
=
(2)用matlab 编制FFT 函数对上述信号进行频谱分析,并画出上述信号谱图。