找回密码
 注册
关于网站域名变更的通知
查看: 286|回复: 3
打印 上一主题 下一主题

[仿真讨论] 用matlab进行时频信号分析器设计

[复制链接]

该用户从未签到

跳转到指定楼层
1#
发表于 2022-8-11 11:20 | 只看该作者 回帖奖励 |倒序浏览 |阅读模式

EDA365欢迎您登录!

您需要 登录 才可以下载或查看,没有帐号?注册

x
产生如下4种信号:频率为5HZ,采样频率为100HZ的正弦信号x;, q7 o& O8 k, D; f' k4 P# c
将信号x混有随机噪声后的带噪信号;; W* w) B0 h& V8 Z) o) j; p, \
x的频谱信号;( k  r, K1 W( @1 ^
带噪信号的频谱信号。
! e' W. j$ J: b" l' [4 p8 X: F设计良好的人机界面,每个按钮对应一种波形6 [; B, h5 n5 `2 G5 o5 ?

" j. x6 ?0 b9 r' l, V, M
  • TA的每日心情
    开心
    2023-6-2 15:15
  • 签到天数: 1 天

    [LV.1]初来乍到

    2#
    发表于 2022-8-11 13:06 | 只看该作者
    wp=0.2*pi;ws=0.3*pi;%你没有给出指标,可以自己改) b: i; z2 V4 t( D" k
    Fs=20000;T=1/Fs;%采样频率可以自己定( g" h( Q& [- d: J/ v: L2 x  b8 N
    OmegaP=(2/T)*tan(wp/2);0 D5 K1 Q0 G' K( \$ T/ G, x3 f+ J, G
    OmegaS=(2/T)*tan(ws/2);
    3 r2 d0 T, x1 X2 ?' h' ^+ grp=1;rs=15;" a0 _6 D+ D7 x+ U! |
    ripple=10^(-rp/20);attn=10^(-rs/20);: `& @4 H6 i* Z; H2 J+ L+ O* G
    [n,wn]=buttord(OmegaP,OmegaS,rp,rs,'s');8 S( z; F$ j6 j0 C0 s
    [z,p,k]=Buttap(n);! ~3 ~3 m5 l/ V) |* ^- r* V# J
    [b,a]=zp2tf(z,p,k);
    ' W7 c  q7 E- @, v% B/ c7 j[bt,at]=lp2lp(b,a,wn);9 h* {6 m! v: R- v1 W0 s+ f
    [b,a]=bilinear(bt,at,Fs);
    0 f( R) ]4 _8 |6 ~' y0 Z7 l$ l[db,mag,pha,grd,w]=freqz_m(b,a);
    ; i5 m! L( |  K* ~%7 ?* }! U4 y1 H" B8 W2 b& I% S; l
    %下面绘出各条曲线4 D; s& H) V6 v1 ^% Q  u
    subplot(2,2,1);plot(w/pi,mag);title('幅频特性');
    2 N2 m! O1 G& mxlabel('w(/pi)');ylabel('|H(jw)|');) C7 b6 U0 y* [  R
    axis([0,1,0,1.1]);, i$ e7 q7 _- @" J4 o  V8 H
    set(gca,'XTickMode','manual','XTick',[0 0.2 0.3 1]);4 P/ @/ {. n1 ^2 F9 f* J/ a5 t
    set(gca,'YTickMode','manual','YTick',[0 attn ripple 1]);grid' H( z) v4 n  a2 p9 V- m
    subplot(2,2,2);plot(w/pi,db);title('幅频特性(db)');7 A3 O9 `# Y$ q* k- h" K
    xlabel('w(/pi)');ylabel('dB');+ }) T( ]( i: [, S5 ]0 q
    axis([0,1,-30,5]);0 C: r+ L: Q. f0 ]' G
    set(gca,'XTickMode','manual','XTick',[0 0.2 0.3 1]);9 ?6 |& N0 V7 s/ r6 S: r
    set(gca,'YTickMode','manual','YTick',[-60 -as -rp 0]);grid
    # p5 d4 A, S. g( nsubplot(2,2,3);plot(w/pi,pha/pi);title('相频特性');
    4 N$ p/ f3 C6 ~' U/ |& Axlabel('w(/pi)');ylabel('pha(/pi)');
    ' v/ `- h5 T" X# @5 \7 Eaxis([0,1,-1,1]);
    : c* @0 A+ m1 P# X) F% }8 Tset(gca,'XTickMode','manual','XTick',[0 0.2 0.3 1]);grid
    3 Y0 [7 B) p3 Y0 r. tsubplot(2,2,4);plot(w/pi,grd);title('群延时');
    2 U5 b5 d1 O7 C' mxlabel('w(/pi)');ylabel('Sample');! a  i% ~) I, G! G9 V
    axis([0,1,0,15]);/ Y2 F9 }$ i; p# \6 m
    set(gca,'XTickMode','manual','XTick',[0 0.2 0.3 1]);grid
    " G, B( y- M$ c6 B0 F! a# |%b,a是滤波器的系数,系统函数根据b和a写出。
    * {( L3 t- Z" [%若程序运行有问题可以Hi我: h# `7 c. O+ B

    该用户从未签到

    3#
    发表于 2022-8-11 15:30 | 只看该作者
    bootloader是dsp芯片上固化的一段引导程序,dsp程序在ram上掉电易丢,所以启动时运行bootloader将辅存中的程序调入ram中
    ; h& G% G4 c8 D1 l5 s至于dsp的框架结构在ti公司有(有点复杂),开发时通常是先在开发板上仿真然后实践,由于dsp一般用于数字信号处理,离不开matlab软件连调。说起来有点复杂,要熟悉dsp开发流程序最简单的方法就是自己动手实践一下,如做个滤波器,频谱分析器等就可以深刻了解了。
    " ^4 X" X) p$ v' b) A& O3 T

    该用户从未签到

    4#
    发表于 2022-8-11 16:04 | 只看该作者
    楼上说的很不错的
    + P! h; I& Z. H3 O0 |2 e" {0 I: @  U& |
    您需要登录后才可以回帖 登录 | 注册

    本版积分规则

    关闭

    推荐内容上一条 /1 下一条

    EDA365公众号

    关于我们|手机版|EDA365电子论坛网 ( 粤ICP备18020198号-1 )

    GMT+8, 2025-7-11 19:22 , Processed in 0.109375 second(s), 23 queries , Gzip On.

    深圳市墨知创新科技有限公司

    地址:深圳市南山区科技生态园2栋A座805 电话:19926409050

    快速回复 返回顶部 返回列表