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

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

[复制链接]

该用户从未签到

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

EDA365欢迎您登录!

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

x
产生如下4种信号:频率为5HZ,采样频率为100HZ的正弦信号x;
  U! s, |% a' }8 |" ^6 _将信号x混有随机噪声后的带噪信号;! q: J( d, G+ k$ ?5 {  u
x的频谱信号;
, m! c' \' ]% d& l带噪信号的频谱信号。" e# i! d# B4 Z0 X9 Y
设计良好的人机界面,每个按钮对应一种波形1 a  M* O5 [3 q6 _% K% D

# Z( A& f( q+ f- A
  • TA的每日心情
    开心
    2023-6-2 15:15
  • 签到天数: 1 天

    [LV.1]初来乍到

    2#
    发表于 2022-8-11 13:06 | 只看该作者
    wp=0.2*pi;ws=0.3*pi;%你没有给出指标,可以自己改
    5 ~% F: \4 n. e. t: pFs=20000;T=1/Fs;%采样频率可以自己定
    ; Z0 k7 e8 k" ~# v0 ]' \, n6 ~OmegaP=(2/T)*tan(wp/2);
    4 |! U% B% J  n. D% `OmegaS=(2/T)*tan(ws/2);
    . h) E, X% l. s9 Q* grp=1;rs=15;
    - E3 H) @! C5 n2 u6 yripple=10^(-rp/20);attn=10^(-rs/20);
    / V9 E  E8 t/ ?8 N' Q+ A[n,wn]=buttord(OmegaP,OmegaS,rp,rs,'s');5 J: P& N0 Q" s( \0 O4 r% J! E& q  [
    [z,p,k]=Buttap(n);2 V4 H7 k! s3 a, F6 X
    [b,a]=zp2tf(z,p,k);
    . A, |+ w: L1 W[bt,at]=lp2lp(b,a,wn);" M9 a- R9 f& w# }1 j  P3 @4 ^
    [b,a]=bilinear(bt,at,Fs);
    1 l! |7 d% a: c$ l, |7 ]3 `[db,mag,pha,grd,w]=freqz_m(b,a);
    7 X" v5 n# @/ u- U  f%9 u4 v. W, F) u% l. s
    %下面绘出各条曲线2 y' `8 [# R9 N; D! t8 C4 O
    subplot(2,2,1);plot(w/pi,mag);title('幅频特性');- b% [9 ?  O! A) t5 M6 y' Z
    xlabel('w(/pi)');ylabel('|H(jw)|');
    5 z+ w& z9 W) @. @axis([0,1,0,1.1]);; A5 z7 _1 ^1 D! d
    set(gca,'XTickMode','manual','XTick',[0 0.2 0.3 1]);
    3 M; U" \0 L! e! W$ I/ nset(gca,'YTickMode','manual','YTick',[0 attn ripple 1]);grid
    9 s! O$ Y# [/ k5 w( m, osubplot(2,2,2);plot(w/pi,db);title('幅频特性(db)');
    # p- e" a+ Z2 B: a( N# j( i* ?xlabel('w(/pi)');ylabel('dB');0 |) w# x6 `  }3 f" ?
    axis([0,1,-30,5]);3 B  G. ]5 i; E9 d: X& Q6 G
    set(gca,'XTickMode','manual','XTick',[0 0.2 0.3 1]);
    , M6 F* \8 Z, ^2 N9 f) l7 J* R6 ~0 z. gset(gca,'YTickMode','manual','YTick',[-60 -as -rp 0]);grid
    1 v6 Q$ M& Z4 _$ G$ \subplot(2,2,3);plot(w/pi,pha/pi);title('相频特性');8 T, Q, V8 |7 d. W' B" u% i! W$ r
    xlabel('w(/pi)');ylabel('pha(/pi)');. P" U) K# e, j+ u+ M+ q
    axis([0,1,-1,1]);
    ( m! i% Z* T8 V9 Q: Bset(gca,'XTickMode','manual','XTick',[0 0.2 0.3 1]);grid* r0 F9 N% T3 Y
    subplot(2,2,4);plot(w/pi,grd);title('群延时');
    ! ?8 w9 M% |' Xxlabel('w(/pi)');ylabel('Sample');
    , M; `! x" q3 s4 O  p2 U4 jaxis([0,1,0,15]);
    ; x3 B8 o) p+ n1 u# lset(gca,'XTickMode','manual','XTick',[0 0.2 0.3 1]);grid
    7 ^8 [# Q0 Y" d) m* [%b,a是滤波器的系数,系统函数根据b和a写出。! U$ c, \" `' u: j2 ^. Z
    %若程序运行有问题可以Hi我
    $ w. l  g/ p! N& ^

    该用户从未签到

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

    该用户从未签到

    4#
    发表于 2022-8-11 16:04 | 只看该作者
    楼上说的很不错的. \- P: ^2 C8 Q
    您需要登录后才可以回帖 登录 | 注册

    本版积分规则

    关闭

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

    EDA365公众号

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

    GMT+8, 2025-6-29 04:16 , Processed in 0.078125 second(s), 23 queries , Gzip On.

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

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

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