- 在线时间
- 57 小时
- 最后登录
- 2018-4-12
- 注册时间
- 2012-2-6
- 听众数
- 4
- 收听数
- 0
- 能力
- 0 分
- 体力
- 4729 点
- 威望
- 0 点
- 阅读权限
- 60
- 积分
- 1673
- 相册
- 0
- 日志
- 0
- 记录
- 1
- 帖子
- 512
- 主题
- 9
- 精华
- 0
- 分享
- 0
- 好友
- 0
升级 67.3% TA的每日心情 | 奋斗 2018-4-11 22:36 |
---|
签到天数: 490 天 [LV.9]以坛为家II
|
2#
发表于 2014-5-28 13:04
|只看该作者
|
|邮箱已经成功绑定
先附上匀加速情况下的代码:# K3 e& v3 s; B% i
clear all; close all;
1 n7 W+ [$ z9 w; H# R8 o6 o%% 参数设置
6 q# w" Q" ` Bdistance = 20.0e3; % 目标初始点到雷达的距离20000米5 F- h- S" @, U' @. K
f0 = 340e9; % 载波频率
+ e9 o; m; l1 g5 H* i+ o5 V3 F%fs = 1000*40;* M4 n# k- |3 u/ Z, E; I2 H4 C, _
%fs = 1000*60; % 采样频率
* b+ l ?2 {" l2 |% ufs = 1000*50;! i! i3 Z1 o& d
% f0 = 10e9; % 载波频率; s9 H9 u9 a0 |
% fs = 1000*16; % 采样频率/ {9 x, a6 I f0 K
t_dur=1;
F/ c6 M* `1 g% B9 M( ]' a%t_dur = 0.2; % 观测时间长度: C* N" h/ Q8 v) _
t_start = 0; % 观测时间起点
5 l+ l$ y5 ^! R6 jM = t_dur*fs; % 时间点数2 @8 V* |8 F' k
t=linspace(t_start+1/fs,t_start+t_dur,M)'; % 各采样时刻
~$ J5 B; J6 N% i/ ~accelerate_1 = 9.6; % 加速度
5 z- A- L. j+ _0 p B% ?0 |2 I) _) ^accelerate_2 = 9.2;: \0 J# K4 `8 B% }" i. G, p
c = 3.0e8; % 电磁波传播速度(m/s), f. j& ~. c6 W7 u* e& u9 c
wave_length = c/f0; % 定义波长4 M5 J$ o( X x' y" t$ Q, m
N=200;
- s1 x+ {# A# o$ E Gv = 0; % 一阶运动速度; T9 N* M- B/ @4 ]* s V
echo_acc = exp(j*4*pi/wave_length*(distance+accelerate_1*t.^2))+exp(j*4*pi/wave_length*(distance+accelerate_2*t.^2));%匀加速目标回波信号
' F! y) l( ?" t% v8 v%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%3 Z' q! C! E' s5 C
%% Gabor算法$ ]8 i9 s6 g |# `+ @6 b
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
' K( T4 ^' J5 b9 |" @n = 4;
3 x5 O- `1 h% @7 c f" QF = N/n; % Gabor变换频率点数为fs/n
/ x) L* W+ e; s! u" h2 O[tfr_acc_gab,t_gab,f] = tfrgabor(echo_acc,N,F);" |/ O, g4 r5 t/ I4 a: W! x9 E2 Y
tfr_acc_gab = fftshift(tfr_acc_gab,1); S& v: K8 {4 V) d% {/ Z- c- s1 w
figure, d. g8 T6 n. }% e; a6 n9 i
colormap(gray(256));
* g( d0 K, w9 B7 S+ W# himagesc(t,linspace(0,fs/2-n,F),255-abs(tfr_acc_gab));
1 Q2 @+ u' ^+ G r/ F% daxis([t_start t_start+t_dur 0 fs/2]);
& X" a D6 z: [* x( t; r" Lset(gca,'YDir','normal'): j4 T, X* F# {9 Y4 W! e" W; o6 C
grid on;
7 D' B/ Q, u& v2 o! exlabel('时间');ylabel('频率'); title('匀加速运动目标微多普勒时频像');" W' g- u# ^3 ~9 E" Y
大家帮我看一下是怎么回事,不报错,有几个warning,说什么不完全匹配,然后即一直卡住,半个小时之后才有结果。 |
|