《尺度和频率的对应可以用命令scal2frq来查》由会员分享,可在线阅读,更多相关《尺度和频率的对应可以用命令scal2frq来查(3页珍藏版)》请在金锄头文库上搜索。
1、尺度和频率的对应可以用命令 scal2frq 来查 eg: load data.txt wname = morl; x=data(1:1000); A = 0; B = 12.5; P =1001; %时间长度 采样点数 t = linspace(A,B,P); delta = (B-A)/(P-1); %计算频率对应的尺度 tab_FREQ = 2.83 10.3; %tab_FREQ 为你要查的频率 scales = 1:0.01:200; tab_PF = scal2frq(scales,wname,delta); for k=1:2dummy,ind = min(abs(tab_PF-
2、tab_FREQ(k);PF_app(k) = tab_PF(ind);SC_app(k) = scales(ind); end %进行连续分解并画图 str1 = 1001 samples of x on 1,200,实际频率(in Hz):2.83 10.3; str2 = pseudo-frequency 矩阵和尺度:; str3 = num2str(tab_PF,scales,2); str4 = Pseudo-frequencies = num2str(PF_app,2); str5 = 相应尺度 = num2str(SC_app,2); figure; w = cwt(x,scal
3、es,wname,plot);ax = gca; colorbar axTITL = get(ax,title); axXLAB = get(ax,xlabel); set(axTITL,String,str1) set(axXLAB,String,strvcat(str4,str5) clc;disp(strvcat(,str1,str2,str3,str4,str5)% 下载信号 load vonkoch vonkoch=vonkoch(1:510);% 尺度 1-32 的连续小波变换 S_Min=1;S_Max=32;index=0; for scale=S_Max:-1:S_Min;i
4、ndex=index+1; cwt_coef(index,:)=CWT_FUNCTION(scale,32*(scale),vonkoch); end% 小波系数取模 cwtcoef_abs=abs(cwt_coef);% 显示 for index=S_Min:S_Max max_coef=max(cwtcoef_abs(index,:); % 系数模最大 min_coef=min(cwtcoef_abs(index,:); % 系数模最小 ext=max_coef-min_coef; % 系数模跨度 cwtcoef_abs(index,:)=64*(cwtcoef_abs(index,:)-
5、min_coef)/ext; % 系数大小变换 endfigure(1)subplot(2,1,1); plot(vonkoch); xlabel(时间) ylabel(幅度) title(分形信号) axis(1 510 0 0.02)subplot(2,1,2) colormap(pink(64); image(cwtcoef_abs) set(gca,YTick,2:3:32) set(gca,YTickLabel,32:-3:2) title(连续小波变换时间尺度图) xlabel(时间) ylabel(尺度)% 某个尺度的连续小波变换的 M 函数% delta 小波变换的尺度 % N
6、 小波函数的长度 % s 原始信号 % g 原始信号某个尺度下的小波变换系数function g=CWT_FUNCTION(delta,N,s);% 原始信号长度 n=length(s);% 构造墨西哥帽子小波函数 for index_x=1:N; x=index_x-(N+1)/2; phi_x(index_x)=(pi(-1/4)*(2/sqrt(3)*(1-x.*x/(delta2)*exp(-(x.*x)/(2*delta2); end; phi_x=phi_x/norm(phi_x); % 能量归一化% 对信号做卷积 g=conv(s,phi_x); % 卷积 g=wkeep(g,n); % 保持信号长度