数字信号处理西电

数字信号处理上机第一次实验

实验一：

设给定模拟信号()1000t a x t e -=，的单位是ms 。

(1) 利用MATLAB 绘制出其时域波形和频谱图（傅里叶变换），估计其等效带宽（忽略谱分

量降低到峰值的3%以下的频谱）。

(2) 用两个不同的采样频率对给定的进行采样。

○1

。

○2 。

比较两种采样率下的信号频谱，并解释。

实验一MATLAB 程序：

（1）

○

1 clc; fs=5000;

ts=1/fs;

N=1000;

t=(-N:N)*ts;

s=exp(-abs(t));

plot(t,s,'linewidth',1.5)

xlabel('时间')

ylabel('幅度')

set(gca,'fontweight','b','fontsize',12)

SPL=N*100;

figure

sp=fftshift(fft(s,SPL));

sp=sp/max(sp)*100;

freqb=-fs/2:fs/SPL:fs/2-fs/SPL;

plot(freqb,abs(sp))

xlabel('频率')

ylabel('频谱幅度')

set(gca,'fontweight','b','fontsize',12)

yy=abs(abs(sp)-3);

[aa,freqind]=min(yy);

(freqind-SPL/2)*fs/SPL

t ()a x t ()()15000s a f x t x n =以样本秒采样得到。()()

11j x n X e ω画出及其频谱()()11000s a f x t x n =以样本得到。()()11j x n X e ω画出及其频谱

○2 clc;

fs=1000;

ts=1/fs;

N=1000;

t=(-N:N)*ts;

s=exp(-abs(t));

plot(t,s,'linewidth',1.5)

xlabel('时间')

ylabel('幅度')

set(gca,'fontweight','b','fontsize',12) SPL=N*100;

figure

sp=fftshift(fft(s,SPL));

sp=sp/max(sp)*100;

freqb=-fs/2:fs/SPL:fs/2-fs/SPL;

plot(freqb,abs(sp))

xlabel('频率')

ylabel('频谱幅度')

set(gca,'fontweight','b','fontsize',12)

yy=abs(abs(sp)-3);

[aa,freqind]=min(yy);

(freqind-SPL/2)*fs/SPL

实验三：

设，，编写MATLAB 程序，计算：

（１） ５点圆周卷积；

（２） ６点圆周卷积；

（３） 线性卷积；

（４） 画出的，和时间轴对齐。

a = [1,2,2];

b = [1,2,3,4];

y1 = cconv(a,b,5)

y2 = cconv(a,b,6)

y3 = conv(a,b)

figure(1);

subplot(311)

stem(y1);

grid on

title('五点圆周卷积y1(n)');

xlabel('n'),ylabel('y1(n)');axis([0 6 0 15])

subplot(312)

(){}11,2,2x n =(){}21,2,3,4x n =()1y n ()2y n ()3y n ()1y n ()2y n ()3y n

grid on

title('六点圆周卷积y2(n)');

xlabel('n'),ylabel('y2(n)');axis([0 6 0 15])

subplot(313)

stem(y3);

grid on

title('线性卷积y3(n)');

xlabel('n'),ylabel('y3(n)');axis([0 6 0 15])

给定因果系统：

（１） 求系统函数并画出零极点示意图。

（２） 画出系统的幅频特性和相频特性。 （３） 求脉冲响应并画序列图。

提示：在ＭＡＴＬＡＢ中，zplane(b,a) 函数可画零极点图；Freqz(b,a,N)可给出范围内均匀间隔的点频率响应的复振幅；Impz(b,a,N)可求的逆变换（即脉冲响应）。 clc

a = [1,0]

b = [1,-0.9]

figure(1)

zplane(b,a);

title('零极点分布图')

w=[-3*pi:0.01:3*pi];

[h,phi]=freqz(b,a,w);

figure(2);

subplot(2,1,1);

plot(w, abs(h));

grid on;

title('幅频特性');

xlabel('f/Hz'),ylabel('H(f)');

subplot(2,1,2);

plot(w, phi);

grid on;

title('相频特性');

xlabel('f/Hz'),ylabel('W(f)');

()()()0.91y n y n x n =-+()H z ()j H e ω

()ϕω()h n []0,πN ()H z