zhuotuhanshu.m

来自「解调时考查频率偏移量的一种鉴频算法 采用了双滤波器结构」· M 代码 · 共 14 行

clear
figure(1)
fs=10000;Ts=1/fs;              %设置抽样频率，因为在fft中运算取2的阶乘比较好，但这里为了与实际频率联系起来就取10000Hz
t=-10:Ts:10;                   %设置抽样的时间范围在－10秒和10秒之间
k1=(t+10)/10;
Nt=randn(1,length(t));         %产生（2fs×10＋1）长的AWGN序列
%Nt=impseq(0,-10*fs,10*fs);     产生冲击序列以备检测
Nk=fft(Nt);                    %进行FFT运算，观察噪声频谱
subplot(3,3,1);plot(k1,abs(Nk));title('Nk');xlabel('k');%频谱为频率归一化的频谱，2对应抽样频率fs
for of=-500:100:500             %对于同一信号，在不同的频偏条件下求出对应的能量差（以100Hz为步长）
    s((of+600)/100)=qiugonglvcha(of,Nt);
end
of=-500:100:500;
figure(1);subplot(3,3,9);plot(of,s);xlabel('f');ylabel('deltP');      %根据上述计算结果画出关系曲线