小波分析函数频谱分析实验报告

实验一函数频谱分析实验报告

一、设函数为f(t)=4sin(10πt)+sin(2πt)+3sin(30πt)+sin(4πt)+2sin(100πt)，对这个函数做频谱分析。

二、程序如下：

clear

t=0:0.01:1;

>> x=4*sin(2*pi*5*t)+sin(2*pi*10*t)+3*sin(2*pi*15*t)+sin(2*pi*20*t)+2*sin(2*pi*50*t); >> plot(x(1:50))

>> y=fft(x,512);

>> pyy=y.*conj(y)/512;

>> f=1000*(0:155)/512;

>> figure,plot(f,pyy(1:156))

xxx=x+1; %加噪声1

>> plot(x(1:50))

>> y=fft(xxx,512);

pyy=y.*conj(y)/512;

f=1000*(0:155)/512;

figure,plot(f,pyy(1:156))

z=randn(size(t)); %加随机信号

xx=x+z;

plot(x(1:50))

y=fft(xx,512);

pyy=y.*conj(y)/512;

f=1000*(0:155)/512;

figure,plot(f,pyy(1:156))

三、运行结果

（1）原始函数的图像：

频谱分析：

（2）加噪声1后的频谱分析：

加随机数后的频谱分析：

应用FFT对信号进行频谱分析实验报告

实验 应用FFT 对信号进行频谱分析 一、实验目的 1、在理论学习的基础上，通过本次实验，加深对快速傅里叶变换的理解，熟悉FFT 算法及其程序的编写。 2、熟悉应用FFT 对典型信号进行频谱分析的方法。 3、了解应用FFT 进行新红啊频谱分析过程中可呢个出现的问题，以便在实际中正确应用FFT 。 二、实验原理 一个连续信号()a x t 的频谱可以用它的傅里叶变换表示为： ()()j t a a X j x t e dt +∞ -Ω-∞Ω=? （2-1） 如果对信号进行理想采样，可以得到离散傅里叶变换： ()()j n X e x n z ω +∞ --∞=∑ （2-2） 在各种信号序列中，有限长序列在数字信号处理中占有很重要的。无限长的序列往往可以用有限长序列来逼近。对于有限长的序列我们可以使用离散傅里叶变换（DFT ），这一序列可以很好的反应序列的频域特性，并且容易利用快速算法在计算机上实现当序列的长度是N 时，我们定义离散傅里叶变换为： 1 0()[()]()N kn N n X k DFT x n x n W -===∑ （2-3） DFT 是对序列傅里叶变换的灯具采样，因此可以用于序列的频谱分析。在利用DFT 进行频谱分析的时候可能有三种误差： （1）混叠现象 序列的频谱是采样信号频谱的周期延拓，周期是2/T π，因此当采样频率不满足奈奎斯特定理，即采样频率1/s f T =小于两倍的信号频率时，经过采样就会发生频谱混叠。这导致采样后的信号序列不能真实的反映原信号的频谱。 （2）泄漏现象 泄漏是不能和混叠完全分开的，因为泄漏导致频谱的扩展，从而造成混淆。为了减小混淆的影响，可以选择适当的窗函数使频谱的扩散减到最小。 （3）栅栏效应 因为DFT 是对单位圆上Z 变换的均匀采样，所以它不可能将频谱视为一个连续的函数。这样就产生了栅栏效应。减小栅栏效应的一个方法是在源序列的末端补一些零值，从而变动DFT 的点数。 三、实验内容和结果 1、观察高斯序列的时域和频域特性 （1）固定高斯序列()a x n 中的参数p=8，当q 为2，4，8时其时域和幅频特性分别如图 2.1，图2.2所示：

Erdas实验报告

E RDAS实验报告 图像融合实验 数据来源 采用Erdas中examples文件内的2000年Atlanta多光谱TM数据和高清全色Pan数据。两图为同一地区不同坐标影像，故使用前需预处理从而得到实验区域。 目的 多光谱TM数据分辨率较低但包含多波段色彩，而全色Pan数据只包含一层高清影像，为了得到研究区域的高清彩色影像，我们将TM和Pan数据在Erdas2014中进行融合以达到实验目的。 方法 在遥感领域运用较多的融合方法有主成分变换法、比值变换法、小波变换法和HIS变换法。本实验则运用HIS变换法。IHS属于色度空间变换，从多光谱彩色合成影像上分离出代表信息的明度（I）和代表光谱信息的色调（H）、饱和度（S）等3个分量，并采用相同区域的高分辨率全色波段数据代替明度（I）进行空间信息融合。 步骤 1.几何校正 因原始图像空间坐标不同，需选取控制点进行几何校正。本实验校正方法为多项式法，选取6个控制点进行校正，其校正叠加截图如下：

2.叠加剪切 由校正结果可知两图像只有部分区域重合，所以建立AOI对重合区域进行剪切，以得到研究区域，截图如下： 3.重采样 因多光谱图像分辨率较低，像元点较大，若要与全色图融合出高清影像需进行重采样来调整像元大小，以达到与高清图一致。 4.二次剪切 因图为栅格，统一像元后，边缘区必然会有一定的扩展(如下图)，虽说扩展的范围较小，但在科研应用方面不符合要求，故须二次剪切。 5.RGB转HIS

TM图像选取前三层再分别赋予蓝、绿、红三色，转化为HIS格式，如下图： 6.直方图匹配 将高清图像直方图以标准图像的直方图为标准作变换,使全色光图和HIS图中I层两图像的直方图相同和近似,从而使两幅图像具有类似的色调和反差，以便作进一步的运算。 7.图像叠加 运用Layer stack功能将全色光高清图和H、S图层进行叠加即所谓的图像融合。它将多波段图层组合到了一起，从而得到新的包含多个有助于研究者使用的多波段影像。 8.IHS转RGB

