数学建模社区-数学中国

标题: 用matlab时频工具箱进行不同频段下的微多普勒时频分析，作图并比较，为什么这么卡？ [打印本页]
作者: tiansifei1    时间: 2014-5-28 13:02
标题: 用matlab时频工具箱进行不同频段下的微多普勒时频分析，作图并比较，为什么这么卡？
这两天赶毕业设计，要用matlab时频工具箱进行不同频段下的微多普勒时频分析，针对匀加速、振动、转动等情况作出对应目标的微多普勒时频像。弄了几个才三四十行的小程序，结果matlab一运行就卡住了，至少半个小时才能跑出结果的图像。这是怎么回事啊？是我的程序有问题？还是操作系统、matlab软件有问题？我matlab用的是7.0（以前想装成2010版的，结果安装过程中要报错，又换回7.0的了）。操作系统用的是XP，内存2G.
' }# C+ R; R. j1 t 系统属性.jpg
作者: tiansifei1    时间: 2014-5-28 13:04
先附上匀加速情况下的代码：
" ^2 G6 j1 t0 |" L& @clear all; close all;. ]) Y0 x) n' v# r' I& T' P2 M
%% 参数设置
/ y1 @" O, {6 f7 r6 R  Fdistance = 20.0e3;                                                          % 目标初始点到雷达的距离20000米
7 w& x% R- ^- v# M9 Hf0 = 340e9;                                                                % 载波频率
% N+ i* `3 o7 k# X# Q: Q%fs = 1000*40;
5 g0 e4 k. o/ Z* M$ ^%fs = 1000*60;                                                               % 采样频率
  [* z' m/ `# c: j! l9 n) Cfs = 1000*50;
$ d7 _$ n! s" R3 R% ?% f0 = 10e9;                                                                % 载波频率
% H- N. s% @) i/ G$ Y( s4 U6 o% fs = 1000*16;                                                               % 采样频率
" }! d, \; }% C1 N  rt_dur=1;, N4 q0 ^' v9 f# `- H6 C' n
%t_dur = 0.2;                                                                 % 观测时间长度
, }/ v3 f; L/ J$ @t_start = 0;                                                               % 观测时间起点9 }4 D; _. Z5 L  B8 W
M = t_dur*fs;                                                              % 时间点数# R' Q6 I' L0 H: e" u. ^
t=linspace(t_start+1/fs,t_start+t_dur,M)';                                  % 各采样时刻
" f$ }6 [) h% b, g. O0 Jaccelerate_1 = 9.6;                                                       % 加速度
7 @( Z4 {( c" D9 |5 d' K7 Caccelerate_2 = 9.2;' S; L( V) ]: X, v) m
c = 3.0e8;                                                                 % 电磁波传播速度(m/s)
* ]6 [8 W, X7 F8 Wwave_length = c/f0;                                                        % 定义波长
& W9 Y2 \1 A+ ^4 oN=200;
- ~$ x) o! y: u! W1 Uv = 0;                                                                   % 一阶运动速度
, N5 U& ], P+ a- ]echo_acc = exp(j*4*pi/wave_length*(distance+accelerate_1*t.^2))+exp(j*4*pi/wave_length*(distance+accelerate_2*t.^2));%匀加速目标回波信号
2 J: u1 G& H( N' }8 Q%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
9 i: P; A. ^9 k8 x* ~2 g# F$ b! u%% Gabor算法
) ^7 a6 B/ j; Y) `5 G! A) W%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
' f& z9 t) Z7 H* M: D( ln = 4;0 Q* c; ^, R5 _
F = N/n;                                                                   % Gabor变换频率点数为fs/n
8 y8 t2 N7 I* l' b[tfr_acc_gab,t_gab,f] = tfrgabor(echo_acc,N,F);' J8 t, ?$ \4 U, `
tfr_acc_gab = fftshift(tfr_acc_gab,1);8 Q& R6 F+ j: P) y8 a; l2 i, S! ~. l
figure  x" h7 D- A3 j; I7 Z! c
colormap(gray(256));) X; r+ Z& j0 e
imagesc(t,linspace(0,fs/2-n,F),255-abs(tfr_acc_gab));
7 G: g3 g% i; j/ `2 U3 J! Maxis([t_start t_start+t_dur 0 fs/2]);
! g7 M! f$ J2 M- f3 E9 Z$ ^6 Bset(gca,'YDir','normal')
7 v4 L+ K3 T/ G# s1 B! V. \grid on;& `7 H7 d" q& f; @
xlabel('时间');ylabel('频率'); title('匀加速运动目标微多普勒时频像');
' P6 N, ?  V, x8 k大家帮我看一下是怎么回事，不报错，有几个warning，说什么不完全匹配，然后即一直卡住，半个小时之后才有结果。
作者: tiansifei1    时间: 2014-5-28 13:21
帮我跑一下这个程序也行，只希望快点出图，不要让我再等半个小时了！！



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