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

tiansifei1

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

tiansifei1

先附上匀加速情况下的代码：
clear all; close all;
%% 参数设置
6 m3 v* k$ D5 j  odistance = 20.0e3;                                                          % 目标初始点到雷达的距离20000米
) C$ E9 B1 k+ k& Hf0 = 340e9;                                                                % 载波频率
) M: h; V$ \5 O%fs = 1000*40;
. ~7 p# h8 K2 q( p0 D9 ~5 H3 Z%fs = 1000*60;                                                               % 采样频率
fs = 1000*50;
% f0 = 10e9;                                                                % 载波频率
% fs = 1000*16;                                                               % 采样频率
4 L& F/ X- c7 ]5 u* o& Wt_dur=1;
6 W: m4 I) n3 e% n- g; u; G% P%t_dur = 0.2;                                                                 % 观测时间长度
- R$ g+ ^5 I- }: d; H; x) v7 yt_start = 0;                                                               % 观测时间起点
M = t_dur*fs;                                                              % 时间点数
t=linspace(t_start+1/fs,t_start+t_dur,M)';                                  % 各采样时刻
7 w5 Y5 W  _6 x$ k. u% baccelerate_1 = 9.6;                                                       % 加速度
  J6 o2 T4 t0 n% [3 v0 Eaccelerate_2 = 9.2;
* i+ A$ O% H# C9 `# _! Kc = 3.0e8;                                                                 % 电磁波传播速度(m/s)
- }# K- z7 w8 q* x% Hwave_length = c/f0;                                                        % 定义波长
7 ?% I8 g0 n% N3 }3 c& h$ X  F* |N=200;
8 U: g# m4 Y+ o& rv = 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));%匀加速目标回波信号
, B( K" i6 Y& R7 L3 }" n%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
5 a& W9 E0 ]6 ^' a& d2 N%% Gabor算法
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
n = 4;
2 j: O, l5 G  \. x2 ]) d  S) mF = N/n;                                                                   % Gabor变换频率点数为fs/n
2 A4 B8 T' m3 g0 d[tfr_acc_gab,t_gab,f] = tfrgabor(echo_acc,N,F);
tfr_acc_gab = fftshift(tfr_acc_gab,1);
figure
" x4 M. ]1 W- }" t3 v. y% scolormap(gray(256));
imagesc(t,linspace(0,fs/2-n,F),255-abs(tfr_acc_gab));
8 E0 ^/ h6 m% }: ?2 ^9 zaxis([t_start t_start+t_dur 0 fs/2]);
; f) f& K* h% t% i% Xset(gca,'YDir','normal')
( d2 L7 o. Q7 F1 C$ o- Qgrid on;
xlabel('时间');ylabel('频率'); title('匀加速运动目标微多普勒时频像');
( I9 K0 T  W) D2 W6 G( J1 Q大家帮我看一下是怎么回事，不报错，有几个warning，说什么不完全匹配，然后即一直卡住，半个小时之后才有结果。

tiansifei1

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