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

tiansifei1

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

tiansifei1

先附上匀加速情况下的代码：
clear all; close all;
%% 参数设置
distance = 20.0e3;                                                          % 目标初始点到雷达的距离20000米
f0 = 340e9;                                                                % 载波频率
%fs = 1000*40;
6 n2 Q2 s( u* W  Z6 I( I1 ~%fs = 1000*60;                                                               % 采样频率
fs = 1000*50;
. y# _* B- A; x- ^7 n9 k% f0 = 10e9;                                                                % 载波频率
0 q8 Z8 C, m9 w+ c- Y, h1 f( E% fs = 1000*16;                                                               % 采样频率
t_dur=1;
- s% O$ e3 j7 ^" h2 {, B! v%t_dur = 0.2;                                                                 % 观测时间长度
, ^4 R1 A: |7 {, @+ Qt_start = 0;                                                               % 观测时间起点
M = t_dur*fs;                                                              % 时间点数
t=linspace(t_start+1/fs,t_start+t_dur,M)';                                  % 各采样时刻
5 c) W& }' p) W1 maccelerate_1 = 9.6;                                                       % 加速度
7 l' ], g+ n& @$ e: h5 y) E) Kaccelerate_2 = 9.2;
c = 3.0e8;                                                                 % 电磁波传播速度(m/s)
wave_length = c/f0;                                                        % 定义波长
) ^/ K; A+ j5 Q% l' \9 t$ kN=200;
v = 0;                                                                   % 一阶运动速度
6 z0 S. S2 u8 A7 Z' k- _( Pecho_acc = exp(j*4*pi/wave_length*(distance+accelerate_1*t.^2))+exp(j*4*pi/wave_length*(distance+accelerate_2*t.^2));%匀加速目标回波信号
' O$ g4 d8 r( P5 E%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
* `. Q. w7 G, n# w%% Gabor算法
* C! k; p- q! T%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
/ V5 Y* b0 }" Fn = 4;
7 b: G; T6 }* r( a, f, s+ U* |F = N/n;                                                                   % Gabor变换频率点数为fs/n
& x0 U$ d/ P3 O0 n/ X' \9 [[tfr_acc_gab,t_gab,f] = tfrgabor(echo_acc,N,F);
( p0 y# V7 O* y& H% ^# f' H, u: H8 R1 H, Jtfr_acc_gab = fftshift(tfr_acc_gab,1);
figure
colormap(gray(256));
* d2 k- Z4 ]! U1 W# p; n; @+ Himagesc(t,linspace(0,fs/2-n,F),255-abs(tfr_acc_gab));
axis([t_start t_start+t_dur 0 fs/2]);
' M* w  o( f2 |% z/ o; X+ vset(gca,'YDir','normal')
9 O) T/ I0 R- ]5 ?9 Lgrid on;
6 w1 b0 v. ^1 L) ]( Gxlabel('时间');ylabel('频率'); title('匀加速运动目标微多普勒时频像');
) w" L" D9 o% Q) c$ `大家帮我看一下是怎么回事，不报错，有几个warning，说什么不完全匹配，然后即一直卡住，半个小时之后才有结果。

tiansifei1

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