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

tiansifei1

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

tiansifei1

先附上匀加速情况下的代码：
6 }/ _1 U6 ]- O9 uclear all; close all;
5 O$ D1 V" ]# o1 n. |0 f( \6 E%% 参数设置
distance = 20.0e3;                                                          % 目标初始点到雷达的距离20000米
4 _2 }3 c+ e; Q6 R5 `; J8 ~$ w, ?f0 = 340e9;                                                                % 载波频率
' s6 Y+ A& ^1 ^' a1 b%fs = 1000*40;
0 q0 ~) c0 N/ d) r' w7 J; n% @! T%fs = 1000*60;                                                               % 采样频率
; L$ w  Z) K$ O0 B3 z0 w+ tfs = 1000*50;
6 K+ I8 w) b1 Z; g% f0 = 10e9;                                                                % 载波频率
% fs = 1000*16;                                                               % 采样频率
7 U" I7 ]9 `/ J6 [0 ~2 W" o) ut_dur=1;
& L- m; V. b- \6 Y7 _4 H%t_dur = 0.2;                                                                 % 观测时间长度
t_start = 0;                                                               % 观测时间起点
  C0 n" z/ w% f5 h5 ?M = t_dur*fs;                                                              % 时间点数
t=linspace(t_start+1/fs,t_start+t_dur,M)';                                  % 各采样时刻
1 X3 B, R' l; q3 C& S* N, Xaccelerate_1 = 9.6;                                                       % 加速度
  {" p1 ]. w, V0 [( Z4 D2 ]accelerate_2 = 9.2;
c = 3.0e8;                                                                 % 电磁波传播速度(m/s)
0 V% {: W  k7 Z- g, `1 P* ]% mwave_length = c/f0;                                                        % 定义波长
9 T, g! r2 F  E$ ?, ]N=200;
- [4 P( z4 A- q( [) g, `! h+ }3 cv = 0;                                                                   % 一阶运动速度
echo_acc = exp(j*4*pi/wave_length*(distance+accelerate_1*t.^2))+exp(j*4*pi/wave_length*(distance+accelerate_2*t.^2));%匀加速目标回波信号
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% Gabor算法
* m0 n$ `2 ~' N$ q6 {6 f4 P%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
" U3 I  ]6 i" a$ B5 rn = 4;
+ j( H- g, _6 n- j* K# ?7 \F = N/n;                                                                   % Gabor变换频率点数为fs/n
[tfr_acc_gab,t_gab,f] = tfrgabor(echo_acc,N,F);
' v. v/ J/ A3 e8 ztfr_acc_gab = fftshift(tfr_acc_gab,1);
figure
colormap(gray(256));
/ y  g& F1 j/ r8 T# v$ C" n& T* y. Z9 Oimagesc(t,linspace(0,fs/2-n,F),255-abs(tfr_acc_gab));
axis([t_start t_start+t_dur 0 fs/2]);
set(gca,'YDir','normal')
8 u" @/ h+ {" h' y9 B$ O% egrid on;
xlabel('时间');ylabel('频率'); title('匀加速运动目标微多普勒时频像');
大家帮我看一下是怎么回事，不报错，有几个warning，说什么不完全匹配，然后即一直卡住，半个小时之后才有结果。

tiansifei1

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