用matlab进行时频信号分析器设计

lahhse

产生如下4种信号：频率为5HZ，采样频率为100HZ的正弦信号x；
将信号x混有随机噪声后的带噪信号；
% J; _2 N) n/ t* Sx的频谱信号；
, k& O$ Q8 W9 g' K' ^  V  I带噪信号的频谱信号。
/ [$ N$ k1 k% G/ V2 P/ J2 U设计良好的人机界面，每个按钮对应一种波形
; \  f* f' @% I/ ^  q
5 L9 p& e( P% h

Blah

wp=0.2*pi;ws=0.3*pi;%你没有给出指标，可以自己改
Fs=20000;T=1/Fs;%采样频率可以自己定
OmegaP=(2/T)*tan(wp/2);
$ b9 E( t0 m; W$ J/ WOmegaS=(2/T)*tan(ws/2);
rp=1;rs=15;
* Q  p- c: E3 J4 @ripple=10^(-rp/20);attn=10^(-rs/20);
[n,wn]=buttord(OmegaP,OmegaS,rp,rs,'s');
) p! p4 x# k* e1 y. N' w' M% Z9 C6 y[z,p,k]=Buttap(n);
, P2 W% E; J. {+ Q6 l& R5 `[b,a]=zp2tf(z,p,k);
[bt,at]=lp2lp(b,a,wn);
[b,a]=bilinear(bt,at,Fs);
! i3 h3 j% r, A' u, y* M[db,mag,pha,grd,w]=freqz_m(b,a);
8 A( N/ j/ Y) F: K%
%下面绘出各条曲线
subplot(2,2,1);plot(w/pi,mag);title('幅频特性');
xlabel('w(/pi)');ylabel('|H(jw)|');
8 E) q( C# O! g8 J* jaxis([0,1,0,1.1]);
& t; w; @* a/ |$ Wset(gca,'XTickMode','manual','XTick',[0 0.2 0.3 1]);
* _1 O7 {1 q7 o3 G% A7 vset(gca,'YTickMode','manual','YTick',[0 attn ripple 1]);grid
subplot(2,2,2);plot(w/pi,db);title('幅频特性(db)');
xlabel('w(/pi)');ylabel('dB');
( j3 C* }( `# X$ R8 @9 caxis([0,1,-30,5]);
9 q1 y$ o  q/ m( p* D; Q0 ?set(gca,'XTickMode','manual','XTick',[0 0.2 0.3 1]);
set(gca,'YTickMode','manual','YTick',[-60 -as -rp 0]);grid
subplot(2,2,3);plot(w/pi,pha/pi);title('相频特性');
  X$ Y, ], l$ T1 w/ d& \7 wxlabel('w(/pi)');ylabel('pha(/pi)');
$ P) o6 q2 B" H& waxis([0,1,-1,1]);
1 v4 Y4 C5 n! j6 l2 U2 Q1 b" N5 z) Rset(gca,'XTickMode','manual','XTick',[0 0.2 0.3 1]);grid
subplot(2,2,4);plot(w/pi,grd);title('群延时');
: l3 R9 h8 A0 }. d  b- v+ _xlabel('w(/pi)');ylabel('Sample');
axis([0,1,0,15]);
) O6 [) `( p- m1 Kset(gca,'XTickMode','manual','XTick',[0 0.2 0.3 1]);grid
%b,a是滤波器的系数，系统函数根据b和a写出。
%若程序运行有问题可以Hi我

land

bootloader是dsp芯片上固化的一段引导程序，dsp程序在ram上掉电易丢，所以启动时运行bootloader将辅存中的程序调入ram中
' p5 o7 A' C, D2 g0 G2 K至于dsp的框架结构在ti公司有（有点复杂），开发时通常是先在开发板上仿真然后实践，由于dsp一般用于数字信号处理，离不开matlab软件连调。说起来有点复杂，要熟悉dsp开发流程序最简单的方法就是自己动手实践一下，如做个滤波器，频谱分析器等就可以深刻了解了。

理论的

楼上说的很不错的