Untitled2.zip_HAAR小波基_haar变换_三维小波_三维小波变换_小波连续变换

%Fourier分析 clear fs=100;N=1024;% 采样频率和数据点数 n=0:N-1;t=n/fs;%时间序列 x=sin(2*pi*30*t)+cos(2*pi*50*t)+rand(size(t));%信号源加随即噪声 y=fft(x,N);%对信号进行快速Fourier变换 mag=abs(y);%求Fourier变换后的振幅 f=n*fs/N; %频域序列 subplot(1,2,1),plot(f,mag); xlabel('频率/Hz'); ylabel('振幅'); subplot(1,2,2);plot(f(1:N/2),mag(1:N/2));%绘出Nyquist频率之前随频率变化的振幅 xlabel('频率/Hz'); ylabel('振幅'); %haar clear clc t=1:100; x=sin(2*pi*30*t)+cos(2*pi*50*t)+rand(size(t));%信号源加随即噪声 figure coef1=cwt(x,[1:0.01:10],'haar','plot');%对信号进行快速连续小波变换 title('对不同尺度的haar小波连续变换的系数值'); xlabel('时间'); ylabel('尺度'); figure coef2=cwt(x,[1:0.01:10],'haar','3Dplot'); title('对不同尺度的haar小波连续变换的系数值的三维图'); xlabel('时间'); ylabel('尺度'); %db % 连续小波变换 clear clc t=1:100; x=sin(2*pi*30*t)+cos(2*pi*50*t)+rand(size(t));%信号源加随即噪声 figure % 用db3小波作母小波函数（如下图形），尺度a分别为1, 1.2, 1.4, 1.6, …, 3. coefs=cwt(x,[1:0.01:10],'db3','plot'); title('对不同的尺度的db3小波连续变换的系数值'); Ylabel('尺度'); Xlabel('时间'); figure % 连续小波变换的三维图形 coefs=cwt(x,[1:0.01:10],'db3','3Dplot'); title('对不同的尺度的db3小波连续变换的系数值'); Ylabel('尺度'); Xlabel('时间');