实验5 抽样定理
一、实验目的：
1、了解用MA TLAB 语言进行时域、频域抽样及信号重建的方法。
2、进一步加深对时域、频域抽样定理的基本原理的理解。
3、观察信号抽样与恢复的图形，掌握采样频率的确定方法和内插公式的编程方法。 二、实验原理：
1、时域抽样与信号的重建 （1）对连续信号进行采样
例5-1 已知一个连续时间信号sin sin(),1Hz 3
ππ=0001f(t)=(2f t)+6f t f ，取最高有限带宽频率f m =5f 0，分别显示原连续时间信号波形和F s >2f m 、F s =2f m 、F s <2f m 三情况下抽样信号的波形。
程序清单如下：
%分别取Fs=fm ，Fs=2fm ，Fs=3fm 来研究问题 dt=0.1; f0=1; T0=1/f0; m=5*f0; Tm=1/fm; t=-2:dt:2;
f=sin(2*pi*f0*t)+1/3*sin(6*pi*f0*t); subplot(4,1,1); plot(t,f);
axis([min(t),max(t),1.1*min(f),1.1*max(f)]); title('原连续信号和抽样信号'); for i=1:3;
fs=i*fm;Ts=1/fs; n=-2:Ts:2;
f=sin(2*pi*f0*n)+1/3*sin(6*pi*f0*n); subplot(4,1,i+1);stem(n,f,'filled');
axis([min(n),max(n),1.1*min(f),1.1*max(f)]); end
程序运行结果如图5-1所示：
原连续信号和抽样信号
图5-1
（2）连续信号和抽样信号的频谱
由理论分析可知，信号的频谱图可以很直观地反映出抽样信号能否恢复原模拟信号。因此，我们对上述三种情况下的时域信号求幅度谱，来进一步分析和验证时域抽样定理。
例5-2编程求解例5-1中连续信号及其三种抽样频率（F s>2f m、F s=2f m、F s<2f m）下的抽样信号的幅度谱。
程序清单如下：
dt=0.1;f0=1;T0=1/f0;fm=5*f0;Tm=1/fm;
t=-2:dt:2;N=length(t);
f=sin(2*pi*f0*t)+1/3*sin(6*pi*f0*t);
wm=2*pi*fm;k=0:N-1;w1=k*wm/N;
F1=f*exp(-j*t'*w1)*dt;subplot(4,1,1);plot(w1/(2*pi),abs(F1));
axis([0,max(4*fm),1.1*min(abs(F1)),1.1*max(abs(F1))]);
for i=1:3;
if i<=2 c=0;else c=1;end
fs=(i+c)*fm;Ts=1/fs;
n=-2:Ts:2;N=length(n);
f=sin(2*pi*f0*n)+1/3*sin(6*pi*f0*n);
wm=2*pi*fs;k=0:N-1;
w=k*wm/N;F=f*exp(-j*n'*w)*Ts;
subplot(4,1,i+1);plot(w/(2*pi),abs(F));
axis([0,max(4*fm),1.1*min(abs(F)),1.1*max(abs(F))]);
end
程序运行结果如图5-2所示。
由图可见，当满足F s≥2f m条件时，抽样信号的频谱没有混叠现象；当不满足F s≥2f m 条件时，抽样信号的频谱发生了混叠，即图5-2的第二行F s＜2f m的频谱图，，在f m=5f0的范围内，频谱出现了镜像对称的部分。
02
4
6
8
10
12
14
16
18
20
02
4
6
8
10
12
14
16
18
20
02
4
6
8
10
12
14
16
18
20
2
4
6
8
10
12
14
16
18
20
图5-2
（3）由内插公式重建信号 信号重建一般采用两种方法：一是用时域信号与理想滤波器系统的单位冲激响应进行卷积积分；二是用低通滤波器对信号进行滤波。本实验只讨论第一种方法。
由理论分析可知，理想低通滤波器的单位冲激响应为
j Ωt πt sin()1T h(t)=H(j Ω)e d Ω=
πt 2πT
∞-∞
? 抽样信号a ?x
(t)通过滤波器输出，其结果应为a ?x (t)与h(t)的卷积积分： sin[()/]
??()()()()()()()
()/a a a a a n t nT T y t x t x
t h t x h t d x nT t nT T πτττπ∞
∞
-∞
=-∞
-==*=-=-∑
?
该式称为内插公式。由式可见，x a (t)信号可以由其抽样值x a (nT)及内插函数重构。MATLAB 中提供了sinc 函数，可以很方便地使用内插公式。
例5-3 用上面推导出的内插公式重建例5-1给定的信号。 程序清单如下：
dt=0.01;f0=1;T0=1/f0;fm=5*f0;Tm=1/fm;
t=0:dt:3*T0;x=sin(2*pi*f0*t)+1/3*sin(6*pi*f0*t);
subplot(4,1,1);plot(t,x);axis([min(t),max(t),1.1*min(x),1.1*max(x)]); title('用时域卷积重建抽样信号'); for i=1:3;
fs=i*fm;Ts=1/fs; n=0:(3*T0)/Ts; t1=0:Ts:3*T0;
x1=sin(2*pi*f0*n/fs)+1/3*sin(6*pi*f0*n/fs);
T_N=ones(length(n),1)*t1-n'*Ts*ones(1,length(t1)); xa=x1*sinc(fs*pi*T_N); subplot(4,1,i+1);plot(t1,xa);
axis([min(t1),max(t1),1.1*min(xa),1.1*max(xa)]); end
程序运行结果如图5-3所示：
用时域卷积重建抽样信号
00.5
1
1.5
2
2.5
3
图5-3
2、频域抽样与信号恢复 （1）频域抽样定理
从理论学习可知，在单位圆上对任意序列的Z 变换等间隔采样N 点得到：
2πj k N
2πj
nk N
z=e n=-X(k)=X (z)
=x(n)e
∞
∞
∑ k=0,1,…,N-1
该式实现了序列在频域的抽样。
那么由频域的抽样得到的频谱的序列能否不失真地恢复原时域信号呢？ 由理论学习又知，频域抽样定理由下列公式表述：
r=-x(n)=
x(n+rN)∞
∞
∑% 表明对一个频谱采样后经IDFT 生成的周期序列x(n)
%是原非周期序列x(n)的周期延拓序列，其时域周期等于频域抽样点数N 。
假定有限长序列x(n)的长度为M ，频域抽样点数为N ，原时域信号不失真地由频域抽样恢复的条件如下：
① 如果x(n)不是有限长序列，则必然造成混叠现象，产生误差；
② 如果x(n)是有限长序列，且频域抽样点数N 小于序列长度M （即N N 为周期进行延拓也将造成混叠，从x(n) %中不能无失真地恢复出原信号x(n)。 ③ 如果x(n)是有限长序列，且频域抽样点数N 大于或等于序列长度M （即N ≥M ）， 则从x(n) %中能无失真地恢复出原信号x(n)，即 N N N N r=-x (n)=x (n)R (n)=x(n+rN)R (n)=x(n)∞ ∞ ∑% （2）从频谱抽样恢复离散时间序列 例5-4 已知一个时间序列的频谱为 j ω -j ωn -j ω-j2ω-j3ω-j4ωn=-X(e )=x(n)e =3+2e +e +2e +3e ∞ ∞ ∑ 用IFFT 计算并求出其时间序列x(n)，并绘图显示时间序列。 分析：该题使用了数字频率，没有给出采样周期，则默认Ts=1S,另外，从j ω X(e )的解析式可以直接看出时域序列xn=[3,2,1,2,3]。但为说明问题，仍编写程序求解如下： 程序清单如下： Ts=1;N0=[3,5,10]; for r=1:3; N=N0(r); D=2*pi/(Ts*N); kn=floor(-(N-1)/2:-1/2); kp=floor(0:(N-1)/2); w=[kp,kn]*D; X=3+2*exp(-j*w)+1*exp(-j*2*w)+2*exp(-j*3*w)+3*exp(-j*4*w); n=0:N-1; x=ifft(X,N) subplot(1,3,r);stem(n*Ts,abs(x)); box end 程序运行结果如图5-4所示： 012 02 4 图5-4 注意：程序中数字频率的排序进行了处理，这是因为j ω X(e )的排列顺序是从0开始，而不是从-(N-1)/2开始。 程序运行后将显示数据： x=5.0000 5.0000 1.0000 x=3.0000 2.0000 1.0000 2.0000 3.0000 x=3.0000 - 0.0000i 2.0000 + 0.0000i 1.0000 - 0.0000i 2.0000 + 0.0000i 3.0000 - 0.0000 -0.0000 + 0.0000i 0 - 0.0000i -0.0000 + 0.0000i 0.0000 - 0.0000i -0.0000 + 0.0000i 由jω X(e)的频谱表达式可知，有限长时间序列x(n)的长度M=5，现分别取频域抽样点 数为N=3,5,10，由图5-4显示的结果可以验证： ①当N=5和N=10时，N≥M，能够不失真地恢复出原信号x(n)； ②当N=3时，N＜M，时间序列有泄漏，形成了混叠，不能无失真地恢复出原信号x(n)。混叠的原因是上一周期的后2点与本周期的前两点发生重叠，如下所示： 3 2 1 2 3 3 2 1 2 3 例5-5已知一个频率范围在[-62.8,62.8]rad/s间的频谱 sin0.275Ω X(jΩ)= sin0.025Ω ，用IFFT计 算并求出时间序列x(n)，用图形显示时间序列。 分析：本题给出了模拟频率Ω，其中Ωm=62.8，需将其归一化为数字频率。根据奈奎斯特定理可知，（1/Ts）=Fs≥(2Ωm/2π),可以推导出Ts≤（π/Ωm）,取Ts=0.05s，即采样频率Fs为20Hz或40π。 程序清单如下： wm=62.8;Ts=pi/wm; N0=[8,20]; for r=1:2 N=N0(r); D=2*pi/(Ts*N); k=[0:N-1]+eps; omg=k*D; X=sin(0.275*omg)./sin(0.025*omg); n=0:N-1; x=abs(ifft(X,N)); subplot(1,2,r);stem(n*Ts,abs(x)); box end 程序运行结果如图5-5所示： 图5-5 由N=20的结果可知，时间序列x(n)是一个矩形窗。根据DFT的循环移位性质可知，非零数据存在于n=-5:5的区域，有限长序列的长度为11。而N=8小于有限长序列的长度，其结果发生了混叠，不能无失真地恢复出原信号x(n)。 3、从频谱恢复连续时间信号 实际应用中，离散信号往往来源于对连续信号的采样，因此，这里要讨论从频谱如何计算连续时间信号。 从本质上讲，计算机处理的都是离散信号。当使用计算机处理连续信号时，实际上是用采样周期极小的序列信号来近似为连续信号。因此在处理时，原来对于离散序列处理的理论依然有效。 （1）选择一个符合奈奎斯特定理的很小的采样周期T ，将主要的模拟频谱限制在奈奎斯特频率范围内，a X (Ω)=0 Ωπ/T ≥当； （2）在[-π/T,π/T ]的频率区间取N 个频率点Ωk ，求出对应的数字频谱： a k k X (Ω) X(Ω)= T （3）对k X(Ω)作IDFT ，求x a (t)。假定没有发生时间混叠，则 a k a t=nT X (Ω) x (t)x(n)=IDFT[ ]T ≈ （4）作图。用plot 自动进行插值，获得连续信号。 例5-6 已知图5-6所示的理想低通滤波器的模拟频率Ωc =3，在|Ω|≤Ωc 范围内幅度为1，|Ω|＞Ωc 时幅度为0。要求计算连续脉冲响应x a (t)。 Ω -Ωc Ωc 图5-6 分析：由奈奎斯特定理可知采样频率Ωs ≥2Ωc ,即采样周期T s ≤πΩc 恢复原信号时不会发生混叠。选的再小一些可以增加样点数，因此可以选T s =0.1π/Ωc =0.1047s 。同时，为使时间信号尽量接近连续信号，需提高N 点的个数。可以由模拟频率的分辨率公式推导： D=2π/NT s ≤0.1Ωc ，使频率分辨率小于有效带宽的1/10，得到：N ≥20π/Ωc T s 程序清单如下： wc=3;Tmax=0.1*pi/wc; Ts=input('(Ts N=input('(N>Nmin)N='); D=2*pi/(Ts*N); M=floor(wc/D); Xa=[ones(1,M+1),zeros(1,N-2*M-1),ones(1,M)]; n=-(N-1)/2:(N-1)/2; xa=abs(fftshift(ifft(Xa/Ts))); plot(n*Ts,xa); 程序执行过程中，在MATLAB 命令窗口将给出提示：输入Ts 和N 的值，再给出绘图结果。图5-7是分别输入Ts=0.1s ，N=300和Ts=0.1s ，N=1000两组数据的运行结果。 -20 -10 10 20 00.20.40.60.8 1 图5-7 三、实验内容： 1、阅读并输入实验原理中介绍的例题程序，观察输出的数据和图形，结合基本原理理解每一条语句的含义。 2、已知一个连续时间信号f(t)=sinc(t)，取最高有限带宽频率f m =1Hz 。 （1）分别显示原连续信号波形和F s =f m 、F s =2f m 、F s =3f m 三种情况下抽样信号的波形； dt=0.1; fm=1; Tm=1/fm; t=-6:dt:6; f=sinc(t); subplot(4,1,1); plot(t,f); axis([min(t),max(t),1.1*min(f),1.1*max(f)]); title('?-á?D?D?o?oí3é?ùD?o?'); for i=1:3; fs=i*fm;Ts=1/fs; n=-6:Ts:6; f=sinc(n); subplot(4,1,i+1);stem(n,f,'filled'); axis([min(n),max(n),1.1*min(f),1.1*max(f)]); end -50 50 00.20.40.60.8 1 dt=0.1;f0=1;T0=1/f0;fm=1;Tm=1/fm; t=-2:dt:2;N=length(t); f=sinc(f0*t); wm=2*pi*fm;k=0:N-1;w1=k*wm/N; F1=f*exp(-j*t'*w1)*dt;subplot(4,1,1);plot(w1/(2*pi),abs(F1)); axis([0,max(4*fm),1.1*min(abs(F1)),1.1*max(abs(F1))]); for i=1:3; if i<=2 c=0;else c=1;end fs=(i+c)*fm;Ts=1/fs; n=-2:Ts:2;N=length(n); f=sinc(f0*n); wm=2*pi*fs;k=0:N-1; w=k*wm/N;F=f*exp(-j*n'*w)*Ts; subplot(4,1,i+1);plot(w/(2*pi),abs(F)); axis([0,max(4*fm),1.1*min(abs(F)),1.1*max(abs(F))]); end dt=0.01;f0=1;T0=1/f0;fm=1;Tm=1/fm; t=-7*T0:dt:7*T0;x=sinc(t); subplot(4,1,1);plot(t,x);axis([min(t),max(t),1.1*min(x),1.1*max(x)]); title('ó?ê±óò?í?y???¨3é?ùD?o?'); for i=1:3; fs=i*fm;Ts=1/fs; n=0:(3*T0)/Ts; t1=-7*T0:Ts:7*T0; x1=sinc(f0*n/fs); T_N=ones(length(n),1)*t1-n'*Ts*ones(1,length(t1)); xa=x1*sinc(fs*pi*T_N); subplot(4,1,i+1);plot(t1,xa); axis([min(t1),max(t1),1.1*min(xa),1.1*max(xa)]); end 3、已知一个时间序列的频谱为：j ω -j ωn -j ω-j2ω-j3ω-j4ωn=-X(e )= x(n)e =2+4e +6e +4e +2e ∞ ∞ ∑ 分别取频域抽样点数N 为3、5和10，用IFFT 计算并求出其时间序列x(n)，绘图显示个时间序列。由此讨论由频域抽样不失真地恢复原时域信号的条件。 Ts=1;N0=[3,5,10]; for r=1:3; N=N0(r); D=2*pi/(Ts*N); kn=floor(-(N-1)/2:-1/2); kp=floor(0:(N-1)/2); w=[kp,kn]*D; X=2+4*exp(-j*w)+6*exp(-j*2*w)+4*exp(-j*3*w)+2*exp(-j*4*w); n=0:N-1; x=ifft(X,N) subplot(1,3,r);stem(n*Ts,abs(x)); box end c 它范围幅度为0.计算其连续信号x a(t)，并绘图显示信号曲线。 wc=3.14;Tmax=0.1*pi/wc; Ts=input('(Ts Nmin=20*pi/wc/Ts; N=input('(N>Nmin)N='); D=2*pi/(Ts*N); M=floor(wc/D); Xa=[zeros(1,M),1,zeros(1,N-2*M-1),1,zeros(1,M-1)]; n=-(N-1)/2:(N-1)/2; xa=abs(fftshift(ifft(Xa/Ts))); plot(n*Ts,xa); 输入Ts=0.1 N=100 1、认真阅读实验原理部分，明确实验目的，读懂例题程序，了解实验方法。 2、根据实验任务预先编写实验程序。 3、预习思考题： ①什么是内插公式？在MATLAB中内插公式可用什么函数来编写？ ②从频域抽样序列不失真地恢复离散时域信号的条件是什么？ 五、实验报告： 1、列写调试通过的实验程序，打印实验程序产生的曲线图形。 2、给出预习思考题答案。 3、思考题： ①试归纳用IFFT数值计算方法从频谱恢复离散时间序列的方法和步骤。 ②从频谱恢复连续时间信号与恢复离散时间序列有何不同？ 实验二 IIR数字滤波器设计及软件实现 1．实验目的 （1）熟悉用双线性变换法设计IIR数字滤波器的原理与方法； （2）学会调用MATLAB信号处理工具箱中滤波器设计函数（或滤波器设计分析工具fdatool）设计各种IIR数字滤波器，学会根据滤波需求确定滤波器指标参数。 （3）掌握IIR数字滤波器的MATLAB实现方法。 （3）通过观察滤波器输入输出信号的时域波形及其频谱，建立数字滤波的概念。 2．实验原理 设计IIR数字滤波器一般采用间接法（脉冲响应不变法和双线性变换法），应用最广泛的是双线性变换法。基本设计过程是：①先将给定的数字滤波器的指标转换成过渡模拟滤波器的指标；②设计过渡模拟滤波器；③将过渡模拟滤波器系统函数转换成数字滤波器的系统函数。MATLAB信号处理工具箱中的各种IIR数字滤波器设计函数都是采用双线性变换法。第六章介绍的滤波器设计函数butter、cheby1 、cheby2 和ellip可以分别被调用来直接设计巴特沃斯、切比雪夫1、切比雪夫2和椭圆模拟和数字滤波器。本实验要求读者调用如上函数直接设计IIR数字滤波器。 本实验的数字滤波器的MATLAB实现是指调用MATLAB信号处理工具箱函数filter对给定的输入信号x(n)进行滤波，得到滤波后的输出信号y(n）。 3. 实验内容及步骤 （1）调用信号产生函数mstg产生由三路抑制载波调幅信号相加构成的复合信号st，该函数还会自动绘图显示st的时域波形和幅频特性曲线，如图1所示。由图可见，三路信号时域混叠无法在时域分离。但频域是分离的，所以可以通过滤波的方法在频域分离，这就是本实验的目的。 图1 三路调幅信号st的时域波形和幅频特性曲线 （2）要求将st中三路调幅信号分离，通过观察st的幅频特性曲线，分别确定可以分离st中三路抑制载波单频调幅信号的三个滤波器（低通滤波器、带通滤波器、高通滤波器）的通带截止频率和阻带截止频率。要求滤波器的通带最大衰减为0.1dB,阻带最小衰减为 实验一MATLAB语言的基本使用方法 实验类别：基础性实验 实验目的： （1）了解MATLAB程序设计语言的基本方法，熟悉MATLAB软件运行环境。 （2）掌握创建、保存、打开m文件的方法，掌握设置文件路径的方法。 （3）掌握变量、函数等有关概念，具备初步的将一般数学问题转化为对应计算机模型并进行处理的能力。 （4）掌握二维平面图形的绘制方法，能够使用这些方法进行常用的数据可视化处理。 实验内容和步骤： 1、打开MATLAB，熟悉MATLAB环境。 2、在命令窗口中分别产生3*3全零矩阵，单位矩阵，全1矩阵。 3、学习m文件的建立、保存、打开、运行方法。 4、设有一模拟信号f(t)=1.5sin60πt,取?t=0.001,n=0,1,2,…,N-1进行抽样，得到 序列f(n),编写一个m文件sy1_1.m,分别用stem,plot,subplot等命令绘制32 点序列f(n)(N=32)的图形，给图形加入标注，图注，图例。 5、学习如何利用MATLAB帮助信息。 实验结果及分析： 1）全零矩阵 >> A=zeros(3,3) A = 0 0 0 0 0 0 0 0 0 2）单位矩阵 >> B=eye(3) B = 1 0 0 0 1 0 0 0 1 3）全1矩阵 >> C=ones(3) C = 1 1 1 1 1 1 1 1 1 4)sy1_1.m N=32; n=0:N-1; dt=0.001; t=n*dt; y=1.5*sin(60*pi*t); subplot(2,1,1), plot(t,y); xlabel('t'); ylabel('y=1.5*sin(60*pi*t)'); legend('正弦函数'); title('二维图形'); subplot(2,1,2), stem(t,y) xlabel('t'); ylabel('y=1.5*sin(60*pi*t)'); legend('序列函数'); title('条状图形'); 00.0050.010.0150.020.0250.030.035 t y = 1 . 5 * s i n ( 6 * p i * t ) 二维图形 00.0050.010.0150.020.0250.030.035 t y = 1 . 5 * s i n ( 6 * p i * t ) 条状图形 《数字信号处理》实验指导书 光电工程学院二○○九年十月 实验一离散时间信号分析 一、实验目的 1．掌握各种常用的序列，理解其数学表达式和波形表示。 2．掌握在计算机中生成及绘制数字信号波形的方法。 3．掌握序列的相加、相乘、移位、反转等基本运算及计算机实现与作用。 4．掌握线性卷积软件实现的方法。 5．掌握计算机的使用方法和常用系统软件及应用软件的使用。 6．通过编程，上机调试程序，进一步增强使用计算机解决问题的能力。 二、实验原理 1．序列的基本概念 离散时间信号在数学上可用时间序列来表示，其中代表序列的第n个数字，n代表时间的序列，n的取值范围为的整数，n取其它值没有意义。离散时间信号可以是由模拟信号通过采样得到，例如对模拟信号进行等间隔采样，采样间隔为T，得到一个有序的数字序列就是离散时间信号，简称序列。 2．常用序列 常用序列有：单位脉冲序列（单位抽样）、单位阶跃序列、矩形序列、实指数序列、复指数序列、正弦型序列等。 3．序列的基本运算 序列的运算包括移位、反转、和、积、标乘、累加、差分运算等。 4．序列的卷积运算 上式的运算关系称为卷积运算，式中代表两个序列卷积运算。两个序列的卷积是一个序列与另一个序列反褶后逐次移位乘积之和，故称为离散卷积，也称两序列的线性卷积。其计算的过程包括以下4个步骤。 （1）反褶：先将和的变量换成，变成和，再将以纵轴为对称轴反褶成。 （2）移位：将移位，得。当为正数时，右移位；当为负数时，左 移位。 （3）相乘：将和的对应点值相乘。 （4）求和：将以上所有对应点的乘积累加起来，即得。 三、主要实验仪器及材料 微型计算机、Matlab软件6.5或更高版本。 四、实验内容 1．知识准备 认真复习以上基础理论，理解本实验所用到的实验原理。 2．离散时间信号（序列）的产生 利用MATLAB或C语言编程产生和绘制下列有限长序列： （1）单位脉冲序列 （2）单位阶跃序列 （3）矩形序列 （4）正弦型序列 （5）任意序列 3．序列的运算 利用MATLAB编程完成上述两序列的移位、反转、加法、乘法等运算，并绘制运算后序列的波形。 4．卷积运算 利用MATLAB编制一个计算两个序列线性卷积的通用程序，计算上述两序列，并绘制卷积后序列的波形。 5．上机调试并打印或记录实验结果。 6．完成实验报告。 五、实验报告要求 1. 简述实验原理及目的。 2. 给出上述序列的实验结果。 3. 列出计算卷积的公式，画出程序框图，并列出实验程序清单 （可略）（包括必要的程序说明）。 4. 记录调试运行情况及所遇问题的解决方法。 5. 给出实验结果，并对结果做出分析。 6. 简要回答思考题。 1 如何产生方波信号序列和锯齿波信号序列? 2 实验中所产生的正弦序列的频率是多少?是否是周期序列? 2.1 clc close all; n=0:15; p=8;q=2; x=exp(-(n-p.^2/q; figure(1; subplot(3,1,1; stem(n,x; title('exp(-(n-p^2/q,p=8,q=2'; xk1=fft(x,16; q=4; x=exp(-(n-p.^2/q; subplot(3,1,2; xk2=fft(x,16; stem(n,x; title('exp(-(n-p^2/q,p=8,q=4'; q=8; x=exp(-(n-p.^2/q; xk3=fft(x,16; subplot(3,1,3; stem(n,x; title('exp(-(n-p^2/q,p=8,q=8';%时域特性figure(2; subplot(3,1,1; stem(n,abs(xk1; title('exp(-(n-p^2/q,p=8,q=2'; subplot(3,1,2; stem(n,abs(xk2; title('exp(-(n-p^2/q,p=8,q=4'; subplot(3,1,3; stem(n,abs(xk3; title('exp(-(n-p^2/q,p=8,q=8';%频域特性%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%% p=8;q=8; figure(3; subplot(3,1,1; stem(n,x; title('exp(-(n-p^2/q,p=8,q=8'; xk1=fft(x,16; p=13; x=exp(-(n-p.^2/q; subplot(3,1,2; xk2=fft(x,16; stem(n,x; title('exp(-(n-p^2/q,p=13,q=8'; p=14; x=exp(-(n-p.^2/q; xk3=fft(x,16; subplot(3,1,3; stem(n,x; title('exp(-(n-p^2/q,p=14,q=8';%时域特性figure(4; subplot(3,1,1; stem(n,abs(xk1; title('exp(-(n-p^2/q,p=8,q=8'; subplot(3,1,2; stem(n,abs(xk2; title('exp(-(n-p^2/q,p=13,q=8'; subplot(3,1,3; 实验6 数字滤波器的网络结构 一、实验目的： 1、加深对数字滤波器分类与结构的了解。 2、明确数字滤波器的基本结构及其相互间的转换方法。 3、掌握用MA TLAB 语言进行数字滤波器结构间相互转换的子函数及程序编写方法。 二、实验原理： 1、数字滤波器的分类 离散LSI 系统对信号的响应过程实际上就是对信号进行滤波的过程。因此，离散LSI 系统又称为数字滤波器。 数字滤波器从滤波功能上可以分为低通、高通、带通、带阻以及全通滤波器；根据单位脉冲响应的特性，又可以分为有限长单位脉冲响应滤波器（FIR ）和无限长单位脉冲响应滤波器（IIR ）。 一个离散LSI 系统可以用系统函数来表示： M -m -1-2-m m m=0 012m N -1-2-k -k 12k k k=1 b z b +b z +b z ++b z Y(z)b(z)H(z)=== =X(z)a(z) 1+a z +a z ++a z 1+a z ∑∑ 也可以用差分方程来表示： N M k m k=1 m=0 y(n)+a y(n-k)=b x(n-m)∑∑ 以上两个公式中，当a k 至少有一个不为0时，则在有限Z 平面上存在极点，表达的是以一个IIR 数字滤波器；当a k 全都为0时，系统不存在极点，表达的是一个FIR 数字滤波器。FIR 数字滤波器可以看成是IIR 数字滤波器的a k 全都为0时的一个特例。 IIR 数字滤波器的基本结构分为直接Ⅰ型、直接Ⅱ型、直接Ⅲ型、级联型和并联型。 FIR 数字滤波器的基本结构分为横截型（又称直接型或卷积型）、级联型、线性相位型及频率采样型等。本实验对线性相位型及频率采样型不做讨论，见实验10、12。 另外，滤波器的一种新型结构——格型结构也逐步投入应用，有全零点FIR 系统格型结构、全极点IIR 系统格型结构以及全零极点IIR 系统格型结构。 2、IIR 数字滤波器的基本结构与实现 （1）直接型与级联型、并联型的转换 例6-1 已知一个系统的传递函数为 -1-2-3 -1-2-3 8-4z +11z -2z H(z)=1-1.25z +0.75z -0.125z 将其从直接型（其信号流图如图6-1所示）转换为级联型和并联型。 一、实验目的 1. 通过本次实验回忆并熟悉MATLAB这个软件。 2. 通过本次实验学会如何利用MATLAB进行序列的简单运算。 3. 通过本次实验深刻理解理论课上的数字信号处理的一个常见方法——对时刻n的样本附近的一些样本求平均，产生所需的输出信号。 3. 通过振幅调制信号的产生来理解载波信号与调制信号之间的关系。 二、实验内容 1. 编写程序在MATLAB中实现从被加性噪声污染的信号中移除噪声的算法，本次试验采用三点滑动平均算法，可直接输入程序P1.5。 2. 通过运行程序得出的结果回答习题Q1.31-Q1.33的问题，加深对算法思想的理解。 3. 编写程序在MATLAB中实现振幅调制信号产生的算法，可直接输入程序P1.6。 4. 通过运行程序得出的结果回答习题Q1.34-Q1.35的问题，加深对算法思想的理解。 三、主要算法与程序 1. 三点滑动平均算法的核心程序： %程序P1.5 %通过平均的信号平滑 clf; R=51; d=0.8*(rand(R,1)-0.5);%产生随噪声 m=0:R-1; s=2*m.*(0.9.^m);%产生为污染的信号 x=s+d';%产生被噪音污染的信号 subplot(2,1,1); plot(m,d','r-',m,s,'g--',m,x,'b-.'); xlabel('时间序号n');ylabel('振幅'); legend('d[n]','s[n]','x[n]'); x1=[0 0 x];x2=[0 x 0];x3=[x 0 0]; y=(x1+x2+x3)/3; subplot(2,1,2); plot(m,y(2:R+1),'r-',m,s,'g--'); legend('y[n]','s[n]'); xlabel('时间序号n');ylabel('振幅'); 2. 振幅调制信号的产生核心程序：（由于要几个结果，因此利用subplot函数画图） %程序P1.6 %振幅调制信号的产生 n=0:100; m=0.1;fH=0.1;fL=0.01; m1=0.3;fH1=0.3;fL1=0.03; xH=sin(2*pi*fH*n); xL=sin(2*pi*fL*n); y=(1+m*xL).*xH; xH1=sin(2*pi*fH1*n); xL1=sin(2*pi*fL1*n); y1=(1+m1*xL).*xH; y2=(1+m*xL).*xH1; y3=(1+m*xL1).*xH; subplot(2,2,1); stem(n,y); grid; xlabel('时间序号n');ylabel('振幅');title('m=0.1;fH=0.1;fL=0.01;'); subplot(2,2,2); stem(n,y1); grid; xlabel('时间序号n');ylabel('振幅');title('m=0.3;fH=0.1;fL=0.01;'); subplot(2,2,3); stem(n,y2); grid; xlabel('时间序号n');ylabel('振幅');title('m=0.3;fH=0.3;fL=0.01;'); subplot(2,2,4); stem(n,y3); grid; 武汉工程大学 数字信号处理实验报告 姓名：周权 学号：1204140228 班级：通信工程02 一、实验设备 计算机，MATLAB语言环境。 二、实验基础理论 1.序列的相关概念 2.常见序列 3.序列的基本运算 4.离散傅里叶变换的相关概念 5.Z变换的相关概念 三、实验内容与步骤 1.离散时间信号（序列）的产生 利用MATLAB语言编程产生和绘制单位样值信号、单位阶跃序列、指数序列、正弦序列及随机离散信号的波形表示。 四实验目的 认识常用的各种信号，理解其数字表达式和波形表示，掌握在计算机中生成及绘制数字信号波形的方法，掌握序列的简单运算及计算机实现与作用，理解离散时间傅里叶变换，Z变换及它们的性质和信号的频域分 实验一离散时间信号（序列）的产生 代码一 单位样值 x=2; y=1; stem(x,y); title('单位样值 ') 单位阶跃序列 n0=0; n1=-10; n2=10; n=[n1:n2]; x=[(n-n0)>=0]; stem(n,x); xlabel('n'); ylabel('x{n}'); title('单位阶跃序列'); 实指数序列 n=[0:10]; x=(0.5).^n; stem(n,x); xlabel('n'); ylabel('x{n}'); title('实指数序列'); 正弦序列 n=[-100:100]; x=2*sin(0.05*pi*n); stem(n,x); xlabel('n'); ylabel('x{n}'); title('正弦序列'); 随机序列 n=[1:10]; x=rand(1,10); subplot(221); stem(n,x); xlabel('n'); ylabel('x{n}'); title('随机序列'); 实验5 抽样定理 一、实验目的： 1、了解用MA TLAB 语言进行时域、频域抽样及信号重建的方法。 2、进一步加深对时域、频域抽样定理的基本原理的理解。 3、观察信号抽样与恢复的图形，掌握采样频率的确定方法和插公式的编程方法。 二、实验原理： 1、时域抽样与信号的重建 （1）对连续信号进行采样 例5-1 已知一个连续时间信号sin sin(),1Hz 3 ππ=0001f(t)=(2f t)+6f t f ，取最高有限带宽频率f m =5f 0，分别显示原连续时间信号波形和F s >2f m 、F s =2f m 、F s <2f m 三情况下抽样信号的波形。 程序清单如下： %分别取Fs=fm ，Fs=2fm ，Fs=3fm 来研究问题 dt=0.1; f0=1; T0=1/f0; m=5*f0; Tm=1/fm; t=-2:dt:2; f=sin(2*pi*f0*t)+1/3*sin(6*pi*f0*t); subplot(4,1,1); plot(t,f); axis([min(t),max(t),1.1*min(f),1.1*max(f)]); title('原连续信号和抽样信号'); for i=1:3; fs=i*fm;Ts=1/fs; n=-2:Ts:2; f=sin(2*pi*f0*n)+1/3*sin(6*pi*f0*n); subplot(4,1,i+1);stem(n,f,'filled'); axis([min(n),max(n),1.1*min(f),1.1*max(f)]); end 程序运行结果如图5-1所示： 原连续信号和抽样信号 图5-1 （2）连续信号和抽样信号的频谱 由理论分析可知，信号的频谱图可以很直观地反映出抽样信号能否恢复原模拟信号。因此，我们对上述三种情况下的时域信号求幅度谱，来进一步分析和验证时域抽样定理。 例5-2编程求解例5-1中连续信号及其三种抽样频率（F s>2f m、F s=2f m、F s<2f m）下的抽样信号的幅度谱。 程序清单如下： dt=0.1;f0=1;T0=1/f0;fm=5*f0;Tm=1/fm; t=-2:dt:2;N=length(t); f=sin(2*pi*f0*t)+1/3*sin(6*pi*f0*t); wm=2*pi*fm;k=0:N-1;w1=k*wm/N; F1=f*exp(-j*t'*w1)*dt;subplot(4,1,1);plot(w1/(2*pi),abs(F1)); axis([0,max(4*fm),1.1*min(abs(F1)),1.1*max(abs(F1))]); for i=1:3; if i<=2 c=0;else c=1;end fs=(i+c)*fm;Ts=1/fs; n=-2:Ts:2;N=length(n); f=sin(2*pi*f0*n)+1/3*sin(6*pi*f0*n); wm=2*pi*fs;k=0:N-1; w=k*wm/N;F=f*exp(-j*n'*w)*Ts; subplot(4,1,i+1);plot(w/(2*pi),abs(F)); axis([0,max(4*fm),1.1*min(abs(F)),1.1*max(abs(F))]); end 程序运行结果如图5-2所示。 由图可见，当满足F s≥2f m条件时，抽样信号的频谱没有混叠现象；当不满足F s≥2f m 条件时，抽样信号的频谱发生了混叠，即图5-2的第二行F s＜2f m的频谱图，，在f m=5f0的围，频谱出现了镜像对称的部分。 实验一 MATLAB 仿真软件的基本操作命令和使用方法 实验容 1、帮助命令 使用 help 命令，查找 sqrt （开方）函数的使用方法； 2、MATLAB 命令窗口 （1）在MATLAB 命令窗口直接输入命令行计算3 1)5.0sin(21+=πy 的值； （2）求多项式 p(x) = x3 + 2x+ 4的根； 3、矩阵运算 （1）矩阵的乘法 已知 A=[1 2;3 4]， B=[5 5;7 8]，求 A^2*B （2）矩阵的行列式 已知A=[1 2 3;4 5 6;7 8 9]，求A （3）矩阵的转置及共轭转置 已知A=[1 2 3;4 5 6;7 8 9]，求A' 已知B=[5+i,2-i,1;6*i,4,9-i], 求B.' , B' （4）特征值、特征向量、特征多项式 已知A=[1.2 3 5 0.9;5 1.7 5 6;3 9 0 1;1 2 3 4] ，求矩阵A的特征值、特征向量、特征多项式； （5）使用冒号选出指定元素 已知：A=[1 2 3;4 5 6;7 8 9]；求A 中第3 列前2 个元素；A 中所有列第2，3 行的元素； 4、Matlab 基本编程方法 （1）编写命令文件：计算1+2+…+n<2000 时的最大n 值； （2）编写函数文件：分别用for 和while 循环结构编写程序，求 2 的0 到15 次幂的和。 5、MATLAB基本绘图命令 （1）绘制余弦曲线 y=cos(t)，t∈[0，2π] （2）在同一坐标系中绘制余弦曲线 y=cos(t-0.25)和正弦曲线 y=sin(t-0.5)， t∈[0，2π] （3）绘制[0，4π]区间上的 x1=10sint 曲线，并要求： （a）线形为点划线、颜色为红色、数据点标记为加号； （b）坐标轴控制：显示围、刻度线、比例、网络线 （c）标注控制：坐标轴名称、标题、相应文本； >> clear; 实验1认识实验-MATLAB语言上机操作实践 一、实验目的 ㈠了解MATLAB语言的主要特点、作用。 ㈡学会MATLAB主界面简单的操作使用方法。 ㈢学习简单的数组赋值、运算、绘图、流程控制编程。 二、实验原理 ㈠简单的数组赋值方法 MATLAB中的变量和常量都可以是数组（或矩阵）,且每个元素都可以是复数。 在MATLAB指令窗口输入数组A=[1 2 3；4 5 6；7 8 9]，观察输出结果。然后，键入：A（4，2）= 11 键入：A (5，：) = [-13 -14 -15] 键入：A（4，3）= abs (A(5，1)) 键入：A ([2，5]，：) = [ ] 键入：A/2 键入：A (4，：) = [sqrt(3) (4+5)/6*2 –7] 观察以上各输出结果。将A式中分号改为空格或逗号，情况又如何？请在每式的后面标注其含义。 2．在MATLAB指令窗口输入B=[1+2i，3+4i；5+6i ，7+8i], 观察输出结果。 键入：C=[1，3；5，7]+[2，4；6，8]*i，观察输出结果。 如果C式中i前的*号省略，结果如何？ 键入：D = sqrt (2+3i) 键入：D*D 键入：E = C’, F = conj(C), G = conj(C)’ 观察以上各输出结果, 请在每式的后面标注其含义。 3．在MATLAB指令窗口输入H1=ones(3,2)，H2=zeros(2,3)，H3=eye(4)，观察输出结果。 ㈡、数组的基本运算 1．输入A=[1 3 5]，B= [2 4 6]，求C=A+B，D=A-2，E=B-A 2．求F1=A*3，F2=A.*B，F3=A./B，F4=A.\B, F5=B.\A, F6=B.^A, F7=2./B, F8=B.\2 *3．求B'，Z1=A*B’，Z2=B’*A 观察以上各输出结果,比较各种运算的区别，理解其含义。 ㈢、常用函数及相应的信号波形显示 例1：显示曲线f(t)=2sin(2πt)，（t>0） ⅰ点击空白文档图标（New M-file），打开文本编辑器。 ⅱ键入：t=0:0.01:3; (1) f=2*sin(2*pi*t); (2) plot(t,f); title(‘f（t）-t曲线’)； xlabel（‘t’），ylabel（‘f（t）’）； 数字信号处理实验答案 HEN system office room 【HEN16H-HENS2AHENS8Q8-HENH1688】 实验一熟悉Matlab环境 一、实验目的 1.熟悉MATLAB的主要操作命令。 2.学会简单的矩阵输入和数据读写。 3.掌握简单的绘图命令。 4.用MATLAB编程并学会创建函数。 5.观察离散系统的频率响应。 二、实验内容 认真阅读本章附录，在MATLAB环境下重新做一遍附录中的例子，体会各条命令的含义。在熟悉了MATLAB基本命令的基础上，完成以下实验。 上机实验内容： （1）数组的加、减、乘、除和乘方运算。输入A=[1 2 3 4]，B=[3 4 5 6]，求C=A+B，D=A-B，E=A.*B，F=A./B，G=A.^B并用stem语句画出A、B、C、D、E、F、G。 clear all; a=[1 2 3 4]; b=[3 4 5 6]; c=a+b; d=a-b; e=a.*b; f=a./b; g=a.^b; n=1:4; subplot(4,2,1);stem(n,a); xlabel('n');xlim([0 5]);ylabel('A'); subplot(4,2,2);stem(n,b); xlabel('n');xlim([0 5]);ylabel('B'); subplot(4,2,3);stem(n,c); xlabel('n');xlim([0 5]);ylabel('C'); subplot(4,2,4);stem(n,d); xlabel('n');xlim([0 5]);ylabel('D'); subplot(4,2,5);stem(n,e); xlabel('n');xlim([0 5]);ylabel('E'); subplot(4,2,6);stem(n,f); xlabel('n');xlim([0 5]);ylabel('F'); subplot(4,2,7);stem(n,g); xlabel('n');xlim([0 5]);ylabel('G'); （2）用MATLAB实现下列序列： a) x(n)= 0≤n≤15 b) x(n)=e+3j)n 0≤n≤15 c) x(n)=3cosπn+π)+2sinπn+π) 0≤n≤15 d) 将c)中的x(n)扩展为以16为周期的函数x(n)=x(n+16),绘出四个周期。 数字信号处理实验实验一离散时间信号与系统及MA TLAB实现 1.单位冲激信号： n = -5:5; x = (n==0); subplot(122); stem(n, x); 2.单位阶跃信号： x=zeros(1,11); n0=0; n1=-5; n2=5; n = n1:n2; x(:,n+6) = ((n-n0)>=0); stem(n,x); 3.正弦序列： n = 0:1/3200:1/100; x=3*sin(200*pi*n+1.2); stem(n,x); 4.指数序列 n = 0:1/2:10; x1= 3*(0.7.^n); x2=3*exp((0.7+j*314)*n); subplot(221); stem(n,x1); subplot(222); stem(n,x2); 5.信号延迟 n=0:20; Y1=sin(100*n); Y2=sin(100*(n-3)); subplot(221); stem(n,Y1); subplot(222); stem(n,Y2); 6.信号相加 X1=[2 0.5 0.9 1 0 0 0 0]; X2=[0 0.1 0.2 0.3 0.4 0.5 0.6 0.7]; X=X1+X2; stem(X); 7.信号翻转 X1=[2 0.5 0.9 1]; n=1:4; X2=X1(5-n); subplot(221); stem(n,X1); subplot(222); stem(n,X2); 8.用MATLAB计算序列{-2 0 1 –1 3}和序列{1 2 0 -1}的离散卷积。a=[-2 0 1 -1 3]; b=[1 2 0 -1]; c=conv(a,b); M=length(c)-1; n=0:1:M; stem(n,c); xlabel('n'); ylabel('幅度'); 9．用MA TLAB计算差分方程 当输入序列为时的输出结果。 N=41; a=[0.8 -0.44 0.36 0.22]; b=[1 0.7 -0.45 -0.6]; x=[1 zeros(1,N-1)]; k=0:1:N-1; y=filter(a,b,x); stem(k,y) xlabel('n'); ylabel('幅度') 10.冲激响应impz N=64; a=[0.8 -0.44 0.36 0.22]; 实验一 基于Matlab 的数字信号处理基本操作 一、 实验目的：学会运用MA TLAB 表示的常用离散时间信号；学会运用MA TLAB 实现离 散时间信号的基本运算。 二、 实验仪器：电脑一台，MATLAB6.5或更高级版本软件一套。 三、 实验内容： （一） 离散时间信号在MATLAB 中的表示 离散时间信号是指在离散时刻才有定义的信号，简称离散信号，或者序列。离散序列通常用)(n x 来表示，自变量必须是整数。 离散时间信号的波形绘制在MATLAB 中一般用stem 函数。stem 函数的基本用法和plot 函数一样，它绘制的波形图的每个样本点上有一个小圆圈，默认是空心的。如果要实心，需使用参数“fill ”、“filled ”，或者参数“.”。由于MATLAB 中矩阵元素的个数有限，所以MA TLAB 只能表示一定时间范围内有限长度的序列；而对于无限序列，也只能在一定时间范围内表示出来。类似于连续时间信号，离散时间信号也有一些典型的离散时间信号。 1. 单位取样序列 单位取样序列)(n δ，也称为单位冲激序列，定义为 ) 0() 0(0 1)(≠=?? ?=n n n δ 要注意，单位冲激序列不是单位冲激函数的简单离散抽样，它在n =0处是取确定的值1。在MATLAB 中，冲激序列可以通过编写以下的impDT .m 文件来实现，即 function y=impDT(n) y=(n==0); %当参数为0时冲激为1,否则为0 调用该函数时n 必须为整数或整数向量。 【实例1-1】 利用MATLAB 的impDT 函数绘出单位冲激序列的波形图。 解：MATLAB 源程序为 >>n=-3:3; >>x=impDT(n); >>stem(n,x,'fill'),xlabel('n'),grid on >>title('单位冲激序列') >>axis([-3 3 -0.1 1.1]) 程序运行结果如图1-1所示。 图1-1 单位冲激序列 第二章上机作业 1、ljdt(A,B)函数定义 function ljdt(A,B) p=roots(A); q=roots(B); p=p'; q=q'; x=max(abs([p q 1])); x=x+0.1; y=x; clf hold on axis([-x x -y y]) w=0:pi/300:2*pi; t=exp(i*w); plot(t) axis('square') plot([-x x],[0 0]) plot([0 0],[-y y]) text(0.1,x,'jIm[z]') text(y,1/10,'Re[z]') plot(real(p),imag(p),'x') plot(ral(q),imag(q),'o') title('pole-zero diagram for discrete system') hold off 例2.26 a=[3 -1 0 0 0 1]; b=[1 1]; ljdt(a,b) p=roots(a) q=roots(b) pa=abs(p) 程序运行结果如下： P= 0.7255+0.4633i 0.7255+0.4633i -0.1861+0.7541i -0.1861-0.7541i -0.7455 q= -1 pa= 0.8608 0.8608 0.7768 0.7768 0.7455 例2.27 b=[0 1 2 1];a=[1 -0.5 -0.005 0.3]; subplot 311 zplane(b,a);xlabel('实部');ylabel('虚部'); num=[0 1 2 1];den=[1 -0.5 -0.005 0.3]; h=impz(num,den); subplot 312 本科生实验报告 实验课程数字信号处理基础 学院名称地球物理学院 专业名称地球物理学 学生姓名 学生学号 指导教师王山山 实验地点5417 实验成绩 二〇一四年十一月二〇一四年十二月 填写说明 1、适用于本科生所有的实验报告（印制实验报告册除外）； 2、专业填写为专业全称，有专业方向的用小括号标明； 3、格式要求： ①用A4纸双面打印（封面双面打印）或在A4大小纸上用蓝黑色水笔书写。 ②打印排版：正文用宋体小四号，1.5倍行距，页边距采取默认形式（上下2.54cm， 左右2.54cm，页眉1.5cm，页脚1.75cm）。字符间距为默认值（缩放100%，间距：标准）；页码用小五号字底端居中。 ③具体要求： 题目（二号黑体居中）； 摘要（“摘要”二字用小二号黑体居中，隔行书写摘要的文字部分，小4号宋体）； 关键词（隔行顶格书写“关键词”三字，提炼3-5个关键词，用分号隔开，小4号黑体)； 正文部分采用三级标题； 第1章××(小二号黑体居中，段前0.5行) 1.1 ×××××小三号黑体×××××（段前、段后0.5行） 1.1.1小四号黑体（段前、段后0.5行） 参考文献（黑体小二号居中，段前0.5行），参考文献用五号宋体，参照《参考文献著录规则（GB/T 7714－2005）》。 实验一生成离散信号并计算其振幅谱 并将信号进行奇偶分解 一、实验原理 单位脉冲响应h(t)=exp(-a*t*t)*sin(2*3.14*f*t)进行离散抽样，分别得到t=0.002s，0.009s，0.011s采样的结果。用Excel软件绘图显示计算结果。并将信号进行奇偶分解，分别得到奇对称信号h(n)-h(-n)与偶对称信号h(n)+h(-n)。用Excel 软件绘图显示计算结果。 二、实验程序代码 （1）离散抽样 double a,t; a=2*f*f*log(m); int i; for(i=0;i 文件名：tstem.m（实验一、二需要） 程序： f unction tstem(xn,yn) %时域序列绘图函数 %xn:被绘图的信号数据序列，yn:绘图信号的纵坐标名称（字符串）n=0:length(xn)-1; stem(n,xn,'.'); xlabel('n');ylabel('yn'); axis([0,n(end),min(xn),1.2*max(xn)]); 文件名：tplot.m（实验一、四需要） 程序： function tplot(xn,T,yn) %时域序列连续曲线绘图函数 %xn:信号数据序列,yn:绘图信号的纵坐标名称（字符串） %T为采样间隔 n=0;length(xn)-1;t=n*T; plot(t,xn); xlabel('t/s');ylabel(yn); axis([0,t(end),min(xn),1.2*max(xn)]); 文件名：myplot.m（实验一、四需要） %(1)myplot;计算时域离散系统损耗函数并绘制曲线图。function myplot(B,A) %B为系统函数分子多项式系数向量 %A为系统函数分母多项式系数向量 [H,W]=freqz(B,A,1000) m=abs(H); plot(W/pi,20*log10(m/max(m)));grid on; xlabel('\omega/\pi');ylabel('幅度(dB)') axis([0,1,-80,5]);title('损耗函数曲线'); 文件名：mstem.m（实验一、三需要） 程序： function mstem(Xk) %mstem(Xk)绘制频域采样序列向量Xk的幅频特性图 M=length(Xk); k=0:M-1;wk=2*k/M;%产生M点DFT对应的采样点频率（关于pi归一化值） stem(wk,abs(Xk),'.');box on;%绘制M点DFT的幅频特性图xlabel('w/\pi');ylabel('幅度'); axis([0,2,0,1.2*max(abs(Xk))]); 文件名：mpplot.m（实验一需要） 数字信号处理作业实验题报告 第一章16.(1) 实验目的： 求解差分方程所描述的系统的单位脉冲响应和单位阶跃响应。 实验要求： 运用matlab求出y(n)=0.6y(n-1)-0.08y(n-2)+x(n)的单位脉冲响应和单位阶跃响应的示意图。 源程序： B1=1;A1=[1, -0.6, 0.08]; ys=2; %设差分方程 xn=[1, zeros(1, 20)]; %xn=单位脉冲序列，长度N=31 xi=filtic(B1, A1, ys); hn1=filter(B1, A1, xn, xi); %求系统输出信号hn1 n=0:length(hn1)-1; subplot(2, 1, 1);stem(n, hn1, '.') title('单位脉冲响应'); xlabel('n');ylabel('h(n)') xn=ones(1, 20); sn1=filter(B1, A1, xn, xi); %求系统输出信号sn1 n=0:length(sn1)-1; Subplot(2, 1, 2); stem(n, sn1, '.') title('单位阶跃响应'); xlabel('n'); ylabel('s(n)') 运行结果： 实验分析： 单位脉冲响应逐渐趋于0，阶跃响应保持不变，由此可见，是个稳定系统。 第二章31题 实验目的： 用matlab判断系统是否稳定。 实验要求： 用matlab画出系统的极，零点分布图，输入单位阶跃序列u(n)检查系统是否稳定。 源程序： A=[2, -2.98, 0.17, 2.3418, -1.5147]; B=[0, 0, 1, 5, -50]; subplot(2,1,1); zplane(B,A); %求H(z)的极点 p=roots(A); %求H(z)的模 pm=abs(p); if max(pm)<1 disp('系统因果稳定'), else,disp('系统因果不稳定'),end un=ones(1,800); sn=filter(B, A, un); n=0:length(sn)-1; subplot(2, 1, 2);plot(n, sn) xlabel('n');ylabel('s(n)') 成都理工大学 《信号处理基础》实验 开设时间：2013—2014学年第2学期 题目1：信号的产生和显示 一、实验目的： 认识基本信号 通过使用MATLAB 设计简单程序, 掌握对MATLAB 的基本使用方法 二、实验原理： 找出下列表达式的信号与：正弦信号、最小相位信号、最大相位信号、零相位信号的对应关系。 1、sin60t 2、e-60t sin60t 3、(1- e-60t)sin60t 4、e60t sin60t 三、实验内容： 产生上述信号的信号并显示 （1）t=[-pi/30:0.001:pi/30]; f=sin(60*t); plot(t,f) 产生图形如下： （2）t=[0:0.001:pi/30]; f=exp(-60*t).*sin(60*t); plot(t,f) 产生图形如下: （3）t=[-5*pi/30:0.001:5*pi/30]; f=(1-exp(-60*t)).*sin(60*t); plot(t,f) 产生图形如下： (4) t=[-pi/30:0.001:pi/30]; f=exp(6*t).*sin(60*t); plot(t,f) 产生如下波形: 四、实验结果与讨论： 讨论上述信号的特点 从第一个波形图可以看出，它的波形与正弦函数sin（t）的相像，只是相位上有改变，是一个正弦信号。最大相位信号的能量集中在后面，最小相位能量集中在前面，所以第二个是一个最小相位，第四个是一个最大相位信号。第三个由于波形在t>0时没有，所以是一个零相位信号。 题目2：频谱分析与显示 一、实验目的 初步认识频谱分析 子程序： function myplot(B,A) %myplot(B,A) %时域离散系统损耗函数绘图 %B为系统函数分子多项式系数向量 %A为系统函数分母多项式系数向量 [H,W]=freqz(B,A,1000); m=abs(H); plot(W/pi,20*log10(m/max(m)));grid on; xlabel('\omega/\pi');ylabel('幅度(dB)') axis([0,1,-80,5]);title('损耗函数曲线'); function tplot(xn,T,yn) %时域序列连续曲线绘图函数 % xn:信号数据序列，yn:绘图信号的纵坐标名称（字符串） % T为采样间隔 n=0:length(xn)-1;t=n*T; plot(t,xn); xlabel('t/s');ylabel(yn); axis([0,t(end),min(xn),1.2*max(xn)]) 程序： %实验4程序exp4.m % IIR数字滤波器设计及软件实现 clear all;close all Fs=10000;T=1/Fs; %采样频率 %调用信号产生函数mstg产生由三路抑制载波调幅信号相加构成的复合信号st st=mstg; %低通滤波器设计与实现========================================= fp=280;fs=450; wp=2*fp/Fs;ws=2*fs/Fs;rp=0.1;rs=60; %DF指标（低通滤波器的通、阻带边界频）[N,wp]=ellipord(wp,ws,rp,rs); %调用ellipord计算椭圆DF阶数N和通带截止频率wp [B,A]=ellip(N,rp,rs,wp); %调用ellip计算椭圆带通DF系统函数系数向量B和A y1t=filter(B,A,st); %滤波器软件实现 % 低通滤波器设计与实现绘图部分 figure(2);subplot(3,1,1); myplot(B,A); %调用绘图函数myplot绘制损耗函数曲线 yt='y_1(t)'; subplot(3,1,2);tplot(y1t,T,yt); %调用绘图函数tplot绘制滤波器输出波形 %带通滤波器设计与实现==================================================== fpl=440;fpu=560;fsl=275;fsu=900; wp=[2*fpl/Fs,2*fpu/Fs];ws=[2*fsl/Fs,2*fsu/Fs];rp=0.1;rs=60;数字信号处理实验二报告
数字信号处理实验报告
数字信号处理基础实验指导书
数字信号处理实验程序2.
数字信号处理实验作业
数字信号处理实验一
数字信号处理实验报告一
数字信号处理实验作业
数字信号处理实验报告(实验1_4)
数字信号处理实验1认识实验
数字信号处理实验答案完整版
数字信号处理实验及参考程序
实验一 基于Matlab的数字信号处理基本
数字信号处理第二章上机作业
数字信号处理基础实验报告_
数字信号处理上机实验代码
数字信号处理第二章上机题作业
数字信号处理基础实验报告 (2)
数字信号处理实验
相关主题
文本预览