完善资料让更多小伙伴认识你,还能领取20积分哦, 立即完善>
本帖最后由 bee825 于 2016-6-15 10:52 编辑
这是一个matlab的czt函数实现24GHZ微波雷达频谱局部细化、提高分辨率、测距精度的实例讲解(包含代码) 希望能对大家有所参考帮助,感兴趣的朋友可以加我企鹅号274001657,一起交流更多关于雷达技术的事情。 fs=256;%采样频率 N=512;%采样点数 nfft=512; n=0:1:N-1;%时间序列号 %n/fs:采样频率下对应的时间序列值 n1=fs*(0:nfft/2-1)/nfft;%F F T对应的频率序列 x=3*cos(2*pi*100*n/fs)+3*cos(2*pi*101.45*n/fs)+2*cos(2*pi*102.3*n/fs)+4*cos(2*pi*103.8*n/fs)+5*cos(2*pi*104.5*n/fs); figure; plot(n,x); xlabel('时间t'); ylabel('value'); title('信号的时域波形'); %------- XK=fft(x,nfft);%单边幅值谱 figure; subplot(211);stem(n1,abs(XK(1:(nfft/2))));%用杆状来画FFT的图,也可以 axis([95,110,0,1500]); title('直接利用FFT变换后的频谱'); subplot(212);plot(n1,abs(XK(1:(N/2)))); axis([95,110,0,1500]); title('直接利用FFT变换后的频谱'); %----------- f1=100;%细化频率段起点 f2=110;%细化频率段终点 M=256;%细化频段的频点数,(这里其实就是细化精度) w=exp(-j*2*pi*(f2-f1)/(fs*M));%细化频段的跨度(步长) a=exp(j*2*pi*f1/fs);%细化频段的起始点,这里需要运算一下才能代入czt函数 xk=czt(x,M,w,a); h=0:1:M-1;%细化频点序列 f0=(f2-f1)/M*h+100;%细化的频率值 figure; subplot(211);stem(f0,abs(xk)); xlabel('f'); ylabel('value'); title('利用CZT变换后的细化频谱'); subplot(212);plot(f0,abs(xk)); xlabel('f'); ylabel('value'); title('利用CZT变换后的细化频谱');
|
|
相关推荐
|
|
小黑屋| 手机版| Archiver| 电子发烧友 ( 湘ICP备2023018690号 )
GMT+8, 2024-11-11 19:38 , Processed in 0.656213 second(s), Total 70, Slave 52 queries .
Powered by 电子发烧友网
© 2015 bbs.elecfans.com
关注我们的微信
下载发烧友APP
电子发烧友观察
版权所有 © 湖南华秋数字科技有限公司
电子发烧友 (电路图) 湘公网安备 43011202000918 号 电信与信息服务业务经营许可证:合字B2-20210191 工商网监 湘ICP备2023018690号