用matlab时频工具箱进行不同频段下的微多普勒时频分析,作图并比较,为什么这么卡?
这两天赶毕业设计,要用matlab时频工具箱进行不同频段下的微多普勒时频分析,针对匀加速、振动、转动等情况作出对应目标的微多普勒时频像。弄了几个才三四十行的小程序,结果matlab一运行就卡住了,至少半个小时才能跑出结果的图像。这是怎么回事啊?是我的程序有问题?还是操作系统、matlab软件有问题?我matlab用的是7.0(以前想装成2010版的,结果安装过程中要报错,又换回7.0的了)。操作系统用的是XP,内存2G.先附上匀加速情况下的代码:
clear all; close all;
%% 参数设置
distance = 20.0e3; % 目标初始点到雷达的距离20000米
f0 = 340e9; % 载波频率
%fs = 1000*40;
%fs = 1000*60; % 采样频率
fs = 1000*50;
% f0 = 10e9; % 载波频率
% fs = 1000*16; % 采样频率
t_dur=1;
%t_dur = 0.2; % 观测时间长度
t_start = 0; % 观测时间起点
M = t_dur*fs; % 时间点数
t=linspace(t_start+1/fs,t_start+t_dur,M)'; % 各采样时刻
accelerate_1 = 9.6; % 加速度
accelerate_2 = 9.2;
c = 3.0e8; % 电磁波传播速度(m/s)
wave_length = c/f0; % 定义波长
N=200;
v = 0; % 一阶运动速度
echo_acc = exp(j*4*pi/wave_length*(distance+accelerate_1*t.^2))+exp(j*4*pi/wave_length*(distance+accelerate_2*t.^2));%匀加速目标回波信号
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% Gabor算法
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
n = 4;
F = N/n; % Gabor变换频率点数为fs/n
= tfrgabor(echo_acc,N,F);
tfr_acc_gab = fftshift(tfr_acc_gab,1);
figure
colormap(gray(256));
imagesc(t,linspace(0,fs/2-n,F),255-abs(tfr_acc_gab));
axis();
set(gca,'YDir','normal')
grid on;
xlabel('时间');ylabel('频率'); title('匀加速运动目标微多普勒时频像');
大家帮我看一下是怎么回事,不报错,有几个warning,说什么不完全匹配,然后即一直卡住,半个小时之后才有结果。 帮我跑一下这个程序也行,只希望快点出图,不要让我再等半个小时了!!
页:
[1]