用matlab时频工具箱进行不同频段下的微多普勒时频分析，作图并比较，为什么这么卡？

tiansifei1

这两天赶毕业设计，要用matlab时频工具箱进行不同频段下的微多普勒时频分析，针对匀加速、振动、转动等情况作出对应目标的微多普勒时频像。弄了几个才三四十行的小程序，结果matlab一运行就卡住了，至少半个小时才能跑出结果的图像。这是怎么回事啊？是我的程序有问题？还是操作系统、matlab软件有问题？我matlab用的是7.0（以前想装成2010版的，结果安装过程中要报错，又换回7.0的了）。操作系统用的是XP，内存2G.
  n: t- L/ J, A3 d: z$ l/ k9 z

tiansifei1

先附上匀加速情况下的代码：
  ^( @+ K, ^% C: d5 Dclear all; close all;
%% 参数设置
% Z& t) `2 k8 _. E- V$ f# Y: ~0 adistance = 20.0e3;                                                          % 目标初始点到雷达的距离20000米
f0 = 340e9;                                                                % 载波频率
%fs = 1000*40;
, ^* a% c, V' v% w0 u) ~%fs = 1000*60;                                                               % 采样频率
0 O# ]( D; X4 H& \4 ~' w2 K" Q& ufs = 1000*50;
% f0 = 10e9;                                                                % 载波频率
* |  L9 q( H( V" b' t# Z% fs = 1000*16;                                                               % 采样频率
t_dur=1;
0 u2 \$ M. y' [%t_dur = 0.2;                                                                 % 观测时间长度
0 Z9 z' q7 _9 i3 \- m( E' W5 \9 i: j8 Ft_start = 0;                                                               % 观测时间起点
* V) _3 ^7 u# X; V1 KM = t_dur*fs;                                                              % 时间点数
t=linspace(t_start+1/fs,t_start+t_dur,M)';                                  % 各采样时刻
accelerate_1 = 9.6;                                                       % 加速度
accelerate_2 = 9.2;
6 ]: Y7 o3 _1 J: Gc = 3.0e8;                                                                 % 电磁波传播速度(m/s)
* s+ r# J4 S2 R" hwave_length = c/f0;                                                        % 定义波长
1 _+ H# D$ b, F7 C/ MN=200;
+ }& `# K7 y1 `/ {v = 0;                                                                   % 一阶运动速度
! f8 i! {6 w) c  N1 \echo_acc = exp(j*4*pi/wave_length*(distance+accelerate_1*t.^2))+exp(j*4*pi/wave_length*(distance+accelerate_2*t.^2));%匀加速目标回波信号
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
; E$ [. L% y6 n! b+ v2 Q# g8 U8 M%% Gabor算法
) @% S, }4 ^  F* \! }) u$ r( y$ j8 H%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
9 W* h. G  a, I' q$ E' jn = 4;
F = N/n;                                                                   % Gabor变换频率点数为fs/n
/ _, C0 g5 I" C6 ?[tfr_acc_gab,t_gab,f] = tfrgabor(echo_acc,N,F);
tfr_acc_gab = fftshift(tfr_acc_gab,1);
) D7 H, K: ^8 Q! Ofigure
colormap(gray(256));
imagesc(t,linspace(0,fs/2-n,F),255-abs(tfr_acc_gab));
: y( x- \$ G2 taxis([t_start t_start+t_dur 0 fs/2]);
' m" W- }5 M; o. i7 Yset(gca,'YDir','normal')
grid on;
xlabel('时间');ylabel('频率'); title('匀加速运动目标微多普勒时频像');
0 v# f* E  ^& g5 B2 }1 `& g* M# F大家帮我看一下是怎么回事，不报错，有几个warning，说什么不完全匹配，然后即一直卡住，半个小时之后才有结果。

tiansifei1

帮我跑一下这个程序也行，只希望快点出图，不要让我再等半个小时了！！