数学建模社区-数学中国

标题: 用matlab时频工具箱进行不同频段下的微多普勒时频分析，作图并比较，为什么这么卡？ [打印本页]
作者: tiansifei1    时间: 2014-5-28 13:02
标题: 用matlab时频工具箱进行不同频段下的微多普勒时频分析，作图并比较，为什么这么卡？
这两天赶毕业设计，要用matlab时频工具箱进行不同频段下的微多普勒时频分析，针对匀加速、振动、转动等情况作出对应目标的微多普勒时频像。弄了几个才三四十行的小程序，结果matlab一运行就卡住了，至少半个小时才能跑出结果的图像。这是怎么回事啊？是我的程序有问题？还是操作系统、matlab软件有问题？我matlab用的是7.0（以前想装成2010版的，结果安装过程中要报错，又换回7.0的了）。操作系统用的是XP，内存2G.' ~5 w! [8 }, N+ K
系统属性.jpg
作者: tiansifei1    时间: 2014-5-28 13:04
先附上匀加速情况下的代码：1 i3 }2 J% y9 R2 z( Q: j
clear all; close all;! n! {& ?2 {* Y8 G" w3 F
%% 参数设置
  h" S' F, K3 p" ddistance = 20.0e3;                                                          % 目标初始点到雷达的距离20000米+ F9 b1 z( v/ H% _& T/ l( O- @
f0 = 340e9;                                                                % 载波频率$ T- `' B* Z. \" P7 @, \4 N9 v
%fs = 1000*40;7 I& x6 i7 A9 G0 i3 f. K2 y5 M
%fs = 1000*60;                                                               % 采样频率
# y. g: l7 N, M  ^, Qfs = 1000*50;
, i4 r& r% x1 ]; R( u" |% f0 = 10e9;                                                                % 载波频率
0 C5 }  w- m1 d- G: H+ P, ^, |% fs = 1000*16;                                                               % 采样频率6 M8 @1 }$ x1 }8 T; {
t_dur=1;9 x( v! T; \/ S$ K4 i- V) [# t) P
%t_dur = 0.2;                                                                 % 观测时间长度2 P% s5 ]) K7 a/ ~- V7 @7 N! }
t_start = 0;                                                               % 观测时间起点
% d' B2 y' ~: Z# M4 b- QM = t_dur*fs;                                                              % 时间点数
; m6 g3 V4 d4 d* At=linspace(t_start+1/fs,t_start+t_dur,M)';                                  % 各采样时刻4 j) e, @+ R- e* K" U% X
accelerate_1 = 9.6;                                                       % 加速度) Y. j/ [% e5 I$ `7 p0 H! D
accelerate_2 = 9.2;
6 V* P5 J& I9 v+ zc = 3.0e8;                                                                 % 电磁波传播速度(m/s)7 \8 n8 P" ?! y5 o& ]3 U
wave_length = c/f0;                                                        % 定义波长
' u4 E+ Z2 K6 {6 O. V7 d6 V) U; `N=200;1 c/ N* i1 {2 }9 O- z3 Q
v = 0;                                                                   % 一阶运动速度
1 x; {) H6 @$ w) ^" X0 decho_acc = exp(j*4*pi/wave_length*(distance+accelerate_1*t.^2))+exp(j*4*pi/wave_length*(distance+accelerate_2*t.^2));%匀加速目标回波信号7 Y2 w; ~$ V3 A! z' S) F
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
! A" \9 w; d. s, L& z6 Z3 m: N%% Gabor算法
2 x- b" `( _5 F. @4 o%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
4 H  l$ @+ ?$ g" N; x% Fn = 4;- ?. t. [6 s! e
F = N/n;                                                                   % Gabor变换频率点数为fs/n
" n4 X- Z* K. M[tfr_acc_gab,t_gab,f] = tfrgabor(echo_acc,N,F);
4 J7 O% w& p+ P/ w9 |+ Y7 F" ^/ ?tfr_acc_gab = fftshift(tfr_acc_gab,1);8 i% U* Z' w5 z6 J; W( c
figure
$ M7 a3 F/ u2 Kcolormap(gray(256));6 H- W$ c' P# q9 d
imagesc(t,linspace(0,fs/2-n,F),255-abs(tfr_acc_gab));
5 v0 a) N6 R/ ?! ?axis([t_start t_start+t_dur 0 fs/2]);
+ \0 m  }% E, N1 m0 H5 mset(gca,'YDir','normal')7 z! b6 _4 Z, a5 }- ^8 t" ^
grid on;
2 j- a; s3 A4 s) Yxlabel('时间');ylabel('频率'); title('匀加速运动目标微多普勒时频像');
; J' _2 G. J6 f& i. j: v: h  p7 O大家帮我看一下是怎么回事，不报错，有几个warning，说什么不完全匹配，然后即一直卡住，半个小时之后才有结果。
作者: tiansifei1    时间: 2014-5-28 13:21
帮我跑一下这个程序也行，只希望快点出图，不要让我再等半个小时了！！



欢迎光临 数学建模社区-数学中国 (http://www.madio.net/) Powered by Discuz! X2.5