数字信号处理(北航)实验二报告
数字信号处理实验二
信号的分析与处理综合实验
38152111 张艾一、实验目的
综合运用数字信号处理的理论知识进行信号的采样，重构，频谱分析和滤波器的设计，通过理论推导得出相应结论，再利用Matlab作为编程工具进行计算机实现，从而加深对所学知识的理解，建立概念。
二、基本要求
1．掌握数字信号处理的基本概念、基本理论和基本方法；
2．学会MATLAB的使用，掌握MATLAB的程序设计方法；
3．掌握用MATLAB设计简单实验验证采样定理的方法；
4．掌握在Windows环境下语音信号采集的方法；
5．学会用MATLAB对信号进行频谱分析；
6．掌握MATLAB设计FIR和IIR数字滤波器的方法；
三、实验内容
1．利用简单正弦信号设计实验验证采样定理：
（1）Matlab产生离散信号的方法，作图的方法，以及基本运算操作
（2）对连续正弦信号以不同的采样频率作采样
（3）对采样前后信号进行傅立叶变换，并画频谱图
（4）分析采样前后频谱的有变化，验证采样定理。
掌握画频谱图的方法，深刻理解采样频率，信号频率，采样点数，频率分辨率等概念2．真实语音信号的采样重构：录制一段自己的语音信号，并对录制的信号进行采样；画出采样前后语音信号的时域波形和频谱图；对降采样后的信号进行插值重构，滤波，恢复原信号。
（1）语音信号的采集
（2）降采样的实现（改变了信号的采样率）
（3）以不同采样率采样后，语音信号的频谱分析
（4）采样前后声音的变化
（5）对降采样后的信号进行插值重构，滤波，恢复原信号
3．带噪声语音信号的频谱分析
（1）设计一频率已知的噪声信号，与实验2中原始语音信号相加，构造带噪声信号（2）画出原始语音信号和加噪声后信号，以及它们的频谱图
（3）利用频谱图分析噪声信号和原语音信号的不同特性
4．对带噪声语音信号滤波去噪：给定滤波器性能指标，采样窗函数法或双线性变换设计滤波器，并画出滤波器的频率响应；然后用自己设计的滤波器对采样的语音信号进行滤波，画出滤波后信号的时域波形和频谱，并对滤波前后的信号进行对比，分析信号的变化；
回放语音信号；
（1）分析带噪声信号频谱，找出噪声所在的频率段
（2）利用matlab中已有的滤波器滤波
（3）根据语音信号特点，自己设计滤波器滤波
（4）比较各种滤波器性能（至少四种），选择一种合适的滤波器将噪声信号滤除
（5）回放语音信号，比较滤波前后声音的变化
四、主要实验仪器及材料
微型计算机、Matlab。
五、实验步骤及结果分析
1、设计一简单正弦信号，通过改变采样率观察采样前后的信号变化。
选取正弦信号sin(2100.3)
π??+，
t
原始信号：0.01s的采样间隔（采样频率为100Hz）表示。
分别以5Hz，10Hz，20Hz，40Hz，80Hz，200Hz对原始信号进行采样，画出采样前后的信号，并画出其频谱图，对比前后的变化，验证采样定理。
结果图如下
fs=5Hz fs=10Hz
fs=20Hz fs=40Hz
Fs=80Hz fs=200Hz
对比不同采用那个频率下信号的变化可验证：
只有在fs（采样频率）≥2f（信号频率）时，信号才能保持原有时域、频域特性，不失真。
2、对真实语音信号的采样、重构
⑴读取样本声音文件“hello.wav”，并对其进行1/2倍，1/4倍，1/20倍，1/50倍，1/100倍的降采样，画出降采样前后信号的波形和频谱；
⑵对采样后的语音信号进行插值重构，滤波，恢复原始信号。画出插值前后信号的波形以及频谱图。
结果图如下：原始信号
1/2采样
1/4采样
1/20 采样
1/50 采样
1/100 采样
将重构后的信号与原信号进行比较，采样率为1/2 和1/4 时信号还能基本保持原样，采样率为1/20 时已经丢失了一部分细节，但大体上还保持原信息，而采样率提高到1/50 和
1/100时，基本上看不出信号原样。随着采样率的提高，重构后的信号听起来也是越来越模糊，1/20 时就听不清楚了。
3、对原始语音信号加噪声
对原始信号“hello.wav”加上幅值为0.01，频率为5000的正弦波噪声信号，
时域图及频域图如下：
4、设计数字滤波器
*用窗函数法设计FIR高通，低通，带通，带阻滤波器(fs=22050)
1）低通滤波器：fp=4500Hz(0.41)，fc=6500(0.59)Hz，Rs=30dB，Rp=1dB。
（fp：通带截至频率；fc：阻带截至频率；Rs：通带波纹；Rp：阻带波纹）由题意，阻带衰减不小于30dB，根据FIR滤波器各种窗函数的基本参数，选择Hanning窗。在窗函数设计法中，要求设计的频率归一化到0~π之间，Nyquist频率对应于π，因此通带和阻带边界频率为0.41π和0.59π。
低通滤波器的幅频和相频特性如下：
2）高通滤波器：fc=4500Hz(0.41)，fp=6500Hz(0.59)，Rs=30dB，Rp=1dB。
（fp：通带截至频率；fc：阻带截至频率；Rs：通带波纹；Rp：阻带波纹）选择Hanning窗。通带和阻带边界频率为0.59π和0.41π。
高通滤波器的幅频和相频特性如下：
3）带阻滤波器：fp1=4800Hz(0.44)，fp2=5200Hz(0.47)，fc1=4600 Hz(0.42)，fc2=5400 Hz(0.49)，Rs=30dB，Rp=1dB。（[fp1 fp2]：阻带截至频率；[fc1 fc2]：通带截至频率）
选择Hanning窗。通带边界频率为0.42π和0.49π，阻带边界频率为0.44π和0.47π。
带阻滤波器的幅频和相频特性如下：
4）带通滤波器：fc1=4800 Hz(0.44)，fc2=5200 Hz(0.47)，fp1=4600 Hz(0.42)，fp2=5400 Hz(0.49)，Rs=30dB，Rp=1dB。（[fp1 fp2]：阻带截至频率；[fc1 fc2]：通带截至频率）
选择Hanning窗，阻带边界频率为0.42π、0.49π，通带边界频率为0.44π、0.47π。
带通滤波器的幅频和相频特性如下：
*用完全设计函数设计IIR滤波器(fs=22050)
1）低通滤波器性能指标，fp=4500Hz(0.41)，fc=6500Hz(0.59)，Rs=100，Rp=1。
（fp：通带截至频率；fc：阻带截至频率；Rs：通带波纹；Rp：阻带波纹）低通滤波器的幅频和相频特性如下：
2）高通滤波器性能指标，fc=4500Hz(0.41)，fp=6500Hz(0.59)，Rs=100，Rp=1。（fp：通带截至频率；fc：阻带截至频率；Rs：通带波纹；Rp：阻带波纹）高通滤波器的幅频和相频特性如下：
3）带通滤波器性能指标，fc1=4800Hz(0.44)，fc2=5200Hz(0.47)，fp1=4600Hz(0.42)，fp2=5400Hz(0.49)，Rs=30dB，Rp=1dB。（[fp1 fp2]：阻带截至频率；[fc1 fc2]：通带截至频率）
带通滤波器的幅频和相频特性如下：
4）带阻滤波器性能指标，fp1=4800Hz(0.44)，fp2=5200Hz(0.47)，fc1=4600Hz(0.42)，fc2=5400Hz(0.49)，Rs=30dB，Rp=1dB。（[fp1 fp2]：阻带截至频率；[fc1 fc2]：通带截至频率）
带通滤波器的幅频和相频特性如下：
5、用滤波器对信号进行滤波
*用FIR低通滤波器：fp=4000Hz(0.36)，fc=4500(0.41)Hz，Rs=30dB，Rp=1 dB对“hello.wav”加噪声信号进行滤波。
结果图如下：
* 用IIR低通滤波器：fp=4000Hz(0.36)，fc=4500(0.41)Hz，Rs=30dB，Rp=1 dB对“hello.wav”加噪声信号进行滤波。
结果图如下：
六、思考题
1、IIR与FIR 设计方法的各自特点是什么？
IIR的特点：先按指标设计模拟滤波器，再转换为数字滤波器；
FIR的特点：直接按指标设计数字滤波器。
2、IIR与FIR各自优缺点是什么？
IIR：优点：相同的性能下阶次低；零极点可同时起作用；
缺点：相位是非线性的；不一定稳定；运算误差比较大，对频率分量的选择性不好。
FIR：优点：相位一定是线性的；系统一定是稳定的；运算误差比较小；对频率分量的选择性好。
缺点：阶次高。
3、为什么有这么多的设计方法？
为了满足各种方面的需要。
当滤波器类型简单，参数固定时宜采用窗函数或者脉冲响应不变法，已达到最为精确的滤波效果；当滤波器类型较复杂，有多个通带阻带时用窗函数法很繁琐，宜采用最优化设计方法进行快速设计；而且还要根据滤波器的不同类型选择方法，比如椭圆滤波器就和使用双线性变换法设计。
4、有没有一种滤波器在所有情况下都是最佳的？
没有。任何一种滤波器都有性能上的优缺点，满足了一方面的需求，就难以满足另一方面的需求，比如说阶次和线性相位的矛盾就是这样。我们在设计的时候应该根据具体情况的要求，优先满足可以达到较好总体效果的标准。
七、收获和总结
这是数字信号处理的第二次实验课，主要进行了语音信号的采样、滤波、插值、重构练习，学习并实践了IIR、FIR数字滤波器的设计。
语音信号的处理相对简单，在练习过程中，我在滤波器设计上花了比较多的时间。对FIR滤波器的原理和设计过程本来是比较清楚的，但在编程时出现了很多细小的错误，如：变量定义不正确，标点符号中英文的错误等等。应该在以后的编程中养成一种良好的习惯，争取一次正确，不要寄希望于在执行出错后再回头检查。IIR的设计用完全设计函数设计中
没有出现问题，但是在脉冲响应法中设计的低通和带通滤波器都没有问题，高通和带阻明显不正确，图如下：
源代码：
clear all;
fs=22050;
wp1=0.41*pi;ws1=0.59*pi;rp1=1,rs1=100;%数字滤波器截止频率、通带波纹和阻带衰减
T=1/fs;Nn=128;%采样间隔
Wp1=wp1/T;Ws1=ws1/T;%得到模拟滤波器的频率—采用脉冲响应不变法的频率转换形式
[N1,Wn1]=cheb1ord(Wp1,Ws1,rp1,rs1,'s');%计算模拟滤波器的最小阶数
[z1,p1,k1]=cheb1ap(N1,rp1);%设计低通原型数字滤波器
[Bap1,Aap1]=zp2tf(z1,p1,k1); %零点极点增益形式转换为传递函数形式
[b1,a1]=lp2lp(Bap1,Aap1,Wn1);%低通滤波器频率转换
[bz1,az1]=impinvar(b1,a1,1/T);%脉冲响应不变法设计数字滤波器传递函数
figure(1)
[H1,f1]=freqz(bz1,az1,Nn,1/T);%输出幅频响应和相频响应
subplot(2,1,1);plot(f1,20*log10(abs(H1)));
xlabel('频率/Hz');ylabel('振幅/dB');title('切比雪夫1型低通IIR');grid on;
subplot(2,1,2);plot(f1,180/pi*unwrap(angle(H1)));
xlabel('频率/Hz');ylabel('相位/^o');grid on;
wp=0.59*pi;ws=0.41*pi;rp=1,rs=100;%数字滤波器截止频率、通带波纹和阻带衰减T=1/fs;Nn=128;%采样间隔
Wp=wp/T;Ws=ws/T;%得到模拟滤波器的频率—采用脉冲响应不变法的频率转换形式[N,Wn]=cheb1ord(Wp,Ws,rp,rs,'s');%计算模拟滤波器的最小阶数
[z,p,k]=cheb1ap(N,rp);%设计高通原型数字滤波器
[Bap,Aap]=zp2tf(z,p,k); %零点极点增益形式转换为传递函数形式
[b,a]=lp2hp(Bap,Aap,Wn);%高通滤波器频率转换
[bz,az]=impinvar(b,a,1/T);%脉冲响应不变法设计数字滤波器传递函数
figure(2)
[H,f]=freqz(bz,az,Nn,1/T);%输出幅频响应和相频响应
subplot(2,1,1);plot(f,20*log10(abs(H)));
xlabel('频率/Hz');ylabel('振幅/dB');title('切比雪夫1型高通IIR');grid on;
subplot(2,1,2);plot(f,180/pi*unwrap(angle(H)));
xlabel('频率/Hz');ylabel('相位/^o');grid on;
wp21=0.42*pi;wp22=0.49*pi;ws21=0.44*pi;ws22=0.47*pi;rp2=1,rs2=30;%数字滤波器截止频率、通带波纹和阻带衰减
T=1/fs;Nn=128;%采样间隔
Wp21=wp21/T;Wp22=wp22/T;Ws21=ws21/T;Ws22=ws22/T;%得到模拟滤波器的频率—采用脉冲响应不变法的频率转换形式
Wo = sqrt(Ws21*Ws22);
Bw = Ws22-Ws21;
[N2,Wn2]=cheb1ord([Wp21 Wp22],[Ws21 Ws22],rp2,rs2,'s');%计算模拟滤波器的最小阶数
[z2,p2,k2]=cheb1ap(N2,rp2);%设计带通原型数字滤波器
[Bap2,Aap2]=zp2tf(z2,p2,k2); %零点极点增益形式转换为传递函数形式
[b2,a2]=lp2bp(Bap2,Aap2,Wo,Bw);%带通滤波器频率转换
[bz2,az2]=impinvar(b2,a2,1/T);%脉冲响应不变法设计数字滤波器传递函数
figure(3)
[H2,f2]=freqz(bz2,az2,Nn,1/T);%输出幅频响应和相频响应
subplot(2,1,1);plot(f2,20*log10(abs(H2)));
xlabel('频率/Hz');ylabel('振幅/dB');title('切比雪夫1型带通IIR');grid on;
subplot(2,1,2);plot(f2,180/pi*unwrap(angle(H2)));
xlabel('频率/Hz');ylabel('相位/^o');grid on;
wp31=0.44*pi;wp32=0.47*pi;ws31=0.42*pi;ws32=0.49*pi;rp3=1,rs3=30;%数字滤波器截止频率、通带波纹和阻带衰减
T=1/fs;Nn=128;%采样间隔
Wp31=wp31/T;Wp32=wp32/T;Ws31=ws31/T;Ws32=ws32/T;%得到模拟滤波器的频率—采用脉冲响应不变法的频率转换形式
Wo1 = sqrt(Wp31*Wp32);
Bw1 = Wp32-Wp31;
[N3,Wn3]=cheb1ord([Wp31 Wp32],[Ws31 Ws32],rp3,rs3,'s');%计算模拟滤波器的最小阶数
[z3,p3,k3]=cheb1ap(N3,rp3);%设计带阻原型数字滤波器
[Bap3,Aap3]=zp2tf(z3,p3,k3); %零点极点增益形式转换为传递函数形式
[b3,a3]=lp2bs(Bap3,Aap3,Wo1,Bw1);%带阻滤波器频率转换
[bz3,az3]=impinvar(b3,a3,1/T);%脉冲响应不变法设计数字滤波器传递函数
figure(4)
[H3,f3]=freqz(bz3,az3,Nn,1/T);%输出幅频响应和相频响应
subplot(2,1,1);plot(f3,20*log10(abs(H3)));
xlabel('频率/Hz');ylabel('振幅/dB');title('切比雪夫1型带阻IIR');grid on;
subplot(2,1,2);plot(f3,180/pi*unwrap(angle(H3)));
xlabel('频率/Hz');ylabel('相位/^o');grid on;
我分析了很久，看书、网络上查也没找到原因，还请老师检查指正。
信号与系统实验5
信号与系统实验（五） 班级：通信5班姓名：刘贺洋学号：11081515 班级：通信5班姓名：章仕波学号：11081522 1.符号函数的傅里叶变换 (1)下面参考程序和运行结果是信号||2 f- t =的傅里叶变换，分析程序，判 e ) (t 断运行结果正确与否。 syms t; %时间符号 f=exp(-2*abs(t)); %符号函数 F=fourier(exp(-2*abs(t))); subplot(1,2,1); ezplot(f); subplot(1,2,2); ezplot(F); 1（1）图
(2)参考上述程序试画出信号)(32 )(3t u e t f t -=的波形及其幅频特性曲线。 1（2）源程序： syms t ; %时间符号 f=2/3*exp(-3*t)*heaviside(t); %符号函数 F=fourier(f); subplot(1,2,1); ezplot(f); subplot(1,2,2); ezplot(abs(F)); 1（2）图： 2.符号函数的傅里叶变换
(1)下面参考程序是求信号211)(ωω+=j F 的逆傅里叶变换，分析程 序，比较运行结果。 源程序2（1） syms t w; F=1/(1+w^2); f=ifourier(F,t); subplot(1,2,1); ezplot(F); subplot(1,2,2); ezplot(f); 2（1）图： （2）求信号ωωωsin 2)(=j F 的逆傅里叶变换，并用程序验证。
源程序2（2） syms t w; F=2*sin(w)/w; f=ifourier(F,t); subplot(1,2,1); ezplot(F); subplot(1,2,2); ezplot(f); 图2（2）：
自动控制原理实验报告
第一章Matlab 基本运算 [范例1-2] 建立矩阵A={7 8 9}，B={7 8 9} >> A=[7,8,9] A = 7 8 9 >> B=A' B = 7 8 9 （2） >> B=[1 1 2 ; 3 5 8 ; 10 12 15] B= 1 1 2 3 5 8 10 12 15 （3） >> a=1:1:10 a = 1 2 3 4 5 6 7 8 9 10 >> t=10:-1:1
t = 10 9 8 7 6 5 4 3 2 1 [范例1-3]求多项式D（S）=（5S^2+3）(S+1)(S-1)的展开式 >> D=conv([5 0 3],conv([1 1],[1 -2])) D = 5 -5 -7 -3 -6 [范例1-4]求多项式P（X）=2X^4-5X^3-X+9 （1） >> P=[2 -5 6 -1 9] P = 2 -5 6 -1 9 >> x=roots(P) x = 1.6024 + 1.2709i 1.6024 - 1.2709i -0.3524 + 0.9755i -0.3524 - 0.9755i 第二章控制系统的数学模型 [范例2-1]已知系统传递函数G(S)= s + 3/ s^3 + 2 s^2 + 2 s + 1 >> num=[0 1 3]; >> den=[1 2 2 1]; >> printsys(num,den) num/den = s + 3 --------------------- s^3 + 2 s^2 + 2 s + 1 [范例2-2]已知系统传递函数G(S)=【5*（S+2）^2(S^2+6S+7)】/S(S+1)^3(S^3+2S+1)],试
北航实验报告实验实验
实验三UC-OS移植实验 一、实验目的 在内核移植了uCOS-II 的处理器上创建任务。 二、实验内容 1．运行实验十，在超级终端上观察四个任务的切换。 2. 任务1~3，每个控制“红”、“绿”、“蓝”一种颜色的显示，适当增加OSTimeDly()的时间，且优先级高的任务延时时间加长，以便看清三种颜色。 3.引入一个全局变量BOOLEAN ac_key，解决完整刷屏问题。 4. #define rUTRSTAT0 (*(volatile unsigned *)0x50000010) #define RdURXH0（）(*(volatile unsigned char *)0x50000024) 当键盘有输入时在超级终端上显示相应的字符。 三、实验设备 硬件：ARM嵌入式开发平台、用于ARM920T的JTAG仿真器、PC机Pentium100以上。 软件：PC机操作系统Win2000或WinXP、ARM ADS1.2集成开发环境、仿真器驱动程序、超级终端通讯程序。 四、实验原理 所谓移植，指的是一个操作系统可以在某个微处理器或者微控制器上运行。虽然uCOS-II的大部分源代码是用C语言写成的，仍需要用C语言和汇编语言完成一些与处理器相关的代码。比如:uCOS-II在读写处理器、寄存器时只能通过汇编语言来实现。因为uCOS-II在设计的时候就己经充分考虑了可移植性，所以，uCOS-II的移植还是比较容易的。 要使uCOS一工工可以正常工作，处理器必须满足以下要求: 1）处理器的C编译器能产生可重入代码。 2）在程序中可以打开或者关闭中断。 3）处理器支持中断，并A能产生定时中断(通常在10Hz}1000Hz之间)。 4）处理器支持能够容纳一定量数据的硬件堆栈。 5）处理器有将堆栈指针和其它CPU寄存器存储和读出到堆栈(或者内存)的指
2014信号与系统实验五(2)
Matlab拓展实验1：周期方波的傅立叶级数分析 一、周期方波信号的产生 1、周期方波函数：square(t,duty) （1）其中t为时间向量。 （2）duty为占空比，此参数可缺省，缺省时，duty=50。 注：周期方波的占空比为x%，是指一个周期内，高电平持续时间占整个周期时长的x%。但在square函数中，参数duty是一个0~100的数值。例如，若需产生占空比为20%的方波，则应设置duty=20。 2、功能描述：square函数产生一个周期为2π、高低电平分别为±1的周期方波，通过适当的编程可以调整为任意周期，任意幅度，任意中心值。 【例1】在时间范围[?4π,4π]产生周期的2π的周期方波。 clear;clf;clc; step=pi/200; %设置时间步长 t=-4*pi:step:4*pi; %设置时间范围 x1=square(t); %产生占空比为50%的周期方波 duty=20;x2=square(t,duty); %产生占空比为20%的周期方波 subplot(2,1,1);plot(t,x1); axis([-4*pi,4*pi,-1.5,1.5]);title('占空比为0.5的周期方波'); subplot(2,1,2);plot(t,x2); axis([-4*pi,4*pi,-1.5,1.5]);title('占空比为0.2的周期方波'); 上述程序的运行结果如下：
【思考】如将例1程序中的作图部分修改如下（修改部分以红色标出），这种处理方式称之为横坐标对pi归一化，图形结果会有何不同？这时横坐标的数值标注含义有何变化？ subplot(2,1,1);plot(t/pi,x1); axis([-4,4,-1.5,1.5]);title('占空比为0.5的周期方波'); subplot(2,1,2);plot(t/pi,x2); axis([-4,4,-1.5,1.5]);title('占空比为0.2的周期方波'); 【实践1】产生周期为1，占空比为30%，高电平为1，低电平为0的周期方波，时间范围为[?2,2]，即能观察到4个完整的周期。（时间变量的尺度变换） 二、傅立叶级数分解的计算 1、周期信号傅立叶级数展开的理论分析 周期信号的傅立叶级数有若干种形式（如三角形式、指数形式），这里以指数形式的傅立叶级数为例进行说明。 信号x(t)满足，x t=x(t+T)，其中T称为周期（或基波周期，单位s），F=1T称为频率（或基波频率，单位Hz），ω=2πT=2πF称为基波角频率，单位rad s。 （1）指数形式的傅立叶级数：
北航自动控制元件复习提纲
自控元件复习提纲 一、关于考试（从学长处获得信息，仅供大家参考）： 1、简答：eg：1）、为什么空载实验可以测r m ，x m？ 2）、变压器中是否必须有无功功率？ 2、论述：eg：1）、变压器工作原理。 2）、电机基本特点。 论述时必要时要作图说明。如要论述两相伺服电机为何有无自转特性时应该作出单相工作时两相伺服电机的机械特性曲线来说明。 3、计算：两个计算，一个直流，一个交流。 二、各章重点内容概述（参考往届笔记）： 第1章直流磁路及其计算 重点章节：1-1，1-2，1-3 磁路总是闭合的；磁路计算的正反两类任务；等效磁路的画法第2章直流电磁铁及其典型应用 重点章节：2-1，2-2，另外第三节中的继电器的主要技术指标大家也要注意第3章直流电机的一般问题 重点章节：3-1，3-2，3-4，3-5，3-6 发电机、电动机的识别，电角的概念，电枢反应的概念，直流电机的电枢电动 势和电磁转矩的计算方法，直流电机的电势平衡关系、转矩平衡关系、功率平 衡关系，电枢绕组的具体原理如果不明白大家不必深究，只要知道概念就可以， 这里不是重点。 PS：本章是比较重要的一章，计算题可能出在这一章。 第4章直流测速发电机和直流伺服电动机 重点章节：4-1，4-2 第三节要掌握直流力矩电动机的特点，知道其应用场合以 及为什么要用在这些场合。 准确理解直流伺服电动机的工作原理、四种工作状态，准确掌握直流测速发电 机和直流伺服电动机的工作原理及特性（输入、输出等） 第5章变压器 重点章节：5-2，5-3，5-4，5-5，5-6 空载电压平衡式和等效电路及相量图，负载运行的电压平衡式和磁动势平衡式 及相量图和等效电路P101的图5-12大家仔细看。参数、额定数据和特性。磁 场问题转化到电路问题。 第6章异步电动机 重点章节：6-1，6-3，6-4，6-5 转差率S，功率传递，相电动势，第三节可能考计算，注意一下习题中的最后 两道计算题6.12，6.15 另外，第六章和第五章联系很紧密，大家可以结合起来复习效果更好。 第7章两相电机 重点章节：7-1，7-2 分解磁场，无自转的分析，稳定运行范围大，第一节的第（四）部分其它大家 可以不用管。异步测速发电机的输出电压的特点，原理，频率与转速无关，只 与电源有关。 第8章同步电动机
北航计软实验报告一
北航计软实验报告一
————————————————————————————————作者: ————————————————————————————————日期：
计算机软件基础上机实验报告（一) ＸXXXＸX班XXX １.实验目的 掌握线性表在顺序分配下的插入与删除运算；掌握线性表的链式存储结构;掌握插入排序的方法;并掌握一种产生随机数的方法。 2.实验内容 1．产生１000个0至999间的随机整数，并以产生的次序存入一个数据文件中。 2．编制一个程序,依次实现以下功能： （1)?定义一个有序(非递减）线性表,其最大容量为１0０0,初始时为空。 (2)?从由１产生的数据文件中依次取前N个随机整数，陆续插入到此线性表中，并要求在每次插入后保持线性表的有序性。最后将此有序线性表打印输出。 (３)?在由（2）产生的线性表中，依在1中产生的次序逐个将元素删除，直至表空为止。 3. 以N＝１0０及N＝400分别运行2的程序,并比较它们的运行时间。 ４. 编写一个程序，用插入排序依次将1中产生的１000个随机整数链接成有序链表(不改变原随机数在存储空间中的顺序）。 ３.源代码与运行结果 #incｌｕdｅ<ｓｔdｉo.h> #ｉnclude<ｓｔｄliｂ.h＞ ＃include
数字信号处理(北航)实验二报告
数字信号处理实验二 信号的分析与处理综合实验 38152111 张艾一、实验目的 综合运用数字信号处理的理论知识进行信号的采样，重构，频谱分析和滤波器的设计，通过理论推导得出相应结论，再利用Matlab作为编程工具进行计算机实现，从而加深对所学知识的理解，建立概念。 二、基本要求 1．掌握数字信号处理的基本概念、基本理论和基本方法； 2．学会MATLAB的使用，掌握MATLAB的程序设计方法； 3．掌握用MATLAB设计简单实验验证采样定理的方法； 4．掌握在Windows环境下语音信号采集的方法； 5．学会用MATLAB对信号进行频谱分析； 6．掌握MATLAB设计FIR和IIR数字滤波器的方法； 三、实验内容 1．利用简单正弦信号设计实验验证采样定理： （1）Matlab产生离散信号的方法，作图的方法，以及基本运算操作 （2）对连续正弦信号以不同的采样频率作采样 （3）对采样前后信号进行傅立叶变换，并画频谱图 （4）分析采样前后频谱的有变化，验证采样定理。
掌握画频谱图的方法，深刻理解采样频率，信号频率，采样点数，频率分辨率等概念2．真实语音信号的采样重构：录制一段自己的语音信号，并对录制的信号进行采样；画出采样前后语音信号的时域波形和频谱图；对降采样后的信号进行插值重构，滤波，恢复原信号。 （1）语音信号的采集 （2）降采样的实现（改变了信号的采样率） （3）以不同采样率采样后，语音信号的频谱分析 （4）采样前后声音的变化 （5）对降采样后的信号进行插值重构，滤波，恢复原信号 3．带噪声语音信号的频谱分析 （1）设计一频率已知的噪声信号，与实验2中原始语音信号相加，构造带噪声信号（2）画出原始语音信号和加噪声后信号，以及它们的频谱图 （3）利用频谱图分析噪声信号和原语音信号的不同特性 4．对带噪声语音信号滤波去噪：给定滤波器性能指标，采样窗函数法或双线性变换设计滤波器，并画出滤波器的频率响应；然后用自己设计的滤波器对采样的语音信号进行滤波，画出滤波后信号的时域波形和频谱，并对滤波前后的信号进行对比，分析信号的变化； 回放语音信号； （1）分析带噪声信号频谱，找出噪声所在的频率段 （2）利用matlab中已有的滤波器滤波 （3）根据语音信号特点，自己设计滤波器滤波 （4）比较各种滤波器性能（至少四种），选择一种合适的滤波器将噪声信号滤除 （5）回放语音信号，比较滤波前后声音的变化
信号与系统-实验报告-实验五
实验五 连续信号与系统的S 域分析 学院 班级 姓名 学号 一、实验目的 1. 熟悉拉普拉斯变换的原理及性质 2. 熟悉常见信号的拉氏变换 3. 了解正/反拉氏变换的MATLAB 实现方法和利用MATLAB 绘制三维曲面图的方法 4. 了解信号的零极点分布对信号拉氏变换曲面图的影响及续信号的拉氏变换与傅氏变换的关系 二、 实验原理 拉普拉斯变换是分析连续时间信号的重要手段。对于当t ∞时信号的幅值不衰减的时间信号，即在f(t)不满足绝对可积的条件时，其傅里叶变换可能不存在，但此时可以用拉氏变换法来分析它们。连续时间信号f(t)的单边拉普拉斯变换F(s)的定义为： 拉氏反变换的定义为： 显然，上式中F(s)是复变量s 的复变函数，为了便于理解和分析F(s)随s 的变化规律，我们将F(s)写成模及相位的形式：()()()j s F s F s e ?=。其中，|F(s)|为复信号F(s)的模，而()s ?为F(s)的相位。由于复变量s=σ+jω，如果以σ为横坐标（实轴），jω为纵坐标（虚轴），这样，复变量s 就成为一个复平面，我们称之为s 平面。从三维几何空间的角度来看，|()|F s 和()s ?分别对应着复平面上的两个曲面，如果绘出它们的三维曲面图，就可以直观地分析连续信号的拉氏变换F(s)随复变量s 的变化情况，在MATLAB 语言中有专门对信号进行正反拉氏变换的函数，并且利用 MATLAB 的三维绘图功能很容易画出漂亮的三维曲面图。 ①在MATLAB 中实现拉氏变换的函数为： F=laplace( f ) 对f(t)进行拉氏变换，其结果为F(s) F=laplace (f,v) 对f(t)进行拉氏变换，其结果为F(v) F=laplace ( f,u,v) 对f(u)进行拉氏变换，其结果为F(v) ②拉氏反变换 f=ilaplace ( F ) 对F(s)进行拉氏反变换，其结果为f(t) f=ilaplace(F,u) 对F(w)进行拉氏反变换，其结果为f(u) f=ilaplace(F,v,u ) 对F(v)进行拉氏反变换，其结果为f(u) 注意： 在调用函数laplace( )及ilaplace( )之前，要用syms 命令对所有需要用到的变量（如t,u,v,w ）等进行说明，即要将这些变量说明成符号变量。对laplace( )中的f 及ilaplace( )中的F 也要用符号定义符sym 将其说明为符号表达式。具体方法参见第一部分第四章第三节。 例①：求出连续时间信号 ()sin()()f t t t ε=的拉氏变换式，并画出图形 求函数拉氏变换程序如下： syms t s %定义符号变量 ft=sym('sin(t)*Heaviside(t)'); %定义时间函数f(t)的表达式
电气工程专业考研专业课初试科目及复试内容汇总
《电气工程专业考研专业课初试科目及复试内容汇总》 自动化专业的考研方向 自动化专业方向很广，考的时候还分双控，模式，电力电机等等方面，你可以参看学校是否在这个方面有无国家重点实验室，是不是国家重点学科来比较。 1. 清华， 2.中科院， 3.上海交大， 4.浙大，5华工，北航，东南，东北大学，西安交大，哈 尔滨工业大学，中国科技大学，华北电力，天津大学，东南大学，华中科技，武汉大学天津大学自动化 一般说来，初试的分数是最重要的，特别是考外校。当然，你的动手能力也是很重要的，还有你的英语口语，考研复试都是要考虑的。例如上海交大的复试，双控353的复试线，有380的被刷下来，就是英语口语已经专业课不是很扎实的。考外校的话依据学校而定是否要找导师 动手能力强，参加电子设计大赛都是作为你考研复试的参考，还是好好的准备初试的考试吧，毕竟它是个门槛。 【电气工程及其自动化】 北京工业大学 421自动控制原理 复试：1、电子技术2、计算机原理 北京航空航天大学 [双控] 432控制理论综合或433控制工程综合 [检测] 433控制工程综合或436检测技术综合 [系统] 431自动控制原理或451材料力学或841概率与数理统计 [模式] （自动化学院）433控制工程综合或436检测技术综合、（宇航学院）423信息类专业综合或431自动控制原理或461计算机专业综合 [导航] （自动化学院）432控制理论综合或433控制工程综合、（宇航学院）431自动控制原理 复试：无笔试。1) 外语口语与听力考核；2) 专业基础理论与知识考核；3) 大学阶段学习成绩、科研活动以及工作业绩考核；4) 综合素质与能力考核 北京化工大学 440电路原理 复试：综合1(含自动控制原理和过程控制系统及工程)、综合2(含自动检测技术装置和传感器原理及应用)、综合3(含信号与系统和数字信号处理) 注：数学可选择301数学一或666数学（单） 北京交通大学 [双控/检测]404控制理论 [模式]405通信系统原理或409数字信号处理 复试： [电子信息工程学院双控]常微分方程 [机械与电子控制工程学院检测]综合复试（单片机、自动控制原理） [计算机与信息技术学院模式] 信号与系统或操作系统
北航自动控制原理实验报告(完整版)
自动控制原理实验报告 一、实验名称：一、二阶系统的电子模拟及时域响应的动态测试 二、实验目的 1、了解一、二阶系统阶跃响应及其性能指标与系统参数之间的关系 2、学习在电子模拟机上建立典型环节系统模型的方法 3、学习阶跃响应的测试方法 三、实验内容 1、建立一阶系统的电子模型，观测并记录在不同时间常数T时的响应曲线，测定过渡过程时间T s 2、建立二阶系统电子模型，观测并记录不同阻尼比的响应曲线，并测定超调量及过渡过程时间T s 四、实验原理及实验数据 一阶系统 系统传递函数： 由电路图可得，取则K=1，T分别取：0.25, 0.5, 1 T 0.25 0.50 1.00 R2 0.25MΩ0.5M Ω1MΩ C 1μ1μ1μ T S 实测0.7930 1.5160 3.1050 T S 理论0.7473 1.4962 2.9927 阶跃响应曲线图1.1 图1.2 图1.3 误差计算与分析 （1）当T=0.25时，误差==6.12%； （2）当T=0.5时，误差==1.32%； （3）当T=1时，误差==3.58% 误差分析：由于T决定响应参数，而,在实验中R、C的取值上可能存在一定误差，另外，导线的连接上也存在一些误差以及干扰，使实验结果与理论值之间存在一定误差。但是本实验误差在较小范围内，响应曲线也反映了预期要求，所以本实验基本得到了预期结果。 实验结果说明 由本实验结果可看出，一阶系统阶跃响应是单调上升的指数曲线，特征有T确定，T越小，过度过程进行得越快，系统的快速性越好。 二阶系统 图1.1 图1.2 图1.3
系统传递函数： 令 二阶系统模拟线路 0.25 0.50 1.00 R4 210.5 C2 111 实测45.8% 16.9% 0.6% 理论44.5% 16.3% 0% T S实测13.9860 5.4895 4.8480 T S理论14.0065 5.3066 4.8243 阶跃响应曲线图2.1 图2.2 图2.3 注：T s理论根据matlab命令[os,ts,tr]=stepspecs(time,output,output(end),5)得出，否则误差较大。 误差计算及分析 1）当ξ=0.25时，超调量的相对误差= 调节时间的相对误差= 2）当ξ=0.5时，超调量的相对误差==3.7% 调节时间的相对误差==3.4% 4）当ξ=1时，超调量的绝对误差= 调节时间的相对误差==3.46% 误差分析：由于本试验中，用的参量比较多，有R1,R2,R3,R4;C1,C2;在它们的取值的实际调节中不免出现一些误差，误差再累加，导致最终结果出现了比较大的误差，另外，此实验用的导线要多一点，干扰和导线的传到误差也给实验结果造成了一定误差。但是在观察响应曲线方面，这些误差并不影响，这些曲线仍旧体现了它们本身应具有的特点，通过比较它们完全能够了解阶跃响应及其性能指标与系统参数之间的关系，不影响预期的效果。 实验结果说明 由本实验可以看出，当ωn一定时，超调量随着ξ的增加而减小，直到ξ达到某个值时没有了超调；而调节时间随ξ的增大，先减小，直到ξ达到某个值后又增大了。 经理论计算可知，当ξ=0.707时，调节时间最短，而此时的超调量也小于5%，此时的ξ为最佳阻尼比。此实验的ξ分布在0.707两侧，体现了超调量和调节时间随ξ的变化而变化的过程，达到了预期的效果。 图2.2 图2.1 图2.3
信号与系统实验实验报告
信号与系统实验实验报 告 文档编制序号：[KKIDT-LLE0828-LLETD298-POI08]
实验五连续系统分析一、实验目的 深刻理解连续时间系统的系统函数在分析连续系统的时域特性、频域特性及稳定性中的重要作用及意义，掌握根据系统函数的零极点设计简单的滤波器的方法。掌握利用MATLAB分析连续系统的时域响应、频响特性和零极点的基本方法。 二、实验原理 MATLAB提供了许多可用于分析线性时不变连续系统的函数，主要包含有系统函数、系统时域响应、系统频域响应等分析函数。 三、实验内容 1.已知描述连续系统的微分方程为，输入，初始状态 ，计算该系统的响应，并与理论结果比较，列出系统响应分析的步骤。 实验代码： a=[1 10]; b=[2]; [A B C D]=tf2ss(b,a); sys=ss(A,B,C,D); t=0: :5; xt=t>0; sta=[1]; y=lsim(sys,xt,t,sta); subplot(3,1,1); plot(t,y); xlabel('t'); title('系统完全响应 y(t)'); subplot(3,1,2); plot(t,y,'-b'); hold on yt=4/5*exp(-10*t)+1/5; plot(t,yt,' : r'); legend('数值计算','理论计算'); hold off xlabel('t'); subplot(3, 1 ,3); k=y'-yt; plot(t,k); k(1) title('误差');
实验结果： 结果分析： 理论值 y(t)=0. 8*exp(-10t)*u(t)+ 程序运行出的结果与理论预期结果相差较大误差随时间增大而变小，初始值相差最大，而后两曲线基本吻合，表明该算法的系统响应在终值附近有很高的契合度，而在初值附近有较大的误差。 2.已知连续时间系统的系统函数为，求输入分别为，， 时，系统地输出，并与理论结果比较。 a=[1,3,2,0]; b=[4,1]; sys=tf(b,a); t=0: :5; x1=t>0; x2=(sin(t)).*(t>0); x3=(exp(-t)).*(t>0); y1=lsim(sys,x1,t); y2=lsim(sys,x2,t); y3=lsim(sys,x3,t); subplot(3,1,1); plot(t,y1); xlabel('t'); title('X(t)=u(t)'); subplot(3,1,2); plot(t,y2); xlabel('t'); title('X(t)=sint*u(t)'); subplot(3, 1 ,3); plot(t,y3); xlabel('t'); title('X(t)=exp(-t)u(t)'); 实验结果： 结果分析： a=[1,3,2,0]; b=[4,1]; sys=tf(b,a); t=0: :5; x1=t>0; x2=(sin(t)).*(t>0); x3=(exp(-t)).*(t>0); y1=lsim(sys,x1,t); y2=lsim(sys,x2,t); y3=lsim(sys,x3,t); subplot(3,1,1); plot(t,y1,'-b');
北航自动控制原理实验报告- 一、二阶系统的电子模拟及时域响应的动态测试
成绩 北京航空航天大学 自动控制原理实验报告 学院机械工程及自动化学院 专业方向机械工程及自动化 班级 学号 学生姓名刘帆 自动控制与测试教学实验中心
实验一 一、二阶系统的电子模拟及时域响应的动态测试 实验时间2014年11月15日 实验编号 同组同学 一、实验目的 1、 了解一、二阶系统阶跃响应及其性能指标与系统参数之间的关系。 2、 学习在电子模拟机上建立典型环节系统模型的方法。 3、 学习阶跃响应的测试方法。 二、实验内容 1、 建立一阶系统的电子模型,观测并记录在不同时间常数T 时的跃响应曲线,并测定其过渡过程时间T s 。 2、 建立二阶系统的电子模型,观测并记录在不同阻尼比ζ时的跃响应曲线,并测定其超调量σ%及过渡过程时间T s 。 三、实验原理 1、一阶系统阶跃响应性能指标的测试 系统的传递函数为:()s ()1 C s K R s Ts φ=+（）= 模拟运算电路如下图 : 其中2 1 R K R = ,2T R C =;在实验中，始终保持21,R R =即1K =，通过调节2R 和C 的不同取值，使得T 的值分别为0.2,0.51,1.0。记录实验数据，测量过度过程的性能指标，其中取正负5%误差带，按照经验公式取3s t T =
2、二阶系统阶跃响应性能指标的测试 系 统 传递函数为： 令ωn=1弧度/秒，则系统结构如下图： 二阶系统的 模拟电路图如下： 在实验过程中，取22321,1R C R C ==，则 442312R R C R ζ==，即42 12R C ζ=；在实验当中取123121,1R R R M C C F μ===Ω==，通过调整4R 取不同的值，使得ζ分别为0.25,0.5,0.707,1；记录所测得的实验数据以及其性能指标，取正负5%误差 带，其中当ζ<1时经验公式为2 1 3.5 %100%,s n e t ζσζω- -=?= ，当ζ=1时经验公式 为n 4.75 ts ω= 四、试验设备： 1、HHMN-1型电子模拟机一台。 2、PC 机一台。 3、数字万用表一块。 4、导线若干。
北航数字图象处理实验报告
数字图像处理实验报告 实验二图像变换实验 1．实验目的 学会对图像进行傅立叶等变换，在频谱上对图像进行分析，增进对图像频域上的感性认识，并用图像变换进行压缩。 2．实验内容 对Lena或cameraman图像进行傅立叶、离散余弦、哈达玛变换。在频域，对比他们的变换后系数矩阵的频谱情况，进一步，通过逆变换观察不同变换下的图像重建质量情况。 3. 实验要求 实验采用获取的图像，为灰度图像，该图像每象素由8比特表示。具体要求如下： （1）输入图像采用实验1所获取的图像（Lena、Cameraman）； （2）对图像进行傅立叶变换、获得变换后的系数矩阵； （3）将傅立叶变换后系数矩阵的频谱用图像输出，观察频谱； （4）通过设定门限，将系数矩阵中95%的（小值）系数置为0，对图像进行反变换，获得逆变换后图像； （5）观察逆变换后图像质量，并比较原始图像与逆变后的峰值信噪比（PSNR）。 （6）对输入图像进行离散余弦、哈达玛变换，重复步骤1-5； （7）比较三种变换的频谱情况、以及逆变换后图像的质量（PSNR）。 4. 实验结果 1. DFT的源程序及结果 J=imread('10021033.bmp'); P=fft2(J); for i=0:size(P,1)-1 for j=1:size(P,2) G(i*size(P,2)+j)=P(i+1,j); end end Q=sort(G); for i=1:size(Q,2) if (i 数字信号处理实验实验一离散时间信号与系统及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]; 北航《信号与系统》在线作业二 一、单选题（共10 道试题，共30 分。） 1. 信号〔ε(t)－ε(t－2)〕的拉氏变换的收敛域为________。 A. Re[s]>0 B. Re[s]>2 C. 全S平面 D. 不存在 -----------------选择：C 2. 信号的时宽与信号的频宽之间呈________。 A. 正比关系 B. 反比关系 C. 平方关系 D. 没有关系 -----------------选择：B 3. If f1(t) ←→F1(jω)，f2(t) ←→F2(jω) Then________。 A. [a f1(t) + b f2(t) ] ←→ [a F1(jω) *b F2(jω) ] B. [a f1(t) + b f2(t) ] ←→ [a F1(jω) - b F2(jω) ] C. [a f1(t) + b f2(t) ] ←→ [a F1(jω) + b F2(jω) ] D. [a f1(t) + b f2(t) ] ←→ [a F1(jω) /b F2(jω) ] -----------------选择：D 4. 某信号的频谱是周期的离散谱，则对应的时域信号为________。 A. 连续的周期信号 B. 连续的非周期信号 C. 离散的非周期信号 D. 离散的周期信号 -----------------选择：D 5. 信号在时域拥有的总能量，________其频谱在频域内能量的总和。 A. 大于 B. 等于 C. 小于 D. 不等于 -----------------选择：B 6. 理想低通滤波器是________。 A. 因果系统 B. 物理可实现系统 C. 非因果系统 D. 响应不超前于激励发生的系统 -----------------选择：C 7. 连续周期信号的傅氏变换是________。 A. 连续的 B. 周期性的 C. 离散的 D. 与单周期的相同 -----------------选择：C 电涡流传感器实验报告 38030414蔡达 一、实验目的 1.了解电涡流传感器原理; 2.了解不同被测材料对电涡流传感器的影响。 二、实验仪器 电涡流传感器实验模块，示波器:DS5062CE，微机电源:WD990型，士12V，万用表:VC9804A型，电源连接电缆，螺旋测微仪 三、实验原理 电涡流传感器由平面线圈和金属涡流片组成，当线圈中通以高频交变电流后，在与其平行的金属片上会感应产生电涡流，电涡流的大小影响线圈的阻抗Z，而涡流的大小与金属涡流片的电阻率、导磁率、厚度、温度以及与线圈的距离X有关，当平面线圈、被测体(涡流片)、激励源确定，并保持环境温度不变，阻抗Z只与距离X有关，将阻抗变化转为电压信号V输出，则输出电压是距离X的单值函数。 四． 实验数据及处理 1.铁片 0.5 1 1.52 2.5 3 3.5 电涡流传感器电压位移曲线—铁片 电压/V 位移/mm 0.5 1 1.5 2 2.53 3.5 电涡流传感器电压位移拟合曲线—铁片 电压/V 位移/mm 其线性工作区为0.6——3.4，对该段利用polyfit 进行函数拟合，可得V=-1.0488X-1.2465 2.铜片 电涡流传感器电压位移曲线—铜片 电压/V 位移/mm 2.2 2.4 2.6 2.83 3.2 3.4 3.6 -6-5.95-5.9-5.85 -5.8-5.75-5.7 -5.65-5.6-5.55-5.5电涡流传感器电压位移拟合曲线—铜片 电压/V 位移/mm 其线性工作区为2.4——3.4，对该段利用polyfit 进行函数拟合，可得V= -0.4500X -4.4667 实验5 连续时间系统的复频域分析 一、实验目的 1．掌握拉普拉斯变换及其反变换的定义，并掌握MATLAB实现方法。 2．学习和掌握连续时间系统函数的定义及复频域分析方法。 3．掌握系统零极点的定义，加深理解系统零极点分布与系统特性的关系。 二、实验原理与方法 1.拉普拉斯变换 连续时间信号x(t)的拉普拉斯变换、反变换定义为 MATLAB中，可以采用符号数学工具箱的laplace函数和ilaplace函数进行拉氏变换和拉氏反变换。 L=laplace(F)符号表达式F的拉氏变换，F中时间变量为t，返回变量为s的结果表达式。 L=laplace(F,t)用t替换结果中的变量t。 F=ilaplace(L)以s为变量的符号表达式L的拉氏反变换，返回时间变量为t的结果表达式。 F=ilaplace(L,x)用x替换结果中的变量t。 2.连续时间系统的系统函数 连续时间系统的系统函数是系统单位冲击响应的拉氏变换 连续时间系统的系统函数还可以由输入和输出信号的拉氏变换之比得到 单位冲击响应h(t)反映了系统的固有性质，而H(s)从复频域反映了系统的固有性质。由上式描述的连续时间系统，其系统系数为s的有理函数 3.连续时间系统的零极点分析 系统的零点指使上式的分子多项式为零的点，极点指使分母多项式为零的点，零点使系统的值为零，极点使系统的值为无穷大。通常将系统函数的零极点绘在s平面上，零点用○表示，极点用×表示，得到零极点分布图。 由零点定义可知，零点和极点分别指上式分子和分母多项式的根，利用MATLAB求多项式的根可以通过roots来实现，该函数调用格式为： r=roots(c) c为多项式的系数向量，返回值r为多项式的根向量。 分别对上式的分子、分母多项式求根即可求得零极点。 此外，在MATLAB中还提供了pzmap函数来求取零极点和绘制系统函数的零极点分布图，该函数的调用格式为： Pzmap(sys) 绘出系统模型sys描述的系统的零极点分布图。 [p,z]=pzmap(sys) 这种调用方法返回极点和零点，而不绘出零极点分布图。其中sys为系统传函模型，由t命令sys=tf(b,a)实现，b、a为传递函数的分子和分母多项式的系数向量。 MATLAB还提供了两个专用函数tf2zp和zp2tf来实现系统传递函数模型和灵机电网增益模型的转换，其调用格式为 [z,p,k]=th2zp(b,a) [b,a]=th2zp(z,p,k) 其中b、a为传递函数的分子和分母多项式的系数向量，返回值z为零点列向量，p为极点列向量，k为系统函数零极点形式的增益。 三、实验内容 1.已知系统的冲激响应h(t)u(t)u(t2) =，试采用复频域的 =--，输入信号x(t)u(t) 方法求解系统的响应，编写MATLAB程序实现。 解：MATLAB程序如下： h=sym('heaviside(t)-heaviside(t-2)'); H=laplace(f) x=sym('heaviside(t)'); X=laplace(x) Y=H*X 自动控制原理 实验报告 实验一一、二阶系统的电子模拟及时域响应的动态测试 实验目的 1.了解一、二阶系统阶跃响应及其性能指标与系统参数之间的关系。 2.学习在电子模拟机上建立典型环节系统模型的方法。 3.学习阶跃响应的测试方法。 二、实验内容 1.立一阶系统的电子模型,观测并记录在不同时间常数T时的跃响应曲线, 并测定其过渡过程时间TS。 2.立二阶系统的电子模型,观测并记录在不同阻尼比ζ时的跃响应曲线, 并测定其超调量σ%及过渡过程时间TS。 三、实验原理 1．一阶系统： 系统传递函数为：错误！未找到引用源。 模拟运算电路如图1-1所示： 图1-1 由图得： 在实验当中始终取错误！未找到引用源。, 则错误！未找到引用源。, 错误！未找到引用源。 取不同的时间常数T分别为： 0.25、 0.5、1。 记录不同时间常数下阶跃响应曲线，测量纪录其过渡过程时 ts。（取错误！ 未找到引用源。误差带） 2．二阶系统: 其传递函数为： 错误！未找到引用源。 令错误！未找到引用源。，则系统结构如图1-２所示： 图1-2 根据结构图，建立的二阶系统模拟线路如图1-３所示： 图1-3 取错误！未找到引用源。，错误！未找到引用源。，则错误！未找到引用源。及错误！未找到引用源。 错误！未找到引用源。取不同的值错误！未找到引用源。 , 错误！未找到引用源。, ,观察并记录阶跃响应曲线，测量超调量σ%（取错误！未找到引用源。误差带）,计算过渡过程时间Ts。 四、实验设备 1.HHMN-1型电子模拟机一台。 2.PC 机一台。 3.数字式万用表一块。 4.导线若干。 五、实验步骤 1.熟悉HHMN-1型电子模拟机的使用方法，将各运算放大器接成比例器，通电调零。 2.断开电源，按照实验说明书上的条件和要求，计算电阻和电容的取值，按照模拟线路图搭接线路，不用的运算放大器接成比例器。 3.将D/A1与系统输入端Ui连接，将A/D1与系统输出端UO连接（此处连接必须谨慎，不可接错）。线路接好后，经教师检查后再通电。 4.在Windows XP桌面用鼠标双击MATLAB图标后进入，在命令行处键入autolab 进入实验软件系统。 5.在系统菜单中选择实验项目，选择实验一，在窗口左侧选择实验模型，其它步骤察看概述3.2节内容。 6.观测实验结果，记录实验数据，绘制实验结果图形，填写实验数据表格，完成实验报告。 7.研究性实验方法。实验者可自行确定典型环节传递函数，并建立系统的SIMULINK模型，验证自动控制理论相关的理论知识。实现步骤可察看概述3.3节内容。 目录 实验一 (2) 实验二 (8) 实验三 (14) 实验四 (27) 实验一 实验目的：熟悉硬件开发流程，掌握Modelsim设计与仿真环境，学会简单组合逻辑电路、简单时序逻辑电路设计，不要求掌握综合和综合后仿真。 实验内容：必做实验：练习一、简单的组合逻辑设计 练习二、简单分频时序逻辑电路的设计 选做实验：选做一、练习一的练习题 选做二、7段数码管译码电路 练习一、简单的组合逻辑设计 描述一个可综合的数据比较器，比较数据a 、b的大小，若相同，则给出结果1，否则给出结果0。 实验代码： 模块源代码： module compare(equal,a,b); input a,b; output equal; assign equal=(a==b)?1:0; endmodule 测试模块源代码： `timescale 1ns/1ns `include "./compare.v" module t; reg a,b; wire equal; initial begin a=0; b=0; #100 a=0;b=1; #100 a=1;b=1; #100 a=1;b=0; #100 a=0;b=0; #100 $stop; end compare m(.equal(equal),.a(a),.b(b)); endmodule 实验波形 练习二、简单分频时序逻辑电路的设计 用always块和@(posedge clk)或@(negedge clk)的结构表述一个1/2分频器的可综合模型，观察时序仿真结果。 实验代码： 模块源代码： module halfclk(reset,clkin,clkout); input clkin,reset; output clkout; reg clkout; always@(posedge clkin) begin if(!reset) clkout=0; else clkout=~clkout; end endmodule 测试模块源代码： `timescale 1ns/100ps `define clkcycle 50 module tt; reg clkin,reset; wire clkout; always#`clkcycle clkin=~clkin; initial begin clkin=0; reset=1; #10 reset=0; #110 reset=1; #100000 $stop; end halfclk m0(.reset(reset),.clkin(clkin),.clkout(clkout)); endmodule数字信号处理实验及参考程序
15秋北航《信号与系统》在线作业二100分答案
北航电涡流传感器实验报告
北理工信号与系统实验(5)
自动控制原理实验报告
北航verilog实验报告(全)