数学建模社区-数学中国

标题: 用matlab时频工具箱进行不同频段下的微多普勒时频分析，作图并比较，为什么这么卡？ [打印本页]
作者: tiansifei1    时间: 2014-5-28 13:02
标题: 用matlab时频工具箱进行不同频段下的微多普勒时频分析，作图并比较，为什么这么卡？
这两天赶毕业设计，要用matlab时频工具箱进行不同频段下的微多普勒时频分析，针对匀加速、振动、转动等情况作出对应目标的微多普勒时频像。弄了几个才三四十行的小程序，结果matlab一运行就卡住了，至少半个小时才能跑出结果的图像。这是怎么回事啊？是我的程序有问题？还是操作系统、matlab软件有问题？我matlab用的是7.0（以前想装成2010版的，结果安装过程中要报错，又换回7.0的了）。操作系统用的是XP，内存2G.
  J4 l+ y' G( m. t1 ~  Q% X 系统属性.jpg
作者: tiansifei1    时间: 2014-5-28 13:04
先附上匀加速情况下的代码：, J% v. T' f  V7 X( Q+ c0 I, U
clear all; close all;
2 Z. P& R9 x* X: H  T0 \+ M7 E%% 参数设置& v' [! ?6 f/ E( ~6 _' }
distance = 20.0e3;                                                          % 目标初始点到雷达的距离20000米. @" o0 |% B  S* t4 ^8 e$ |
f0 = 340e9;                                                                % 载波频率8 `  X9 {! |; k! r6 t4 K
%fs = 1000*40;
% j- V$ p. W2 D3 y%fs = 1000*60;                                                               % 采样频率
3 l4 z. N2 b" ?; E8 ~fs = 1000*50;
5 P4 o7 Q  I5 H9 \7 i% f0 = 10e9;                                                                % 载波频率
. I' B& M3 T0 g$ r% fs = 1000*16;                                                               % 采样频率5 y" l) c" J* ~7 ~2 i" E
t_dur=1;
( j" U5 O4 F  l%t_dur = 0.2;                                                                 % 观测时间长度' r3 R& k; M: D4 o' t9 \* W' ?/ X0 i
t_start = 0;                                                               % 观测时间起点# _# u2 q# x. a1 c+ q
M = t_dur*fs;                                                              % 时间点数( W* n8 w# c0 g1 ]. b
t=linspace(t_start+1/fs,t_start+t_dur,M)';                                  % 各采样时刻' w* P2 V0 _, Y6 p1 o/ Y
accelerate_1 = 9.6;                                                       % 加速度' a) A# Q' ?8 I( B
accelerate_2 = 9.2;
7 ]' t0 I: z: ?' j; @0 z: Q0 Nc = 3.0e8;                                                                 % 电磁波传播速度(m/s)0 q0 q+ E1 x3 ?3 x- j3 V
wave_length = c/f0;                                                        % 定义波长3 k9 y- z! `3 j, `- q2 ^( m4 Q
N=200;4 E, w7 I( t9 J3 b; O; W3 U
v = 0;                                                                   % 一阶运动速度
5 l. h5 N( q% I; @echo_acc = exp(j*4*pi/wave_length*(distance+accelerate_1*t.^2))+exp(j*4*pi/wave_length*(distance+accelerate_2*t.^2));%匀加速目标回波信号
) v& C( w, U; q$ z. C+ e%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%" ]( n) G0 H. \1 @, F! D) w
%% Gabor算法3 k! F, J) h8 J
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
( h6 N, C8 J/ i' @+ Z' e; gn = 4;5 O! i7 u+ ^$ L# a6 r% H$ z
F = N/n;                                                                   % Gabor变换频率点数为fs/n
8 }3 X$ ^  Y* L9 F[tfr_acc_gab,t_gab,f] = tfrgabor(echo_acc,N,F);
, Q4 x+ F; P+ j7 v; ktfr_acc_gab = fftshift(tfr_acc_gab,1);, O$ U. D! j1 b# W3 N
figure
0 @! ?: t& q* F# i3 h0 Jcolormap(gray(256));
% P5 E3 {! ?) pimagesc(t,linspace(0,fs/2-n,F),255-abs(tfr_acc_gab)); 9 Z1 r+ H. ?1 y2 X: o" _
axis([t_start t_start+t_dur 0 fs/2]);
' E! G3 c. z: P, g& `; u) `set(gca,'YDir','normal')
5 p5 N! X1 b5 N) }4 g% K9 g: r( Egrid on;
+ U/ I! f7 \# U) C  o" dxlabel('时间');ylabel('频率'); title('匀加速运动目标微多普勒时频像');  W" w' e3 A; R$ i
大家帮我看一下是怎么回事，不报错，有几个warning，说什么不完全匹配，然后即一直卡住，半个小时之后才有结果。
作者: tiansifei1    时间: 2014-5-28 13:21
帮我跑一下这个程序也行，只希望快点出图，不要让我再等半个小时了！！



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