OBLF光谱仪操作规程-1

4. 4. 1 氩气供应要求 气体型号： 只有纯氩气（含A r99. 995%）及更高纯度的氩气可用。减 压阀门要求： 两级减压阀门要保证纯净的氩气的初始压强为3bar ，流量最低要达到800l / h。 5、样品的准备 ·样品激发之前必须将其磨光。对于钢和铸铁（颗粒大小，对钢是60 到80 之间， 对铸铁是30 到60 之间），建议用适合于本仪器样品的砂轮干磨。 ·砂轮一定不要被高合金含量的样品污染（例如：应该先用于低合金样品，再用于高合金样品）。 ·必须确保样品在磨的过程中没有过热。如果需要，样品应该用水冷却，再干燥，用尽可能短的时间干磨。 ·软质样品（例如：铜、铝、锌、铅）必须车削表面或用酒精湿磨。 ·磨好的表面一定不要玷污（例如用手触摸）。 6、操作 6. 1 仪器的开机与关机 用仪器内部电源接线板主开关开机与关机。 一旦光谱仪的开关置于“OFF”，仪器的所有组件都与电源断开。警 告： 在维修和打开设备面板之前，主开关必须处于“OFF”位置。只有专业人员除外。 在预先知道要停电的情况下，一定要关机。 仪器程序一定要在关机之前退出，以保证所有的硬盘操作全部完成，否则可能造成数据丢失。 关机前注意： 跟所有用电设备一样，这台仪器在打开时也消耗能量，也有一个消耗大小的问题。这样，就该确定一下多长时间开关一次，对于这一问题，取决于设备操作的各个条件， 并要考虑以下这些： ·所有O B L F光谱分析仪能长期开着而无危险。 ·本仪器是复杂的测量系统，在打开之后到达设备的稳定状况和最大测量精度，自然需要一个时间。 ·长期关闭之后，例如，几小时或几天，必然会使真空室的压强增加，自动恒温控制的光学系统的温度降低。这意味着在光谱仪重新开机之后，真空泵和加热系统需要一段时间才能达到分析样品所需要的稳定状态。 ·作为一种防止不必要的能量消耗和物质损耗的有效措施，O B L F为所有的光谱 仪安装了一种辅助功能，即在启动之后的低功耗功能，它只让那些保证分析仪的稳定性和自动检测系统而必需的组件处于工作状态。 与此同时，光谱仪只有在长期不用的情况下（例如：长期的工作假期期间），在可预见的停电，或者在要进行修理时才关闭。在光谱仪没有真正使用期间，为了合理的利用能源，辅助功能应被激活，可采用手动模式，或者是自动模式。

荧光光谱分析仪工作原理

X 荧光光谱分析仪工作原理 用x 射线照射试样时,试样可以被激发出各种波长得荧光x 射线，需要把混合得x 射线 按波长（或能量）分开，分别测量不同波长（或能虽:）得X 射线得强度，以进行左性与定疑 分析，为此使用得仪器叫X 射线荧光光谱仪。由于X 光具有一泄波长，同时又有一立能量, 因此，X 射线荧光光谱仪有两种基本类型:波长色散型与能量色散型。下图就是这两类仪器 得原理图. 用X 射线照射试样时，试样可以被激发出各种波长得荧光X 射线，需要把混合得X 射 线按波长（或能疑）分开，分别测量不同波长（或能量）得X 射线得强度，以进行定性与左疑 分析，为此使用得仪器叫X 射线荧光光谱仪。由于X 光具有一左波长，同时又有一左能量， 因此，X 射线荧光光谱仪有两种基本类型：波长色散型与能量色散型。下图就是这两类仪器 得原理图。 （a ）波长色散谱仪 （b ）能虽色散谱仪 波长色散型和能量色散型谱仪原理图 现将两种类型X 射线光谱仪得主要部件及工作原理叙述如下: X 射线管 酥高分析器 分光晶体 计算机 再陋电源

丝电源 灯丝 电了悚 X则线 BeiV 輪窗型X射线管结构示意图 两种类型得X射线荧光光谱仪都需要用X射线管作为激发光源?上图就是X射线管得结构示意图。灯丝与靶极密封在抽成貞?空得金属罩内，灯丝与靶极之间加高压（一般为4OKV）, 灯丝发射得电子经高压电场加速撞击在靶极上，产生X射线。X射线管产生得一次X射线, 作为激发X射线荧光得辐射源.只有当一次X射线得波长稍短于受激元素吸收限Imi n时，才能有效得激发出X射线荧光?笥？SPAN Ian g =EN-U S >lmin得一次X射线其能量不足以使受激元素激发。 X射线管得靶材与管工作电压决立了能有效激发受激元素得那部分一次X射线得强度。管 工作电压升高，短波长一次X射线比例增加，故产生得荧光X射线得强度也增强。但并不就是说管工作电压越髙越好，因为入射X射线得荧光激发效率与苴波长有关，越靠近被测元素吸收限波长，激发效率越髙。A X射线管产生得X射线透过彼窗入射到样品上, 激发岀样品元素得特征X射线，正常工作时，X射线管所消耗功率得0、2%左右转变为X 射线辐射，其余均变为热能使X射线管升温，因此必须不断得通冷却水冷却靶电极。 2、分光系统 第?准讥器 平面晶体反射X线示意图 分光系统得主要部件就是晶体分光器，它得作用就是通过晶体衍射现彖把不同波长得X射线分开.根据布拉格衍射左律2d S in 0 =n X ,当波长为X得X射线以0角射到晶体，如果晶面间距为d,则在出射角为0得方向，可以观测到波长为X =2dsi n 0得一级衍射及波长为X/2, X /3 ------ ―等髙级衍射。改变（）角，可以观测到另外波长得X

