title('信源频谱')

% 乘以相干载波后的频谱

mmf=fftshift(abs(fft(nst))); subplot(413) plot(f,mmf)

title('乘以相干载波后的频谱') %经过低通滤波后的频谱

dmf=fftshift(abs(fft(dst))); subplot(414) plot(f,dmf)

title('经过低通滤波后的频谱');

function[out]= tongji(d,M)

%将输入序列扩展成间隔15个0的序列 N = length(d);

out = zeros(1,M*N); % 生产全零矩阵 for i=0:N-1

out(i*M+1)=d(i+1);

%out = reshape(out,1,M*N); % 改变矩阵行列 1行和M*N列 end;

function [f,sf]=FFT_SHIFT(t,st) %dt=t(2)-t(1); T=t(end); df=1/T;

N=length(st);

%f=-N/2*df:df:N/2*df-df; %频率划分间距df f=[-N/2:N/2-1]*df;

sf=fft(st);% 傅里叶变换 sf=T/N*fftshift(sf);