QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 3345|回复: 2
打印 上一主题 下一主题

用matlab时频工具箱进行不同频段下的微多普勒时频分析,作图并比较,为什么这么卡?

[复制链接]
字体大小: 正常 放大

9

主题

4

听众

1673

积分

升级  67.3%

  • TA的每日心情
    奋斗
    2018-4-11 22:36
  • 签到天数: 490 天

    [LV.9]以坛为家II

    邮箱绑定达人

    跳转到指定楼层
    1#
    发表于 2014-5-28 13:02 |只看该作者 |倒序浏览
    |招呼Ta 关注Ta |邮箱已经成功绑定
    这两天赶毕业设计,要用matlab时频工具箱进行不同频段下的微多普勒时频分析,针对匀加速、振动、转动等情况作出对应目标的微多普勒时频像。弄了几个才三四十行的小程序,结果matlab一运行就卡住了,至少半个小时才能跑出结果的图像。这是怎么回事啊?是我的程序有问题?还是操作系统、matlab软件有问题?我matlab用的是7.0(以前想装成2010版的,结果安装过程中要报错,又换回7.0的了)。操作系统用的是XP,内存2G.9 ^1 A) z5 ^6 v0 i2 ^
    系统属性.jpg
    zan
    转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信
    暂无签名。

    9

    主题

    4

    听众

    1673

    积分

    升级  67.3%

  • TA的每日心情
    奋斗
    2018-4-11 22:36
  • 签到天数: 490 天

    [LV.9]以坛为家II

    邮箱绑定达人

    先附上匀加速情况下的代码:
    8 u# P( D6 ^2 D! I( ?- W, K  S. bclear all; close all;( t5 W4 E2 R  T. k4 [. C( ?3 s
    %% 参数设置4 |% C; j; q$ V  X
    distance = 20.0e3;                                                          % 目标初始点到雷达的距离20000米
    . F" a$ w* T5 K. lf0 = 340e9;                                                                % 载波频率
    , M$ t: o% O+ W. `" A%fs = 1000*40;
      Y* U% {/ L6 i%fs = 1000*60;                                                               % 采样频率
    3 a; M$ n0 h8 B7 Xfs = 1000*50;/ C0 G5 p. N' x
    % f0 = 10e9;                                                                % 载波频率1 o! B" n5 C. z
    % fs = 1000*16;                                                               % 采样频率
    ' j& U# v! H3 w3 E& c) M- nt_dur=1;
    . D% K% R. M) H6 R9 g0 G%t_dur = 0.2;                                                                 % 观测时间长度
    : }7 ]* b3 |; Ot_start = 0;                                                               % 观测时间起点7 I$ p( m- d. e
    M = t_dur*fs;                                                              % 时间点数
      U8 a5 X& J* D& E1 Y+ _t=linspace(t_start+1/fs,t_start+t_dur,M)';                                  % 各采样时刻
      L# R5 ~$ c; ?5 saccelerate_1 = 9.6;                                                       % 加速度
    8 Q! G3 a: G$ Eaccelerate_2 = 9.2;
    8 d0 G: g4 E) L3 s( C  a$ Y! Wc = 3.0e8;                                                                 % 电磁波传播速度(m/s)/ ~8 Q3 I, n5 T2 @) v, ^) e, f1 G
    wave_length = c/f0;                                                        % 定义波长
    8 L4 S3 M! a) _* I+ D0 ^1 wN=200;
    + f: m! k  G, s, e; t& tv = 0;                                                                   % 一阶运动速度
    ; x" R+ T5 c* Q! K& |echo_acc = exp(j*4*pi/wave_length*(distance+accelerate_1*t.^2))+exp(j*4*pi/wave_length*(distance+accelerate_2*t.^2));%匀加速目标回波信号$ p* ^  o0 Z; l6 o1 u
    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    0 T* `7 d. g4 \2 v# D; ~%% Gabor算法
    : \1 C* `& R: U/ Z. ^' x%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    6 ?) P* y: y1 w  _6 tn = 4;
    ' A; ^8 g% S" uF = N/n;                                                                   % Gabor变换频率点数为fs/n
    6 d: @) c% E: s9 Y5 n. N/ k[tfr_acc_gab,t_gab,f] = tfrgabor(echo_acc,N,F);# X0 f$ a4 o7 R* ^0 l
    tfr_acc_gab = fftshift(tfr_acc_gab,1);
    8 }. S4 i% ?* z* g* p4 a3 s, V! ^figure
    9 n4 c1 M( k, J  s- O! Ucolormap(gray(256));9 t# \2 s; q6 g4 y
    imagesc(t,linspace(0,fs/2-n,F),255-abs(tfr_acc_gab));   D! C% e* ~+ `' j' Z5 P( `
    axis([t_start t_start+t_dur 0 fs/2]);
    : W( ]$ q0 j" k/ W# f9 o. ^set(gca,'YDir','normal')
    8 \* t- }8 D# E# \% D) A, ogrid on;
    $ A) T$ q5 j; {) Dxlabel('时间');ylabel('频率'); title('匀加速运动目标微多普勒时频像');
    " {# i' @, ^, r大家帮我看一下是怎么回事,不报错,有几个warning,说什么不完全匹配,然后即一直卡住,半个小时之后才有结果。
    暂无签名。
    回复

    使用道具 举报

    9

    主题

    4

    听众

    1673

    积分

    升级  67.3%

  • TA的每日心情
    奋斗
    2018-4-11 22:36
  • 签到天数: 490 天

    [LV.9]以坛为家II

    邮箱绑定达人

    回复

    使用道具 举报

    您需要登录后才可以回帖 登录 | 注册地址

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

    关于我们| 联系我们| 诚征英才| 对外合作| 产品服务| QQ

    手机版|Archiver| |繁體中文 手机客户端  

    蒙公网安备 15010502000194号

    Powered by Discuz! X2.5   © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 )     论坛法律顾问:王兆丰

    GMT+8, 2026-7-9 10:08 , Processed in 0.546046 second(s), 64 queries .

    回顶部