matlab fx函数图像,五种常用小波基含MATLAB实现 五种基本函数图像
五种基本函数图像
1.给出五种常用小波基的时域和频域波形图。 与标准的傅里叶变换相比,小波分析中使用到的小波函数具有不唯一性,即小波函数(t) 具有多样性。小波分析在工程应用中,一个十分重要的问题就是最优小波基的选择问题,因为用不同的小波基分析同一个问题会产生不同的结果。目前我们主要是通过用小波分析方法处理信号的结果与理论结果的误差来判定小波基的好坏,由此决定小波基。常用小波基有Haar小波、Daubechies(dbN)小波、Mexican Hat(mexh)小波、Morlet小波、Meyer小波等5种。
(1)Haar小波 Haar函数是小波分析中最早用到的一个具有紧支撑的正交小波函数,也是最简答的一个小波函数,它是支撑域在t[0,1]范围内的单个矩形波。 Haar10t112函数的定义如下:(t)-12t1 0其他 Haar小波在时域上是不连续的,所以作为基本小波性能不是特别好。但它也有自己的优点,如: 计算简单; (t)不但与(2jt)[jz]正交,而且与自己的整数位移正交。
因此,在a族。 (t)的傅里叶变换是: 2j的多分辨率系统中Haar小波构成一组最简单的正交归一的小波()=j4sin2()ej/2 a
Haar小波的时域和频域波形图 haar 时域1.510.50-0.5-1-1.500.511.5t i=20; wav = 'haar'; [phi,g1,xval] = wavefun(wav,i); subplot(1,2,1); plot(xval,g1,'-r','LineWidth',1.5); xlabel('t') title('haar 时域'); g2=fft(g1); g3=abs(g2); subplot(1,2,2);plot(g3); xlabel('f') title('haar 频域') 5haar 频域7x 106543210051015fx 105
(2)Daubechies(dbN)小波 Daubechies小波是世界著名的小波分析学者Inrid Daubechies构造的小波函数,简写为dbN,N是小波的阶数。
小波(t)和尺度函数(t)中的支撑区为除N1外,dbN不具有对称性(即非线性相位)。2N1,(t)的消失矩为N。dbN没有明确的表达式(除N1外),但转换函数h的平方模是明确的。 Daubechies小波系是由法国学者Daubechies提出的一系列二进制小波的总称,在Matlab中记为dbN,N为小波的序号,N值取2,3,…,10。该小波没有明确的解析表达式,小波函数φ与尺度函数Φ的有效支撑长度为2N-1.当N取1时便成为Haar小波。 令p(y)Ck0N1N-1kkyk,其中Ck2N-1k为二项式的系数,则有 m0()(cos2-jkk)p(sin) 2221式中,m()he22N1k00。
Daubechies小波具有以下特点:
(1)在时域是有限支撑的,即(t)长度有限。
(2)在频域()在=0处有N阶零点。
k(t)(t-k)dt。
(4)小波函数(t)可以由所谓“尺度函数”(t)求出来。
尺度函数(t)为
(3)(t)和它的整数位移正交归一,即低通函数,长度有限,支撑域在t=0~(2N-1)范围内。
Daubechies小波的时域和频域波形图 db4 时域1.510.50-0.5-102468t i=10; wname = 'db4'; [phi,g1,xval] = wavefun(wname,i); subplot(1,2,1); plot(xval,g1,'-r','LineWidth',1.5); xlabel('t') title('db4 时域'); g2=fft(g1); g3=abs(g2); subplot(1,2,2);plot(g3,'-r','LineWidth',1.5); xlabel('f') title('db4 频域')db4 频域1000900800700600500400300200100002000400060008000f
注意 Daubechies小波常用来分解和重构信号,作为滤波器使用。波形如下: 分解低通滤波器0.80.60.400.20-0.202468-0.5-110.5分解高通滤波器02468重构低通滤波器10.510.500-0.502468-10重构高通滤波器-0.52468 wname = 'db4'; % 计算该小波的4个滤波器 [Lo_D,Hi_D,Lo_R,Hi_R] = wfilters(wname); subplot(2,2,1); stem(Lo_D); title('分解低通滤波器'); subplot(2,2,2); stem(Hi_D); title('分解高通滤波器'); subplot(2,2,3); stem(Lo_R); title('重构低通滤波器'); subplot(2,2,4); stem(Hi_R); title('重构高通滤波器');
(3)Mexican Hat(mexh)小波 Mexican Hat函数为Gauss函数的二阶导数:(t)(1-t)e2-t22 ()2e2-22 因为它的形状像墨西哥帽的截面,所以也称为墨西哥帽函数。 Mexican Hat(mexh)小波的时域和频域波形图 Mexihat 时域115mexihat 频域0.51005-0.5-10-50t5100050f100 d=-6; h=6; n=100; [g1,x]=mexihat(d,h,n); subplot(2,2,1);plot(x,g1,'-r','LineWidth',1.5); xlabel('t') title('Mexihat 时域'); g2=fft(g1); g3=(abs(g2)); subplot(2,2,2);plot(g3); xlabel('f') title('mexihat 频域'); 特点:墨西哥帽函数在时间域与频率域都有很好的局部化,并且满足(t)dt0。由于它不存在尺度函数,所以小波函数不具有正交性。 R
(4)Morlet小波 它是高斯包络下的单频率副正弦函数: (t)Cecos(5x )2-t2C是重构时的归一化常数。Morlet小波没有尺度函数(t),而且是非正交分解。 Morlet 频域Morlet小波的时域波形图和频域波形图 Morlet 时域10.50-0.5-1-10d=-6; h=6; n=100; [g1,x]=morlet(d,h,n); subplot(2,2,1);plot(x,g1,'-r','LineWidth',1.5); xlabel('t') title('Morlet 时域'); g2=fft(g1); g3=abs(g2); subplot(2,2,2);plot(g3); xlabel('f') title('Morlet 频域') 15105-50t5100050f100
(5)Meyer小波 Meyer小波的小波函数和尺度函数都是在频率域中进行定义的,其定义为: 324(2)esin(2v(2-1))33348()(2)ecos(v(-1))2233280,33-12jw2-12jw2 其中,v(a)为构造Meyer小波的辅助函数,具有 v(a)a(35-84a70a-20a) a0,1 423(2) 324()(2)cos(v(-1))2233043--121223 Meyer小波不是紧支撑的,但它收敛的速度很快: (t)C(1t)n2-n (t)无限可微。 Meyer小波的时域和频域波形图 Meyer 时域210-1-6-4-20tmeyer 频域2461510500204060f80100120140
d = -6; h = 6; n = 128; [psi,x] = meyer(d,h,n,'psi'); subplot(2,1,1), plot(x,psi,'-r','LineWidth',1.5) xlabel('t') title('Meyer 时域'); PSI=fft(psi); PSII=abs(PSI); subplot(2,1,2), plot(PSII); xlabel('f') title('meyer 频域')
2、在信号 x(t)=sin(2π*30t)+cos(2π*50t)加上噪音后分别进行FFT和CWT变换。 解:引入随机噪声randn(1,N) 原信号x(t)波形图250x(t)的fft变换图x(t)0-20050tx(t)加噪声后fft变换图100050fmorlet100500642120406080time (or space) b尺度为2100050f尺度为110010-1scales a05010020-2050100 N=100;fs=1000; n=
0:N-1; t=n/fs; x=sin(60*pi*t)+cos(100*pi*t); %原信号 subplot(3,2,1); plot(x,'-r','LineWidth',1.5); xlabel('t') ylabel('x(t)') title('原信号x(t)波形图') F1=fft(x); m1=abs(F1); subplot(3,2,2);plot(m1);
xlabel('f') title('x(t)的fft变换图') x1=randn(1,N); %加入噪声 x2=x+x1; F2=fft(x2); m2=abs(F2); subplot(3,2,3);plot(m2); xlabel('f') title('x(t)加噪声后fft变换图') scale=[1 2 4 6]; %设置尺度 subplot(3,2,4);x3=cwt(x2,scale,'morl','plot'); title('morlet'); %加噪声后CWT变换结果图 subplot(3,2,5);plot(x3(1,:));title('尺度为1'); subplot(3,2,6);plot(x3(2,:));title('尺度为2');
五种基本函数图像。
Tags:
matlab fx函数图像,五种常用小波基含MATLAB实现 五种基本函数图像相关推荐
- 【RRT三维路径规划】基于matlab RRT算法无人机三维路径规划【含Matlab源码 155期】
一.获取代码方式 获取代码方式1: 通过订阅紫极神光博客付费专栏,凭支付凭证,私信博主,可获得此代码. 获取代码方式2: 完整代码已上传我的资源:[三维路径规划]基于matlab RRT算法无人机三维 ...
- 【RRT三维路径规划】基于matlab RRT算法无人机三维路径规划【含Matlab源码 1363期】
一.获取代码方式 获取代码方式1: 通过订阅紫极神光博客付费专栏,凭支付凭证,私信博主,可获得此代码. 获取代码方式2: 完整代码已上传我的资源:[三维路径规划]基于matlab RRT算法无人机三维 ...
- 【电力预测】基于matlab GUI灰色模型电力负荷预测【含Matlab源码 769期】
一.获取代码方式 获取代码方式1: 完整代码已上传我的资源: [电力负荷预测]基于matlab GUI灰色模型电力负荷预测[含Matlab源码 769期] 获取代码方式2: 通过订阅紫极神光博客付费专 ...
- 【BFS三维路径规划】基于matlab广度优先搜索算法无人机三维路径规划【含Matlab源码 270期】
一.获取代码方式 获取代码方式1: 通过订阅紫极神光博客付费专栏,凭支付凭证,私信博主,可获得此代码. 获取代码方式2: 完整代码已上传我的资源:[三维路径规划]基于matlab广度优先搜索算法无人机 ...
- 【单目标优化求解】基于matlab黑猩猩算法求解单目标问题【含Matlab源码 1413期】
一.获取代码方式 获取代码方式1: 完整代码已上传我的资源:[单目标优化求解]基于matlab黑猩猩算法求解单目标问题[含Matlab源码 1413期] 点击上面蓝色字体,直接付费下载,即可. 获取代 ...
- 【芯片识别】基于matlab GUI形态学PCB板芯片识别【含Matlab源码 1820期】
⛄一.获取代码方式 获取代码方式1: 完整代码已上传我的资源:[芯片识别]基于matlab GUI形态学 PCB板芯片识别[含Matlab源码 1820期] 点击上面蓝色字体,直接付费下载,即可. 获 ...
- 【大米粒计数】基于matlab GUI形态学大米粒颗粒识别【含Matlab源码 915期】
⛄一.获取代码方式 获取代码方式1: 完整代码已上传我的资源:[大米粒计数]基于matlab GUI形态学大米粒颗粒识别[含Matlab源码 915期] (https://download.csdn. ...
- 【雷达干扰】基于matlab速度聚类欺骗式干扰仿真【含Matlab源码 2221期】
⛄一.获取代码方式 获取代码方式1: 完整代码已上传我的资源:[雷达干扰]基于matlab速度聚类欺骗式干扰仿真[含Matlab源码 2221期] 点击上面蓝色字体,直接付费下载,即可. 获取代码方式 ...
- 【TS TSP】基于matlab禁忌搜索算法求解31城市旅行商问题【含Matlab源码 1143期】
一.获取代码方式 获取代码方式1: 完整代码已上传我的资源:[TSP]基于matlab禁忌搜索算法求解31城市旅行商问题[含Matlab源码 1143期] 点击上面蓝色字体,直接付费下载,即可. 获取 ...
最新文章
- 学完 CompletionService,可以做时间管理大师?
- 【python图像处理】几何图形的绘制与文字的绘制(ImageDraw类详解)
- VTK:vtkCubeAxesActor用法实战
- 使用Sublime text 3打造一个小巧但强大的Go语言开发IDE
- python更新到什么版本_Linux更新Python版本及修改python默认版本的方法
- 年末巨制:知识图谱嵌入方法研究总结
- 阿里云云效发布研发协同工具,以新的产研协同工作方式助力实现BizDevOps
- git 删除tag_Git常用命令
- 输出矩阵的左下半三角
- .net中序列化读写xml方法的总结
- 蓝桥杯2017年第八届C/C++省赛C组第二题-兴趣小组
- .java和.class的区别_Java中class与Class有什么区别
- 什么叫误差的反向传播,反向传播误差怎么算的
- redis压测工具redis-benchmark
- urt-8转成GBK 之多种方法
- 2019数据结构考研复习指导习题代码(王道论坛)
- 游戏建模小白必看的游戏人物模型贴图制作方法及制作小技巧
- 详细说明如何在pycharm不联网的情况下,离线安装第三方库及依赖包(如sklearn)
- Redis安装(完整版)
- 华师计算机学院新上任书记,华师大4位校领导履新 任友群任常务副书记