哈工大威海信号系统实验报告完整版
- 格式:docx
- 大小:4.50 MB
- 文档页数:27
《信号与系统》实验报告(完整版)长江大学电工电子实验中心电路与系统(2)实验报告姓名高文昌班级电信10909班序号06指导教师黄金平老师成绩实验名称:连续信号的绘制一、实验目的1.掌握用Matlab 绘制波形图的方法,学会常见波形的绘制。
2.掌握用Matlab 编写函数的方法。
3.周期信号与非周期信号的观察。
加深对周期信号的理解。
二、实验内容1、用MATLAB 画出下列信号的波形。
(a) ][cos )(1t t f ε=; (b) )]2()2([2||)(2--+=t t t t f εε; (c) )]2()([sin )(3t t t t f ---=εεπ; (d) )sgn()()(24t t G t f =; (e) )2()(265-=t Q t G f ; (f) )sin(|)|2()(6t t t f πε-= (a )t=linspace(-10,10,400);f1=u(cos(t));figure(1),myplot(t,f1)xlabel('Time(sec)'),ylabel('f1(t)')(b)t=linspace(-4,4,400);f2=abs(t)/2.*(u(t+2)-u(t-2)); figure(2),myplot(t,f2)xlabel('Time(sec)'),ylabel('f2(t)');(c)t=linspace(-1,3,400);f3=sin(pi*t).*(u(-t)-u(2-t)); figure(3),myplot(t,f3)xlabel('Time(sec)'),ylabel('f3(t)')(d)t=linspace(-2,2,400); f4=sign(t).*rectpuls(t,2); figure(4),myplot(t,f4)xlabel('Time(sec)'),ylabel('f3(t)')(e)t=linspace(-1,4,400);f5=rectpuls(t,6).*tripuls(t-2,4); figure(5),myplot(t,f5)xlabel('Time(sec)'),ylabel('f5(t)')(f)t=linspace(-4,4,400); f6=u(2-abs(t)).*sin(pi*t) figure(6),myplot(t,f6)xlabel('Time(sec)'),ylabel('f6(t)')2、用基本信号画出图2.1-10中的信号。
哈尔滨⼯业⼤学威海校区_《数字信号处理》实验⼀数字信号处理实验报告实验名称:实验⼀离散傅⾥叶变换的性质实验⽇期:2011.11.16姓名:尤伟学号:090240328哈尔滨⼯业⼤学(威海)实验⼀离散傅⾥叶变换的性质⼀、实验⽬的1、掌握离散傅⾥叶变换的性质,包括线性特性、时移特性、频移特性、对称性和循环卷积等性质;2、通过编程验证傅⾥叶变换的性质,加强对傅⾥叶变换性质的认识。
⼆、实验原理和⽅法1.线性特性1212D FT [()()]()()ax n bx n aX k bX k +=+ 2.时移特性DFT[()]()DFT[()]()km kmx n m W X k x n m WX k -+=-=3.频移特性()()nlN IDFT X k l IDFT X k W +=4. 对称性设由x(n)开拓成的周期序列为 ()p x n 则()()()p pe po x n x n x n =+ 偶序列()()()*12pe p p x n x n x N n ??=+-?奇序列()()()*12pop p x n x n x N n ??=--?? 将()pe x n 和()po x n 截取主周期,分别得()()()pet pe N x n x n R n = ()()()p o tp oN x n x n R n =则()()()()()p N pet pot x n x n R n x n x n ==+ x(n)序列的实部和虚部的离散⽴叶变换(){}()R e petD FT x n X k = (){}()Im potj x n Xk =[][]()()()()()()()()()()()arg ()arg ()R R R I I I X k X k X N k X k X k X N k X k X k X N k X k X N k X k X k * =-=-=-=-=--=--=-=-- 5.循环卷积()3123121()()()()()x n x n x n X k X k X k N=?=有限长序列线性卷积与循环卷积的关系 X1(n)和x2(n)的线性卷积:11312120()()()()()N m m x n x m x n m x m x n m -∞=-∞==-=-∑∑112()()N m x m xn m -==-∑将X1(n)和x2(n)开拓成以N 为周期的周期序列11()()p r x n x n rN ∞=-∞=+∑22()()p q x n x n qN ∞=-∞=+∑则它们的周期卷积为14120()()()N p p p m x n xm x n m -==-∑12()()N p m x m xn m -==-∑1120()()N m q x m x n m qN -∞==-∞=-+∑∑1120()()N q m x m x n qN m ∞-=-∞=??=+-∑∑ 3()q x n qN ∞=-∞=+∑X1(n)和x2(n)周期开拓后的周期卷积等于他们的线性卷积的的周期开拓。
信号与系统实验实验报告一、实验目的本次信号与系统实验的主要目的是通过实际操作和观察,深入理解信号与系统的基本概念、原理和分析方法。
具体而言,包括以下几个方面:1、掌握常见信号的产生和表示方法,如正弦信号、方波信号、脉冲信号等。
2、熟悉线性时不变系统的特性,如叠加性、时不变性等,并通过实验进行验证。
3、学会使用基本的信号处理工具和仪器,如示波器、信号发生器等,进行信号的观测和分析。
4、理解卷积运算在信号处理中的作用,并通过实验计算和观察卷积结果。
二、实验设备1、信号发生器:用于产生各种类型的信号,如正弦波、方波、脉冲等。
2、示波器:用于观测输入和输出信号的波形、幅度、频率等参数。
3、计算机及相关软件:用于进行数据处理和分析。
三、实验原理1、信号的分类信号可以分为连续时间信号和离散时间信号。
连续时间信号在时间上是连续的,其数学表示通常为函数形式;离散时间信号在时间上是离散的,通常用序列来表示。
常见的信号类型包括正弦信号、方波信号、脉冲信号等。
2、线性时不变系统线性时不变系统具有叠加性和时不变性。
叠加性意味着多个输入信号的线性组合产生的输出等于各个输入单独作用产生的输出的线性组合;时不变性表示系统的特性不随时间变化,即输入信号的时移对应输出信号的相同时移。
3、卷积运算卷积是信号处理中一种重要的运算,用于描述线性时不变系统对输入信号的作用。
对于两个信号 f(t) 和 g(t),它们的卷积定义为:\(f g)(t) =\int_{\infty}^{\infty} f(\tau) g(t \tau) d\tau \在离散时间情况下,卷积运算为:\(f g)n =\sum_{m =\infty}^{\infty} fm gn m \四、实验内容及步骤实验一:常见信号的产生与观测1、连接信号发生器和示波器。
2、设置信号发生器分别产生正弦波、方波和脉冲信号,调整频率、幅度和占空比等参数。
3、在示波器上观察并记录不同信号的波形、频率和幅度。
信号与系统实验报告一、实验目的(1) 理解周期信号的傅里叶分解,掌握傅里叶系数的计算方法;(2)深刻理解和掌握非周期信号的傅里叶变换及其计算方法;(3) 熟悉傅里叶变换的性质,并能应用其性质实现信号的幅度调制;(4) 理解连续时间系统的频域分析原理和方法,掌握连续系统的频率响应求解方法,并画出相应的幅频、相频响应曲线。
二、实验原理、原理图及电路图(1) 周期信号的傅里叶分解设有连续时间周期信号()f t ,它的周期为T ,角频率22fT,且满足狄里赫利条件,则该周期信号可以展开成傅里叶级数,即可表示为一系列不同频率的正弦或复指数信号之和。
傅里叶级数有三角形式和指数形式两种。
1)三角形式的傅里叶级数:01212011()cos()cos(2)sin()sin(2)2cos()sin()2n n n n a f t a t a t b t b t a a n t b n t 式中系数n a ,n b 称为傅里叶系数,可由下式求得:222222()cos(),()sin()T T T T nna f t n t dtb f t n t dtTT2)指数形式的傅里叶级数:()jn tn nf t F e式中系数n F 称为傅里叶复系数,可由下式求得:221()T jn tT nF f t edtT周期信号的傅里叶分解用Matlab进行计算时,本质上是对信号进行数值积分运算。
Matlab中进行数值积分运算的函数有quad函数和int函数。
其中int函数主要用于符号运算,而quad函数(包括quad8,quadl)可以直接对信号进行积分运算。
因此利用Matlab进行周期信号的傅里叶分解可以直接对信号进行运算,也可以采用符号运算方法。
quadl函数(quad系)的调用形式为:y=quadl(‘func’,a,b)或y=quadl(@myfun,a,b)。
其中func是一个字符串,表示被积函数的.m文件名(函数名);a、b分别表示定积分的下限和上限。
数字信号处理实验报告实验题目:实验一用FFT做谱分析院系:班级:姓名:学号:哈尔滨工业大学实验一: 用FFT 作谱分析一、实验目的(1) 进一步加深DFT 算法原理和基本性质的理解(因为FFT 只是DFT 的一种快速算法, 所以FFT 的运算结果必然满足DFT 的基本性质)。
(2) 熟悉FFT 算法原理和FFT 子程序的应用。
(3) 学习用FFT 对连续信号和时域离散信号进行谱分析的方法,了解可能出现的分析误差及其原因,以便在实际中正确应用FFT 。
二、实验原理1.DFT 的运算量:一次完整的DFT 运算总共需要2N 次复数乘法和(1)N N -复数加法运算,因而直接计算DFT 时,乘法次数和加法次数都和2N 成正比,当N 很大时,运算量很客观的。
例如,当N=8时,DFT 运算需64位复数乘法,当N=1024时,DFT 运算需1048576次复数乘法。
而N 的取值可能会很大,因而寻找运算量的途径是很必要的。
2.FFT 算法原理:大多数减少离散傅里叶变换运算次数的方法都是基于nkN W 的对称性和周期性。
(1)对称性()*()k N n kn knN N NW W W --== (2)周期性()(mod`)()()kn N kn n N k n k N N N N NW W W W ++===由此可得()()/2(/2)1n N k N n k nk N N N N N k N k N N W W W W W W ---+⎧==⎪=-⎨⎪=-⎩这样:1.利用第三个方程的这些特性,DFT 运算中有些项可以合并;2.利用nkN W 的对称性和周期性,可以将长序列的DFT 分解为短序列的DFT 。
前面已经说过,DFT 的运算量是与2N 成正比的,所以N 越小对计算越有利,因而小点数序列的DFT 比大点数序列的DFT 运算量要小。
快速傅里叶变换算法正是基于这样的基本思路而发展起来的,她的算法基本上可分成两大类,即按时间抽取法和按频率抽取法。
信号与系统实验报告
实验名称:信号与系统实验
一、实验目的:
1.了解信号与系统的基本概念
2.掌握信号的时域和频域表示方法
3.熟悉常见信号的特性及其对系统的影响
二、实验内容:
1.利用函数发生器产生不同频率的正弦信号,并通过示波器观察其时域和频域表示。
2.通过软件工具绘制不同信号的时域和频域图像。
3.利用滤波器对正弦信号进行滤波操作,并通过示波器观察滤波前后信号的变化。
三、实验结果分析:
1.通过实验仪器观察正弦信号的时域表示,可以看出信号的振幅、频率和相位信息。
2.通过实验仪器观察正弦信号的频域表示,可以看出信号的频率成分和幅度。
3.利用软件工具绘制信号的时域和频域图像,可以更直观地分析信号的特性。
4.经过滤波器处理的信号,可以通过示波器观察到滤波前后的信号波形和频谱的差异。
四、实验总结:
通过本次实验,我对信号与系统的概念有了更深入的理解,掌
握了信号的时域和频域表示方法。
通过观察实验仪器和绘制图像,我能够分析信号的特性及其对系统的影响。
此外,通过滤波器的处理,我也了解了滤波对信号的影响。
通过实验,我对信号与系统的理论知识有了更加直观的了解和应用。
《随机信号分析》实验一班级学号姓名实验一实验内容:1 . 熟悉并练习使用下列Matlab的函数,给出各个函数的功能说明和内部参数的意义,并给出至少一个使用例子和运行结果:(1)randn()产生随机数数组或矩阵,其元素服从均值为0,方差为1的正态分布(1)Y = randn 产生一个伪随机数(2)Y = randn(n) 产生n×n的矩阵,其元素服从均值为0,方差为1的正态分布(3)Y = randn(m,n) 产生m×n的矩阵,其元素服从均值为0,方差为1的正态分布(4)Y= randn([m n]) 产生m×n的矩阵,其元素服从均值为0,方差为1的正态分布例:以(2)为例Y = randn(4)结果为:Y =-0.1941 -1.0722 -1.9609 0.8252-2.1384 0.9610 -0.1977 1.3790-0.8396 0.1240 -1.2078 -1.05821.3546 1.43672.9080 -0.4686(2)rand()(1)Y = rand(n) 生成n×n 随机矩阵,其元素在(0,1)例:以(2)为例Y = rand(3,4)内(2)Y = rand(m,n) 生成m×n 随机矩阵(3)Y = rand([m n]) 生成m×n 随机矩阵(4)Y = rand(m,n,p,…) 生成m×n×p×…随机矩阵或数组(5)Y = rand([m n p…]) 生成m×n×p×…随机矩阵或数组(6)Y = rand(size(A)) 生成与矩阵A 相同大小的随机矩阵结果为:Y =0.5797 0.8530 0.5132 0.23990.5499 0.6221 0.4018 0.12330.1450 0.3510 0.0760 0.1839(3)normrnd()产生服从正态分布的随机数(1)Y = normrnd(mu,sigma) 产生服从均值为mu,标准差为sigma的随机数,mu和sigma可以为向量、矩阵、或多维数组。
哈尔滨工业大学(威海)操作系统实验报告说明:本实验报告实验答案,是本人在上实验时的测试数据,由于操作系统实验中后面实验与当时所做实验的计算机的配置有关,因此本实验报的数据仅供参考。
实验1进程的描述与控制Windows 2000编程(实验估计时间:100分钟)1.1 背景知识Windows 2000 可以识别的应用程序包括控制台应用程序、GUI应用程序和服务应用程序。
控制台应用程序可以创建GUI,GUI应用程序可以作为服务来运行,服务也可以向标准的输出流写入数据。
不同类型应用程序间的惟一重要区别是其启动方法。
Windows 2000是以NT技术构建的,它提供了创建控制台应用程序的能力,使用户可以利用标准的C++工具,如iostream库中的cout和cin对象,来创建小型应用程序。
当系统运行时,Windows 2000的服务通常要向系统用户提供所需功能。
服务应用程序类型需要ServiceMail()函数,由服务控制管理器(SCM)加以调用。
SCM是操作系统的集成部分,负责响应系统启动以开始服务、指导用户控制或从另一个服务中来的请求。
其本身负责使应用程序的行为像一个服务,通常,服务登录到特殊的LocalSystem账号下,此账号具有与开发人员创建的服务不同的权限。
当C++编译器创建可执行程序时,编译器将源代码编译成OBJ文件,然后将其与标准库相链接。
产生的EXE文件是装载器指令、机器指令和应用程序的数据的集合。
装载器指令告诉系统从哪里装载机器代码。
另一个装载器指令告诉系统从哪里开始执行进程的主线程。
在进行某些设置后,进入开发者提供的main()、Servicemain()或WinMain()函数的低级入口点。
机器代码中包括控制逻辑,它所做的事包括跳转到Windows API函数,进行计算或向磁盘写入数据等。
Windows允许开发人员将大型应用程序分为较小的、互相有关系的服务模块,即动态链接库(DLL)代码块,在其中包含应用程序所使用的机器代码和应用程序的数据。
《信号与系统》实验报告姓名:学号:同组人:无指导教师:成绩:实验一典型连续时间信号描述及运算实验报告要求:(1)仿照单边指数信号的示例程序,按要求完成三种典型连续信号,即:正弦信号、衰减正弦信号、钟型信号的波形绘制。
(要求:要附上程序代码,以下均如此,不再说明)(2)根据《信号与系统》教材第一章的习题(1,3,5,8)函数形式绘制波形。
(3)完成三种奇异信号,即:符号函数、阶跃信号、单位冲激信号的波形绘制。
(4)完成实验一中信号的运算:三、6 实验内容中的(1)(2)(3)(4)。
(5)求解信号的直流/交流分量,按第四部分的要求完成。
正文:(1)<1>正弦信号:代码:>> t=-250:1:250;>> f1=150*sin(2*pi*t/100);>> f2=150*sin(2*pi*t/200);>> f3=150*sin(2*pi*t/200+pi/5);>> plot(t,f1,'-',t,f2,'--',t,f3,'-.')<2>衰减正弦信号<3>代码:>> t=-250:1:250;>> f1=400*exp(-1.*t.*t./10000);>> f1=400*exp(-1.*t.*t./22500);>> f1=400*exp(-1.*t.*t./62500);>> plot(t,f1,'-',t,f2,'--',t,f3,'-.')(2)习题1,3,5,8<1>代码:t=0:1:10;f=t;plot(t,f)<3>代码:t=1:1:10;f=t;plot(t,f)<5>代码:t=0:1:10;f=2-exp(-1.*t.);plot(t,f)<8>代码:t=1::2;f=exp(-1.*t.)*cos(10*pi*t);plot(t,f)(3)三种奇异函数<1>符号函数代码: t=-5::5;f=sign(t);plot(t,f)<2>阶跃信号代码:>> t=-5::5;>> f=u(t);>> plot(t,f)<3>单位冲激信号代码:function chongji(t1,t2,t0)dt=;t=t1:dt:t2;n=length(t);x=zeros(1,n);x(1,(-t0-t1)/dt+1)=1/dt;stairs(t,x);axis([t1,t2,0,dt]) title('单位冲激信号δ(t) ')(4)实验三1234<1>syms tf1=sym('(-t+4)*(u(t)-u(t-4))'); subplot(1,2,1);ezplot(f1);y1=subs(f1,t,-t);f3=f1+y1;subplot(1,2,2);ezplot(f3);function f=u(t) f=(t>0);<2>4、function f=u(t)f=(t>0)syms tf1=sym('(-t+4)*(u(t)-u(t-4))'); subplot(1,3,1);ezplot(f1);f2=sym('sin(2*pi*t)'); subplot(1,3,2);ezplot(f2);f6=f1.*f2;subplot(1,3,3);ezplot(f6);5、function f=u(t)f=(t>0)syms tf1=sym('(-t+4)*(u(t)-u(t-4))'); f2=sym('sin(2*pi*t)'); subplot(1,3,1);ezplot(f2);f6=f1.*f2;y6=subs(f6,t,t-2);subplot(1,3,2);ezplot(y6);f7=y6+f2;subplot(1,3,3);ezplot(f7);四、t=0::500;f=100.*abs(sin(2.*pi.*t./50));plot(t,f,t,fD,t,fA)调用子程序:function fD=fDC(f)fD=mean(f);function fA=fAC(f,fD)fA=f-fD;(5)求解信号的交直流分量代码:function fD=fDC(f)fD=mean(f);function fA=fAC(f,fD)fA=f-fD;t=0::500;f(t)=100|sin(2*PI*t/50)|;plot(t,fD,t,fA)实验二线性系统时域分析实验报告要求:(1)求解下面两个信号的卷积积分。
要求:1) 在实验报告中推导出这两个信号卷积积分运算表达式;(手写)2) 利用MA TLAB 进行求解验证,附程序代码和波形。
(2)已知描述系统的微分方程和激励信号如下r??(t ) ??3r?(t ) ?2r(t ) ??e?(t)+3 e(t) ,e(t)=u(t) 。
要求:1) 用解析法求系统的零状态响应r(t );(手写)2)利用MATLAB绘出系统零状态响应的时域仿真波形,并验证1)的结果是否正确,附程序代码和波形;3)利用MA TLAB绘出系统的冲激响应和阶跃响应波形,附程序代码和波形。
正文:(1)求解下面两个信号的卷积积分。
<1>在实验报告中推导出这两个信号卷积积分运算表达式;(手写)<2>利用MATLAB 进行求解验证,附程序代码和波形。
代码:p=;k1=-1/2:p:1;f1=1;k2=0:p:2;f2=*k2;[f,k]=sconv(f1,f2,k1,k2,p)function [f,k]=sconv(f1,f2,k1,k2,p)f=conv(f1,f2);f=f*p;k0=k1(1)+k2(1);k3=length(f1)+length(f2)-2;k=k0:p:(k3*p+k0);subplot(2,2,1)plot(k1,f1)xlabel('t')ylabel('f1(t)')subplot(2,2,2)plot(k2,f2)ylabel('f2(t)')subplot(2,2,3)plot(k,f);h=get(gca,'position');h(3)=*h(3);set(gca,'position',h)title('f(t)=f1(t)*f2(t)')xlabel('t')ylabel('f(t)')(2)已知描述系统的微分方程和激励信号如下r??(t ) ??3r?(t ) ?2r(t ) ??e?(t)+3 e(t) ,e(t)=u(t) 。
要求:<1> 用解析法求系统的零状态响应r(t );(手写)<2>利用MATLAB绘出系统零状态响应的时域仿真波形,并验证<1>的结果是否正确,附程序代码和波形;代码:>> a=[1,3,2];>> b=[0,1,3];>> impulse(a,b)<3>利用MATLAB绘出系统的冲激响应和阶跃响应波形,附程序代码和波形。
代码:a=[1,3,2];b=[0,1,3];step(b,a)代码;a=[1,3,2];b=[0,1,3];>> p=;>> t=0:p:5;>> x=exp(-t);>> lsim(b,a,x,t)实验三:实验报告内容:对所给音频信号,进行时域压缩和扩展,画出时域波形与幅度谱,使其满足以下要求。
(1)将的音调变低a倍(<a<1),及变高a倍(1<a<),比较变换前后的10%带宽(最大值10%处的带宽)的变化情况。
(2)将的音调变化到与儿童和男声相似的声音,分析变换前后的10%带宽的变化情况(变宽或变窄多少)。
(1)①a=[ ];a1=30;b=11025;[x0 b]=wavread('F:\信号与系统\');for m=1:3;a2=a1*a(m);x=resample(x0,a1,a2) ;wavplay(x ,b);t=0:1/b:(length(x)-1)/b;subplot(3,3,3*m-2);plot(t,x);axis([0 6 -inf inf]);title(['a=',num2str(a(m)),'的时域图']);N=(length(resample(x0,a1,a1*max(a))));f=(-N/2:N/2-1)/N*b;X= fftshift(fft(x,N));AMP=abs(X);AMP0=AMP/max(AMP);subplot(3,3,3*m-1);plot(f,AMP0);title(['a=',num2str(a(m)),'的频域图']);subplot(3,3,3*m);plot(f,AMP0);axis([-inf inf 0 ])title(['a=',num2str(a(m)),'的10%带宽图']);end②a=[1 ];str=['原' '男' ‘童'];a1=30;b=11025;[x0 b]=wavread('F:\信号实验_2018\');for m=1:3;a2=a1*a(m);x=resample(x0,a1,a2);wavplay(x,b);t=0:1/b:(length(x)-1)/b;subplot(3,3,3*m-2); plot(t,x);axis([0 3 -inf inf])title([str(m),'声的时域图 ']);N=(length(resample(x0,a1,a1*max(a)))); f=(-N/2:N/2-1)/N*b;X=fftshift(fft(x,N));AMP=abs(X); AMP0=AMP/max(AMP);subplot(3,3,3*m-1); plot(f,AMP0);title([str(m),'声的频域图 ']);subplot(3,3,3*m); plot(f,AMP0);axis([-inf inf 0 ])title([str(m),'声的10%带宽图 '])end分析:扩展倍,变换前后的10%带宽由变为,变化率:变窄%;压缩倍,变换前后的10%带宽由变为,变化率:变宽%;(2)男声:扩展倍,变换前后的10%带宽由变为,变化率:变窄%;童声:压缩倍,变换前后的10%带宽由变为,变化率:变宽%;实验七连续系统零极点分析实验报告内容:(1)根据例7-1、7-2、7-3的要求和提示完成。