信号处理与系统析分课程设计内容.docVIP

信号处理与系统分析课程设计内容 ※三个系统的系统函数分别为h1(s)=5/[(5s+1)]; h2(s)=(4s+5)/[s(5s+1)];h3(s)= (4s+5)/[5s2+5s+5)]用simulink来仿真三个系统的阶跃响应，并分析系统的稳定性。 单边带信号的产生 （1）用simulink来仿真单边带信号的产生过程，假设载波频率为wc=100。X(t)=sin(pi*t/4)其原理如图： x(t) sinwct y(t) 相移π/2 - y(t)= sin(pi*t/4)sinwct- cos(pi*t/4)coswct （2）将信号的频谱分量均相移，在求和点处，取两个信号相加，得到下边带信号；取两个信号相减，得到上边带信号。试分析信号的频谱。 频谱分析 设一序列中含有两种频率成分，f1=2hz，f2=2.05hz，采样频率为fs=10Hz，即x(n)=sin(2πf1n/fs)+ sin(2πf2n/fs)，分析其频谱。 ※IIR滤波器的设计 脉冲响应不变法、双线性变换法设计IIR数字巴特沃斯低通数字滤波器。例如：通带截止频率0.2*pi,阻带截止频率0.3*pi,通带波动1dB；在频率0.3π到π之间的阻带衰减大于10dB。 五、设计一个高通数字滤波器，要求通带截止频率ωp=0.8πrad,通带最大衰减αp=3db,阻带截止频率ωs=0.5πrad，阻带最小衰减αs=18db。 六、※FIR滤波器的设计 用窗函数法设计FIR数字低通滤波器。假设N=11，wc=0.2*πrad。（P210） 七、使用MATLAB进行功率谱的噪声分析（（功率谱：信号频谱幅度|Xk|2/N），注意a,n的取值情况） 十、※设信号： x(t)=s(t)+u(t),其中u(t)是干扰信号，s(t)+u(t)的频谱不混叠，其幅度谱如图所示，要求设计一数字滤波器，将干扰滤除掉，指标是允许|s(f)|在0≤f≤15Khz频率范围中幅度失真为±2%（δ1=0.02）；f20Khz,衰减大于40db(δ1=0.01）,希望分别用FIR 和IIR两种滤波器进行滤除干扰，最后进行比较。 |X(f)| |s(f)| |v(f)| 0 20 35 f/khz 十一、ap=3db,as=35db,fp=50hz,fs=1000hz,试利用极、零点对消法设计一简单整系数低通滤波器，求系统函数，并画出幅频特性曲线。 1.参考资料 （1）数字信号处理??丁玉美??西安电子科技大学出版社 （2）应用MATLAB实现信号分析和处理?? 张明照? ?科学出版社 （3）数字信号处理及MATLAB实现?? 余成波 ?清华大学出版社 （4）MATLAB7.0在数字信号处理中的应用? ?罗军辉??机械工业出版社 （5）MATLAB信号处理? 刘波?? 电子工业出版社 （6）Matlab信号处理与应用?? 董长虹?? 国防工业出版社 （7）数字信号处理原理及其MATLAB实现? 从玉良 ??电子工业出版社 （8）数字信号处理基础及MATLAB实现 ?周辉 ?中国林业出版社 （10）MATLAB6.x信号处理邹鲲等清华大学出版社出版社 2.3 设计思想 （1）利用单片机产生方波、正弦波、三角波和锯齿波等信号波形，信号的频率和幅度可变。 （2）将一个周期的信号分离成256个点（按X轴等分），每两点之间的时间间隔为T，用单片机的定时器产生，其表示式为：T=T/256。 如果单片机的晶振为12MHz，采用定时器方式0，则定时器的初值为： X=213—T/Tmec 定时时间常数为： TL =（8192—T）/MOD(2.2) TH=(8192T)/256 (2.3) MOD32表示除32取余数 （3）正弦波的模拟信号是D/A转换器的模拟量输出，其计算公式为： Y=（A/2sint）+A/2(其中A=VREF)t=NT (N=1~256) (2.5) 那么对应着存放在计算机里的这一点的数据为：