当前位置：文档之家 › 实验五 谱分析

实验五 谱分析

  • 格式：doc
  • 大小：1.20 MB
  • 文档页数：17

实验五 谱分析 一 实验目的 1 研究不同类型的窗函数，研究一些不同的方法来测试窗的性能； 2 专注于有关窄带信号的几个不同的情形。

二 实验原理 信号是无限长的，而在进行信号处理时只能采用有限长信号，所以需要将窄带信号“截断”。在信号处理中，“截断”被看成是用一个有限长的“窗口”看无限长的信倍号，或者从分析的角度是无限长的信号x(t)乘以有限长的窗函数w(t),由傅立叶变换性质可知： x(t)w(t)==1/2πX(jw)W(jw) 如果x(t)是频带有限信号，而w(t)是频带无限函数，截断后的信号也必是频带无限信号，从而产生所谓的频谱泄露。频谱泄露是不可避免的，但是尽量减小，因此设计了不同的窗函数满足不同的要求。从能量的角度，频谱泄露也是能量泄露，因为加窗后，是原来的信号集中在宅频带内的能量分散到无限的频带范围。

三 实验内容、程序及结果 1&&2：用MATLAB编程绘制各种窗函数的形状及其幅度响应。 %矩形窗 N=20; figure(1) w1=boxcar(N); n=[0:N-1]; [H,W]=dtft(w1,1024); subplot(121); stem(n,w1),title('矩形窗形状'); xlabel('n'),ylabel('w(n)'); subplot(122); plot(W/2/pi,abs(H)),title('矩形窗幅度响应'); xlabel('f'),ylabel('|W(w)|') %汉宁窗 figure(2) w2=hanning(N); n=[0:N-1]; [H,W]=dtft(w2,1024); subplot(121); stem(n,w2),title('汉宁窗形状'); xlabel('n'),ylabel('w(n)'); subplot(122); plot(W/2/pi,abs(H)),title('汉宁窗幅度响应'); xlabel('f'),ylabel('|W(w)|') %汉明窗 figure(3) w3=hamming(N); n=[0:N-1]; [H,W]=dtft(w3,1024); subplot(121); stem(n,w3),title('汉明窗形状'); xlabel('n'),ylabel('w(n)'); subplot(122); plot(W/2/pi,abs(H)),title('汉明窗幅度响应'); xlabel('f'),ylabel('|W(w)|') %巴特利特窗 figure(4) w4=bartlett(N); n=[0:N-1]; [H,W]=dtft(w4,1024); subplot(121); stem(n,w4),title('巴特利特窗形状'); xlabel('n'),ylabel('w(n)'); subplot(122); plot(W/2/pi,abs(H)),title('巴特利特窗幅度响应'); xlabel('f'),ylabel('|W(w)|') %布莱克曼窗 figure(5) w5=blackman(N); n=[0:N-1]; [H,W]=dtft(w5,1024); subplot(121); stem(n,w5),title('布莱克曼窗形状'); xlabel('n'),ylabel('w(n)'); subplot(122); plot(W/2/pi,abs(H)),title('布莱克曼窗幅度响应'); xlabel('f'),ylabel('|W(w)|') %Triang窗 figure(6) w6=triang(N); n=[0:N-1]; [H,W]=dtft(w6,1024); subplot(121); stem(n,w6),title('Triang窗形状'); xlabel('n'),ylabel('w(n)'); subplot(122); plot(W/2/pi,abs(H)),title('Triang窗幅度响应'); xlabel('f'),ylabel('|W(w)|') %Kaiser窗 figure(7) w7=kaiser(N,2); n=[0:N-1]; [H,W]=dtft(w7,1024); subplot(121); stem(n,w7),title('Kaiser窗形状'); xlabel('n'),ylabel('w(n)'); subplot(122); plot(W/2/pi,abs(H)),title('Kaiser窗幅度响应'); xlabel('f'),ylabel('|W(w)|') %切比雪夫窗 figure(8) w8=chebwin(N,16); n=[0:N-1]; [H,W]=dtft(w8,1024); subplot(121); stem(n,w8),title('切比雪夫窗形状'); xlabel('n'),ylabel('w(n)'); subplot(122); plot(W/2/pi,abs(H)),title('切比雪夫窗幅度响应'); xlabel('f'),ylabel('|W(w)|') 附图： 3.绘制矩形窗的幅频响应，窗长度分别为：N=10，N=20，N=50，N=100。 N=10; w=boxcar(N); n=[0:N-1]; [H,W]=dtft(w,1024); subplot(411) plot(W/2/pi,abs(H)); title('N=10') xlabel('f'),ylabel('|W(w)|')

