❶ 三级提升小波变换怎么用matlab实现

[YC,YS]=wavedec2(Y,2,'db1');
Y为要分解的图像矩阵，2为分解的层数，‘db1'为采用的小波基

返回两个矩阵YC和YS。Yh2=detcoef2('h',YC,YS,2);这是提取出图像2层分解后的水平分量，h改v是垂直分量，h该d是对角分量。细节分量用另外一个方法提取。

❷ 求提取小波系数的matlab程序

% 装载leleccum信号
load leleccum;
s = leleccum(1:3920);
% 用db1小波函数对信号进行三尺度小波分解
[C,L]=wavedec(s,2,'db1');
figure(1);
plot(s);
title('leleccum原始信号');
% 提取尺度1的低频系数
cA1 = appcoef(C,L,'db1',1); %用小波分解框架[C.L]计算1层低频系数的近似值，小波基为db1
% 提取尺度2的低频系数
cA2 = appcoef(C,L,'db1',2);
figure(2);
subplot(2,1,1);
plot(cA1);
title('尺度1的低频系数');
subplot(2,1,2);
plot(cA2);
title('尺度2的低频系数');
% 提取尺度1的高频系数
cD1 = detcoef(C,L,1); %用小波分解框架[C.L]计算1层高频系数的近似值，小波基为db1
% 提取尺度2的高频系数
cD2 = detcoef(C,L,2);
figure(3);
subplot(2,1,1);
plot(cD1);
title('尺度1的高频系数');
subplot(2,1,2);
plot(cD2);
title('尺度2的高频系数');

❸ 用MATLAB 进行整数小波分解的程序，最好有程序详细介绍的，谢谢！

举个例子，希望有所帮助。clc;clearall;closeall;X=rgb2gray(imread('lena.jpg'));x=X;x=double(x);line=250;column=250;J=6;Tu=7;smooth_area_num=20;edge_area_num=40;fori=1:linecoef1(i,:)=cwt(x(i,:),1,'bior4.4');coef2(i,:)=cwt(x(i,:),2,'bior4.4');coef3(i,:)=cwt(x(i,:),J,'bior4.4');endmax=0;max_loc=0;max_num=15;fori=1:linej=1;whilej0&coef3(i,j+1)0)%找到左位置jfork=j+1:column-1ifcoef3(i,k)*coef3(i,j)>=0%找到右位置kform=j:kifabs(coef3(i,m))>maxmax=abs(coef3(i,m));max_loc=m;endendifmax>max_numform(i,max_loc)=1;endj=k;max=0;break;endend%forkifk==column-1j=column-1;endelsej=j+1;end%ifend%whileend%对列作小波变换x=x';fori=1:linecoef1(i,:)=cwt(x(i,:),1,'bior4.4');coef2(i,:)=cwt(x(i,:),2,'bior4.4');coef3(i,:)=cwt(x(i,:),J,'bior4.4');end%正值1，负值0的矩阵temp表示(突变点左位置）max=0;max_loc=0;fori=1:linej=1;whilej0&coef3(i,j+1)0)%找到左位置jfork=j+1:column-1ifcoef3(i,k)*coef3(i,j)>=0%找到右位置kform=j:kifabs(coef3(i,m))>maxmax=abs(coef3(i,m));max_loc=m;endendifmax>max_numform(max_loc,i)=1;endj=k;max=0;break;endend%forkifk==column-1j=column-1;endelsej=j+1;end%ifend%whileendfori=2:line-2forj=2:column-2%链长为1ifform(i,j-1)==0&form(i,j+1)==0&form(i-1,j-1)==0&form(i-1,j)==0&form(i-1,j+1)==0&form(i+1,j-1)==0&form(i+1,j)==0&form(i+1,j+1)==0form(i,j)=0;%endifform(i,j)==1&form(i,j+1)==1&sum(sum(form(i-1:i+1,j-1:j+2)))==2form(i,j)=0;form(i,j+1)=0;endifform(i,j)==1&form(i+1,j)==1&sum(sum(form(i-1:i+2,j-1:j+1)))==2form(i,j)=0;form(i+1,j)=0;endifform(i,j)==1&form(i+1,j+1)==1&sum(sum(form(i:i+1,j:j+1)))==2ifsum(sum(form(i-1:i+2,j-1:j+2)))==3&(form(i-1,j+2)==1|form(i+2,j-1)==1)form(i,j)=0;form(i+1,j+1)=0;endifsum(sum(form(i-1:i+2,j-1:j+2)))==4&form(i-1,j+2)==1&form(i+2,j-1)==1form(i,j)=0;form(i+1,j+1)=0;endendifform(i,j+1)==1&form(i+1,j)==1&sum(sum(form(i:i+1,j:j+1)))==2ifsum(sum(form(i-1:i+2,j-1:j+2)))==3&(form(i-1,j-1)==1|form(i+2,j+2)==1)form(i,j+1)=0;form(i+1,j)=0;endifsum(sum(form(i-1:i+2,j-1:j+2)))==4&form(i-1,j-1)==1&form(i+2,j+2)==1form(i,j+1)=0;form(i+1,j)=0;endendendendmax=0;fori=3:line-3forj=3:column-3s=sum(sum(form(i-2:i+2,j-2:j+2)));form=i-2:i+2forn=j-2:j+2ifs>=7&abs(coef3(m,n))<=s*3form(m,n)=0;endendendendendfigure;imshow(form,[]);

