信号进行频谱分析和滤波处理的程序
信号分析处理流程
信号分析处理流程下载温馨提示:该文档是我店铺精心编制而成，希望大家下载以后，能够帮助大家解决实际的问题。
文档下载后可定制随意修改，请根据实际需要进行相应的调整和使用，谢谢!Download Tip: This document has been carefully written by the editor. I hope that after you download, they can help you solve practical problems. After downloading, the document can be customized and modified. Please adjust and use it according to actual needs. Thank you!信号分析处理流程：①信号采集：利用传感器、接收器等设备捕捉待分析的物理量或数据流，转换为电信号或数字信号形式。
②预处理：对采集的原始信号进行滤波（如去噪、低通、高通滤波），消除干扰成分，以及放大、偏移调整等，改善信号质量。
③信号分割：根据信号特性或分析需求，将连续信号分割成有意义的片段，便于后续分析，如语音识别中的单词分割。
④特征提取：从信号中提取有助于分析的关键特征，如频谱分析、时域统计特征、波形形状参数等，为后续处理提供基础。
⑤数据分析：运用统计学、信号处理算法分析特征，识别模式、趋势或异常，如频谱分析识别频率成分，时序分析预测趋势。
⑥信号重构与压缩：根据分析目的，可能需要对信号进行重构，如降噪后的信号重建，或为存储传输需求进行信号压缩。
⑦模式识别与分类：利用机器学习、深度学习等技术，基于提取的特征对信号进行分类或识别，如语音识别、图像分类。
⑧结果验证与优化：对比分析结果与实际情况或预期目标，评估处理效果，根据误差反馈调整算法参数，优化处理流程。
⑨报告生成与应用：整理分析结果，生成报告或可视化展示，为决策支持、故障诊断、系统优化等实际应用提供依据。
工程信号处理实验
工程信号处理实验1. 实验目的本实验旨在通过处理工程信号，掌握信号处理的基本原理和方法，提高对信号处理算法的理解与应用能力。
2. 实验背景工程信号处理是一门研究如何对信号进行采集、传输、存储、处理和分析的学科。
在实际工程应用中，信号处理技术被广泛应用于通信、图像处理、音频处理、生物医学工程等领域。
通过本实验，我们将学习信号采集、滤波、频谱分析等基本信号处理技术。
3. 实验设备与材料- 信号发生器- 示波器- 模拟滤波器- 数据采集卡- 个人电脑- 实验电路板- 信号处理软件4. 实验步骤4.1 信号发生器设置将信号发生器与示波器连接，设置信号发生器的频率、幅度和波形类型，例如正弦波、方波或三角波。
4.2 信号采集将信号发生器输出的信号通过模拟滤波器，然后连接到数据采集卡。
使用数据采集卡将信号转换为数字信号，并传输到个人电脑。
4.3 信号滤波在个人电脑上使用信号处理软件，对采集到的信号进行滤波处理。
选择适当的滤波器类型（低通、高通、带通或带阻）和截止频率，根据需要对信号进行滤波操作。
4.4 频谱分析利用信号处理软件进行频谱分析，将滤波后的信号进行频谱转换，得到信号的频谱图。
观察频谱图，分析信号的频谱特征。
4.5 实验结果分析根据实验得到的数据和频谱图，进行结果分析。
对滤波效果进行评估，比较不同滤波器类型和截止频率对信号的影响。
5. 实验注意事项- 操作仪器时应注意安全，避免电流过大或触电等意外情况发生。
- 操作示波器时，应先调整示波器的触发、增益和时间基准等参数，确保信号能够正常显示。
- 在进行信号滤波和频谱分析时，应根据实际需要选择合适的滤波器类型和截止频率。
- 实验结束后，应将仪器设备恢复到初始状态，保持实验环境整洁。
6. 实验结果及讨论根据实验数据和频谱图，我们可以看到滤波器对信号的影响。
不同类型的滤波器可以对信号的频率进行选择性的增强或衰减，从而实现对信号的滤波处理。
通过频谱分析，我们可以观察到信号的频谱特征，了解信号的频率分布情况。
用FFT对信号作频谱分析
用FFT对信号作频谱分析快速傅立叶变换（FFT）是一种在信号处理中常用于频谱分析的方法。
它是傅立叶变换的一种快速算法，通过将信号从时间域转换到频域，可以提取信号的频率信息。
FFT算法的原理是将信号分解为不同频率的正弦波成分，并计算每个频率成分的幅度和相位。
具体而言，FFT将信号划分为一系列时间窗口，每个窗口内的信号被认为是一个周期性信号，然后对每个窗口内的信号进行傅立叶变换。
使用FFT进行频谱分析可以得到信号的频率分布情况。
频谱可以显示信号中各个频率成分的强度。
通过分析频谱可以识别信号中的主要频率成分，判断信号中是否存在特定频率的干扰或噪声。
常见的应用包括音频信号处理、图像处理、通信系统中的滤波和解调等。
使用FFT进行频谱分析的步骤如下：1.首先，获取待分析的信号，并确保信号是离散的，即采样频率与信号中的最高频率成分满足奈奎斯特采样定理。
2.对信号进行预处理，包括去除直流分量和任何不需要的干扰信号。
3.对信号进行分段，分段后的每个窗口长度在FFT算法中通常为2的幂次方。
常见的窗口函数包括矩形窗、汉明窗等。
4.对每个窗口内的信号应用FFT算法，将信号从时间域转换到频域，并计算每个频率成分的幅度和相位。
5.对所有窗口得到的频谱进行平均处理，以得到最终的频谱分布。
在使用FFT进行频谱分析时需要注意的问题有：1.噪声的影响：FFT对噪声敏感，噪声会引入幅度偏差和频率漂移。
可以通过加窗等方法来减小噪声的影响。
2.分辨率的选择：分辨率是指在频谱中能够分辨的最小频率间隔。
分辨率与信号长度和采样频率有关，需要根据需求进行选择。
3.漏泄效应：当信号中的频率不是FFT长度的整数倍时，会出现漏泄效应。
可以通过零填充等方法来减小漏泄效应。
4.能量泄露：FFT将信号限定在一个周期内进行计算，如果信号过长，则可能导致部分频率成分的能量泄露到其他频率上。
总之，FFT作为信号处理中常用的频谱分析方法，能够提取信号中的频率信息，广泛应用于多个领域。
通信系统中的频谱分析和信号调整
通信系统中的频谱分析和信号调整随着科技的飞速发展，通信系统在我们的日常生活中扮演着越来越重要的角色。
而频谱分析和信号调整是保证通信系统有效运行的关键步骤。
下面将详细介绍频谱分析和信号调整的步骤。
一、频谱分析1. 理解频谱：频谱是指信号在频率域上的分布情况，是衡量信号特性的重要指标。
通过频谱分析，可以了解信号的频率分布和能量分布等信息。
2. 具体步骤：(1) 采集信号样本：利用合适的设备和传感器，采集到待分析的信号样本。
可以采用天线接收无线电波、音频设备记录声音等方式。
(2) 信号预处理：对采集到的信号进行预处理，包括滤波、降噪等操作，以提高分析准确度。
(3) 时域转频域：将预处理后的信号从时域转换为频域，通过傅里叶变换等数学方法，得到信号的频谱。
(4) 分析与解读：对得到的频谱进行分析和解读，包括观察峰值、能量分布等信息，以了解信号的特性和属性。
(5) 结果报告：将分析结果进行整理和报告，记录信号的频谱特征，并进行归档和备份。
二、信号调整1. 理解信号调整：信号调整是指在通信系统中对信号进行处理和调整，以满足通信需求和提高信号的质量和可靠性。
2. 具体步骤：(1) 信号检测：利用相应的设备检测信号的强度、误码率等参数，获取关于信号的基本信息。
(2) 信号分析：对检测到的信号进行分析，包括频率分析、时域分析等，以了解信号的特性和问题。
(3) 问题定位：根据信号分析的结果，确定信号出现的问题和原因，如信号干扰、衰减等。
(4) 信号优化：根据问题定位，采取相应的优化措施，如增加信号发射功率、优化天线布局等，以提高信号的质量和可靠性。
(5) 重复检测与调整：在进行信号调整后，再次进行信号检测和分析，以验证调整效果，并根据需要进行进一步的调整。
(6) 结果验证与记录：对信号调整的结果进行验证，对调整后的信号进行监测和记录，以备后续分析和参考。
以上就是通信系统中的频谱分析和信号调整的详细步骤。
通过频谱分析，可以了解信号的频率分布和能量分布等特性；而信号调整则是为了满足通信需求和提高信号质量与可靠性。
使用MATLAB进行信号处理与滤波
使用MATLAB进行信号处理与滤波信号处理与滤波是数字信号处理领域中的重要技术，而MATLAB是一种广泛应用于信号处理的工具。
本文将介绍如何使用MATLAB进行信号处理与滤波，包括信号采样、信号重构、频谱分析以及常用的滤波器设计和应用。
首先，我们先了解一下信号处理的基本概念。
信号处理是对信号进行采样、重构、滤波、增强、压缩等操作的过程。
信号可以是连续的或离散的，常常通过采样将连续信号转换为离散信号进行处理。
在MATLAB中，可以使用`plot`函数来绘制信号的波形图。
假设有一个正弦信号，我们可以通过以下代码绘制其波形图：```matlabfs = 1000; % 采样率为1000Hzt = 0:1/fs:1; % 时间向量，从0到1sf = 10; % 正弦信号的频率为10Hzx = sin(2*pi*f*t); % 构造正弦信号plot(t, x); % 绘制波形图xlabel('Time (s)'); % x轴标签ylabel('Amplitude'); % y轴标签title('Sinusoidal Signal'); % 图片标题```这段代码中，首先定义了采样率`fs`、时间向量`t`和信号频率`f`，然后使用`sin`函数构造了正弦信号`x`，最后通过`plot`函数绘制出信号的波形图。
在进行信号处理时，经常需要进行频谱分析来研究信号的频率特性。
MATLAB 提供了多种函数来计算信号的频谱，其中最常用的是`fft`函数。
以下代码演示了如何计算信号的频谱，并绘制频谱图：```matlabFs = 1000; % 采样率为1000HzT = 1/Fs; % 采样间隔L = 1000; % 信号长度为1000t = (0:L-1)*T; % 时间向量x = sin(2*pi*50*t) + 0.5*sin(2*pi*120*t); % 构造含有两个频率成分的信号Y = fft(x); % 对信号进行傅里叶变换P2 = abs(Y/L); % 计算双边频谱P1 = P2(1:L/2+1); % 取单边频谱P1(2:end-1) = 2*P1(2:end-1); % 幅度归一化f = Fs*(0:(L/2))/L; % 频率向量plot(f,P1); % 绘制频谱图title('Single-Sided Amplitude Spectrum of x(t)'); % 图片标题xlabel('Frequency (Hz)'); % x轴标签ylabel('Amplitude'); % y轴标签```这段代码中，首先定义了采样率`Fs`、采样间隔`T`、信号长度`L`和时间向量`t`，然后使用两个正弦信号相加的方式构造了含有两个频率成分的信号`x`，接着使用`fft`函数对信号进行傅里叶变换得到频谱`Y`，最后根据频谱进行幅度归一化并绘制频谱图。
信号分析与处理实验报告
实验一图像信号频谱分析及滤波一：实验原理FFT不是一种新的变化，而是DFT的快速算法。
快速傅里叶变换能减少运算量的根本原因在于它不断地把长序列的离散傅里叶变换变为短序列的离散傅里叶变换，在利用的对称性和周期性使DFT运算中的有些项加以合并，达到减少运算工作量的效果。
为了消除或减弱噪声，提取有用信号，必须进行滤波，能实现滤波功能的系统成为滤波器。
按信号可分为模拟滤波器和数字滤波器两大类。
数字滤波器的关键是如何根据给定的技术指标来得到可以实现的系统函数。
从模拟到数字的转换方法很多，常用的有双线性变换法和冲击响应不变法，本实验主要采用双线性变换法。
双线性变换法是一种由s平面到z平面的映射过程，其变换式定义为：数字域频率与模拟频率之间的关系是非线性关系。
双线性变换的频率标度的非线性失真是可以通过预畸变的方法去补偿的。
变换公式有Ωp=2/T*tan(wp/2)Ωs=2/T*tan(ws/2)二：实验内容1．图像信号的采集和显示选择一副不同彩色图片，利用Windows下的画图工具，设置成200*200像素格式。
然后在Matlab软件平台下，利用相关函数读取数据和显示图像。
要求显示出原始灰度图像、加入噪声信号后的灰度图像、滤波后的灰度图像。
2．图像信号的频谱分析要求分析和画出原始灰度图像、加入噪声信号后灰度图像、滤波后灰度图像信号的频谱特性。
3.数字滤波器设计给出数字低通滤波器性能指标:通带截止频率fp＝10000 Hz，阻带截止频率fs＝15000 Hz，阻带最小衰减Rs＝50 dB，通带最大衰减Rp＝3 dB，采样频率40000Hz。
三：实验程序clear allx=imread('D:\lan.jpg');%原始彩色图像的数据读取x1=rgb2gray(x);%彩色图像值转化为灰度图像值[M,N]=size(x1);%数据x1的长度,用来求矩阵的大小x2=im2double(x1);%unit8转化为double型x3=numel(x2);%计算x2长度figure(1);subplot(1,3,1);imshow(x2);title('原始灰度图')z1=reshape(x2,1,x3);%将二维数据转化成一维数据g=fft(z1);%对图像进行二维傅里叶变换mag=fftshift(abs(g));%fftshift是针对频域的，将FFT的DC分量移到频谱中心K=40000;Fs=40000;dt=1/Fs;n=0:K-1;f1=18000;z=0.1*sin(2*pi*f1*n*dt);x4=z1+z;%加入正弦噪声f=n*Fs/K;y=fft(x4,K);z2=reshape(x4,M,N);%将一维图转换为二维图subplot(1,3,2);imshow(z2);title('加入噪声后')g1=fft(x4);mag1=fftshift(abs(g1));%设计滤波器ws=0.75*pi;wp=0.5*pi;fs=10000;wp1=2*fs*tan(wp/2);ws1=2*fs*tan(ws/2);rs=50;rp=3;% [n,wn]=buttord(wp/pi,ws/pi,rp,rs);% [bz,az]=butter(n,wn);[n,wn]=buttord(wp1,ws1,rp,rs,'s');[z,p,k]=buttap(n);[b,a]=zp2tf(z,p,k);[B,A]=lp2lp(b,a,wn);[bz,az]=bilinear(B,A,fs);[h,w]=freqz(bz,az,128,fs);L=numel(z2);z3=reshape(z2,1,L);x6=filter(bz,az,double(z3));x7=reshape(x6,M,N);subplot(1,3,3);imshow(x7);g2=fft(x6);mag2=fftshift(abs(g2));title('滤波后')%建立频谱图figure(2);subplot(1,3,1);plot(mag);title('原始Magnitude')subplot(1,3,2);plot(mag1);title('加噪声Magnitude')subplot(1,3,3);plot(mag2);title('滤波后Magnitude')figure(3);subplot(1,2,1)plot(w,abs(h));xlabel('f');ylabel('h');title('滤波器幅谱');subplot(1,2,2);plot(w,angle(h));title('滤波器相谱');四：实验结果与分析图一图二分析：由图二可以知道加入噪声后的幅值谱和原始图的幅值谱明显多了两条幅值线，而这两条幅值线就是我们对原始灰度图加入的正弦噪声，而相应的图一中的加噪声后的图与原始图相比，出现了明显的变化。
无线电频谱分析仪的工作原理与应用
无线电频谱分析仪的工作原理与应用无线电频谱分析仪是一种用于测量和分析无线电频谱的仪器。
它可以实时显示频谱，帮助工程师了解无线电信号的特征及其在各个频率范围内的分布情况。
本文将介绍无线电频谱分析仪的工作原理以及在不同领域的应用。
一、工作原理无线电频谱分析仪的工作原理可以简单概括为以下几个步骤：1. 信号接收：无线电频谱分析仪通过内置或外接天线接收到要分析的无线电信号。
2. 信号放大：接收到的信号经过前置放大电路进行信号放大，以提高信号的幅度和灵敏度。
3. 信号混频：经过放大后的信号和本地振荡器产生的中频信号进行混频操作，得到中频信号。
4. 信号滤波：对混频得到的中频信号进行滤波，去除不需要的频率成分，以便进行后续的频谱分析。
5. 信号解调：对滤波后的中频信号进行解调，恢复信号的原始调制方式，如调幅、调频等。
6. 信号转换：将解调后的信号转换为数字信号，以便进行数字信号处理和显示。
7. 数字信号处理：使用数字信号处理技术对信号进行频谱分析、频谱显示和信号参数计算等操作。
8. 频谱显示：将处理后的信号转换为频谱图形并显示在仪器的显示屏上，供用户查看和分析。
二、应用领域无线电频谱分析仪在多个领域有着广泛的应用，以下将介绍其中几个主要的应用领域。
1. 电信领域：无线电频谱分析仪在电信领域中被广泛应用于对无线电信号进行调制解调、频谱分析、调频定位、无线电干扰监测等工作。
它可以帮助工程师更好地分析和监测无线电信号的质量以及各种干扰情况，从而保证通信系统的正常运行。
2. 广播电视领域：广播电视频谱分析是保障广播电视信号质量的重要手段之一。
无线电频谱分析仪可以帮助广播电视工程师进行频谱监测、频谱规划以及无线电干扰分析等工作，从而提高广播电视信号的传输质量和覆盖范围。
3. 电子设备测试领域：在电子设备测试领域中，无线电频谱分析仪可以用于对设备的射频性能进行测试和分析。
通过对设备发出的无线电信号进行频谱分析，工程师可以了解到设备的发射功率、频率稳定性、谐波等参数，从而评估设备的性能和合格性。
使用LabVIEW进行信号处理实现信号滤波和频谱分析
使用LabVIEW进行信号处理实现信号滤波和频谱分析信号处理在科学研究和工程应用中扮演着重要的角色。
信号滤波和频谱分析是信号处理的两个基本任务，而LabVIEW是一款功能强大的可视化编程环境，适合用于信号处理的实现。
本文将介绍如何使用LabVIEW进行信号滤波和频谱分析的步骤和方法。
一、信号滤波信号滤波是通过改变信号的频率特性，去除不需要的频率成分，使得信号更加清晰和准确。
LabVIEW提供了多种滤波器模块，可以方便地实现信号滤波的功能。
1. 数据获取首先，需要从外部设备或者文件中获取待处理的信号。
LabVIEW 提供了多种数据采集模块，可以选择合适的模块进行数据获取。
2. 滤波器设计在信号滤波过程中，首先需要设计滤波器。
LabVIEW中的滤波器设计模块可以根据具体需求选择滤波器类型，并进行参数设置。
根据信号的特性和应用要求，可以选择低通、高通、带通或带阻滤波器。
3. 滤波器应用设计好滤波器后，需要将其应用到待处理的信号上。
LabVIEW提供了滤波器模块，可以直接调用已设计好的滤波器进行信号滤波。
4. 数据输出滤波后的信号经过处理后，可以将结果输出到显示模块或者保存到文件中，以便后续分析或应用。
二、频谱分析频谱分析是对信号进行频域分析，得到信号的频率分布和功率谱等信息。
LabVIEW提供了丰富的频谱分析工具和函数，可以方便地进行频谱分析。
1. 数据获取首先，需要获取待分析的信号数据。
可以利用LabVIEW的数据采集模块或者导入外部文件的方式获取数据。
2. 数据预处理在进行频谱分析之前，有时需要对数据进行预处理，例如去除噪声、降低采样率等。
LabVIEW提供了多种数据处理函数和模块，可以方便地进行数据预处理。
3. 频谱分析LabVIEW中的频谱分析模块可以对信号进行快速傅里叶变换（FFT）或者其他频谱分析算法。
可以选择合适的分析模块，并进行参数设置，如分辨率、窗函数等。
4. 结果展示频谱分析完成后，可以将结果以图表、曲线等形式展示出来，使得分析结果更加直观和易于理解。
如何使用Matlab进行信号处理和滤波
如何使用Matlab进行信号处理和滤波信号处理和滤波在工程领域中扮演着重要的角色，它们可以帮助我们从一系列的数据中提取有用的信息，并消除噪声。
Matlab作为一种强大的工具，提供了丰富的函数和工具箱，可以方便地进行信号处理和滤波。
本文将介绍如何使用Matlab进行信号处理和滤波的基本方法，并使用实例进行演示。
一、Matlab的信号处理工具箱Matlab的信号处理工具箱是一个强大的工具集，它包含了许多用于处理各种类型信号的函数和算法。
通过引入信号处理工具箱，我们可以方便地处理音频、图像和视频信号，并进行频域分析、滤波和解调等操作。
在Matlab中，可以使用命令"toolbox"来查看已安装的工具箱，对于信号处理，我们需要确保已经安装了"Signal Processing Toolbox"。
如果没有安装，可以通过访问Matlab官方网站下载并安装。
二、信号处理的基本操作1. 读取和显示信号在进行信号处理之前，首先需要将信号加载进Matlab中。
可以使用函数"audioread"来读取音频信号，例如读取一个.wav格式的音频文件：```[x,Fs] = audioread('audio.wav');```其中，x是音频信号的数据，Fs是信号的采样率。
读取完成后，可以使用函数"soundsc"来播放信号，并使用函数"plot"来绘制信号的波形图：```soundsc(x,Fs);plot(x);```2. 频谱分析频谱分析可以帮助我们了解信号的频率特性。
在Matlab中，可以使用函数"fft"进行快速傅里叶变换（FFT），将信号从时域转换到频域。
例如，对于上文中读取的音频信号x，可以使用以下代码计算其频谱：```X = fft(x);```频谱的结果是一个复数向量，表示信号在不同频率上的幅值和相位。
应用MATLAB对信号进行频谱分析及滤波
应用MATLAB对信号进行频谱分析及滤波频谱分析和滤波是信号处理中常用的技术，可以帮助我们了解信号的频率特性并对信号进行去噪或增强。
MATLAB是一个强大的数学计算和工程仿真软件，提供了各种工具和函数用于频谱分析和滤波。
频谱分析是通过将信号在频域上进行分解来研究信号的频率特性。
MATLAB提供了几种进行频谱分析的函数，包括FFT（快速傅里叶变换）、periodogram和spectrogram等。
下面将以FFT为例，介绍如何使用MATLAB进行频谱分析。
首先，我们需要先生成一个信号用于频谱分析。
可以使用MATLAB提供的随机信号生成函数来生成一个特定频率和幅度的信号。
例如，可以使用以下代码生成一个包含两个频率成分的信号：```MATLABFs=1000;%采样率t=0:1/Fs:1;%时间向量，从0秒到1秒，采样率为Fsf1=10;%第一个频率成分f2=50;%第二个频率成分A1=1;%第一个频率成分的幅度A2=0.5;%第二个频率成分的幅度x = A1*sin(2*pi*f1*t) + A2*sin(2*pi*f2*t);```上述代码生成了一个采样率为1000Hz的信号，包含10Hz和50Hz两个频率的成分。
接下来，我们可以使用MATLAB的FFT函数对信号进行频谱分析，并将频谱绘制出来。
FFT函数将信号从时域转换到频域，并返回频谱幅度和频率信息。
以下是使用FFT函数对上述生成的信号进行频谱分析的代码：```MATLABN = length(x); % 信号长度X = abs(fft(x))/N; % 计算FFTf=(0:N-1)*(Fs/N);%计算频率坐标plot(f,X)xlabel('频率(Hz)')ylabel('幅度')title('信号频谱')```上述代码中，我们首先计算FFT并将结果除以信号长度，以得到正确的幅度值。
然后，我们计算频率坐标，并将频谱幅度与频率绘制出来。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
关于用MATLAB设计对信号进行频谱分析和滤波处理的程序设计出一套完整的系统，对信号进行频谱分析和滤波处理；1．产生一个连续信号，包含低频，中频，高频分量，对其进行采样，进行频谱分析，分别设计三种高通，低通，带通滤波器对信号进行滤波处理，观察滤波后信号的频谱。
2.采集一段含有噪音的语音信号（可以录制含有噪音的信号，或者录制语音后再加进噪音信号），对其进行采样和频谱分析，根据分析结果设计出一合适的滤波器滤除噪音信号。
只要有程序就可以了，谢谢各位了，提供点帮助喽！！没学过MATLAB，自学的，要做很难的这个就有点困难啦！！这两个题目随便哪个都可以.先谢谢啦!!好的答案会追加分的，目前我只有这么多了!完整的程序%写上标题%设计低通滤波器：[N,Wc]=buttord()%估算得到Butterworth低通滤波器的最小阶数N和3dB截止频率Wc[a,b]=butter(N,Wc); %设计Butterworth低通滤波器[h,f]=freqz(); %求数字低通滤波器的频率响应figure(2); % 打开窗口2subplot(221); %图形显示分割窗口plot(f,abs(h)); %绘制Butterworth低通滤波器的幅频响应图title(巴氏低通滤波器'');grid; %绘制带网格的图像sf=filter(a,b,s); %叠加函数S经过低通滤波器以后的新函数subplot(222);plot(t,sf); %绘制叠加函数S经过低通滤波器以后的时域图形xlabel('时间(seconds)');ylabel('时间按幅度');SF=fft(sf,256); %对叠加函数S经过低通滤波器以后的新函数进行256点的基—2快速傅立叶变换w= %新信号角频率subplot(223);plot()); %绘制叠加函数S经过低通滤波器以后的频谱图title('低通滤波后的频谱图');%设计高通滤波器[N,Wc]=buttord()%估算得到Butterworth高通滤波器的最小阶数N和3dB截止频率Wc[a,b]=butter(N,Wc,'high'); %设计Butterworth高通滤波器[h,f]=freqz(); %求数字高通滤波器的频率响应figure(3);subplot(221);plot()); %绘制Butterworth高通滤波器的幅频响应图title('巴氏高通滤波器');grid; %绘制带网格的图像sf=filter(); %叠加函数S经过高通滤波器以后的新函数subplot(222);plot(t,sf); ;%绘制叠加函数S经过高通滤波器以后的时域图形xlabel('Time(seconds)');ylabel('Time waveform');w; %新信号角频率subplot(223);plot()); %绘制叠加函数S经过高通滤波器以后的频谱图title('高通滤波后的频谱图');%设计带通滤波器[N,Wc]=buttord([)%估算得到Butterworth带通滤波器的最小阶数N和3dB截止频率Wc[a,b]=butter(N,Wc); %设计Butterworth带通滤波器[h,f]=freqz(); %求数字带通滤波器的频率响应figure(4);subplot(221);plot(f,abs(h)); %绘制Butterworth带通滤波器的幅频响应图title('butter bandpass filter');grid; %绘制带网格的图像sf=filter(a,b,s); %叠加函数S经过带通滤波器以后的新函数subplot(222);plot(t,sf); %绘制叠加函数S经过带通滤波器以后的时域图形xlabel('Time(seconds)');ylabel('Time waveform');SF=fft(); %对叠加函数S经过带通滤波器以后的新函数进行256点的基—2快速傅立叶变换w=( %新信号角频率subplot(223);plot(')); %绘制叠加函数S经过带通滤波器以后的频谱图title('带通滤波后的频谱图');光盘使用说明B.1光盘文件的结构在光盘上有如下五个文件夹：●matlab_api 专放第12章所涉及的各类文件。
●matlab_c 专放第11章所涉及的各类M、C、DLL、EXE文件。
●mbook 存放着包含本书全部算例的M-book形式DOC文档。
●mfiles 汇集了全书90%以上算例的M或MDL文件。
●PowerMatlab 存放着PowerPoint制作的科技演讲稿简例文件。
B.2光盘对软件环境的要求●需要（包含Word、PowerPoint的）Offiece2000和MATLAB6.5支持。
●假如要运行mbook文件夹上的DOC文件，则需要Notebook环境。
（关于Notebook环境的设置参见第13.1节）●假如要编译产生EXE文件，则需要与MATLAB6.5适配的Borland 、MicroSoft 、Visual 、Watcom C/C++（详见11.2和11.3节）。
●假若要制作API接口，则需要相应的配套软件。
B.3光盘文件的操作准备在运行光盘文件之前，应首先使用MATLAB的路径浏览器把 \mfiles , \ma tlab_c , \matlab_api等三个文件夹设置在MATLAB的搜索路径上。
这种设置可以是“永久”的，即今后再打开MATLAB时，总确认 \mfiles 文件夹等在搜索路径上。
假如在今后运行MATLAB时，光盘没有插入，那么在MATLAB指令窗中会出现警告性提示。
这对光盘以外的MATLAB其他运作没有任何影响。
B.4mbook文件夹上DOC文件的使用该文件夹包含第1章到第13章的全部算例，以及前言、附录等。
章节的编号、名称与印刷版完全一致。
光盘DOC文件都是在MATLAB 6.5的Notebook环境中生成的。
（1）光盘DOC文件的功用●弥补了印刷版丢失的色彩信息在MATLAB中，用M文件编辑器或Notebook编写的指令或文件运行的结果（尤其图形）都采用不同的色彩鲜明地表现对象特征。
但目前印刷版书籍出于价格和技术原因，不得不牺牲色彩信息而采用“黑白”处理。
读者借助本光盘可克服印刷版丢失色彩信息的遗憾。
●提供了与印刷版对应的Notebook演练环境本光盘中DOC文件的章节结构、算例编号与印刷版完全相同。
因此在学习过程中，读者可在本光盘启动的Notebook环境中，或直接运行算例，观察运行结果；或改变若干指令，举一反三地观察运行结果的变化；或通过简单的复制操作，使相应指令在MATLAB指令窗中运行，而避免自己键入的错误。
●提供制作M-book的样板读者若想制作自己的M-book文档，可调用本光盘任何一章的DOC文件为样板进行。
注意：假若调用原英文M-book模版，则需通过若干设置才能在中文状态下正常运行。
（2）DOC文件的开启●所有DOC文档都是在“Word2000中文版 + MATLAB6.1”构成的Notebook环境中生成的。
因此，在相同环境下开启是最佳选择。
此时，文档具有“活性”。
●假如读者的MATLAB与Word联接正确，用鼠标双击光盘上的DOC文件，就能直接进入Notebook 环境。
●若前述启动失败，可以先启动MATLAB，然后再启动光盘DOC文档。
在很多情况下，这样的操作很可能奏效。
●在不具备相同Notebook环境的情况下，DOC文档最好在Word2000环境中开启。
此时电子文档虽然不“活”，但所有指令都能被准确地复制。
●在DOC文件开启后，节次编号有可能显现得不很正确。
这是由Word软件本身问题引起的。
此时只要进行以下操作就可使编号正确：先使光标位于“章名”大标题上；然后点中工具条样式栏右边的下拉菜单按键，引出下拉菜单；在这下拉菜单中再点选“标题1”菜单项，就能自动整理节次编号，使之正确。
（3）光盘DOC文件的使用方法●作为演练环境使用在正常打开的光盘DOC文件中，读者只要把光标放在绿色的输入细胞内，按组合键【Ctrl + Enter】，就可使该输入细胞重新执行计算。
在演练中，读者可以通过对指令的修改、变化和重新运行，观察运算结果的变化，从而达到举一反三的效果。
●作为样板使用先打开光盘DOC文件，然后删去原光盘文件内容，再写入读者自己所需的内容，最后通过菜单项的“另存为”操作保存为自己的文件。
这样获得的文件能正常地在Notebook环境下工作，也就是既可以输入文字、公式，又可以运行MATLAB指令、嵌入数值或图形结果；既拥有Word的所有文字处理能力，又具备MATLAB的运算、表现能力。
B.5mfiles文件夹上的M、MDL文件的使用除不能用M、MDL文件表达的少数算例外，其余算例（占总数的90%以上）都以M文件或MDL文件形式刻录在光盘的 \mfiles 文件夹中。
本光盘M文件应在MATLAB6.5以上版本运行；涉及符号计算的应有Symboli c Math Toobox 2.1.3以上版本适配；MDL文件应有SIMULINK5.0以上版本适配。
对于其他版本，或较低版的工具包，有些文件的运行可能会失败。
但只要对个别指令稍加修改就可，有关这方面的叙述可参见相应印刷版。
（1）光盘M、MDL文件的功用●提供可直接运作的M源码文件只要有MATLAB环境，本光盘上的M文件就可以运行。
它的适用条件比\mbook 文件夹上的DOC文件宽松得多，也就是不管读者是否正确安装Note book，不管文件产生的是动画还是交互操作界面，它们都能在MATLAB环境中正确执行。
每个算例文件都是完整的，可在MATLAB环境中直接运行的，所得结果与印刷版相对应。
但出于运行方式不同的考虑，有些光盘M文件与印刷版文件指令可能会存在少许差别，目的是为把算例特征表现得更充分。
此外，本光盘提供的M文件中，有许多是很通用的，读者只要稍加修改，就可为己所用。
●弥补了印刷版没有SIMULINK模型文件的缺陷由于SIMULINK工作特点的缘故，所以迄今为止所有涉及SIMULINK的印刷版书籍中都没有能直接运行的模型文件。
这给读者带来许多困惑和麻烦：一，读者如想验证书中结论，那就不得不从建模做起；二，仿真模块中的参数设置常使初学者顾此失彼，而造成仿真失败。
本光盘上MDL模型文件都可直接在MATLAB中运行，进行验证。
用户也可以在模型打开后，修改参数，观察变化。