用FFT对信号作频谱分析 实验报告

实验报告 实验三：用FFT 对信号作频谱分析 一、 实验目的与要求 学习用FFT 对连续信号和时域离散信号进行谱分析的方法，了解可能出现的分析误差及其原因，以便正确应用FFT 。 二、 实验原理 用FFT 对信号作频分析是学习数字信号处理的重要内容，经常需要进行分析的信号是模拟信号的时域离散信号。对信号进行谱分析的重要问题是频谱分辨率D 和分析误差。频谱分辨率直接和FFT 的变换区间N 有关，因为FFT 能够实现的频率分辨率是2π/N ，因此要求2π/N 小于等于D 。可以根据此式选择FFT 的变换区间N 。误差主要来自于用FFT 作频谱分析时，得到的是离散谱，而信号（周期信号除外）是连续谱，只有当N 较大时，离散谱的包络才能逼近连续谱，因此N 要适当选择大一些。 三、 实验步骤及内容（含结果分析） （1）对以下序列进行FFT 分析： x 1(n)=R 4(n) x 2(n)= x 3(n)= 选择FFT 的变换区间N 为8和16两种情况进行频谱分析，分别打印出幅频特性曲线，并进行讨论、分析与比较。 【实验结果如下】： n+1 0≤n ≤3 8-n 4≤n ≤7 0 其它n 4-n 0≤n ≤3 n-3 4≤n ≤7 0 其它 n

实验结果图形与理论分析相符。（2）对以下周期序列进行谱分析： x4(n)=cos[(π/4)*n]

x5(n)= cos[(π/4)*n]+ cos[(π/8)*n] 选择FFT的变换区间N为8和16两种情况进行频谱分析，分别打印出幅频特性曲线，并进行讨论、分析与比较。 【实验结果如下】： （3）对模拟周期信号进行频谱分析： x6(n)= cos(8πt)+ cos(16πt)+ cos(20πt) 选择采样频率Fs=64Hz，FFT的变换区间N为16、32、64三种情况进行频谱分析，分别打印出幅频特性曲线，并进行讨论、分析与比较。 【实验结果如下】：

基于小波信号的噪声消除matlab实验报告

南京师范大学物理科学与技术学院 医用电子学论文 论文名称：基于小波变换的心电信号噪声消除 院系：物科院 专业：电路与系统 姓名：聂梦雅 学号： 121002043 指导教师：徐寅林

摘要 以小波变换的多分辨率分析为基础, 通过对体表心电信号(ECG) 及其噪声的分析, 对ECG信号中存在的基线漂移、工频干扰及肌电干扰等几种噪声, 设计了不同的小波消噪算法; 并利用MIT/BIH 国际标准数据库中的ECG 信号和程序模拟所产生的ECG 信号, 分别对算法进行了仿真与实验验证。结果表明, 算法能有效地滤除ECG 信号检测中串入的几类主要噪声, 失真度很小, 可满足临床分析与诊断对ECG 波形的要求。 关键词: ECG 信号, 小波变换, 基线漂移, 工频干扰, 肌电干扰

Abstract We apply the multi-resolution analysis (MRA ) of wavelet transform ( WT ) , which was proposed by Mallat [ 5 ] , to suppress the three main types of noises existing in electrocardiogram ( ECG ) signals : baseline wander, power line interference and electro my ographical interference. We apply Mallat algorithm [ 4 ] to suppress the baseline wander in ECG signals. We apply the sof t-thresholding algorithm, proposed by donohoetal on the basis of MRA of WT , to suppress power line interference in ECG signals. We apply Mallat algorithm and then the algorithm proposed by Donohoetal to suppress the electro my ographical interference in ECG signals ,who sefrequency range varies f rom 5Hz to 2kHz. We performed simulations ,using both ECG signals from MIT/BIH database, and ECG signals generated via computer simulation .The results show that the algorithm can suppress the main no isesexisting in ECG signals efficiently with very little distortion, and can satisfy the requirement s of clinical analysis and diagnosis on ECG waveforms. Key words: ECG (electro cardio gram ) signal, wavelet transform , baseline wander, power line interference , electro my ographical interference

光谱仪操作规程

