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

tiansifei1

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

tiansifei1

先附上匀加速情况下的代码：
/ S# g9 z+ P$ {9 u9 V, ]) ^clear all; close all;
4 \- y' w: f/ ~; s6 ]# v% a%% 参数设置
distance = 20.0e3;                                                          % 目标初始点到雷达的距离20000米
f0 = 340e9;                                                                % 载波频率
%fs = 1000*40;
%fs = 1000*60;                                                               % 采样频率
fs = 1000*50;
: B! E! `2 N" E4 }& R, u/ h3 ?% f0 = 10e9;                                                                % 载波频率
8 H( q- T* m/ N( E6 I% fs = 1000*16;                                                               % 采样频率
; R4 G5 n  p. t2 |) K5 i$ g1 ct_dur=1;
' O" K  ?4 b5 F# r( s%t_dur = 0.2;                                                                 % 观测时间长度
7 b* \4 n# @0 C$ \t_start = 0;                                                               % 观测时间起点
M = t_dur*fs;                                                              % 时间点数
4 a  b- h( H* s# Pt=linspace(t_start+1/fs,t_start+t_dur,M)';                                  % 各采样时刻
( y" i* w% n2 `/ Z# d' @accelerate_1 = 9.6;                                                       % 加速度
# L) u. C, p, q) n: B% X/ ~# O5 vaccelerate_2 = 9.2;
* Q0 T. v* @9 H$ T* x4 F- Uc = 3.0e8;                                                                 % 电磁波传播速度(m/s)
# p! O- e- m* R/ L: T: @0 Bwave_length = c/f0;                                                        % 定义波长
! o0 S" ]) {5 C( T" UN=200;
  i, `$ P- s( ]5 t% W% M7 pv = 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算法
8 X6 W! t* H4 @1 f' r1 ~! O2 A%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
n = 4;
* O* L+ D1 Y9 _# `! e( A0 V& K( g# H; eF = N/n;                                                                   % Gabor变换频率点数为fs/n
4 t0 F# r; g( V& R/ ]; @7 p[tfr_acc_gab,t_gab,f] = tfrgabor(echo_acc,N,F);
tfr_acc_gab = fftshift(tfr_acc_gab,1);
figure
, `/ C$ i' s  ~: v- M) Kcolormap(gray(256));
, \, R5 Y+ u, C$ l8 ?0 z9 z) ?imagesc(t,linspace(0,fs/2-n,F),255-abs(tfr_acc_gab));
! x3 i( h) P% T& X( |, Haxis([t_start t_start+t_dur 0 fs/2]);
- R  {8 U9 A5 w* oset(gca,'YDir','normal')
5 y/ R7 M$ U+ O4 n# B2 j; U. Bgrid on;
xlabel('时间');ylabel('频率'); title('匀加速运动目标微多普勒时频像');
* o5 ]5 Z1 e& N大家帮我看一下是怎么回事，不报错，有几个warning，说什么不完全匹配，然后即一直卡住，半个小时之后才有结果。

tiansifei1

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