大学毕业设计一席谈之三十二 GPS卫星采样信号的捕获(9)补零和内插总结及工程化

继续讲解！经过补零和内插两种方法来测试FFT在相关运算下的性能，我们明显感觉到任意点数的FFT的计算是最好的。其他的方式都会对相关峰造成损伤！原理如下图所示。

纯手工推导

补零过程涉及的原理和其类似，大家可以自行推导。本系列文章的工程实用性非常强，只有这么循序渐进的学，研发出来的产品才具备最优性价比。学的知识能够真正用起来，这才是工程师的价值体现。整天靠编内容发文章的人应该成不了工程师，那这么做是为了什么？也许就是为了那么点利益！弄虚作假出不了真正的产品，只会害人害己。我拒绝和这样的人交流，以前身边有这么一两位这样的同事，现在已经没有这样的人，专心写作直到退休！会不会坚持十年以后，公众 号会扬名万里？只要我坚持传播真正有用的知识，也许会的！加油，为了梦想还是为了出名？我又开始苦恼起来。抑郁就是这么一天一天的加重了。还是安心写作，其他的事情不要去想了！我在企业工作多年，深知企业不需要弄虚作假的人！但这样的人在高校里面也许会活的很滋润，这就是当下的乱像之一。唉，我只能保证我的一亩三分地里面出的都是真货！世态炎凉，但更需要真正的人才！不为名利的人才！2025年3月，本文增加了产品中涉及此方面的代码示例，内容超过七千字，极具借鉴作用，能够极大的增长工程化经验！

知识的提高过程是漫长而艰辛的，这个系列堪称经典（自我评价），需要大家认真学习才能真正掌握！如果你想到企业上班后能快速上手研发产品，这个系列必须看！一篇都不能拉下！！！文章涉及了大量产品化研制的思路，值得借鉴。先看看仿真结果！看看补零或者内插后相关峰的增益变化情况？

% sig_inter = signal_Interpolate1023to1024(signal);

%%% 正弦波内插后测试FFT效果

sig_inter = func_Interp1023to1024(signal);

% figure(1)

% plot(sig_inter,'->');

% hold on;

% plot(round(signal*1024),'-r*');

% 

% figure(2)

% plot(abs(fft(sig_inter)),'-');

% hold on;

% plot(abs(fft(round(signal*1024))),'-r');

%%% 测试FFT效果

% 补零效果比对

signal_fft = fft(signal);

signal_fft_1024 = fft([signal 0]);

aaa = 1;

%%% 对CA码测试FFT效果

for PRN = 1:32

    %- Generate CA code for given PRN --

    caCodesTable(PRN, :) = generateCAcode(PRN); 

end % for PRN = 1:32

satID = 1; %:32           

caCode = caCodesTable(satID, :); % 转换为存储表

%% 补零