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

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

[复制链接]

该用户从未签到

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

EDA365欢迎您登录!

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

x
产生如下4种信号:频率为5HZ,采样频率为100HZ的正弦信号x;
* U) l- u$ Y. X- S将信号x混有随机噪声后的带噪信号;
5 e: B* m) M" v1 D* K! Ex的频谱信号;3 Q9 Q0 L0 a# K; @
带噪信号的频谱信号。
: F( x3 j6 Q9 _! [$ C. U& K设计良好的人机界面,每个按钮对应一种波形/ Z; _7 |+ {. v. s. H

, C: R  `6 G+ ?. K+ U9 \
  • TA的每日心情
    开心
    2023-6-2 15:15
  • 签到天数: 1 天

    [LV.1]初来乍到

    2#
    发表于 2022-8-11 13:06 | 只看该作者
    wp=0.2*pi;ws=0.3*pi;%你没有给出指标,可以自己改
    8 }; r. E, b! i! t/ T+ a  W+ L& r5 iFs=20000;T=1/Fs;%采样频率可以自己定
    3 M2 H+ K. `0 z1 SOmegaP=(2/T)*tan(wp/2);
    6 u* d- Q1 k+ p/ hOmegaS=(2/T)*tan(ws/2);
    . a: P: e$ V  q" {rp=1;rs=15;
    8 V# Z9 Q8 j: L; }. xripple=10^(-rp/20);attn=10^(-rs/20);/ w" {+ b  |  `( T; n+ d! C  Y
    [n,wn]=buttord(OmegaP,OmegaS,rp,rs,'s');
      q/ ]) p0 B5 n4 t; N1 J4 p[z,p,k]=Buttap(n);
    - ~! C# a6 _! d2 ][b,a]=zp2tf(z,p,k);
      Q& y+ i+ `# U6 @  ][bt,at]=lp2lp(b,a,wn);
    ' {! t8 V+ ]5 X$ d) x- i[b,a]=bilinear(bt,at,Fs);
    8 a( d- o/ Q, E3 P[db,mag,pha,grd,w]=freqz_m(b,a);0 [% B4 A) ~4 r( u4 ?9 w- ]
    %
    9 J! }. _! h( V( _5 H%下面绘出各条曲线* l" ?9 O1 _9 o2 J: [/ L
    subplot(2,2,1);plot(w/pi,mag);title('幅频特性');9 `9 ?! b3 D" o! Y/ n
    xlabel('w(/pi)');ylabel('|H(jw)|');, E9 \8 F; w/ m
    axis([0,1,0,1.1]);
    / K& C% ?6 F$ \% Q! X* z3 h4 Qset(gca,'XTickMode','manual','XTick',[0 0.2 0.3 1]);
    & d' @: I( P% [4 @( v, G3 ]set(gca,'YTickMode','manual','YTick',[0 attn ripple 1]);grid
    5 Z: s! v/ h5 rsubplot(2,2,2);plot(w/pi,db);title('幅频特性(db)');
    ) `# k8 y0 Y0 ]# `1 u% f" ]7 k: Vxlabel('w(/pi)');ylabel('dB');- p: i. F5 \) f. _
    axis([0,1,-30,5]);+ M1 g/ _+ C; h: |% u9 Z- P
    set(gca,'XTickMode','manual','XTick',[0 0.2 0.3 1]);
    ' C) M9 Z+ D  j+ B, Aset(gca,'YTickMode','manual','YTick',[-60 -as -rp 0]);grid
    * V7 g# R" o4 w9 O# ksubplot(2,2,3);plot(w/pi,pha/pi);title('相频特性');# q6 w. q# Z) [" Y. d9 W) r, ?  Z* [
    xlabel('w(/pi)');ylabel('pha(/pi)');
    ' S3 ], |! ^3 Saxis([0,1,-1,1]);
    / v/ N$ ^# z; p% z( vset(gca,'XTickMode','manual','XTick',[0 0.2 0.3 1]);grid
    & L6 v- |4 p& \/ t$ ?; Usubplot(2,2,4);plot(w/pi,grd);title('群延时');2 V, ^4 c2 X( N4 R5 M5 L! L
    xlabel('w(/pi)');ylabel('Sample');
    4 X8 i& e+ n: ]' E# E/ Qaxis([0,1,0,15]);
    * N. k4 ?4 O6 Z0 Tset(gca,'XTickMode','manual','XTick',[0 0.2 0.3 1]);grid
    3 m2 F2 C- f' @; x& N%b,a是滤波器的系数,系统函数根据b和a写出。: o6 r) K" e/ ^. b# {! H1 a
    %若程序运行有问题可以Hi我
    2 j% b' p' s! M* g

    该用户从未签到

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

    该用户从未签到

    4#
    发表于 2022-8-11 16:04 | 只看该作者
    楼上说的很不错的
    8 G- e8 z$ g0 E7 b/ J+ @
    您需要登录后才可以回帖 登录 | 注册

    本版积分规则

    关闭

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

    EDA365公众号

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

    GMT+8, 2025-6-25 15:15 , Processed in 0.062500 second(s), 23 queries , Gzip On.

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

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

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