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

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

[复制链接]

该用户从未签到

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

EDA365欢迎您登录!

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

x
产生如下4种信号:频率为5HZ,采样频率为100HZ的正弦信号x;" S" V2 e) R/ i8 [+ o  f* q
将信号x混有随机噪声后的带噪信号;
2 z3 {  g3 ^: F5 Fx的频谱信号;
% ~) c' `* h" p) T; @带噪信号的频谱信号。  l8 k  Q' e% \4 ~9 X
设计良好的人机界面,每个按钮对应一种波形) s, F) H. t8 i' z" M/ s8 P
5 s" S; o  S9 I1 p8 T
  • TA的每日心情
    开心
    2023-6-2 15:15
  • 签到天数: 1 天

    [LV.1]初来乍到

    2#
    发表于 2022-8-11 13:06 | 只看该作者
    wp=0.2*pi;ws=0.3*pi;%你没有给出指标,可以自己改
    . w* x- T) G0 a+ P/ jFs=20000;T=1/Fs;%采样频率可以自己定1 D: |$ ~7 n# {
    OmegaP=(2/T)*tan(wp/2);
    2 T/ X, O7 Y5 L1 a) ~) TOmegaS=(2/T)*tan(ws/2);/ o# E: y, s. J) K, v3 X7 ^
    rp=1;rs=15;
    9 f' _+ W$ u1 V" J3 H7 kripple=10^(-rp/20);attn=10^(-rs/20);# _& j' a" Z- e0 h
    [n,wn]=buttord(OmegaP,OmegaS,rp,rs,'s');! B/ S6 Z" t# m7 N! V$ S$ C
    [z,p,k]=Buttap(n);
    ! A' S0 B# S3 D8 w[b,a]=zp2tf(z,p,k);  C/ W0 u. s$ {' C
    [bt,at]=lp2lp(b,a,wn);
    ! E- L/ K4 v8 c( I[b,a]=bilinear(bt,at,Fs);
    % \3 V! P; C/ J2 b[db,mag,pha,grd,w]=freqz_m(b,a);9 i. @1 l, R: y- F+ t, c9 D
    %4 x/ x3 }+ z8 X% n
    %下面绘出各条曲线
    8 D1 N+ d! U6 o8 M2 t: O" qsubplot(2,2,1);plot(w/pi,mag);title('幅频特性');& p$ ~1 ~* C2 c
    xlabel('w(/pi)');ylabel('|H(jw)|');7 I8 T- r, c# b
    axis([0,1,0,1.1]);
    7 M: E. M; Y. rset(gca,'XTickMode','manual','XTick',[0 0.2 0.3 1]);+ l2 N9 _0 ]: G4 l
    set(gca,'YTickMode','manual','YTick',[0 attn ripple 1]);grid
    % I. c$ A( B, H  l# d- m3 Z# Csubplot(2,2,2);plot(w/pi,db);title('幅频特性(db)');$ z, L9 [* u9 Z. U
    xlabel('w(/pi)');ylabel('dB');
    7 }! Z+ p* Z: i/ @8 J8 J7 \axis([0,1,-30,5]);3 E* n( m1 I" c
    set(gca,'XTickMode','manual','XTick',[0 0.2 0.3 1]);! T* c7 t3 F  S: _/ Z2 L
    set(gca,'YTickMode','manual','YTick',[-60 -as -rp 0]);grid& M1 y6 R) t& N" }: E
    subplot(2,2,3);plot(w/pi,pha/pi);title('相频特性');
    , a; p! k9 h, k5 p; Oxlabel('w(/pi)');ylabel('pha(/pi)');! X( ~5 z4 N- L# Q+ F
    axis([0,1,-1,1]);9 w4 q% K. j# M- d; f7 b' `
    set(gca,'XTickMode','manual','XTick',[0 0.2 0.3 1]);grid
    # E/ x- y9 `6 }( n" Z% K! u( l* }subplot(2,2,4);plot(w/pi,grd);title('群延时');  l9 \" g7 [( n& w
    xlabel('w(/pi)');ylabel('Sample');/ ~4 N6 {- E9 b3 E, S  F  n
    axis([0,1,0,15]);
    - b! _* w8 ]4 U3 Fset(gca,'XTickMode','manual','XTick',[0 0.2 0.3 1]);grid
    : Z  S$ q4 [8 Q- \%b,a是滤波器的系数,系统函数根据b和a写出。
    % }) E& Z4 u: p1 y; r%若程序运行有问题可以Hi我
    ' Q& A# \% r. R7 N7 e; _

    该用户从未签到

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

    该用户从未签到

    4#
    发表于 2022-8-11 16:04 | 只看该作者
    楼上说的很不错的
    2 s4 V# d/ h) O  O
    您需要登录后才可以回帖 登录 | 注册

    本版积分规则

    关闭

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

    EDA365公众号

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

    GMT+8, 2025-11-3 03:11 , Processed in 0.125000 second(s), 23 queries , Gzip On.

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

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

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