新手求助,用matlab来改进hilbert算法,下面是源程序:
function testhilbert()%t = linspace(0,10*pi,1024);
data = sin(t)+cos(5*t)+sin(t/2.0);figure(1);plot(t,data);hdata = hilbert(data);
a = (hdata.*conj(hdata)).^(0.5);p = angle(hdata);f = diff(p);
%figure(1);plot(t,real(a));figure(2);plot(t,p);figure(3);plot(t(1:length(f)),f);
fid = fopen('c:\\data\\fivtanks.dzt','rb')datagpr = fread(fid,[512,400],'uint8');
fclose(fid);
figure(10);imagesc(datagpr);
hdata = hilbert(datagpr);
a = (hdata.*conj(hdata)).^(0.5);p = angle(hdata);f = diff(p);figure(11);imagesc(real(a));figure(12);imagesc(p);figure(13);imagesc(f);
function testhilbert()%t = linspace(0,10*pi,1024);
data = sin(t)+cos(5*t)+sin(t/2.0);figure(1);plot(t,data);hdata = hilbert(data);
a = (hdata.*conj(hdata)).^(0.5);p = angle(hdata);f = diff(p);
%figure(1);plot(t,real(a));figure(2);plot(t,p);figure(3);plot(t(1:length(f)),f);
fid = fopen('c:\\data\\fivtanks.dzt','rb')datagpr = fread(fid,[512,400],'uint8');
fclose(fid);
figure(10);imagesc(datagpr);
hdata = hilbert(datagpr);
a = (hdata.*conj(hdata)).^(0.5);p = angle(hdata);f = diff(p);figure(11);imagesc(real(a));figure(12);imagesc(p);figure(13);imagesc(f);




