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

tiansifei1

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

tiansifei1

先附上匀加速情况下的代码：
- \0 I  V' M% k( {clear all; close all;
9 e0 j0 F( P, I3 ]% C8 ~7 p%% 参数设置
, I$ Z: v* z6 ~1 Z8 n* ~3 C' o8 ldistance = 20.0e3;                                                          % 目标初始点到雷达的距离20000米
1 s  y) D3 X" n  S% F+ x9 mf0 = 340e9;                                                                % 载波频率
%fs = 1000*40;
& d; ?. ]) Y$ d6 ]  k; J%fs = 1000*60;                                                               % 采样频率
; Y7 v3 I6 a. jfs = 1000*50;
% f0 = 10e9;                                                                % 载波频率
% fs = 1000*16;                                                               % 采样频率
7 J9 T& j: Y/ b8 T" Lt_dur=1;
( I: V% p; T6 [8 y%t_dur = 0.2;                                                                 % 观测时间长度
t_start = 0;                                                               % 观测时间起点
5 x2 G9 O) j1 P. J4 LM = t_dur*fs;                                                              % 时间点数
t=linspace(t_start+1/fs,t_start+t_dur,M)';                                  % 各采样时刻
6 d2 @! b2 @5 ?0 o- Kaccelerate_1 = 9.6;                                                       % 加速度
accelerate_2 = 9.2;
1 V; b! v( M# W3 m$ Jc = 3.0e8;                                                                 % 电磁波传播速度(m/s)
2 s/ u5 C4 z+ Mwave_length = c/f0;                                                        % 定义波长
: K) M* o- f! I% k, N9 Y/ N7 aN=200;
( a: k+ ^& s! }, p" |9 Nv = 0;                                                                   % 一阶运动速度
( U# U9 x, z) Q, w+ R- c) ?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算法
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
n = 4;
, d* {  r8 D7 |& L) d( nF = N/n;                                                                   % Gabor变换频率点数为fs/n
[tfr_acc_gab,t_gab,f] = tfrgabor(echo_acc,N,F);
/ T: I: R# U! S: Q/ o) Mtfr_acc_gab = fftshift(tfr_acc_gab,1);
% \; z0 J/ m0 _+ J" a* [( L. j; wfigure
8 v  I6 J/ w! Y3 n9 M7 g3 F5 |8 e/ ~+ ycolormap(gray(256));
imagesc(t,linspace(0,fs/2-n,F),255-abs(tfr_acc_gab));
% F) w1 M% X& a9 e: Caxis([t_start t_start+t_dur 0 fs/2]);
set(gca,'YDir','normal')
7 k& D1 T* @+ H( }4 Cgrid on;
xlabel('时间');ylabel('频率'); title('匀加速运动目标微多普勒时频像');
* }' ?6 Z3 U4 z3 H& @8 h大家帮我看一下是怎么回事，不报错，有几个warning，说什么不完全匹配，然后即一直卡住，半个小时之后才有结果。

tiansifei1

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