標(biāo)題: 線性調(diào)頻LFM信號脈沖壓縮matlab代碼 [打印本頁]
作者: exp(jpai)+1=0    時間: 2022-12-1 15:52
標(biāo)題: 線性調(diào)頻LFM信號脈沖壓縮matlab代碼
親測可用
  1. %%線性調(diào)頻信號經(jīng)過匹配濾波
  2. T=10e-6;                                   %脈沖持續(xù)時間10us
  3. B=30e6;                                    %線性調(diào)頻信號的帶寬30MHz
  4. K=B/T;                                     %調(diào)頻斜率
  5. Fs=10*B;Ts=1/Fs;                     %采樣頻率和采樣間隔
  6. N=T/Ts;
  7. t=linspace(-T/2,T/2,N);
  8. St=exp(j*pi*K*t.^2);                     %線性調(diào)頻信號
  9. Ht=exp(-j*pi*K*t.^2);                    %匹配濾波
  10. Sot=conv(St,Ht);                         %線性調(diào)頻信號經(jīng)過匹配濾波

  12. subplot(411)
  13. L=2*N-1;
  14. t1=linspace(-T,T,L);
  15. Z=abs(Sot);
  16. Z=Z/max(Z);             %歸一化
  17. Z=20*log10(Z+1e-6);
  18. Z1=abs(sinc(B.*t1));                   %辛格函數(shù)
  19. Z1=20*log10(Z1+1e-6);
  20. t1=t1*B;                                         
  21. plot(t1,Z,t1,Z1,'r.');
  22. axis([-15,15,-50,inf]);grid on;
  23. legend('匹配濾波','辛格');
  24. xlabel('時間/s ');
  25. ylabel('幅值/dB');
  26. title('線性調(diào)頻信號經(jīng)過匹配濾波');

  28. subplot(412)                              %展開
  29. N0=3*Fs/B;
  30. t2=-N0*Ts:Ts:N0*Ts;
  31. t2=B*t2;
  32. plot(t2,Z(N-N0:N+N0),t2,Z1(N-N0:N+N0),'r.');
  33. axis([-inf,inf,-50,inf]);grid on;
  34. set(gca,'Ytick',[-13.4,-4,0],'Xtick',[-3,-2,-1,-0.5,0,0.5,1,2,3]);
  35. xlabel('時間/s');
  36. ylabel('幅度/dB');
  37. title('線性調(diào)頻信號經(jīng)過匹配濾波（補(bǔ)零展開之后）');


  40. Z_value = Z(2950:3050)
  41. subplot(413)
  42. plot(Z_value);
  43. grid on;
  44. [value, locmax ] = max( Z_value(1:40));
  45. locmax
  46. value

  48. subplot(414)
  49. plot(Z);
  50. sum = 0;
  51. for mm = 1:length(Z)
  52.     sum = sum +Z(mm)
  53.    
  54. end

  56. sum = -sum;
  57. sum
復(fù)制代碼






歡迎光臨 (http://www.torrancerestoration.com/bbs/) Powered by Discuz! X3.1