利用希尔伯特变换进行的单边带幅度调制matlab程序

function [M,m,df]=fftseq(m,ts,df)                         

%   [M,m,df]=fftseq(m,ts,df)

%   [M,m,df] = fftseq(m,ts)

fs = 1/ts;

if nargin ==2

    n1 =0;

else

    n1 = fs/df;

end

n2 = length(m);

n = 2^(max(nextpow2(n1),nextpow2(n2)));

M = fft(m,n);

m = [m,zeros(1,n-n2)];

df = fs/n;

function [M,m,df1,f]=T2F(m,ts,df,fs)

[M,m,df1]=fftseq(m,ts,df);

f=[0:df1:df1*(length(m)-1)]-fs/2;

M=M/fs;