tfGabor.rar_matlab 绘制曲线_图形绘制

%小波Gabor变换的实现程序 clf; clear all ; close all; fs=100; %采样率 Ts=1/fs; t=0:Ts:10; gass=2^(1/4)*exp(-pi*(t).^2).*cos(5*pi*t); % 生成一个高斯函数 subplot(221),plot(t,gass); title('高斯函数'); xlabel('t'); ylabel('幅度'); T=0:Ts:10; ft=cos(T.^2+2*T)+cos(T.^2); % 生成要变换的信号函数 subplot(222),plot(t,ft); title('信号函数'); xlabel('time'); ylabel('幅度'); y=fft(ft); %信号做FFT变换 amp=abs(y); subplot(223);plot(amp); title('信号的FFT变换'); xlabel('F(Hz)'); ylabel('幅度'); subplot(224),plot(t,imag(hilbert(ft))); shl=100; %高斯窗每次平移点数 shn=(length(t)-1)/shl; %求高斯窗平移总次数 y2=zeros(shn,2001); for k=0:shn-1; gassc=2^(1/4)*exp(-pi*(t-k*shl*Ts).^2).*cos(5*pi*t); %平移后的高斯函数 gassc2=gassc/sum(gassc.^2) %归一化 yl=conv(hilbert(ft),gassc2); %短时傅立叶变换，即对信号与Gauss函数做卷积 y2(k+1,:)=yl; end % 最后的图形显示 [F,T]=size(y2); [F,T]=meshgrid(1:T,1:F); figure(2),mesh(F,T,abs(y2)) title('信号 a-F图'); xlabel('F(Hz)'); ylabel('尺度') zlabel('幅度'); figure(3),contour(F,T,abs(y2)) % 等高线图 % Dinga's Blog　http://www.dinga.cn +++++++