fig9_21.m
上传用户:szahd2008
上传日期:2020-09-25
资源大小:1275k
文件大小:2k
源码类别:

传真(Fax)编程

开发平台:

Matlab

  1. %Reproduce Fig.21 of text
  2. k = .00035/25;
  3. a = 25*k;
  4. b = 30*k;
  5. c = 27*k;
  6. d = 31*k;
  7. v2 = linspace(0,1345,10000);
  8. f2 = (2.*v2)/.0375;
  9. % H1(f)
  10. T1 = exp(-j*2*pi.*f2*a);
  11. X1 = 1/2.*(1 - T1).*conj(1 - T1);
  12. H1 = 10*log10(abs(X1));
  13. % H2(f)
  14. T2 = exp(-j*2*pi.*f2*b);
  15. X2 = 1/2.*(1 - T2).*conj(1 - T2);
  16. H2 = 10*log10(abs(X2));
  17. % H3(f)
  18. T3 = exp(-j*2*pi.*f2*c);
  19. X3 = 1/2.*(1 - T3).*conj(1 - T3);
  20. H3 = 10*log10(abs(X3));
  21. % H4(f)
  22. T4 = exp(-j*2*pi.*f2*d);
  23. X4 = 1/2.*(1 - T4).*conj(1 - T4);
  24. H4 = 10*log10(abs(X4));
  25. % Plot of the four components of H(f)
  26. figure(1)
  27. subplot(2,1,1)
  28. % H(f) Average
  29. ave2 = abs((X1 + X2 + X3 + X4)./4);
  30. Have2 = 10*log10(abs((X1 + X2 + X3 + X4)./4));
  31. plot(v2,Have2);
  32. axis([0 1345 -25 5]);
  33. title('Two pulse MTI stagger ratio 25:30:27:31');
  34. xlabel('Radial Velocity (m/s)');
  35. ylabel('MTI Gain (dB)');
  36. grid on
  37. % %Mean value of H(f)
  38. v4 = v2; 
  39. f4 = (2.*v4)/.0375;
  40. % H1(f)
  41. T1 = exp(-j*2*pi.*f4*a);
  42. T2 = exp(-j*2*pi.*f4*(a + b));
  43. T3 = exp(-j*2*pi.*f4*(a + b + c));
  44. X1 = 1/20.*(1 - 3.*T1 + 3.*T2 - T3).*conj(1 - 3.*T1 + 3.*T2 - T3);
  45. H1 = 10*log10(abs(X1));
  46. % H2(f)
  47. T3 = exp(-j*2*pi.*f4*b);
  48. T4 = exp(-j*2*pi.*f4*(b + c));
  49. T5 = exp(-j*2*pi.*f4*(b + c + d));
  50. X2 = 1/20.*(1 - 3.*T3 + 3.*T4 - T5).*conj(1 - 3.*T3 + 3.*T4 - T5);
  51. H2 = 10*log10(abs(X2));
  52. % H3(f)
  53. T6 = exp(-j*2*pi.*f4*c);
  54. T7 = exp(-j*2*pi.*f4*(c + d));
  55. T8 = exp(-j*2*pi.*f4*(c + d + a));
  56. X3 = 1/20.*(1 - 3.*T6 + 3.*T7 - T8).*conj(1 - 3.*T6 + 3.*T7 - T8);
  57. H3 = 10*log10(abs(X3));
  58. % H4(f)
  59. T9 = exp(-j*2*pi.*f4*d);
  60. T10 = exp(-j*2*pi.*f4*(d + a));
  61. T11 = exp(-j*2*pi.*f4*(d + a + b));
  62. X4 = 1/20.*(1 - 3.*T9 + 3.*T10 - T11).*conj(1 - 3.*T9 + 3.*T10 - T11);
  63. H4 = 10*log10(abs(X4));
  64. % H(f) Average
  65. ave4 = abs((X1 + X2 + X3 + X4)./4);
  66. Have4 = 10*log10(abs((X1 + X2 + X3 + X4)./4));
  67. % Plot of H(f) Average
  68. subplot(2,1,2)
  69. plot(v4,Have4);
  70. axis([0 1345 -25 5]);
  71. title('Four pulse MTI stagger ratio 25:30:27:31');
  72. xlabel('Radial Velocity (m/s)');
  73. ylabel('MTI Gain (dB)');
  74. grid on