N=20; w=boxcar(N); n=[0:N-1]; [H,W]=dtft(w,1024); subplot(412) plot(W/2/pi,abs(H)); title('N=20') xlabel('f'),ylabel('|W(w)|')

N=50; w=boxcar(N); n=[0:N-1]; [H,W]=dtft(w,1024); subplot(413) plot(W/2/pi,abs(H)); title('N=50') xlabel('f'),ylabel('|W(w)|')

N=100; w=boxcar(N); n=[0:N-1]; [H,W]=dtft(w,1024); subplot(414) plot(W/2/pi,abs(H)); title('N=100') xlabel('f'),ylabel('|W(w)|')

4.已知周期信号x(t)=0.75+3.4cos2pift+2.7cos4pift+1.5sin3.5pift+2.5sin7pift,其中f=(25/16)Hz，若截断时间长度分别为信号周期的0.9倍和1.1倍，试绘制和比较采用下面窗函数提取的x(t)的频谱。 经计算，本题中周期T=2.56s %0.9倍周期 f=25/16; dt=0.01; N=230; n=0:(N-1); x=0.75+3.4*cos(2*pi*f*dt*n)+2.7*cos(4*pi*f*dt*n)+1.5*sin(3.5*pi*f*dt*n)+2.5*sin(7*pi*f*dt*n);

figure(1) w1=boxcar(N); y1=x.*w1'; [Y1,W]=dtft(y1,1000); subplot(121),plot(W/2/pi,abs(Y1)); grid,title('信号周期0.9的矩形窗（幅频）'); xlabel('f');ylabel('|Y1|'); subplot(122),plot(W/2/pi,angle(Y1)); grid,title('信号周期0.9的矩形窗（相频）'); xlabel('f');ylabel('

figure(2) w2=hanning(N); y2=x.*w2'; [Y2,W]=dtft(y2,1000); subplot(121),plot(W/2/pi,abs(Y2)); grid,title('信号周期0.9的汉宁窗（幅频）'); xlabel('f');ylabel('|Y2|'); subplot(122),plot(W/2/pi,angle(Y2)); grid,title('信号周期0.9的汉宁窗（相频）'); xlabel('f');ylabel('

figure(3) w3=hamming(N); y3=x.*w3'; [Y3,W]=dtft(y3,1000); subplot(121),plot(W/2/pi,abs(Y3)); grid,title('信号周期0.9的汉明窗（幅频）'); xlabel('f');ylabel('|Y3|'); subplot(122),plot(W/2/pi,angle(Y3)); grid,title('信号周期0.9的汉明窗（相频）'); xlabel('f');ylabel('

figure(4) w4=bartlett(N); y4=x.*w4'; [Y4,W]=dtft(y4,1000); subplot(121),plot(W/2/pi,abs(Y4)); grid,title('信号周期0.9的巴特利特窗（幅频）'); xlabel('f');ylabel('|Y4|'); subplot(122),plot(W/2/pi,angle(Y4)); grid,title('信号周期0.9的巴特利特窗（相频）'); xlabel('f');ylabel('

figure(5) w5=blackman(N); y5=x.*w5'; [Y5,W]=dtft(y5,1000); subplot(121),plot(W/2/pi,abs(Y5)); grid,title('信号周期0.9的布莱克曼窗（幅频）'); xlabel('f');ylabel('|Y5|'); subplot(122),plot(W/2/pi,angle(Y5)); grid,title('信号周期0.9的布莱克曼窗（相频）'); xlabel('f');ylabel('

figure(6) w6=triang(N); y6=x.*w2'; [Y6,W]=dtft(y6,1000);

合集下载

实验报告质谱分析实验

实验报告质谱分析实验

实验报告质谱分析实验实验报告：质谱分析实验一、实验目的质谱分析作为一种强大的分析技术，在化学、生物、医药等领域有着广泛的应用。

本次实验的主要目的是通过实际操作，熟悉质谱分析的基本原理和实验流程，掌握仪器的使用方法，能够对给定的样品进行准确的质谱分析，并根据所得数据进行物质的鉴定和结构解析。

二、实验原理质谱分析是通过将样品分子转化为离子，并根据其质荷比（m/z）的不同在电磁场中进行分离和检测的方法。

