第十一章：
1.
clear all;
close all;
t=-10:0.01:10;
ft=exp(abs(-0.4*t));
plot(t,ft)
2
clear all;
close all;
h=0.001
t=-3:h:3;
%ft的函数表达式%
ft=0.5*exp(-2*t).*stepfun(t,0);
%ft1的函数表达式%
ft1=0.5*exp(-2*(t-1)).*stepfun(t-1,0); %ft的导数表达式%
y1=diff(ft)*1/h;
%产生ft的波形%
subplot(3,1,1)
plot(t,ft);
title('f(t）波形');
%产生ft1的波形%
subplot(3,1,2)
plot(t,ft1);
title('f(t-1)的波形');
%产生ft的微分波形%
subplot(3,1,3)
plot(t(1:length(t)-1),y1);
title('f(t）的微分波形');
第十二章：
1). clear all;
close all;
f=fourier(cos(2*t).*sin(3*t)); ezplot(f);
2). clear all;
close all;
r=0.02;
t=-5:r:5;
N=200;
w=2*pi;
k=-N:N;
w=k*w/N;
f1=sinc(t/pi);
F=r*f1*exp(-j*t'*w);
F1=abs(F);
P1=angle(F);
subplot(3,1,1);
plot(t,f1);
grid on;
xlabel('t');
ylabel('f(t)');
title('f(t)');
subplot(3,1,2)
plot(w,F1);
xlabel('w');
grid on;
ylabel('F(jw)');
subplot(3,1,3);
plot(w,P1*180/pi);
grid;
xlabel('w');
ylabel('相位(度)');
3). clear all;
close all;
syms t w;
ifourier((((sin(w/4))/(w/4))^2),t);ans 4).
clear all;
close all;
r=0.02;
t=-5:r:5;
w=2*pi;
k=-N:N;
w=k*w/N;
f1=sin(2*pi*3*t+eps)./(pi*3*t+eps); F=r*f1*exp(-j*t'*w);
F1=abs(F);
P1=angle(F);
subplot(3,1,1);
plot(t,f1);
grid on;
xlabel('t');
ylabel('f(t)');
title('f(t)');
subplot(3,1,2)
plot(w,F1);
xlabel('w');
grid on;
ylabel('F(jw)');
subplot(3,1,3);
plot(w,P1*180/pi);
grid;
xlabel('w');
ylabel('相位(度)');
clear all;
close all;
r=0.02;
t=-5:r:5;
N=200;
w=2*pi;
k=-N:N;
w=k*w/N;
f1=sin(2*pi*(t-2)+eps)./(pi*(t-2)+eps); F=r*f1*exp(-j*t'*w);
F1=abs(F);
P1=angle(F);
subplot(3,1,1);
plot(t,f1);
grid on;
xlabel('t');
ylabel('f(t)');
title('f(t)');
subplot(3,1,2)
plot(w,F1);
xlabel('w');
grid on;
ylabel('F(jw)');
subplot(3,1,3);
plot(w,P1*180/pi);
grid;
xlabel('w');
ylabel('相位(度)');
clear all;
close all;
r=0.02;
t=-5:r:5;
N=200;
w=2*pi;
k=-N:N;
w=k*w/N;
f1=sin(2*pi*(t/3)+eps)./(pi*(t/3)+eps); F=r*f1*exp(-j*t'*w);
F1=abs(F);
P1=angle(F);
subplot(3,1,1);
plot(t,f1);
grid on;
xlabel('t');
ylabel('f(t)');
title('f(t)');
subplot(3,1,2)
plot(w,F1);
xlabel('w');
grid on;
ylabel('F(jw)');
subplot(3,1,3);
plot(w,P1*180/pi);
grid;
xlabel('w');
ylabel('相位(度)');
第十三章
syms s real;
a=[0,1,8,16];
b=[1,0,2,1];
sn=[s^3,s^2,s,1];
A=a*sn';
B=b*sn';
H=B/A;
e=sym('(1-t)*(Heaviside(t)-Heaviside(t-2))-(Heaviside(t-2)-Heaviside(t-3))'); E=laplace(e);
V os=H.*E;
V o=ilaplace(Vos);
V o=vpa(Vo,4);
ezplot(Vo,[0,10])
axis([0,10,-8,1])
第十四章
clear all;
close all;
for k=1:20;
f1(k)=2^k*k;
f2(k)=(-1)^k-(-2)^k;
end
y=conv(f1,f2);
y0=0;
yf=[y0 y(1:19)];
ny=1:20;
figure(1);
stem(ny-1,yf);
grid on;
xlabel('ny');
ylabel('yf');
title('离散系统的零状态响应');
y0=1;
y(1)=2;
y(2)=-6;
for k=3:20;
y(k)=2^(k-1)*(k-1)-3*y(k-1)-2*y(k-2);
end
yy=[y0 y(1:19)];
k=1:20;
figure(2);
stem(k-1,yy);
grid on;
xlabel('k');
ylabel('y(k)');
title('系统全响应');
k=1:20;
figure(3);
stem(k-1,yy-yf); grid on;
xlabel('k');
ylabel('y(k)');
title('零输入响应');
信号、系统与信号处理实验Ⅱ 实验报告 姓名：王健 学号：14072119 班级：14083413 上课时间：周五-六七八
实验名称：用双线性变换法设计IIR数字滤波器 一、实验目的 熟悉模拟巴特沃兹滤波器设计和用双线性变换法设计IIR数字滤波器的方法 二、实验原理与要求 实验原理 利用双线性变换法设计IIR数字滤波器，首先要设计出满足指标要求的模拟滤波器的传递函数Ha(s),然后由Ha(s)通过双线性变换可得要设计的IIR数字滤波器的系统函数H(z),如果给定的指标为数字滤波器的指标，直接利用模拟滤波器的低通原理，通过式子 到式子 的频率变换关系，可一步完成数字滤波器的设计。式中是低通模拟滤波器的截止频率 实验要求 （1）编写用双线性变换法设计的巴特沃兹低通IIR滤波器的程序，要求通带内频率低于,容许幅度误差在1dB之内，频率在到之间的阻带衰减大于10dB。 （2）用法设计的巴特沃兹低通IIR滤波器，要求使用buttord,butter和biliner函数，滤波器技术指标：取样频率为1Hz;通带内衰减小于1Db; 阻带临界频率0.3Hz,阻带内衰减大于25dB。 （3）以pi/64为取样间隔，在屏幕上打印出数字滤波器的频率区间[0 pi]上的幅频响应特性曲线。 （4）在屏幕上打印出H(z)的分子，分母多项式系数。 三、实验程序与结果 1. 用双线性变换法设计的巴特沃兹低通IIR滤波器的程序，要求通带内频率低于,容许幅度误差在1dB之内，频率在到之间的阻带衰减大于10dB。 clear;clc;close all; Rp=1; Rs=10; Fs=1; Ts=1/Fs
《信号与系统及MATLAB实现》实验指导书
前言 长期以来，《信号与系统》课程一直采用单一理论教学方式，同学们依靠做习题来巩固和理解教学内容，虽然手工演算训练了计算能力和思维方法，但是由于本课程数学公式推导较多，概念抽象，常需画各种波形，作题时难免花费很多时间，现在，我们给同学们介绍一种国际上公认的优秀科技应用软件MA TLAB，借助它我们可以在电脑上轻松地完成许多习题的演算和波形的绘制。 MA TLAB的功能非常强大，我们此处仅用到它的一部分，在后续课程中我们还会用到它，在未来地科学研究和工程设计中有可能继续用它，所以有兴趣的同学，可以对MA TLAB 再多了解一些。 MA TLAB究竟有那些特点呢？ 1．高效的数值计算和符号计算功能，使我们从繁杂的数学运算分析中解脱出来； 2．完备的图形处理功能，实现计算结果和编程的可视化； 3．友好的用户界面及接近数学表达式的自然化语言，易于学习和掌握； 4．功能丰富的应用工具箱，为我们提供了大量方便实用的处理工具； MA TLAB的这些特点，深受大家欢迎，由于个人电脑地普及，目前许多学校已将它做为本科生必须掌握的一种软件。正是基于这些背景，我们编写了这本《信号与系统及MA TLAB实现》指导书，内容包括信号的MA TLAB表示、基本运算、系统的时域分析、频域分析、S域分析、状态变量分析等。通过这些练习，同学们在学习《信号与系统》的同时，掌握MA TLAB的基本应用，学会应用MA TLAB的数值计算和符号计算功能，摆脱烦琐的数学运算，从而更注重于信号与系统的基本分析方法和应用的理解与思考，将课程的重点、难点及部分习题用MA TLAB进行形象、直观的可视化计算机模拟与仿真实现，加深对信号与系统的基本原理、方法及应用的理解，为学习后续课程打好基础。另外同学们在进行实验时，最好事先预习一些MA TLAB的有关知识，以便更好地完成实验，同时实验中也可利用MA TLAB的help命令了解具体语句以及指令的使用方法。
实验5 IIR和FIR滤波器过滤信号的实现及比较：以心电信号为例 一、实验目的 1、探究心电信号的初步分析。心电信号(频率-般在0.05Hz ~100Hz范围)是一种基本的人体生理信号，体表检测人体心电信号中常带有工频干扰(50HZ)、基线漂移(频率低于0.5Hz)和肌电干扰等各种噪声。 2、为了得到不失真的原始心电信号，需要滤波预处理。设计数字低通滤波器、高通滤波器、带阻滤波器，用MATLAB软件对含噪心电信号分别进行高通、带阻和低通滤波等处理，将心电信号中的低频基线漂移、50Hz 工频高频和高频杂波进行滤除。 3、通过观察对含噪心电图信号的滤波作用，获得数字滤波的感性知识。 二、实验要求及内容 实验题目： 给定一组干净心电信号数据，数据文件存于C盘Ecg.txt。采样频率Fs = 500Hz。 1、编写程序读出心电信号，并在屏幕上打印出其波形。 2、产生模拟高斯白噪声信号，与干净心电混合,设计一个IIR低通滤波器和一个FIR 低通滤波器分别滤除心电信号中的白噪声干扰，调整白噪声信噪比大小，对滤波前后的心电信号的频谱进行分析比较。其中数字低通滤波器指标要求，通带截止频率Wp=0.1π，阻带截止频率 Ws=0.16π，阻带衰减不小于15 dB，通带衰减不大于1 dB。 要求:编写一个IIR低通滤波器和一个FIR低通滤波器仿真程序，在屏幕上打印出数字滤波器的频率区间[0, π]上的幅频响应特性由线(H(e^jw)) ;计算其对含噪心电信号的低通滤波响应序列，并在屏幕上打印出干净心电信号波形，含工频干扰的心电信号波形以及IIR低通滤波和FIR低通后的信号波形，并进行比较;同时对滤波前后的心电信号的频谱进行分析比较，并在屏幕上打印出滤波前后的心电信号的频谱，观察其变化。 3、产生模拟工频信号，与干净心电混合，设计一个带阻滤波器(50Hz 陷波器)滤除心电信号中的电源线干扰，调整工频幅度大小，对滤波前后的心电信号的频谱进行分析比较。其中带阻滤波器指标要求，通带下限频率Wp1=0.18π，阻带下截止频率Ws1=0.192 π，阻带上截止频率Ws2=0.208π，通带上限频率Wp2=0.22π，阻带衰减不小于15 dB, 通带衰减不大于1 dB。 要求:编写IIR带阻滤波器仿真程序,在屏幕上打印出数字滤波器的频率区间[0, π]上的幅频响应特性由线(H(e^jw ));计算其对含工频干扰的心电信号的带阻滤波响应序列，并在屏幕上打印出干净心电信号波形，含工频干扰的心电信号波形以及滤波后的信号波形，并进行比较;同时对滤波前后的心电信号的频谱进行分析比较，并在屏幕上打印出滤波前后的心电信号的频谱，观察其变化。 4、产生模拟基线漂移信号，与干净心电信号混合，设计一个高通滤波器滤除心电信号中的基线低频干扰，调整基线的幅度大小,对滤波前后的心电信号的频谱进行分析比较。其中，高通滤波器指标要求，通带截止频率Wp=0.0028π，阻带截止频率Ws=0.0012π，阻带衰减不小于15 dB，通带衰减不大于1 dB。 要求:编写IIR高通滤波器(或FIR高通滤波器)仿真程序，在屏幕上打印出数字滤波器的频率区间[0,π]上的幅频响应特性由线(H(e^jw);计算其对含基线低频干扰的心电信号的高通滤波响应序
信号与系统综合实验项目doc 信号与系统综合实验项目 (竞 实 验 指 导 项目一 用MATLAB 验证时域抽样定理 目的： 通过MATLAB 编程实现对时域抽样定理的验证，加深抽样定理的明白得。同时训练应用运算机分析咨询题的能力。 任务： 连续信号f(t)=cos(8*pi*t)+2*sin(40*pi*t)+cos(24*pi*t)，通过理想抽样后得到抽样信号fs(t)，通过理想低通滤波器后重构信号f(t)。 方法： 1、确定f(t)的最高频率fm 。关于无限带宽信号，确定最高频率fm 的方法：设其频谱的模降到10-5左右时的频率为fm 。 2、确定Nyquist 抽样间隔T N 。选定两个抽样时刻：T S 信号与系统期末考试试题 一、选择题（共10题，每题3分 ，共30分，每题给出四个答案，其中只有一个正确的） 1、 卷积f 1(k+5)*f 2(k-3) 等于 。 （A ）f 1(k)*f 2(k) （B ）f 1(k)*f 2(k-8)（C ）f 1(k)*f 2(k+8)（D ）f 1(k+3)*f 2(k-3) 2、 积分 dt t t ? ∞ ∞ --+)21()2(δ等于 。 （A ）1.25（B ）2.5（C ）3（D ）5 3、 序列f(k)=-u(-k)的z 变换等于 。 （A ） 1-z z （B ）-1-z z （C ）11-z （D ）1 1--z 4、 若y(t)=f(t)*h(t),则f(2t)*h(2t)等于 。 （A ） )2(41t y （B ）)2(21t y （C ）)4(41t y （D ）)4(2 1 t y 5、 已知一个线性时不变系统的阶跃相应g(t)=2e -2t u(t)+)(t δ,当输入f(t)=3e —t u(t)时，系 统的零状态响应y f (t)等于 （A ）(-9e -t +12e -2t )u(t) （B ）(3-9e -t +12e -2t )u(t) （C ）)(t δ+(-6e -t +8e -2t )u(t) （D ）3)(t δ +(-9e -t +12e -2t )u(t) 6、 连续周期信号的频谱具有 （A ） 连续性、周期性 （B ）连续性、收敛性 （C ）离散性、周期性 （D ）离散性、收敛性 7、 周期序列2)455.1(0 +k COS π的 周期N 等于 （A ） 1（B ）2（C ）3（D ）4 8、序列和 ()∑∞ -∞ =-k k 1δ等于 （A ）1 (B) ∞ (C) ()1-k u (D) ()1-k ku 9、单边拉普拉斯变换()s e s s s F 22 12-+= 的愿函数等于 ()()t tu A ()()2-t tu B ()()()t u t C 2- ()()()22--t u t D 10、信号()()23-=-t u te t f t 的单边拉氏变换()s F 等于 ()A ()()()232372+++-s e s s ()() 2 23+-s e B s 实验一 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; 学生实验报告 （理工类） 课程名称：信号与线性系统专业班级：M11通信工程 学生学号：1121413017 学生姓名：王金龙 所属院部：龙蟠学院指导教师：杨娟 20 11 ——20 12 学年第 1 学期 金陵科技学院教务处制 实验报告书写要求 实验报告原则上要求学生手写，要求书写工整。若因课程特点需打印的，要遵照以下字体、字号、间距等的具体要求。纸张一律采用A4的纸张。 实验报告书写说明 实验报告中一至四项内容为必填项，包括实验目的和要求；实验仪器和设备；实验内容与过程；实验结果与分析。各院部可根据学科特点和实验具体要求增加项目。 填写注意事项 （1）细致观察，及时、准确、如实记录。 （2）准确说明，层次清晰。 （3）尽量采用专用术语来说明事物。 （4）外文、符号、公式要准确，应使用统一规定的名词和符号。 （5）应独立完成实验报告的书写，严禁抄袭、复印，一经发现，以零分论处。 实验报告批改说明 实验报告的批改要及时、认真、仔细，一律用红色笔批改。实验报告的批改成绩采用百分制，具体评分标准由各院部自行制定。 实验报告装订要求 实验批改完毕后，任课老师将每门课程的每个实验项目的实验报告以自然班为单位、按学号升序排列，装订成册，并附上一份该门课程的实验大纲。 实验项目名称：常用连续信号的表示 实验学时： 2学时 同组学生姓名： 无 实验地点： A207 实验日期： 11.12.6 实验成绩： 批改教师： 杨娟 批改时间： 一、实验目的和要求 熟悉MATLAB 软件；利用MATLAB 软件，绘制出常用的连续时间信号。 二、实验仪器和设备 586以上计算机，装有MATLAB7.0软件 三、实验过程 1. 绘制正弦信号)t Asin t (f 0?ω+=（），其中A=1，πω2=，6/π?=； 2. 绘制指数信号at Ae t (f =），其中A=1，0.4a -=； 3. 绘制矩形脉冲信号，脉冲宽度为2； 4. 绘制三角波脉冲信号，脉冲宽度为4；斜度为0.5； 5. 对上题三角波脉冲信号进行尺度变换，分别得出)2t (f ，)2t 2(f -； 6. 绘制抽样函数Sa （t ）,t 取值在-3π到+3π之间； 7. 绘制周期矩形脉冲信号，参数自定； 8. 绘制周期三角脉冲信号，参数自定。 四、实验结果与分析 1．制正弦信号)t Asin t (f 0?ω+=（），其中A=1，πω2=，6/π?= 实验代码： A=1; 数字信号处理实验实验一离散时间信号与系统及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]; 广州大学 综合设计性实验 报告册 实验项目选频网络的设计及应用研究 学院物电学院年级专业班电子131 姓名朱大神学号成绩 实验地点电子楼316 指导老师 《综合设计性实验》预习报告 实验项目：选频网络的设计及应用研究 一 引言： 选频网络在信号分解、振荡电路及其收音机等方面有诸多应用。比如，利用选频网络可以挑选出一个周期信号中的基波和高次谐波。选频网络的类型和结构有很多，本实验将通过设计有源带通滤波器实现选频。 二 实验目的： （1）熟悉选频网络特性、结构及其应用，掌握选频网络的特点及其设计方法。 （2）学会使用交流毫伏表和示波器测定选频网络的幅频特性和相频特性。 （3）学会使用Multisim 进行电路仿真。 三 实验原理： 带通滤波器： 这种滤波器的作用是只允许在某一个通频带范围内的信号通过，而比通频带下限频率低和比上限频率高的信号均加以衰减和抑制。 典型的带通滤波器可以从二阶低通滤波器中将其中一级改成高通而成，如图1所示。 电路性能参数可由下面各式求出。 通带增益：CB R R R R A f vp 144+= 其中B 为通频带宽。 中心频率：)1 1(121 3 12 20R R C R f += π 通带宽度：)2 1(14 321R R R R R C B f -+= 品质因数：B f Q 0 = 此电路的优点是，改变f R 和4R 的比值，就可以改变通带宽度B 而不会影响中心频率0f 。 四 实验内容： 设计一个中心频率Hz f 20000=，品质因数5>Q 的带通滤波器。 五 重点问题： （1）确定带通滤波器的中心频率、上限频率及下限频率。 （2）验证滤波器是否能筛选出方波的三次谐波。 六 参考文献： [1]熊伟等.Multisim 7 电路设计及仿真应用.北京：清华大学出版社，2005. [2]吴正光，郑颜.电子技术实验仿真与实践.北京：科学出版社，2008. [4]童诗白等.模拟电子技术基础（第三版）.北京：高等教育出版社， 2001. 图1 二阶带通滤波器 重庆大学信号与线性系统期末考试试题 一、填空题：（30分，每小题3分） 1. =-? ∞ ∞ -dt t t )()5cos 2(δ 。 2. ()dt t e t 12-?+∞ ∞ --δ= 。 3. 已知 f (t )的傅里叶变换为F (j ω), 则f (2t -3)的傅里叶变换为 。 4. 已知 6 51 )(2 +++= s s s s F ，则=+)0(f ; =∞)(f 。 5. 已知 ω ωπδεj t FT 1 )()]([+=，则=)]([t t FT ε 。 6. 已知周期信号 )4sin()2cos()(t t t f +=，其基波频率为 rad/s ； 周期为 s 。 7. 已知 )5(2)2(3)(-+-=n n k f δδ，其Z 变换 =)(Z F ；收敛域为 。 8. 已知连续系统函数1342 3)(2 3+--+= s s s s s H ，试判断系统的稳定性： 。 9．已知离散系统函数1.07.02 )(2+-+=z z z z H ，试判断系统的稳定性： 。 10．如图所示是离散系统的Z 域框图，该系统的系统函数H(z)= 。 二．(15分)如下方程和非零起始条件表示的连续时间因果LTI 系统， ?????==+=++-- 5 )0(',2)0() (52)(4522y y t f dt df t y dt dy dt y d 已知输入 )()(2t e t f t ε-=时，试用拉普拉斯变换的方法求系统的零状态响应 )(t y zs 和零输入响应)(t y zi ，0≥t 以及系统的全响应),(t y 0≥t 。 三．（14分） ① 已知2 36 62)(22++++=s s s s s F ,2]Re[->s ,试求其拉氏逆变换f (t )； ② 已知) 2(2 35)(2>+-=z z z z z X ，试求其逆Z 变换)(n x 。 四 （10分）计算下列卷积： 1. }1,0,6,4,3{}4,1,2,1{)()(21--*=*k f k f ； 2． )(3)(23t e t e t t εε--* 。 中南大学 信号与系统试验报告 姓名: 学号: 专业班级:自动化 实验一 基本信号的生成 1．实验目的 ● 学会使用MATLAB 产生各种常见的连续时间信号与离散时间信号； ● 通过MATLAB 中的绘图工具对产生的信号进行观察，加深对常用信号的 理解； ● 熟悉MATLAB 的基本操作，以及一些基本函数的使用，为以后的实验奠 定基础。 2．实验内容 ⑴ 运行以上九个例子程序，掌握一些常用基本信号的特点及其MATLAB 实现方法；改变有关参数，进一步观察信号波形的变化。 ⑵ 在 k [10:10]=- 范围内产生并画出以下信号： a) 1f [k][k]δ=； b) 2f [k][k+2]δ=； c) 3f [k][k-4]δ=； d) 4f [k]2[k+2][k-4]δδ=-。 源程序： k=-10:10; f1k=[zeros(1,10),1,zeros(1,10)]; subplot(2,2,1) stem(k,f1k) title('f1[k]') f2k=[zeros(1,8),1,zeros(1,12)]; subplot(2,2,2) stem(k,f2k) title('f2[k]') f3k=[zeros(1,14),1,zeros(1,6)]; subplot(2,2,3) stem(k,f3k) title('f3[k]') f4k=2*f2k-f3k; subplot(2,2,4) stem(k,f4k) title('f4[k]') ⑶ 在 k [0:31]=范围内产生并画出以下信号： a) ()()k k 144f [k]sin cos π π=； b) ()2k 24f [k]cos π =； c) ()()k k 348f [k]sin cos π π=。 请问这三个信号的基波周期分别是多少？ 源程序： k=0:31; f1k=sin(pi/4*k).*cos(pi/4*k); subplot(3,1,1) stem(k,f1k) title('f1[k]') f2k=(cos(pi/4*k)).^2; subplot(3,1,2) stem(k,f2k) title('f2[k]') f3k=sin(pi/4*k).*cos(pi/8*k); subplot(3,1,3) stem(k,f3k) title('f3[k]') 其中f1[k]的基波周期是4, f2[k]的基波周期是4, f3[k]的基波周期是16。 2006-2007学年第二学期《信号与线性系统》（课内）试卷A 卷 一、计算题（共45分） 1.（5分）计算积分dt t t t )6 ()sin (π δ- +?+∞ ∞-的值。 2.（5分）绘出函数)1()]1()([-+--t u t u t u t 的波形图。 3.（6分）已知)2()()(),1()()(21--=--=t u t u t f t u t u t f ，求卷积)()(21t f t f *。 4.（6分）若)(t f 的傅里叶变换已知，记为)(ωF ，求)1()1(t f t --对应的傅里叶变换。 5.（6分）如下图所示信号，已知其傅里叶变换，记为)(ωF ， 求： （1）)0(F ； （2）?+∞ ∞ -ωωd F )(。 6.（5分）已知)(t f 对应的拉氏变换为)(s F ，求)/(/a t f e a t -（0>a ）对应的拉氏变换。 7.（6分） 已知)(t f 对应的拉氏变换2 3)(2 +-=-s s e s F s ，求)(t f 8．（6分）线性时不变系统的单位样值响应为)(n h ，输入为)(n x ，且有 )4()()()(--==n u n u n x n h ，求输出)(n y ，并绘图示出)(n y 。 二、综合题（共计55分） 1、（10分）系统如图所示，已知t t x 2000cos )(=，t t t f 2000cos 100cos )(=，理想低通滤波器)300()300()(--+=ωωωu u H ，求滤波器的响应信号)(t y 。 x(t) y(t) f(t) 实验六: 用FFT对信号作频谱分析 一、实验目的 1．了解双音多频信号的产生、检测、包括对双音多频信号进行DFT时的参数选择等。 2．初步了解数字信号处理在是集中的使用方法和重要性。 3．掌握matlab的开发环境。 二、实验原理与方法 1、引言 双音多频（Dual Tone Multi Frequency, DTMF）信号是音频电话中的拨号信号，由美国AT&T贝尔公司实验室研制，并用于电话网络中。这种信号制式具有很高的拨号速度，且容易自动监测识别，很快就代替了原有的用脉冲计数方式的拨号制式。这种双音多频信号制式不仅用在电话网络中，还可以用于传输十进制数据的其它通信系统中，用于电子邮件和银行系统中。这些系统中用户可以用电话发送DTMF信号选择语音菜单进行操作。DTMF信号系统是一个典型的小型信号处理系统，它要用数字方法产生模拟信号并进行传输，其中还用到了D/A变换器；在接收端用A/D变换器将其转换成数字信号，并进行数字信号处理与识别。为了系统的检测速度并降低成本，还开发一种特殊的DFT算法，称为戈泽尔(Goertzel)算法，这种算法既可以用硬件（专用芯片）实现，也可以用软件实现。下面首先介绍双音多频信号的产生方法和检测方法，包括戈泽尔算法，最后进行模拟实验。下面先介绍电话中的DTMF信号的组成。在电话中，数字0~9的中每一个都用两个不同的单音频传输，所用的8个频率分成高频带和低频带两组，低频带有四个频率：679Hz,770Hz,852Hz和941Hz；高频带也有四个频率：1209Hz,1336Hz,1477Hz和1633Hz.。每一个数字均由高、低频带中各一个频率构成，例如1用697Hz和1209Hz两个频率，信号用表示。这样8个频率形成16种不同的双频信号。具体号码以及符号对应的频率如表10.6.1所示。表中最后一列在电话中暂时未用。DTMF信号在电话中有两种作用，一个是用拨号信号去控制交换机接通被叫的用户电话机，另一个作用是控制电话机的各种动作，如播 实验一 维纳滤波 1. 实验内容 设计一个维纳滤波器： （1） 产生三组观测数据，首先根据()(1)()s n as n w n =-+产生信号()s n ，将其加噪，（信噪比分别为20,10,6dB dB dB ），得到观测数据123(),(),()x n x n x n 。 （2） 估计()i x n ，1,2,3i =的AR 模型参数。假设信号长度为L ，AR 模型阶数为N ，分析实验结果，并讨论改变L ，N 对实验结果的影响。 2. 实验原理 滤波目的是从被噪声污染的信号中分离出有用的信号来，最大限度地抑制噪声。对信号进行滤波的实质就是对信号进行估计。滤波问题就是设计一个线性滤波器，使得滤波器的输出信号()y n 是期望响应()s n 的一个估计值。下图就是观测信号的组成和信号滤波的一般模型。 观测信号()()()x n s n v n =+ 信号滤波的一般模型 维纳滤波解决从噪声中提取信号的滤波问题，并以估计的结果与真值之间的误差均方值最小作为最佳准则。它根据()()(),1, ,x n x n x n m --估计信号的当前 值，它的解以系统的系统函数()H z 或单位脉冲()h n 形式给出，这种系统常称为最佳线性滤波器。 维纳滤波器设计的任务就是选择()h n ，使其输出信号()y n 与期望信号()d n 误差的均方值最小。 假设滤波系统()h n 是一个线性时不变系统，它的()h n 和输入信号都是复函数，设 ()()()h n a n jb n =+ 0,1, n = 考虑系统的因果性，可得到滤波器的输出 ()()()()()0 *m y n h n x n h m x n m +∞ ===-∑ 0,1, n = 设期望信号()d n ，误差信号()e n 及其均方误差()2 E e n ???? 分别为 ()()()()()e n d n y n s n y n =-=- ()()()()()()22 2 0m E e n E d n y n E d n h m x n m ∞=?? ????=-=--????? ????? ∑ 要使均方误差为最小，需满足： ()() 2 0E e n h j ?????=? 整理得()()0E x n j e n *??-=??，等价于()()0E x n j e n * ??-=?? 上式说明，均方误差达到最小值的充要条件使误差信号与任一进入估计的输入信号正交，这就是正交性原理。 将()()0E x n j e n * ??-=??展开，得 ()()()()00m E x n k d n h m x m +∞ *** =????--=?? ???? ?∑ 整理得 ()()()0 dx xx m r k h m r m k +∞ *=-=-∑ 0,1,2, k = 等价于()()()()()0 dx xx xx m r k h m r k m h k r k +∞ ==-=*∑ 0,1,2, k = 此式称为维纳-霍夫（Wiener-Holf ）方程。解此方程可得到最优权系数 012,,, h h h ，此式是Wiener 滤波器的一般方程。 定义 信号与系统综合设计实验项目 实 验 指 导 项目一 用MATLAB 验证时域抽样定理 目的： 通过MATLAB 编程实现对时域抽样定理的验证，加深抽样定理的理解。同时训练应用计算机分析问题的能力。 任务： 连续信号f(t)=cos(8*pi*t)+2*sin(40*pi*t)+cos(24*pi*t)，经过理想抽样后得到抽样信号fs(t)，通过理想低通滤波器后重构信号f(t)。 方法： 1、确定f(t)的最高频率fm 。对于无限带宽信号，确定最高频率fm 的方法：设其频谱的模降到10-5左右时的频率为fm 。 2、确定Nyquist 抽样间隔T N 。选定两个抽样时间：T S 实验2 利用快速离散傅里叶变换分析信号的频谱 一、实验目的 1、通过这一实验，能够熟练掌握快速离散傅里叶变换FFT的原理及其用FFT 进行频谱分析的基本方法。 2、在通过计算机上用软件实现FFT 及信号的频谱分析。 3、通过实验对离散傅里叶变换的主要性质及FFT 在数字信号处理中的重要作用有进一步的了解。 二、实验要求及内容 第一题代码： clear all; N=32; n=0:N-1; xn=sin(2*pi*n/32); subplot(2,1,1); plot(n,xn); ylabel('x(n)'); title('x(n) N=32'); XK=fft(xn,N); magXK=abs(XK); %幅频特性 phaXK=angle(XK); %相频特性 subplot(2,1,2); k=0:length(magXK)-1; stem(k,magXK,'.'); %信号幅频特性曲线xlabel('k'); ylabel('|X(k)|'); title('X(k) N=32'); 第二题代码： clear all; N=100; n=0:N-1; xn=0.9*sin(2*pi/N*n)+0.6*cos(2*pi/(N/3)*n); subplot(2,1,1); plot(n,xn); xlabel('n'); ylabel('x(n)'); title('x(n) N=32'); XK=fft(xn,N); magXK=abs(XK); %幅频特性 phaXK=angle(XK); %相频特性 subplot(2,1,2); k=0:length(magXK)-1; stem(k,magXK,'.'); %信号幅频特性曲线xlabel('k'); ylabel('|X(k)|'); title('X(k) N=100'); 第三题代码： clear all; N=8; x=[ones(1,8),zeros(1,N-8)]; n=0:N-1; X=dft(x,N); magX=abs(X); phaX=angle(X)*180/pi; k=(0:length(magX)'-1)*N/length(magX); subplot(2,2,1); stem(n,x); 课程名称 信号与线性系统A 考试学期 08-07 得分 适用专业 微电、物理、 考试形式 闭卷 考试时间 120分钟 姓名 班级 学号 一、选择题（每小题可能有一个或几个正确答案，将正确的题号填入[ ]内） 1．f （5-2t ）是如下运算的结果————————（ C ） （A ）f （-2t ）右移5 （B ）f （-2t ）左移5 （C ）f （-2t ）右移 2 5 （D ）f （-2t ）左移25 2．已知)()(),()(21t u e t f t u t f at -==，可以求得=)(*)(21t f t f —————（ C ） （A ）1-at e - （B ）at e - （C ）)1(1at e a -- （D ）at e a -1 3．线性系统响应满足以下规律————————————（AD ） （A ）若起始状态为零，则零输入响应为零。 （B ）若起始状态为零，则零状态响应为零。 （C ）若系统的零状态响应为零，则强迫响应也为零。 （D ）若激励信号为零，零输入响应就是自由响应。 4．若对f （t ）进行理想取样，其奈奎斯特取样频率为f s ，则对)23 1 (-t f 进行取 样，其奈奎斯特取样频率为————————（B ） （A ）3f s （B ） s f 31 （C ）3（f s -2） （D ）)2(3 1 -s f 5．理想不失真传输系统的传输函数H （jω）是 ————————（B ） （A ）0j t Ke ω- （B ）0 t j Ke ω- （C ）0 t j Ke ω-[]()()c c u u ωωωω+-- （D ）00 j t Ke ω- （00,,,c t k ωω为常数） 6．已知Z 变换Z 1 311 )]([--= z n x ，收敛域3z >，则逆变换x （n ）为——（ A ） （A ）)(3n u n （C ）3(1)n u n - （B ）)(3n u n -- （D ）)1(3----n u n 数字信号处理实验答案 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),绘出四个周期。 信号与系统实验综合设计 《二阶电路系统的设计与测试分析》年级： 班级： 姓名： 学号： 日期： 二阶电路系统的设计与测试分析一、实验原理 二阶电路图如下所示 描述这种电路的微分方程为 y’’ (t) + 2αy’ (t) +ω2y (t) =ω2f (t) ，式中ω2= 1/LC ，α= R/2L。 其特征根为 事实上，针对该电路，可列出以下方程： I (t) = C*Uc’(t) ① R*I (t) + L*I’(t) + Uc(t) = Us(t) ② 将①式代入②式：CR*Uc’(t) + LC*Uc’’(t) + Uc(t) = Us(t) 可化为：Uc’’(t) + (R/L)*Uc’(t) + (1/LC)Uc(t) = (1/LC)*Us(t) 将各元件值代入：Uc’’(t) + 15000 *Uc’(t) + 10^5 *Uc(t) = 10^5 *Us(t) ω2= 10^5 ，α= 7500 可以看出，α>ω，这属于过阻尼情况。选取过阻尼情况是考虑到实测电路时电路输出在示波器上的波形观察起来可以更加明显。 接下来，利用软件SystemView对该电路进行时域、频域、S域分析。 二、Systemview仿真 由电路实际方程Uc’’(t) + 15000 *Uc’(t) + 10^5 *Uc(t) = 10^5 *Us(t) 移项可得； Uc’’(t) = -15000 *Uc’(t) - 10^5 *Uc(t) + 10^5 *Us(t) 据此画出系统框图 其系统函数H（jω）为1/（LCs 2+ RCs + 1） 1、阶跃响应 ①时域仿真波形 2012-2013学年第二学期《信号与线性系统》（课内）试卷A 卷 一、计算题（共45分） 1.（5分）计算积分dt t t t )6 ()sin (π δ- +? +∞ ∞ -的值。 2.（5分）绘出函数)1()]1()([-+--t u t u t u t 的波形图。 3.（6分）已知)2()()(),1()()(21--=--=t u t u t f t u t u t f ，求卷积)()(21t f t f *。 4.（6分）若)(t f 的傅里叶变换已知，记为)(ωF ，求)1()1(t f t --对应的傅里叶变换。 5.（6分）如下图所示信号，已知其傅里叶变换，记为)(ωF ， 求： （1）)0(F ； （2）? +∞ ∞ -ωωd F )(。 6.（5分）已知)(t f 对应的拉氏变换为)(s F ，求)/(/a t f e a t -（0>a ）对应的拉氏变换。 7.（6分） 已知)(t f 对应的拉氏变换2 3)(2+-=-s s e s F s ，求)(t f 8．（6分）线性时不变系统的单位样值响应为)(n h ，输入为)(n x ，且有 )4()()()(--==n u n u n x n h ，求输出)(n y ，并绘图示出)(n y 。 二、综合题（共计55分） 1、（10分）系统如图所示，已知t t x 2000 cos )(=，t t t f 2000cos 100cos )(=，理想低通滤波器)300()300()(--+=ωωωu u H ，求滤波器的响应信号)(t y 。 y(t) f(t)信号与系统期末考试试题(有答案的)
数字信号处理实验报告(实验1_4)
信号与系统实验报告1
数字信号处理实验及参考程序
信号与系统综合实验报告-带通滤波器的设计DOC
信号与系统期末考试试题
信号与系统实验报告
《信号与线性系统》期末试卷
数字信号处理实验
数字信号处理实验——维纳滤波器设计..
信号与系统综合实验项目信号与系统综合实验项目竞
杭电_数字信号处理课程设计_实验2
信号与系统期末试题与答案
数字信号处理实验答案完整版
信号与系统实验综合设计课设
《信号与线性系统》期末试卷要点
相关主题
文本预览