数学建模社区-数学中国

标题: 用matlab时频工具箱进行不同频段下的微多普勒时频分析，作图并比较，为什么这么卡？ [打印本页]
作者: tiansifei1    时间: 2014-5-28 13:02
标题: 用matlab时频工具箱进行不同频段下的微多普勒时频分析，作图并比较，为什么这么卡？
这两天赶毕业设计，要用matlab时频工具箱进行不同频段下的微多普勒时频分析，针对匀加速、振动、转动等情况作出对应目标的微多普勒时频像。弄了几个才三四十行的小程序，结果matlab一运行就卡住了，至少半个小时才能跑出结果的图像。这是怎么回事啊？是我的程序有问题？还是操作系统、matlab软件有问题？我matlab用的是7.0（以前想装成2010版的，结果安装过程中要报错，又换回7.0的了）。操作系统用的是XP，内存2G.
/ w) t" ]' q7 f6 m" L! b! j 系统属性.jpg
作者: tiansifei1    时间: 2014-5-28 13:04
先附上匀加速情况下的代码：$ U! Z& S* m+ {" J
clear all; close all;% {' R# T; S9 I9 s! j# U
%% 参数设置
3 w$ o: J3 |% C  }- p' c5 bdistance = 20.0e3;                                                          % 目标初始点到雷达的距离20000米4 }: N9 }4 }) ]! @
f0 = 340e9;                                                                % 载波频率/ J6 s6 B) S; C  ?
%fs = 1000*40;1 {0 p/ H0 H! e8 n" \# J  _* k
%fs = 1000*60;                                                               % 采样频率( R& T7 d$ h3 z+ i7 A
fs = 1000*50;5 O2 I/ M1 A7 E5 D( j
% f0 = 10e9;                                                                % 载波频率& {, B4 p6 C; l& z9 I& |
% fs = 1000*16;                                                               % 采样频率
, p% {+ n2 k+ _# M; Kt_dur=1;+ [) ]* @5 r4 u# ?
%t_dur = 0.2;                                                                 % 观测时间长度
8 n$ H; k/ P* }/ G8 ut_start = 0;                                                               % 观测时间起点
( r, H% E7 L+ c, @$ C6 hM = t_dur*fs;                                                              % 时间点数
% c& `1 O& z) N  M$ v8 Y; Pt=linspace(t_start+1/fs,t_start+t_dur,M)';                                  % 各采样时刻) [0 E, y" ^; e5 R
accelerate_1 = 9.6;                                                       % 加速度
" q9 p& X- D1 iaccelerate_2 = 9.2;
( p% f: K5 ^- N; I/ Gc = 3.0e8;                                                                 % 电磁波传播速度(m/s), z  P* @: i- h
wave_length = c/f0;                                                        % 定义波长
! K5 g" A$ W1 eN=200;4 N8 r1 j- i9 H
v = 0;                                                                   % 一阶运动速度7 X. J4 f5 F1 h  t
echo_acc = exp(j*4*pi/wave_length*(distance+accelerate_1*t.^2))+exp(j*4*pi/wave_length*(distance+accelerate_2*t.^2));%匀加速目标回波信号$ m+ |: d  D& x$ s- Q# |  i# U
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5 k5 `% X$ P; C0 Q" x, t3 M
%% Gabor算法
7 a4 ~+ b' [2 w* a9 s' Z%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%/ z' E. H7 j, o; z' [8 [; N+ C
n = 4;
- q6 \$ k" r7 B2 I) _& l# W2 }- TF = N/n;                                                                   % Gabor变换频率点数为fs/n! D8 @$ ?" y- K* ]. w
[tfr_acc_gab,t_gab,f] = tfrgabor(echo_acc,N,F);
6 g( b: x8 I2 m  ^+ X# ?  G- Otfr_acc_gab = fftshift(tfr_acc_gab,1);& j  J8 y8 E6 u, [/ ?
figure
# w, x2 [. t2 c& |1 r& S% ecolormap(gray(256));2 M  \& O% c: }& u8 a/ {- p4 a, l( ?
imagesc(t,linspace(0,fs/2-n,F),255-abs(tfr_acc_gab));
/ O( w: y% j; l) [4 Taxis([t_start t_start+t_dur 0 fs/2]);. M7 w* ?& Y* N
set(gca,'YDir','normal')/ L. ^3 P% f. v9 c" ~" A+ i4 c" o
grid on;
3 `- w1 B3 B' ^9 r+ m0 Hxlabel('时间');ylabel('频率'); title('匀加速运动目标微多普勒时频像');$ O; T, n% h2 [: _( K( l; \. w: k
大家帮我看一下是怎么回事，不报错，有几个warning，说什么不完全匹配，然后即一直卡住，半个小时之后才有结果。
作者: tiansifei1    时间: 2014-5-28 13:21
帮我跑一下这个程序也行，只希望快点出图，不要让我再等半个小时了！！



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