样品首先经过离子化源，常见的离子化方式有电子轰击电离（EI）、化学电离（CI）、电喷雾电离（ESI）等。

离子化后的样品离子在加速电场中获得一定的动能，进入质量分析器。

质量分析器的作用是根据离子的质荷比将其分离。

常见的质量分析器有扇形磁场分析器、四极杆分析器、飞行时间分析器（TOF）等。

分离后的离子最终到达检测器，产生电信号，经过放大和处理后得到质谱图。

质谱图是以质荷比为横坐标，离子强度为纵坐标绘制的曲线。

通过对质谱图的分析，可以确定样品的分子量、化学式、结构等信息。

三、实验仪器与试剂1、仪器质谱仪（型号：＿____）进样系统（包括自动进样器、注射器等）数据处理系统（计算机及相关软件）2、试剂标准样品（已知分子量和结构的化合物，如苯、甲苯等）待测试样（未知组成的混合物）溶剂（如甲醇、乙腈等，用于溶解样品）四、实验步骤1、仪器准备开启质谱仪电源，预热至稳定状态。

检查仪器的真空度，确保达到工作要求。

校准仪器，使用标准物质进行质量轴的校准。

2、样品制备准确称取一定量的标准样品或待测试样，用适当的溶剂溶解并配制一定浓度的溶液。

将样品溶液转移至进样瓶中。

3、进样与分析设置进样参数，如进样量、进样速度等。

启动进样程序，将样品引入质谱仪进行分析。

4、数据采集与处理在分析过程中，仪器自动采集质谱数据。

利用数据处理软件对采集到的数据进行处理，如平滑、基线校正、峰识别等。

5、结果分析根据质谱图中的峰位置（质荷比）和峰强度，确定样品中所含物质的分子量和相对含量。

FFT算法分析实验实验报告

FFT算法分析实验实验报告

FFT算法分析实验实验报告一、实验目的快速傅里叶变换（Fast Fourier Transform，FFT）是数字信号处理中一种非常重要的算法。

本次实验的目的在于深入理解 FFT 算法的基本原理、性能特点，并通过实际编程实现和实验数据分析，掌握 FFT 算法在频谱分析中的应用。

二、实验原理FFT 算法是离散傅里叶变换（Discrete Fourier Transform，DFT）的快速计算方法。

DFT 的定义为：对于长度为 N 的序列 x(n)，其 DFT 为X(k) ＝∑n=0 到 N-1 x(n) e^（j 2π k n ／ N) ，其中 j 为虚数单位。

FFT 算法基于分治法的思想，将 N 点 DFT 分解为多个较小规模的DFT，从而大大减少了计算量。

常见的 FFT 算法有基 2 算法、基 4 算法等。

三、实验环境本次实验使用的编程语言为 Python，主要依赖 numpy 库来实现 FFT 计算和相关的数据处理。

四、实验步骤1、生成测试信号首先，生成一个包含不同频率成分的正弦波叠加信号，例如100Hz、200Hz 和 300Hz 的正弦波。

设定采样频率为 1000Hz，采样时间为 1 秒，以获取足够的采样点进行分析。

2、进行 FFT 计算使用 numpy 库中的 fft 函数对生成的测试信号进行 FFT 变换。

3、频谱分析计算 FFT 结果的幅度谱和相位谱。

通过幅度谱确定信号中各个频率成分的强度。

4、误差分析与理论上的频率成分进行对比，计算误差。

五、实验结果与分析1、幅度谱分析观察到在 100Hz、200Hz 和 300Hz 附近出现明显的峰值，对应于生成信号中的频率成分。

峰值的大小反映了相应频率成分的强度。

2、相位谱分析相位谱显示了各个频率成分的相位信息。

3、误差分析计算得到的频率与理论值相比，存在一定的误差，但在可接受范围内。

误差主要来源于采样过程中的量化误差以及 FFT 算法本身的近似处理。

光谱分析实验实验报告(3篇)

光谱分析实验实验报告(3篇)

第1篇一、实验目的和要求通过本次实验，掌握光谱分析的基本原理和方法，了解不同光谱仪（如紫外-可见分光光度计、傅里叶变换红外光谱仪、荧光光谱仪等）的原理和操作步骤。

学会如何通过光谱分析技术来鉴定物质、研究物质的组成和结构，并分析实验过程中可能影响结果的因素。

二、实验原理光谱分析是一种基于物质对电磁辐射吸收、发射或散射特性的分析方法。