光谱仪操作规程 一、实验室条件： 1、实验室要求温度控制在16--30摄氏度之间；开启设备 之前应首先将实验室温度升起。 2、实验室相对湿度要求控制在30—70%。 3、实验室应保持清洁，无灰尘、震动。 4、试样要求厚度超过1mm，覆盖面积控制在15-50mm，测 试面必须光滑。用40-80目砂轮或铣床制作，制作后不能用手或水清理，可以用钢丝刷进行清理。 二、开关机、操作注意事项：（操作必须顺序进行） 1、先检查各部位是否正常： （1）、检查循环水位是否正常，控制在总高度的1/3处，最高不超过总高度4/5. （2）、检查真空泵油位是否正常。 （3）、检查氩气是否充足。 2、开机： （1）、打开稳压器电源开关 （2）、打开光谱仪主电源开关即MAINS 16A.(1#) （3）、打开真空泵电源开关即VACUUMPUMP(2) (4）、打开水泵电源开关即 WATERPUMP(3) (5)、打开电子板电源开关即ELECTRONICS(4) (6)、打开负高压电源开关即HVPS(5)

（7）、打开电脑，用户名输入ARL回车进入，选择程序进入即可，再点登陆用户名输入U（一般操作）或M（管理员，可以更改程序）回车进入；查看设备状态，看数据是否控制在范围内。(所有数据必须没有大于或小于号)一般开机后4至5小时后才能使用。 （8）、打开氩气，压力表控制在0.35MPA，如果氩气净化机需再生时将压力表调至0.5MPA。（现在的表指针+0.35）先用净化机的II塔（按净化机的指示旋转手柄到指定位置） 3、关机： （1）、正常关机按开机的反顺序操作，必须按步骤关机。（2）、遇到突然停电时应先关闭光谱仪主电源开关，防止突然来电损坏设备。（禁忌设备用电急停急送） 4、程序操作： （1）、每天第一次应先做标准样测试：先登录程序OX SAS，选择管理员用户名进入界面。点‘测量类型标准样’点‘任务’再点‘TS’，点‘方法’再点‘FELAST’，‘类型标准’点方法中以创建的方法，再创建名称即可。将标准样放于‘激发台’开始打样，连续打三点点击完成。（如45#、35#类型标准样）看测试的数据相对数据小于5%即可。 （2）、点‘定量分析’开始对工作中试样进行打样，打样后点击完成。 三、维护与保养：

各种仪器分析的基本原理及谱图表示方法!!!紫外吸收光谱UV分析

各种仪器分析的基本原理及谱图表示方法！！！ 紫外吸收光谱UV 分析原理：吸收紫外光能量，引起分子中电子能级的跃迁谱图的表示方法：相对吸收光能量随吸收光波长的变化提供的信息：吸收峰的位置、强度和形状，提供分子中不同电子结构的信息荧光光谱法FS 分析原理：被电磁辐射激发后，从最低单线激发态回到单线基态，发射荧光谱图的表示方法：发射的荧光能量随光波长的变化提供的信息：荧光效率和寿命，提供分子中不同电子结构的信息红外吸收光谱法IR 分析原理：吸收红外光能量，引起具有偶极矩变化的分子的振动、转动能级跃迁谱图的表示方法：相对透射光能量随透射光频率变化提供的信息：峰的位置、强度和形状，提供功能团或化学键的特征振动频率拉曼光谱法Ram 分析原理：吸收光能后，引起具有极化率变化的分子振动，产生拉曼散射谱图的表示方法：散射光能量随拉曼位移的变化提供的信息：峰的位置、强度和形状，提供功能团或化学键的特征振动频率核磁共振波谱法NMR 分析原理：在外磁场中，具有核磁矩的原子核，吸收射频能量，产生核自旋能级的跃迁谱图的表示方法：吸收光能量随化学位移的变化提供的信息：峰的化学位移、强度、裂分数和偶合常数，提供核的数目、所处化学环境和几何构型的信息电子顺磁共振波谱法ESR 分析原理：在外磁场中，分子中未成对电子吸收射频能量，产生电子自旋能级跃迁谱图的表示方法：吸收光能量或微分能量随磁场强度变化提供的信息：谱线位置、强度、裂分数目和超精细分裂常数，提供未成对电子密度、分子键特性及几何构型信息 质谱分析法MS 分析原理：分子在真空中被电子轰击，形成离子，通过电磁场按不同m/e 分离 谱图的表示方法：以棒图形式表示离子的相对峰度随m/e 的变化提供的信息：分子离子及碎片离子的质量数及其相对峰度，提供分子量，元素组成及结构的信息气相色谱法GC 分析原理：样品中各组分在流动相和固定相之间，由于分配系数不同而分离谱图的表示方法：柱后流出物浓度随保留值的变化提供的信息：峰的保留值与组分热力学参数有关，是定性依据；峰面积与组分含量有关反气相色谱法IGC 分析原理：探针分子保留值的变化取决于它和作为固定相的聚合物样品之间的相互作用力谱图的表示方法：探针分子比保留体积的对数值随柱温倒数的变化曲线提供的信息：探针分子保留值与温度的关系提供聚合物的热力学参数裂解气相色谱法PGC 分析原理：高分子材料在一定条件下瞬间裂解，可获得具有一定特征的碎片谱图的表示方法：柱后流出物浓度随保留值的变化提供的信息：谱图的指纹性或特征碎片峰，表征聚合物的化学结构和几何构型凝胶色谱法GPC 分析原理：样品通过凝胶柱时，按分子的流体力学体积不同进行分离，大分子先流出谱图的表示方法：柱后流出物浓度随保留值的变化提供的信息：高聚物的平均分子量及其分布热重法TG 分析原理：在控温环境中，样品重量随温度或时间变化谱图的表示方法：样品的重量分数随温度或时间的变化曲线提供的信息：曲线陡降处为样品失重区，平台区为样品的热稳定区热差分析DTA 分析原理：样品与参比物处于同一控温环境中，由于二者导热系数不同产生温差，记录温度随环境温度或时间的变化 谱图的表示方法：温差随环境温度或时间的变化曲线提供的信息：提供聚合物热转变温度及各种热效应的信息示差扫描量热分析DSC 分析原理：样品与参比物处于同一控温环境中，记录维持温差为零时，所需能量随环境温度或时间的变化 谱图的表示方法：热量或其变化率随环境温度或时间的变化曲线提供的信息：提供聚合物热转变温度及各种热效应的信息静态热―力分析TMA 分析原理：样品在恒力作用下产生的形变随温度或时间变化谱图的表示方法：样品形变值随温度或时间变化曲线提供的信息：热转变温度和力学状态

频域分析实验报告

频域分析实验报告 班级： 学号： 姓名：

一、实验内容： 1利用计算机作出开环系统的波特图； 2、观察记录控制系统的开环频率特性； 3、控制系统的开环频率特性分析。 二、仿真原理： 对数频率特性图（波特图）： 对数频率特性图包括了对数幅频特性图和对数相频特性图。横坐标为频率w，采用对数分度，单位为弧度/秒；纵坐标均匀分度，分别为幅值函数20lgA(w)，以dB表示；相角，以度表示。MATLAB提供了函数bode()来绘制系统的波特图，其用法如下： （1）bode(num，den)：可绘制出以连续时间多项式传递函数表示的系统的波特图。 （2）当带输出变量[mag，pha，w]或[mag，pha]引用函数时，可得到系统波特图相应的幅值mag、相角pha及角频率点w矢量或只是返回幅值与相角。相角以度为单位，幅值可转换为分贝单位：magdb=20×log10(mag) 二、实验验证 1、用Matlab作Bode图。要求:画出对应Bode图。 （1）G(S)=25/S2+4s+25 （7）G(S)=9(s2+0.2s+1)/s(s2+1.2s+9)；

图 1 图 2 (1)G(S)=25/S2+4s+25 可以看成是一个比例环节和一个振荡环节组成，所以k=1，T1=0.04,因为v=0，所以在转折频率之前都为20lgk，因为k=1所以斜率为0，经过转折频率，分段直线斜率的变化量为-40db/dec。

（7）G(S)=9(s2+0.2s+1)/s(s2+1.2s+9)； 可以看成是一个二阶微分环节和一个积分环节和一个振荡环节组成，化常数为1后，v=1，t1=1，t2=1/3，所以我们可以看到，在起始阶段是-20*vdb/dec，所以一开始斜率为-20db/dec。当经过1/3的转折频率之后分段直线的改变量为40db/dec，当经过1的转折频率之后分段直线的改变量为-40db/dec。故图像如图所示。 第二题： 典型二阶系统Gs=Wn2/s2+2ζWns+Wn2，试绘制取不同值时的Bode图。取Wn=8，ζ=0.1，0.2，0.3，,0.5，0.6； 图 3 如图所示。

小波变换

《医学图像处理》实验报告 实验十：小波变换 日期： 2014年05月06日 摘要 本次实验的实验目的及主要内容是: 一维小波变换和反变换 二维小波变换和反变换 二维小波细节置零、去噪

一、技术讨论 1.1实验原理 小波变换的原理:是指一组衰减震动的波形，其振幅正负相间变化为零，是具有一定的带宽和中心频率波组。小波变换是用伸缩和平移小波形成的小波基来分解（变换）或重构（反变换）时变信号的过程。不同的小波具有不同带宽和中心频率，同一小波集中的带宽与中心频率的比是不变的，小波变换是一系列的带通滤波响应。它的数学过程与傅立叶分析是相似的，只是在傅立叶分析中的基函数是单频的调和函数，而小波分析中的基函数是小波，是一可变带宽内调和函数的组合。 小波去噪的原理：利用小波变换把含噪信号分解到多尺度中，小波变换多采用二进型，然后在每一尺度下把属于噪声的小波系数去除，保留并增强属于信号的小波系数，最后重构出小波消噪后的信号。其中关键是用什么准则来去除属于噪声的小波系数，增强属于信号的部分。 1.2实验方法 1）dwt函数（实现1-D离散小波变换） [cA,cD]=dwt(X,’wname’)使用指定的小波基函数‘wname’对信号X进行分解，cA和cD分别是近似分量和细节分量； [cA,cD]=dwt(X,Lo_D,Hi_D)用指定的滤波器组Lo_D,Hi_D对信号进行分解 2）idwt函数（实现1-D离散小波反变换） X=idwt(cA,cD,’wname’) X=idwt(cA,cD,Lo_R,Hi_R) X=idwt(cA,cD,’wname’,L) X=idwt(cA,cD,Lo_R,Hi_R,L) 由近似分量cA和细节分量cD经过小波反变换，选择某小波函数或滤波器组，L为信号X中心附近的几个点 3）dwt2函数（实现2-D离散小波变换） [cA,cH,cV,cD]=dwt2(X,’wname’) [cA,cH,cV,cD]=dwt2(X,’wname’) cA近似分量，cH水平细节分量，cV垂直细节分量，cD对角细节分量 4）idwt2函数（实现2-D离散反小波变换） X=idwt2(cA,cH,cV,cD,’wname’) X=idwt2(cA,cH,cV,cD,Lo_R,Hi_R) X=idwt2(cA,cH,cV,cD,’wname’,S) X=idwt2(cA,cH,cV,cD,Lo_R,Hi_R,S)

光谱仪的操作规程(通用版)

光谱仪的操作规程(通用版) The safety operation procedure is a very detailed operation description of the work content in the form of work flow, and each action is described in words. ( 安全管理 ) 单位：______________________ 姓名：______________________ 日期：______________________ 编号：AQ-SN-0916

光谱仪的操作规程(通用版) 1准备 1.1检查实验间的工作环境是否满足要求，温度：22士 2.8℃，湿度≤60％； 1.2气的纯度是否小于99.996％，检查氩气的流量是否小于 30psi； 1.3样品的制备：保证样品表面清洁无油污，无氧化皮，表面均匀没有沙眼。 2操作 2.1在电压稳定的情况下，打开系统电源，打开真空泵开关，要求真空泵空转五分钟后，打开真空阀，稳定一段时间要求真空度达到1μ之内； 2.2在真空度达到要求的情况下，打开高压和低压电源，检查氩

气流量达到要求后，打开激发电源； 2.3打开微机，进入BAIRD操作系统； 2.4汞灯的校准：打开汞灯开关，调节旋转数字指示器的旋钮，按照计算机屏幕指令操作，并在仪器运转记录本上记录俩个半峰值的读数，将指示器最终调至计算机求出的答案为止，锁紧旋钮； 2.5清洗循环系统：取一空样在试样台上，激发数次，待各数值达到稳定。注意在装夹试样的时候，要将弹簧负载紧固轴轻轻的压到样品表面，不要突然松手，以免压碎台板，并注意氮化硼原盘固定在原来的位置； 2.6标准化：放好样品后，进行全体元素的标准化。依照电脑中的提示，依次对各块标钢进行激发，得到一组数据，通过和标钢原有数据的对比，计算机自动得到一个修正系数。修正系数不小于 3.0时，可以进行试样分析，否则，要对不满足要求的元素进行元素的标准化，使该元素的修正系数达到规定要求； 2.7在进行了标准化后，进行试样的分析。一个样品测量多次，计算机自动给出平均值。通过打印机得到试验数据。试验数据记录

控制系统的频域分析实验报告

实验名称： 控制系统的频域分析 实验类型：________________同组学生姓名：__________ 一、实验目的和要求 用计算机辅助分析的方法，掌握频率分析法的三种方法，即Bode 图、Nyquist 曲线、Nichols 图。 二、实验内容和原理 （一）实验原理 1．Bode(波特)图 设已知系统的传递函数模型： 1 1211121)(+-+-+???+++???++=n n n m m m a s a s a b s b s b s H 则系统的频率响应可直接求出： 1 1211121)()()()()(+-+-+???+++???++=n n n m m m a j a j a b j b j b j H ωωωωω MATLAB 中，可利用bode 和dbode 绘制连续和离散系统的Bode 图。 2．Nyquist(奈奎斯特)曲线 Nyquist 曲线是根据开环频率特性在复平面上绘制幅相轨迹，根据开环的Nyquist 线，可判断闭环系统的稳定性。 反馈控制系统稳定的充要条件是，Nyquist 曲线按逆时针包围临界点(-1，j0)p 圈，为开环传递函数位于右半s 一平面的极点数。在MATLAB 中，可利用函数nyquist 和dnyquist 绘出连续和离散系统的乃氏曲线。 3．Nicho1s(尼柯尔斯)图 根据闭环频率特性的幅值和相位可作出Nichols 图，从而可直接得到闭环系统的频率特性。在 MATLAB 中，可利用函数nichols 和dnichols 绘出连续和离散系统的Nichols 图。 （二）实验内容 1．一系统开环传递函数为 ) 2)(5)(1(50)(-++=s s s s H 绘制系统的bode 图，判断闭环系统的稳定性，并画出闭环系统的单位冲击响应。 2．一多环系统 ) 10625.0)(125.0)(185.0(7.16)(+++=s s s s s G 其结构如图所示 试绘制Nyquist 频率曲线和Nichols 图，并判断稳定性。 （三）实验要求

哈工大小波分析上机实验报告

小波分析上机实验报告 院系：电气工程及自动化学院 学科：仪器科学与技术

实验一小波分析在信号压缩中的应用 一、试验目的 （1）进一步加深对小波分析进行信号压缩的理解； （2）学习Matlab中有关信号压缩的相关函数的用法。 二、相关知识复习 用一个给定的小波基对信号进行压缩后它意味着信号在小波阈的表示相对缺少了一些信息。之所以能对信号进行压缩是因为对于规则的信号可以用很少的低频系数在一个合适的小波层上和一部分高频系数来近似表示。 利用小波变换对信号进行压缩分为以下几个步骤来完成： （1）进行信号的小波分解； （2）将高频系数进行阈值量化处理。对从1 到N 的每一层高频系数都可以选择不同的阈值并且用硬阈值进行系数的量化； （3）对量化后的系数进行小波重构。 三、实验要求 （1）对于某一给定的信号（信号的文件名为leleccum.mat），利用小波分析对信号进行压缩处理。 （2）给出一个图像，即一个二维信号（文件名为wbarb.mat），利用二维小波分析对图像进行压缩。 四、实验结果及程序 （1）load leleccum %将信号装入Matlab工作环境 %设置变量名s和ls，在原始信号中，只取2600-3100个点 s = leleccum(2600:3100); ls = length(s); %用db3对信号进行3级小波分解 [c,l] = wavedec(s, 3, 'db3'); %选用全局阈值进行信号压缩 thr = 35; [xd,cxd,lxd,perf0,perfl2] = wdencmp('gbl',c,l,'db3',3,thr,'h',1); subplot(2,1,1);plot(s); title('原是信号s'); subplot(2,1,2);plot(xd); title('压缩后的信号xd');

光谱仪操作规程

光谱仪操作规程 1.目的 光谱仪属于精密贵重仪器，为了操作人员更加合理、安全操作光谱仪，保证光谱仪能正常长时间的工作，特制定本规程。 2.使用范围 本公司操作光谱人员从事有色金属材质成分检测。 3.工作职责 品保部主管负责培训操作光谱人员从事有色金属材质成分检测。 4.工作流程 4.1 操作流程： 4.1.1 操作之前检查光谱仪无异常，仪器与外部设备接线连接是否正常。 4.1.2 开机顺序：把配电柜开关打开，开启稳压器十几秒后再开光谱仪电源。 4..1.3 将准备好的99.999%高纯氩气连接到氩气净化器,打开高纯氩气瓶开关，分压压力表 气压调至0.25-0.3Mpa之间，连接好后用肥皂水或洗衣粉水泡沫分别涂入氩气瓶及净化器各接入口，检验是否有漏气现象。正常后开启净化器电源，待净化器工作温控稳定260度方可与光谱同步工作。 4.1.4 光谱仪电源开机顺序：开启主电源（Main），间隔5-10秒、再开真空泵（Vacuum pump） 电源、接着开水泵(Water pump)电源、然后开电子柜（Electronis）电源、间隔5-10秒,最后开高压（Hvps）电源。光谱仪电源严格按照此顺序操作，否则容易烧坏光谱仪部件。 4.1.5 计算机电源由光谱仪仪器有专门接出电源，打开计算机电源，双击OXSAS打开操作 软件。 4.1.6 光谱仪开机后需仪器稳定后才能工作，一般仪器稳定为5小时以上（含停电后再开 启仪器）。检查仪器是否稳定，从“工具-操作-读取状态【F7】”,检查各参数值应在最小～最大范围内，方可进行下面检测工作。 4.1.7 将制备好的样品正确放在激发台上，然后在OXSAS软件中单击“分析和数据”选择 “定量分析【F10】”，弹出分析界面。 4.1.8 将弹出分析界面“当前位置”，在“任务”里下拉菜单选择Conc Cu,根据样品种类 在“方法”中选择下拉菜单适合的方法,如：黄铜（CUBAS）、青铜（CUBAON）、铝青铜（CUCUAI） 在结果格式选择“Type standard”，元素格式选择“Cu%”，“值”菜单栏输入样品名称，然后将加工好的样品放入激发台，再分析界面鼠标点击“Sidok+开始”进入试验，同一个样品重复试验，点击“继续”，试验结束点击“完成” 4.1.9 打印结果：点击“分析和数据”选择“查看结果【F8】”弹出界面，选定需要打印 数据一栏，点击“打印分析”弹出窗口“打印机选择”，在“打印机选择”窗口中打印格式默认为Standard不需要改动，勾上“包含每次测量”和“分析组”，打印机栏选定打印机型号，然后点击确认。 4.1.10 光谱仪未发生故障，仪器正常运行，不允许关机，否则对测量数据会产生误差。若 需要关机，开机按照上述工作流程操作，然后需对光谱仪作标准化，标准化参照流程操作。 4.1.11 关机顺序：先关高压电源（HVPS），再关电子柜，接着关水泵，然后关真空泵. 4.2光谱仪标准化操作流程： 4.2.1 点击“分析数据”选择“标准化【F3】”，进入界面中“任务（Task）”中下拉菜单选

光谱仪的工作原理

光谱仪的工作原理-标准化文件发布号：（9556-EUATWK-MWUB-WUNN-INNUL-DDQTY-KII

光谱仪的工作原理元素的原子在激发光源的作用下发射谱线，谱线经光栅分光后形成光谱，每种元素都有自己的特征谱线，谱线的强度可以代表试样中元素的含量，用光电检测器将谱线的辐射能转换成电能。检测输出的信号，经加工处理，在读出装置上显示出来。然后根据相应的标准物质制作的分析曲线，得出分析试样中待测元素的含量。 表面轮廓仪介绍 表面轮廓仪 - 简介 表面轮廓仪LK-200M型表面轮廓仪采用广精精密最新的基于windows版本的测量软件，具有强大卓越的数据处理分析功能。测量时，零件装夹位置即使任意放置，也能得到满意的测量结果；即使需要测量长度为220mm的工件，测量软件也能保证其1μm的采样步长。 LK-200H型表面轮廓仪采用耐用可靠的16位A/D功能板，其极高的分辨率量程比（1/65536），用户即使需要大量程测量，仍能保持极高的测量精度。 LK-200M型表面轮廓仪采用工控计算机处理测量数据及仪器控制操作。其高质量、高可靠性及突出的防尘、防振、防油、防静电能力使广精精密用户将使用维护成本降至最低。 表面轮廓仪 - 原理 表面轮廓仪LK-200M型表面轮廓仪采用直角坐标法，传感器移动式。直线运动导轨采用高精度气浮导轨，作为测量基准； 电器部分由高级计算机组成；测量软件采用基于中文版Windows操作系统平台的系统测量软件，完成数据采集、处理及测量数据管理等工作。 表面轮廓仪 - 功能 角度处理：两直线夹角、直线与Y轴夹角、直线与X轴夹角 点线处理：两直线交点、交点到直线距离、交点到交点距离、交点到圆心距离、交点到点距离 圆处理：圆心距离、圆心到直线的距离、交点到圆心的距离、直线到切点的距离线处理:直线度、凸度、LG凸度、对数曲线 表面轮廓仪 - 技术规格 表面轮廓仪测量长度：≤200mm

信号与系统实验报告实验三 连续时间LTI系统的频域分析

实验三 连续时间LTI 系统的频域分析 一、实验目的 1、掌握系统频率响应特性的概念及其物理意义； 2、掌握系统频率响应特性的计算方法和特性曲线的绘制方法，理解具有不同频率响应特性的滤波器对信号的滤波作用； 3、学习和掌握幅度特性、相位特性以及群延时的物理意义； 4、掌握用MA TLAB 语言进行系统频响特性分析的方法。 基本要求：掌握LTI 连续和离散时间系统的频域数学模型和频域数学模型的MATLAB 描述方法，深刻理解LTI 系统的频率响应特性的物理意义，理解滤波和滤波器的概念，掌握利用MATLAB 计算和绘制LTI 系统频率响应特性曲线中的编程。 二、实验原理及方法 1 连续时间LTI 系统的频率响应 所谓频率特性，也称为频率响应特性，简称频率响应（Frequency response ），是指系统在正弦信号激励下的稳态响应随频率变化的情况，包括响应的幅度随频率的变化情况和响应的相位随频率的变化情况两个方面。 上图中x(t)、y(t)分别为系统的时域激励信号和响应信号，h(t)是系统的单位冲激响应，它们三者之间的关系为：)(*)()(t h t x t y =，由傅里叶变换的时域卷积定理可得到： )()()(ωωωj H j X j Y = 3.1 或者： ) () ()(ωωωj X j Y j H = 3.2 )(ωj H 为系统的频域数学模型，它实际上就是系统的单位冲激响应h(t)的傅里叶变换。即 ? ∞ ∞ --= dt e t h j H t j ωω)()( 3.3

由于H(j ω)实际上是系统单位冲激响应h(t)的傅里叶变换，如果h(t)是收敛的，或者说是绝对可积（Absolutly integrabel ）的话，那么H(j ω)一定存在，而且H(j ω)通常是复数，因此，也可以表示成复数的不同表达形式。在研究系统的频率响应时，更多的是把它表示成极坐标形式： ) ()()(ω?ωωj e j H j H = 3.4 上式中，)j (ωH 称为幅度频率相应（Magnitude response ），反映信号经过系统之后，信号各频率分量的幅度发生变化的情况，)(ω?称为相位特性（Phase response ），反映信号经过系统后，信号各频率分量在相位上发生变换的情况。)(ωj H 和)(ω?都是频率ω的函数。 对于一个系统，其频率响应为H(j ω)，其幅度响应和相位响应分别为)(ωj H 和)(ω?，如果作用于系统的信号为t j e t x 0)(ω=，则其响应信号为 t j e j H t y 0)()(0ωω= t j j e e j H 00)(0)(ωω?ω=))((000)(ω?ωω+=t j e j H 3.5 若输入信号为正弦信号，即x(t) = sin(ω0t )，则系统响应为 ))(sin(|)(|)sin()()(00000ω?ωωωω+==t j H t j H t y 3.6 可见，系统对某一频率分量的影响表现为两个方面，一是信号的幅度要被)(ωj H 加权，二是信号的相位要被)(ω?移相。 由于)(ωj H 和)(ω?都是频率ω的函数，所以，系统对不同频率的频率分量造成的幅度和相位上的影响是不同的。 2 LTI 系统的群延时 从信号频谱的观点看，信号是由无穷多个不同频率的正弦信号的加权和（Weighted sum ）所组成。正如刚才所述，信号经过LTI 系统传输与处理时，系统将会对信号中的所有频率分量造成幅度和相位上的不同影响。从相位上来看，系统对各个频率分量造成一定的相位移（Phase shifting ），相位移实际上就是延时（Time delay ）。群延时（Group delay ）的概念能够较好地反

小波实验报告一维Haar小波2次分解

一、题目：一维Haar 小波2次分解 二、目的：编程实现信号的分解与重构 三、算法及其实现：离散小波变换 离散小波变换是对信号的时－频局部化分析，其定义为：/2200()(,)()(),()()j j Wf j k a f t a t k dt f t L R φ+∞---∞=-∈? 本实验实现对信号的分解与重构： （１）信号分解：用小波工具箱中的dwt 函数来实现离散小波变换，函数dwt 将信号分解为两部分，分别称为逼近系数和细节系数（也称为低频系数和高频系数）,实验中分别记为cA1,cD1，它们的长度均为原始信号的一半，但dwt 只能实现原始信号的单级分解。在本实验中使用小波函数db1来实现单尺度小波分解，即： [cA1,cD1]＝dwt(s,’db1’)，其中s 是原信号;再通过[cA2,cD2]＝dwt(cA1,’db1’)进行第二次分解，长度又为cA2的一半。 （２）信号重构：用小波工具箱中的upcoef 来实现，upcoef 是进行一维小波分解系数的直接重构,即： A1 = upcoef('a',cA1,'db1'); D1 = upcoef('a',cD1,'db1')。 四、实现工具：Matlab 五、程序代码： %装载leleccum 信号 load leleccum; s = leleccum(1:3920); %用小波函数db1对信号进行单尺度小波分解 [cA1,cD1]=dwt(s,'db1'); subplot(3,2,1); plot(s); title('leleccum 原始信号'); %单尺度低频系数cA1向上一步的重构信号 A1 = upcoef('a',cA1,'db1'); %单尺度高频系数cD1向上一步的重构信号 D1 = upcoef('a',cD1,'db1'); subplot(3,2,3); plot(A1); title('单尺度低频系数cA1向上一步的重构信号'); subplot(3,2,5); plot(D1); title('单尺度高频系数cD1向上一步的重构信号'); [cA1,cD1]=dwt(cA1,’db1'); subplot(3,2,2); plot(s); title('leleccum 第一次分解后的cA1信号'); %第二次分解单尺度低频系数cA2向上一步的重构信号 A2= upcoef('a',cA2,'db1',2); %第二次分解单尺度高频系数cD2向上一步的重构信号 D2 = upcoef('a',cD2,'db1',2); subplot(3,2,4); plot(A2);