郑州航院

end

diff_nsfmn = abs(hilbert(diff_nsfm)); %hilbert变换，求绝对值得到瞬时幅度（包络检波）

zero=(max(diff_nsfmn)-min(diff_nsfmn))/2; diff_nsfmn1=diff_nsfmn-zero;

%*****************************************

%·*·*·*·*·*·*·*·*·*·*·*·*·*·*·*· %**************时域到频域转换**************

ts=0.001; %抽样间隔 fs=1/ts; %抽样频率

df=0.25; %所需的频率分辨率，用在求傅里叶变换

%时，它表示FFT的最小频率间隔

%*****对调制信号m(t)求傅里叶变换*****

m=ac*cos(2*pi*fm1*t); %原调信号 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)];

df1=fs/n; %以上程序是对调制后的信号u求傅里变换

M=M/fs; %缩放，便于在频谱图上整体观察

f=[0:df1:df1*(length(m)-1)]-fs/2; %时间向量对应的频率向量

%************对已调信号u求傅里变换********** fs=1/ts;

if nargin==2 n1=0; else

n1=fs/df; end

n2=length(sfm);

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

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

df1=fs/n; %以上是对已调信号u求傅里变换