当物质与电磁波相互作用时，会发生能量转移，从而产生吸收、发射或散射现象。

通过分析这些现象，可以获得有关物质的定量和定性信息。

1. 紫外-可见分光光度计：基于物质对紫外和可见光的吸收特性，通过测量吸光度来定量分析物质的浓度。

2. 傅里叶变换红外光谱仪（FTIR）：基于物质对红外光的吸收特性，通过分析红外光谱中的吸收峰来鉴定物质的结构。

3. 荧光光谱仪：基于物质对紫外光的吸收和荧光发射特性，通过分析荧光光谱来研究物质的性质。

三、主要仪器设备1. 紫外-可见分光光度计2. 傅里叶变换红外光谱仪（FTIR）3. 荧光光谱仪4. 标准样品5. 待测样品6. 空白溶液四、实验内容和原理1. 紫外-可见分光光度计实验- 原理：根据比尔-朗伯定律，吸光度与物质的浓度成正比。

- 步骤：配制标准溶液，测量吸光度，绘制标准曲线，测定待测样品的浓度。

2. 傅里叶变换红外光谱仪（FTIR）实验- 原理：根据红外光谱的吸收峰位置和强度，鉴定物质的结构。

- 步骤：将待测样品制成薄片，进行红外光谱扫描，与标准光谱图进行比对，鉴定物质的结构。

3. 荧光光谱仪实验- 原理：根据物质的荧光发射光谱，研究物质的性质。

- 步骤：将待测样品制成薄片，进行荧光光谱扫描，分析荧光光谱，研究物质的性质。

五、实验数据记录和处理1. 紫外-可见分光光度计实验数据：- 标准溶液浓度：C1, C2, C3, ...- 吸光度：A1, A2, A3, ...- 标准曲线：y = ax + b2. 傅里叶变换红外光谱仪（FTIR）实验数据：- 待测样品红外光谱图3. 荧光光谱仪实验数据：- 待测样品荧光光谱图六、实验结果与分析1. 紫外-可见分光光度计实验结果：- 标准曲线线性良好，相关系数R² > 0.99。

气相色谱法实验报告

气相色谱法实验报告

气相色谱法实验报告色谱和光谱实验实验5-气相色谱实验姓氏:张瑞芳薛浩:12月30日，XXXX第二组色谱和光谱实验色谱和光谱实验气相色谱实验1，实验目的299了解气相色谱仪各部件的功能2.加深对气相色谱原理和应用的理解 3.掌握气相色谱分析的一般实验方法 4.学习使用氢火焰离子化气相色谱分析未知物质2，实验原理1。

气相色谱的基本原理气相色谱的流动方向是惰性气体。

具有大表面积和一定活性的吸附剂在气固色谱中用作固定相。

当多组分混合样品进入色谱柱时，由于吸附剂对各组分的吸附力不同，经过一定时间后，色谱柱中各组分的运行速度不同。

吸附力弱的组分容易解吸，先离开色谱柱进入检测器，而吸附力强的组分最不容易解吸，最后离开色谱柱这样，组分可以在色谱柱中相互分离，然后依次进入检测器进行检测和记录。

气相色谱仪的框图见图1:图1。

气相色谱仪仪器的框图由以下五个系统组成:气路、进样、分离、温度控制、检测和记录系统2。

气相色谱定性和定量分析的原理通常用于描述样品中各成分的浓度换句话说，每个成分的分离光谱带的浓度变化被输入到能量转换装置中，并被转换成电信号的变化然后，电信号的变化被输入到记录器以进行记录，并且获得如图2所示的曲线。

表示成分进入检测器进行色谱检测和光谱实验后，检测器给出的信号随时间的变化规律。

它是柱内组分分离结果的反映，是研究色谱分离过程机理的基础，也是定性和定量的基础。

图2。

典型色谱流动曲线3。

FID原理本试验使用氢火焰离子化检测器(FID)，该检测器使用氢火焰和空气燃烧作为能源，使用含碳有机物在火焰中燃烧产生离子，在外加电场的作用下，离子形成离子流，并根据离子流产生的电信号强度检测色谱柱分离的成分iii。

实验试剂和仪器(1)试剂:甲醇、异丙醇、异丁醇(2)仪器:带氢火焰离子化检测器的气相色谱仪(GC-2014气相色谱仪)；氢-空气发生器(SPH-300氢发生器)，氮气瓶；色谱柱；微量注射器4。

实验步骤1。

打开稳定的电源2。