QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 3350|回复: 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.
    . l' Q+ m4 a! J5 J' j  ?1 S 系统属性.jpg
    zan
    转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信
    暂无签名。

    9

    主题

    4

    听众

    1673

    积分

    升级  67.3%

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

    [LV.9]以坛为家II

    邮箱绑定达人

    先附上匀加速情况下的代码:
    2 U  Q$ [  G, w. m. d, hclear all; close all;+ B( {5 h) u) U
    %% 参数设置
    , s: o5 J$ O& G: [" @% I* S+ i' mdistance = 20.0e3;                                                          % 目标初始点到雷达的距离20000米4 u5 n1 Q3 I+ M1 b
    f0 = 340e9;                                                                % 载波频率
    $ S& l# R3 t5 }%fs = 1000*40;
    ) E9 i5 t1 [/ R%fs = 1000*60;                                                               % 采样频率) m$ C( `. C. Q% k+ s
    fs = 1000*50;5 j8 ]8 _: V/ Y
    % f0 = 10e9;                                                                % 载波频率4 F; k) S5 J+ x& b+ Y0 t
    % fs = 1000*16;                                                               % 采样频率0 B& [! B7 G; _! Y9 t2 U7 `* N
    t_dur=1;, J1 ?& W# M9 ]) }- Z
    %t_dur = 0.2;                                                                 % 观测时间长度' r# a. W9 Y4 F" o
    t_start = 0;                                                               % 观测时间起点0 o% m0 O' W: ]# J  W* `7 C3 H" h
    M = t_dur*fs;                                                              % 时间点数2 s) C; V' w/ f0 Y% i8 _: W$ G
    t=linspace(t_start+1/fs,t_start+t_dur,M)';                                  % 各采样时刻$ S6 `* [- C: X2 J) ]
    accelerate_1 = 9.6;                                                       % 加速度
    ( a, Y( @9 J* N) {. P& haccelerate_2 = 9.2;! j# `% g5 _9 R8 m  n4 P
    c = 3.0e8;                                                                 % 电磁波传播速度(m/s): C3 R0 _" }' W; P& u" G' m
    wave_length = c/f0;                                                        % 定义波长2 ~' E6 {0 \$ N4 y
    N=200;
    , ?: A, q) d) Z/ m0 Xv = 0;                                                                   % 一阶运动速度) Z3 I- m  v9 ~" I! n6 z
    echo_acc = exp(j*4*pi/wave_length*(distance+accelerate_1*t.^2))+exp(j*4*pi/wave_length*(distance+accelerate_2*t.^2));%匀加速目标回波信号
    ( N0 u. }7 E' A9 H5 M3 j, t5 ^$ f%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    * A) }; }/ a/ L: m; l  e. `%% Gabor算法
    0 G! B7 W8 u, Z4 n+ [! h%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%2 G, Y0 s7 b' P$ C+ b& I3 c* b
    n = 4;5 S4 f8 I/ P; I9 y% ^+ q
    F = N/n;                                                                   % Gabor变换频率点数为fs/n# |6 U! U- e, A9 b/ f8 I
    [tfr_acc_gab,t_gab,f] = tfrgabor(echo_acc,N,F);- M3 |% x1 N  t3 F+ j
    tfr_acc_gab = fftshift(tfr_acc_gab,1);$ R$ k- l( C  E0 M$ R- b) A& X
    figure
    0 C" I/ k) d$ w: p9 G. C1 Vcolormap(gray(256));$ P+ P' p& l( n4 @8 t' I) u& m- u6 v( i
    imagesc(t,linspace(0,fs/2-n,F),255-abs(tfr_acc_gab));
    : H0 g. t4 e4 B9 L% eaxis([t_start t_start+t_dur 0 fs/2]);1 P) P; w& x; ?. n) k' F; j
    set(gca,'YDir','normal')
    : a3 x) c9 I8 j5 a2 D. ugrid on;- E. A# R1 T9 ]; E: J
    xlabel('时间');ylabel('频率'); title('匀加速运动目标微多普勒时频像');' Z% Q* M$ W; M- x
    大家帮我看一下是怎么回事,不报错,有几个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 11:03 , Processed in 0.333283 second(s), 63 queries .

    回顶部