資源簡介
小波(Wavelet)這一術語，顧名思義，“小波”就是小區域、長度有限、均值為0的波形。所謂“小”是指它具有衰減性；而稱之為“波”則是指它的波動性，其振幅正負相間的震蕩形式。與Fourier變換相比，小波變換是時間(空間)頻率的局部化分析，它通過伸縮平移運算對信號(函數)逐步進行多尺度細化，最終達到高頻處時間細分，低頻處頻率細分，能自動適應時頻信號分析的要求，從而可聚焦到信號的任意細節，解決了Fourier變換的困難問題，成為繼Fourier變換以來在科學方法上的重大突破。有人把小波變換稱為“數學顯微鏡”。
代碼片段和文件信息
//?DWTHelp.cpp:?implementation?of?the?CDWTHelp?class.
//
//////////////////////////////////////////////////////////////////////
#include?“stdafx.h“
#include?“DWTHelp.h“
#include?“math.h“
#ifdef?_DEBUG
#undef?THIS_FILE
static?char?THIS_FILE[]=__FILE__;
#define?new?DEBUG_NEW
#endif
//////////////////////////////////////////////////////////////////////
//?Construction/Destruction
//////////////////////////////////////////////////////////////////////
CDWTHelp::CDWTHelp（）
{
m_haar_filter=2;
}
CDWTHelp::~CDWTHelp（）
{
}
//全局變量
LONG?hBen?hEnd?h1Ben?h1End?hLen?h1Len;
double??*h?*g?*h1?*g1;
void CDWTHelp::FilterSet_haar（?）
{
hLen?=?2;
h1Len?=?2;
hBen?=?0;
hEnd?=?0;
h1Ben?=?0;
h1End?=?0;
if（h!=NULL）?delete[]?h;
h?=?new?double[hLen];
if（h1!=NULL）?delete[]?h1;
h1?=?new?double[h1Len];
if（g!=NULL）?delete[]?g;
g?=?new?double[h1Len];
if（g1!=NULL）?delete[]?g1;
g1?=?new?double[hLen];
for（int?k=0;k
h[k]=0.7071;
h1[k]=0.7071;
g[k]=0.7071;
g1[k]=0.7071;
}
}
void?CDWTHelp::FilterSet_97（float?t）
{
//好像長度分別為9/7，然后4/3是中點的位置.
hLen?=?9;
h1Len?=?7;
hBen?=?-4;
hEnd?=?4;
h1Ben?=?-3;
h1End?=?3;
if（h!=NULL）?delete[]?h;
h?=?new?double[hLen];
if（h1!=NULL）?delete[]?h1;
h1?=?new?double[h1Len];
if（g!=NULL）?delete[]?g;
g?=?new?double[h1Len];
if（g1!=NULL）?delete[]?g1;
g1?=?new?double[hLen];
????*h?=?（8*t*t*t-6*t*t+3*t）/（1+2*t）*（1/32.）*sqrt（2）;
*（h+1）?=?（-16*t*t*t+20*t*t-12*t+3）/（1+2*t）*（1/32.）*sqrt（2）;
*（h+2）?=?（2*t-3）/（1+2*t）*（1/8.）*sqrt（2）;
*（h+3）?=?（16*t*t*t-20*t*t+28*t+5）/（1+2*t）*（1/32.）*sqrt（2）;
*（h+4）?=?（-8*t*t*t+6*t*t+5*t+20）/（1+2*t）*（1/16.）*sqrt（2）;
*（h+5）?=?*（h+3）;
*（h+6）?=?*（h+2）;
*（h+7）?=?*（h+1）;
*（h+8）?=?*（h+0）;
*h1?=?-t/16.*sqrt（2）;
*（h1+1）?=?（1-2*t）/16.*sqrt（2）;
*（h1+2）?=?（t+4）/16.*sqrt（2）;
*（h1+3）?=?（3+2*t）/8.*sqrt（2）;
*（h1+4）?=?*（h1+2）;
*（h1+5）?=?*（h1+1）;
*（h1+6）?=?*（h1+0）;
*g1?=?*（h+8）;
*（g1+1）?=?-*（h+7）;
*（g1+2）?=??*（h+6）;
*（g1+3）?=?-*（h+5）;
*（g1+4）?=??*（h+4）;
*（g1+5）?=?-*（h+3）;
*（g1+6）?=??*（h+2）;
*（g1+7）?=?-*（h+1）;
*（g1+8）?=??*h;
*g?=?-*（h1+6）;
*（g+1）?=??*（h1+5）;
*（g+2）?=?-*（h1+4）;
*（g+3）?=??*（h1+3）;
*（g+4）?=?-*（h1+2）;
*（g+5）?=??*（h1+1）;
*（g+6）?=?-*h1;
}
//小波逆變換函數
void?CDWTHelp::ICWT2Dhelp（double?*?buffer?LONG?height?LONG?width?LONG?lv）
{
LONG?ijk;
LONG?nheightnwidth;
for?（k=0;?k
nheight=height>>（lv-k-1）;
nwidth=width>>（lv-k-1）;
double?*p;
p?=?new?double?[nheight];
for（j=0;?j
for?（i=0;?i
ICWT1Dhelp（pnheighthBenhEndh1Benh1Endh1g1）;
for（i=0;?i
}
delete?[]p;
double?*p1;
p1?=?new?double?[nwidth];
for?（i=0;?i
for（j=0;?j
*（p1+j）?=?*（buffer+i*width+j）;//-xia
ICWT1Dhelp（p1nwi?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件??????10689??2008-11-19?21:04??dwt\DWTHelp.cpp
?????文件???????1301??2008-07-11?15:36??dwt\DWTHelp.h
?????目錄??????????0??2010-07-21?11:23??dwt
-----------?---------??----------?-----??----
????????????????11990????????????????????3
- 上一篇：某企業ERP系統需求分析報告
- 下一篇：東南大學通信原理課程試卷+答案+復習提綱
評論
共有 條評論
川公網安備 51152502000135號