信处理工具箱最新课件

《信处理工具箱最新课件》由会员分享，可在线阅读，更多相关《信处理工具箱最新课件（55页珍藏版）》请在装配图网上搜索。

1、信号处理工具箱最新课件第八章 信号处理工具箱信号处理工具箱最新课件一、工具箱简介功能型工具箱 通用型 功能型工具箱主要用来扩充matlab的数值计算、符号运算功能、图形建模仿真功能、文字处理功能以及与硬件实时交互功能，能够用于多种学科。信号处理工具箱最新课件领域型工具箱 专用型 领域型工具箱是学科专用工具箱，其专业性很强，比如控制系统工具箱（ Control System Toolbox）；信号处理工具箱（Signal Processing Toolbox）；财政金融工具箱（ Financial Toolbox）等等。只适用于本专业。信号处理工具箱最新课件例如：控制系统工具箱连续系统设计和离散

2、系统设计状态空间和传递函数以及模型转换时域响应（脉冲响应、阶跃响应、斜坡响应）频域响应（Bode图、Nyquist图）根轨迹、极点配置信号处理工具箱最新课件优化工具箱线型规划和二次规划求函数的最大值和最小值多目标优化约束条件下的优化非线型方程求解信号处理工具箱最新课件信号处理工具箱数字和模拟滤波器设计、应用及仿真谱分析和估计FFT、DCT等变换参数化模型信号处理工具箱最新课件学科前沿最新的工具箱 模糊控制逻辑工具箱 友好的交互设计界面，自适应神经-模糊学习、聚类以及Sugeno推理 神经网络工具箱 神经网络系统具有集体运算的能力和自适应的学习能力。具有很强的容错性和鲁棒性，善于联想、综合和推广

3、。信号处理工具箱最新课件工具箱的扩充功能： 用户可以修改工具箱中的函数，更为重要的是用户可以通过编制 m文件来任意地添加工具箱中原来没有的工具函数。此功能充分体现了matlab语言的开发性。信号处理工具箱最新课件Matlab常用工具箱Matlab Main Toolboxmatlab主工具箱Control System Toolbox控制系统工具箱Communication Toolbox通讯工具箱Financial Toolbox财政金融工具箱System Identification Toolbox系统辨识工具箱Fuzzy Logic Toolbox模糊逻辑工具箱Higher-Order

4、Spectral Analysis Toolbox高阶谱分析工具箱信号处理工具箱最新课件Image Processing Toolbox图象处理工具箱LMI Control Toolbox线性矩阵不等式工具箱Model predictive Control Toolbox模型预测控制工具箱-Analysis and Synthesis Toolbox分析工具箱Neural Network Toolbox神经网络工具箱Optimization Toolbox优化工具箱Partial Differential Toolbox偏微分方程工具箱Robust Control Toolbox鲁棒控制工具箱

5、信号处理工具箱最新课件Signal Processing Toolbox信号处理工具箱Spline Toolbox样条工具箱Statistics Toolbox统计工具箱Symbolic Math Toolbox符号数学工具箱Simulink Toolbox动态仿真工具箱System Identification Toolbox系统辨识工具箱Wavele Toolbox小波工具箱每个新出的版本都在增加、更新完善 信号处理工具箱最新课件Matlab主工具箱前面课程所介绍的数值计算、符号运算、绘图以及句柄绘图都是matlab主工具箱的内容，是matlab的基本部分。Matlab主工具箱位于： ma

6、tlabtoolboxmatlab matlab主工具箱是任何版本的matlab都不可缺少的。matlab主工具箱中有许多函数库信号处理工具箱最新课件audiovideo 声频、视频支持函数库datafun 数据分析函数库datatypes 数据类型函数库demos matlab演示函数库elfun 初等数学函数库elmat 初等矩阵和时间函数库funfun 函数功能和数学分析函数库general 通用命令函数库graph2d 二维绘图graph3d 三维绘图graphics 句柄绘图函数库信号处理工具箱最新课件iofun 底层输入输出函数库lang 语言结构函数库matfun 矩阵线性代数函

