作业帮 > 综合 > 作业

谐波分析中的小波变换求低频系数、高频系

来源:学生作业帮 编辑:神马作文网作业帮 分类:综合作业 时间:2024/11/11 20:04:16
谐波分析中的小波变换求低频系数、高频系
小波变换求的低频系数、高频系数、重构的高频和低频信号和别人的不同,别人用的是二进小波变换或者离散小波变换.我的程序如下,
f0=50;
fs=5000;
N=500;%2048/256=8
rad=180/pi;
n=[0:1:N-1]; ;t=n/fs;
s=10*sin(2*pi*f0*t)+2*sin(6*pi*f0*t+pi/2)+1.2*sin(10*pi*f0*t+pi/3)+sin(14*pi*f0*t+pi/4)+0.7*sin(22*pi*f0*t+pi/6)+2*randn(1,length(t))+2*sin(26*pi*f0*(t-0.02)).*exp(-20*t);
[C,L] = wavedec(s,5,'db8');
cA1 = appcoef(C,L,'db8',1);
cA2 = appcoef(C,L,'db8',2);
cA3 = appcoef(C,L,'db8',3);
cA4 = appcoef(C,L,'db8',4);
cA5 = appcoef(C,L,'db8',5);
cD1 = detcoef(C,L,1);
cD2 = detcoef(C,L,2);
cD3 = detcoef(C,L,3);
cD4 = detcoef(C,L,4);
cD5 = detcoef(C,L,5);
t1=linspace(0,0.1,length(cA1));
t2=linspace(0,0.1,length(cA2));
t3=linspace(0,0.1,length(cA3));
t4=linspace(0,0.1,length(cA4));
t5=linspace(0,0.1,length(cA5));
A5 = wrcoef('a',C,L,'db8',5);
D5= wrcoef('d',C,L,'db8',5);
谐波分析中的小波变换求低频系数、高频系
有啥问题?
再问: 波形和别人的不同,同样的信号用同样的小波,分解层数也一样,但是得到的波形不同。
再答: 你的程序中没有绘制波形显示的语句,你的程序同样使用matlab提供的DWT标准函数,也同样属于使用mallat算法的二进小波变换。波形不同只看你显示语句的问题;或是看看你显示的波形是小波系数还是经由小波系数重构的各阶信号,这两者波形当然不同;再或者你的原始信号的点数是否和别人的相同,如果点数不同,同阶小波结果的波形也可能和别人做的不同。