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

tiansifei1

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

tiansifei1

先附上匀加速情况下的代码：
6 h# |4 G# c, S- a- u; J  kclear all; close all;
%% 参数设置
8 f* c9 U3 u! {; C8 ~; c3 O# U4 mdistance = 20.0e3;                                                          % 目标初始点到雷达的距离20000米
f0 = 340e9;                                                                % 载波频率
$ z* ~( u7 f; W& E) c  s%fs = 1000*40;
%fs = 1000*60;                                                               % 采样频率
fs = 1000*50;
% f0 = 10e9;                                                                % 载波频率
( J4 t$ R% t# ~* z* F& S# W$ X% fs = 1000*16;                                                               % 采样频率
9 L4 F3 l; C& h/ i8 d, @& tt_dur=1;
6 n; b: M- [4 F0 V3 R2 K%t_dur = 0.2;                                                                 % 观测时间长度
, _0 s% l8 ?2 I- J0 _# Mt_start = 0;                                                               % 观测时间起点
M = t_dur*fs;                                                              % 时间点数
+ ^- p) O8 j8 I. f3 d7 xt=linspace(t_start+1/fs,t_start+t_dur,M)';                                  % 各采样时刻
7 X9 o* j! y. e8 V( ]accelerate_1 = 9.6;                                                       % 加速度
' t/ l1 c5 y2 R. ?( f! p% uaccelerate_2 = 9.2;
" l. s7 {& z1 O6 x- L6 f) Nc = 3.0e8;                                                                 % 电磁波传播速度(m/s)
wave_length = c/f0;                                                        % 定义波长
3 V9 i% ~4 z. [& d0 v* n# ~N=200;
6 {2 O; d  r5 q9 J( A1 h, K* Fv = 0;                                                                   % 一阶运动速度
- Q, U  B; T3 N5 [+ }: B4 @echo_acc = exp(j*4*pi/wave_length*(distance+accelerate_1*t.^2))+exp(j*4*pi/wave_length*(distance+accelerate_2*t.^2));%匀加速目标回波信号
! x/ o; q/ O0 R) X%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% Gabor算法
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
n = 4;
F = N/n;                                                                   % Gabor变换频率点数为fs/n
[tfr_acc_gab,t_gab,f] = tfrgabor(echo_acc,N,F);
" y1 }% N5 a5 z! @: Y( E7 f+ f0 Atfr_acc_gab = fftshift(tfr_acc_gab,1);
figure
6 Y# k5 a+ F2 f6 ^% a" Z; S% scolormap(gray(256));
; d, e3 B+ E9 f3 B  _! D8 Y& simagesc(t,linspace(0,fs/2-n,F),255-abs(tfr_acc_gab));
+ k4 L* U0 X+ v/ E7 Jaxis([t_start t_start+t_dur 0 fs/2]);
set(gca,'YDir','normal')
grid on;
2 f5 R# z$ q7 I  k+ Fxlabel('时间');ylabel('频率'); title('匀加速运动目标微多普勒时频像');
, |8 z" {7 a( }# x6 Y; S2 w+ `大家帮我看一下是怎么回事，不报错，有几个warning，说什么不完全匹配，然后即一直卡住，半个小时之后才有结果。

tiansifei1

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