7、数库ops 运算符和逻辑函数库polyfun 多项式函数库sparfun 稀疏矩阵函数库strfun 字符串函数库uitools 图形界面函数库信号处理工具箱最新课件1、各函数库或工具箱中的函数可用help 工具箱（函数库）名查询help optim2、具体函数的内容可以用type 函数名方法查看type laplace3、函数文件定位使用whichwhich laplace信号处理工具箱最新课件工具箱的添加很多时候我们要将一个工具箱添加到到系统中来运行。如果是Matlab安装光盘上的工具箱，重新执行安装程序，选中即可。 如果是单独下载的工具箱，一般情况下仅需要： 把新的工具箱解压到某个目录(

8、例如toolbox目录下) 用addpath命令或者pathtool添加工具箱的路径 用which newtoolbox_command.m来检验是否可以访问 如果能够显示新设置的路径，则表明该工具箱可以使用了。 如果是自己编写的工具箱添加方法同上。信号处理工具箱最新课件二、 信号处理工具箱信号处理工具箱最新课件（一）（一） 波形产生波形产生1.sin 正弦波格式: x= sin(t)例：t=0:0.001:1;y=sin(2*pi*t);plot(t,y)z=sin(2*pi*10*t);plot(t,z)设信号的采样频率为F，信号的自变量通常取为t= 0:1/F:n ，n表示信号的时间长度

9、为n秒。则sin(2*pi*t)即为频率为1的正弦波， sin(2*pi*f*t)即为频率为f的正弦波信号处理工具箱最新课件产生一个幅度为产生一个幅度为2，频率为，频率为4Hz,相位为相位为 的正弦信号的正弦信号 A=2;f=4;phi=pi/6;w0=2*pi*f;t=0:0.01:1;x=A*sin(w0*t+phi);plot(t,x);6/p信号处理工具箱最新课件2. sawtooth周期锯齿波或三角波格式: x= sawtooth(t) sawtooth(t,1) x= sawtooth(2*pi*f*t,width)3. square周期方波格式: x= square(t) squ