❹ matlab小波包程序

t = wpdec(x,3,'db1','shannon');改为抄T = wpdec(x,3,'db1','shannon');

plot(t)改为plot(T)；

rcfs = wprcoef(t,[2 1]);改为rcfs = wprcoef(T,[2 1]);

变量名除非用袭于递归，不要前后重复使用。后面的t只是部分覆盖了开头的t所以会出现问题。
希望对你能有所帮助。

❺ 求一个可运行的5/3或者9/7小波提升算法的MATLAB源代码

您好，定义二维5/3小波提升方案的分解函数如下：
int multilevel_decomposition_2D( int * image, int width, int height, int * wavelet_data, int
levels );
定义二维5/3小波提升方案的重构函数如下：
int multilevel_reconstruction_2D ( int * image, int width, int height, int * wavelet_data, int levels );
类似地，对于9/7小波，我们用C++语言分别实现了二维的多层9/7小波提升方案的分解和重构函数如下：
void WaveletForward2D(int *data, int height, int width, int level); void WaveletInvert2D(int *data, int height, int width, int level);
值得注意的是，DSP平台上实现的5/3、9/7小波提升算法的输入和输出都是整型，这十分有利于直接对图像进行变换。另外，9/7小波提升算法中将对每次乘积的结果通过floor函数进行取整，从而实现了整型到整型的无失真变换。DSP上只实现了CDF97小波。
在CCS 2 ('C6000)编译环境下编程，程序编译优化后，加载到DSP工作平台上。由于本算法使用的存储空间较多，程序中大量使用了堆空间，默认的存储分配方式无法满足要求，所以重新编写cmd文件分配堆栈空间大小如下：
-stack 0xf00 -heap 0x20000

❻ 怎样在matlab中自己编写一个小波的程序

怎样在matlab中自己编写一个小波的程序
a1.a2,d2,d1是ca1,ca2,cd1,cd2这些小波系数的重构。ca1,ca2,cd1,cd2是小波系数，它们的数据点内数随分解层次的增容大而减少，这就难以与原始信号对比分析，通常会经过重构变为与原始信号个数相同的a1.a2,d2,d1，从物理意义上讲，只有a1.a2,d2,d1才是有实际量纲的信号，ca1,ca2,cd1,cd2是没有量纲和物理意义的。

❼ 关于小波变换的Matlab编程

这是一个带参数的函数 根据下面的参数校验可以看出你应该是直接点运行了吧 所以报错了
你至少应该产生一个输入信号，比如一个正弦信号来做输入进行变换
。。。谁给你的程序啊 这bug也太多了 很多低级bug 我只是改的能运行了
function [WT,FreqBins,Scales]=CWT_Morlet(Sig,WinLen,nLevel)
%continuous Wavelet Transform using Morlet function
%Sig:信号
%WinLen:小波函数在尺度参数a=1时的长度（默认为10）
%nLevel:频率轴划分区间段（默认为1024）
%
%WT:返回的小波变换计算结果
%FreqBins:返回频率轴划分结果（归一化频率，最高频率为0.5）
%Scales:返回与频率轴划分值相对应的尺度划分（频率0.5对应的尺度为1）
if (nargin==0),
error('At least 1 parameter required!');
end;

if (nargin<3),
nLevel=1024;
end
if (nargin<2),
WinLen=10;
end
Sig=hilbert(real(Sig));%计算信号的解析信号
SigLen=length(Sig);%获取信号的长度
fmax=0.5;%设置最高分析频率
fmin=0.005;%设置最低分析频率
FreqBins=logspace(log10(fmin),log10(0.5),nLevel);%将频率轴在分析范围内等
%对数坐标划分
Scales=fmax*ones(size(FreqBins))./FreqBins;%计算响应的尺度参数
omg0=WinLen/6;
WT=zeros(nLevel,SigLen);%分配计算结果的存储单元
wait=waitbar(0,'Under calculation,please wait');
for m=1:nLevel,
waitbar(m/nLevel,wait);
a=Scales(m);%提取尺度参数
t=-round(a*WinLen):1:round(a*WinLen);
Mor1=pi^(-1/4)*exp(1i*2*pi*0.5*t/a).*exp(-t.^2/2/(2*omg0*a)^2);
%计算当前尺度下的小波函数
temp=conv(Sig,Mor1)/sqrt(a); %计算信号与小波函数的卷积
WT(m,:)=temp(round(a*WinLen)+1:length(temp)-round(a*WinLen));
end;
close(wait);
WT=WT/WinLen;

%%%%%%运行示例%%%%%%%%%
Fs=100;
t=0:1/Fs:2*pi;
sig=sin(t);
[WT,FreqBins,Scales]=CWT_Morlet(sig);