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

用matlab进行时频信号分析器设计 悬赏分：150 - 解决时间：2010-7-1 22:50 产生如下4种信号：频率为5HZ，采样频率为100HZ的正弦信号x； 将信号x混有随机噪声后的带噪信号； x的频谱信号； 带噪信号的频谱信号。 设计良好的人机界面，每个按钮对应一种波形 高分求助，任务紧急，谢谢各位。最好是成型的论文，如果没有的话有程序也行！ 问题补充： 不好意思忘记说了 我要用MATLAB模拟，不用单片机 提问者： 迪卡朋鑫 - 五级 最佳答案 wp=0.2*pi;ws=0.3*pi;%你没有给出指标，可以自己改 Fs=20000;T=1/Fs;%采样频率可以自己定 OmegaP=(2/T)*tan(wp/2); OmegaS=(2/T)*tan(ws/2); rp=1;rs=15; ripple=10^(-rp/20);attn=10^(-rs/20); [n,wn]=buttord(OmegaP,OmegaS,rp,rs,s); [z,p,k]=Buttap(n); [b,a]=zp2tf(z,p,k); [bt,at]=lp2lp(b,a,wn); [b,a]=bilinear(bt,at,Fs); [db,mag,pha,grd,w]=freqz_m(b,a); % %下面绘出各条曲线 subplot(2,2,1);plot(w/pi,mag);title(幅频特性); xlabel(w(/pi));ylabel(|H(jw)|); axis([0,1,0,1.1]); set(gca,XTickMode,manual,XTick,[0 0.2 0.3 1]); set(gca,YTickMode,manual,YTick,[0 attn ripple 1]);grid subplot(2,2,2);plot(w/pi,db);title(幅频特性(db)); xlabel(w(/pi));ylabel(dB); axis([0,1,-30,5]); 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(相频特性); xlabel(w(/pi));ylabel(pha(/pi)); axis([0,1,-1,1]); set(gca,XTickMode,manual,XTick,[0 0.2 0.3 1]);grid subplot(2,2,4);plot(w/pi,grd);title(群延时); xlabel(w(/pi));ylabel(Sample); axis([0,1,0,15]); set(gca,XTickMode,manual,XTick,[0 0.2 0.3 1]);grid %b,a是滤波器的系数，系统函数根据b和a写出。 %若程序运行有问题可以Hi我 Note:Matlab中DFT在连续信号谱分析中的应用 2008-06-04 21:08 ?? 本来是想用Matlab做通信课的模拟调制分析的，结果弄一个时频变换就吭哧了两三天时间，把原来的老底信号系统、DSP的书又翻出来看，总算稀里糊涂画出来了..... 所谓信号的谱分析，就是时频域转换，变成数学知识，就是对信号进行傅立叶变换（FT），连续信号与系统的傅立叶变换不便于用计算机进行计算。而在Matlab中的绘图操作，是通过将一些离散的点连接而成。求连续信号的频谱，应用的是离散傅立叶变换（DFT）。 ??? 下面这段程序是绘制简单的Sa函数时域和频域图像。??? Sa(t)=sin(200t)/200t ????? 对于这个函数我们要说明一点。由傅立叶变换理论知道，若信号持续时间有限长，则其频谱无限宽；若信号频谱有限宽，则其持续时间无限长。因此要对其进行谱分析就必须在时域采样无穷多个点，这又显然不满足DFT变换的要求。从工程实际角度来看，滤除幅度很小的高频成分和截去幅度很小的部分时间信号是允许的，因此这里我们截取（-5≤t≤5）时间段的Sa（t）。下面是程序代码： Sa.m Clear t0=10;??????????????????????? %定义时间长度10s，即（-5≤t≤5）?? ts=0.001; fs=1/ts;??? %定义采样周期、采样频率 df=0.5;???????????????????? %定义频谱分辨率(见说明1)。 t=[-t0/2:ts:t0/2];??