10、are(t,50) x= square(2*pi*f*t,duty)4.sinc 产生sinc或sin(pi*t)/(pi*t) 函数格式: x= sinc(t)信号处理工具箱最新课件sinc(x)N=1000;t=-10:20/N:10;x=sinc(t/pi);plot(t,x);grid on信号处理工具箱最新课件5.随机信号：随机信号： rand产生均匀分布的白噪声，产生均匀分布的白噪声，randn产生高斯分布的产生高斯分布的白噪声白噪声t=0:0.01:1;y=randn(1,length(t);plot(t,y);grid on;信号处理工具箱最新课件 6. 单位脉冲序列单位脉冲序

11、列 和单位阶跃序列和单位阶跃序列u(n-3) )3( nn=-2:30;x=zeros(1,5),1,zeros(1,27);y=zeros(1,5),ones(1,28);subplot(2,1,1);stem(n,x,fill);grid on;subplot(2,1,2)stem(n,y,fill);grid on;00001)(nnnnnn00001)(nnnnnnu信号处理工具箱最新课件7. 非周期信号非周期信号产生一个高度为产生一个高度为1，宽度为，宽度为3，延时为，延时为2秒的矩形脉冲秒的矩形脉冲 信号信号t=0:0.02:6;plot(t,rectpuls(t-2-1.5,3)

12、; axis(0,6,0,1.5);注：注：rectpuls(t-延时延时-width/2,width); 信号处理工具箱最新课件（二）信号的基本运算1 1 信号的相加与相乘信号的相加与相乘 y(n)=x1(n)+x2(n) y(n)=x1(n)x2(n) MATLAB实现：实现：y=x1+x2; y=x1.*x22 2 序列移位与周期延拓运算序列移位与周期延拓运算序列移位：y(n)=x(n-m)。MATLAB实现：实现：y=x; ny=nx-m序列周期延拓：y(n)=x(n)M，MATLAB实现：实现：ny=nxs:nxf；y=x(mod(ny,M)+1)3 3 序列翻褶与序列累加运算序列翻

13、褶与序列累加运算序列翻褶：y(n)=x(-n)。MATLAB可实现: y=fliplr(x)序列累加的数学描述为： nnisixny)()(MATLAB实现：实现：y=cumsum(x)信号处理工具箱最新课件尺度变换、翻转、时移、尺度变换、翻转、时移、 相加、相乘相加、相乘 t=-3:0.001:3; ft1=tripuls(2*t,4,0.5); subplot(2,1,1) plot(t,ft1) title(f(2t) ft2=tripuls(2-2*t),4,0.5); subplot(2,1,2) plot(t,ft2) title(f(2-2t) 信号处理工具箱最新课件4 差分与求

14、和、微分与积分差分与求和、微分与积分 例例:已知三角波已知三角波f(t)，画出其微分与积分的波形，画出其微分与积分的波形%differentiation %integration h=0.001;t= -3:h:3; t= -3:0.1:3; y1=diff(f2_2(t)*1/h; for x=1:length(t) plot(t(1:length(t)-1),y1) y2(x)=quad(f2_2, -3,t(x); end plot(t,y2)差分差分 y=diff(f); 求和求和 y=sum(f(k1:k2); 微分微分 y=diff(f)/h; h为数值计算所取时间间隔为数值计算所

15、取时间间隔 定积分定积分 quad(function_name,a,b); %function_name为被积函数名，为被积函数名，a和和b指定积分区间。指定积分区间。 信号处理工具箱最新课件其中函数其中函数 f2_2(t)为：为： function ft2=f2_2(t) ft2=tripuls(t,4,0.5); 信号处理工具箱最新课件5 两序列的卷积运算两序列卷积运算： mmnxmxnxnxny)()()()()(2121MATLAB实现：实现：y=conv(x1,x2)。序列序列x1(n)和和x2(n)必须长度有限。必须长度有限。 6 两序列的相关运算两序列相关运算： nmnxnxmy

16、)()()(21。MATLAB实现：实现：y=xcorr(x1,x2)。信号处理工具箱最新课件NnnxnxE0*1021NnnxNP7.信号能量数学定义：MATLAB实现实现: E=sum(x.*conj(x); 或或: E=sum(abs(x).2);数学定义：8. 8. 信号功率信号功率MATLAB实现: P=sum(x.*conj(x)/N; 或: P=sum(abs(x).2)/N;信号处理工具箱最新课件（三）傅里叶(Fourier)变换1 1 连续时间、连续频率傅里叶变换连续时间、连续频率傅里叶变换2 2 连续时间、离散频率傅里叶级数连续时间、离散频率傅里叶级数dtetxjXtj)(

17、)(pdejXtxtj)(21)(正变换： 逆变换： NitjitetXjX1)()(teeetxtxtxntjtjtjn21,)(,),(),(212/2/00000)(1)(TTtjkdtetxTjkXktjkejkXtx0)()(0正变换： 逆变换： 信号处理工具箱最新课件3 3 时间离散、连续频率序列傅里叶变换时间离散、连续频率序列傅里叶变换4 4 离散时间、离散频率离散傅里叶级数离散时间、离散频率离散傅里叶级数5 离散时间、离散频率离散傅里叶变换(DFT)nnjjenxeX)()(pppdeeXnxnjj)(21)(正变换： 逆变换： 10)()()(NnnkNWnxnxDFSkX正

18、变换： 逆变换： 10)(1)()(NknkNWkXNnXIDFSnx1, 2, 1, 0Nk1, 2, 1, 0Nn10)()()(NnnkNWnxnxDFTkX1, 2, 1, 0Nk正变换： 10)(1)()(NknkNWkXNkXIDFSnx1, 2, 1, 0Nn逆变换： 信号处理工具箱最新课件1一维快速正傅里叶变换函数一维快速正傅里叶变换函数fft格式：格式：X=fft(x, N)功能：采用FFT算法计算序列向量x的N点DFT变换; 当N缺省时，fft函数自动按x的长度计算DFT;当N为2整数次幂时，fft按基-2算法计算,否则用混合算法。2一维快速逆傅里叶变换函数一维快速逆傅里叶

19、变换函数ifft格式：格式：x=ifft(X, N)功能：采用FFT算法计算序列向量X的N点IDFT变换。, 信号处理工具箱最新课件N=512;F=1000;n=1;t=0:1/F:n;x=sin(2*pi*50*t)+ sin(2*pi*120*t);y=x+1.5*randn(1,length(t); Y=fft(y,N);P=Y.*conj(Y)/N; %计算功率谱密度f=F*(0:N/2-1)/N; plot(f,P(1:N/2)05010015020025030035040045050001020304050607080信号处理工具箱最新课件（四）（四） 统计信号处理统计信号处理1c

20、ov 协方差矩阵格式: c=cov(x)当x为矢量时，cov(x)可求出矢量x的方差(标量)当x为矩阵时，cov(x)可得到协方差矩阵而diag(conv(x)则为由每列数据的方差所构成的矢量信号处理工具箱最新课件2xcov 互协方差(自协方差)函数估计格式: v=xcov(x,y) v=xcov(x) 3xcorr 互相关(自相关)函数估计格式: v=xcorr(x,y) v=xcorr(x) 信号处理工具箱最新课件（五）（五） 滤波滤波IIR滤波器结构：信号处理工具箱最新课件M阶IIR滤波器：差分方程表达式：信号处理工具箱最新课件FIR滤波器结构：信号处理工具箱最新课件M阶FIR滤波器：差

21、分方程表达式：信号处理工具箱最新课件1filter 利用递归滤波器(IIR)或非递归滤波器(FIR)对数据进行数字滤波格式: y=filter(b,a,x) b,a为滤波器系数，x为待滤波的数据 2fftfilt 利用基于FFT的重叠相加法对数据进行滤波，只适用于非递归滤波器(FIR)格式: y=fftfilt(b,x) 信号处理工具箱最新课件3freqz 数字滤波器的频率响应格式: h,w=freqz(b,a) freqz(b,a,n)freqz(b,a,w) 4freqs 模拟滤波器的频率响应格式: h,w= freqs(b,a) freqs(b,a,n) freqs(b,a,w) 信号处

22、理工具箱最新课件（六）（六） IIR滤波器设计滤波器设计1besself 贝塞尔模拟滤波器设计 b,a=besself(n,Wn)n阶截止频率Wn的低通模拟滤波器 b,a=besself(n,W1 W2) W1W22n阶带通模拟滤波器 b,a=besself(n,Wn,high) n阶截止频率Wn的高通模拟滤波器 b,a=besself(n,W1 W2,stop) W1W22n阶带阻模拟滤波器信号处理工具箱最新课件2butter 比特沃思模拟和数字滤波器设计 b,a=butter(n,Wn) 0=Wn=1n阶截止频率Wn的低通数字滤波器 b,a=butter(n,W1 W2) W1W22n阶带

23、通数字滤波器 b,a=butter(n,Wn,high) Wn=1相当于fs/2n阶截止频率Wn的高通数字滤波器 b,a=butter(n,W1 W2,stop) W1W22n阶带阻数字滤波器 b,a=butter(n,Wn,high,s) n阶截止频率Wn的高通模拟滤波器信号处理工具箱最新课件例：t=0:0.001:1;x=sin(t*2*pi*200)+sin(t*2*pi);b,a=butter(10,0.2);figure(1);freqz(b,a)y=filter(b,a,x);figure(2)subplot(2,1,1);plot(x)subplot(2,1,2);plot(y)

24、信号处理工具箱最新课件3yulewalk 递归数字滤波器设计格式: b,a=yulewalk(n,f,m) f为频率点，m为相应的响应幅度例：f=0 0.6 0.65 0.7 1; m=1 1 0.5 0 0; b,a= yulewalk(8,f,m) ; h,w= freqz(b,a,128); plot(f,m,w/pi,abs(h),-)信号处理工具箱最新课件（七）（七） FIR滤波器设计滤波器设计1fir1 基于窗函数的FIR滤波器设计标准频率响应b=fir1(n,Wn) 0=Wn=1, Wn=1相当于fs/2 n阶截止频率Wn的加汉明Hamming窗线形相位低通FIR滤波器 b=fi

25、r1(n,W1 W2,stop) W1W22n阶带阻模拟滤波器 b=fir1(n,Wn,Window) Window长度为n+1采用向量Window中指定的窗函数进行设计的滤波器信号处理工具箱最新课件窗函数窗函数1矩形窗 w=boxcar(n)2三角窗 w=triang(n)3巴特利特窗 w=bartlett(n)4汉明窗 w=hamming(n)5汉宁窗 w=hanning(n)6布莱克曼窗 w=blackman(n)7切比雪夫窗 w=chebwin(n,r)8凯泽窗 w=kaiser(n,beta)信号处理工具箱最新课件例：W=chebwin(35,30);b=fir1(34,0.48,h

26、igh,W);freqz(b,1,512);信号处理工具箱最新课件2fir2 基于窗函数的FIR滤波器设计任意标准频率响应 b=fir2(n,f,m)b=fir2(n,f,m,Window)其中：f为频率点，m为相应的响应幅度，Window为指定的窗函数信号处理工具箱最新课件例：f=0 0.6 0.6 1; m=1 1 0 0; b= fir2(30,f,m); h,w= freqz(b,1,128); figure(1);plot(f,m,w/pi,abs(h),-) t=0:0.001:1; x=sin(t*2*pi*350)+sin(t*2*pi*10); y=fftfilt(b,x);

27、figure(2) subplot(2,1,1);plot(x) subplot(2,1,2);plot(y)信号处理工具箱最新课件3firls 最小二乘线形相位FIR滤波器设计b=firls(n,f,m)3remez 最优FIR滤波器设计b= remez(n,f,m)信号处理工具箱最新课件（八）（八） 其他其他1medfilt1 一维中值滤波y= medfilt1(x,n)信号处理工具箱最新课件2interp 提高采样速率（内插）y= interp(x,r)结果：y的长度是x的r倍3decimate 降低采样速率（内插）y= decimate(x,r)结果：y的长度是x的1/r信号处理工具箱最新课件例：t=0:0.00025:1;x=sin(2*pi*30*t)+sin(2*pi*60*t);y= decimate(x,4);subplot(2,1,1);stem(x(1:120);subplot(2,1,2);stem(y (1:30);