文档详情

DSP数字信号处理课程设计单片机红绿灯课程设计.doc

发布:2018-10-07约5.23千字共5页下载文档
文本预览下载声明
DSP数字信号处理课程设计单片机红绿灯-课程设计 DSP数字信号处理课程设计单片机红绿灯|DSP课程设计及总结|数字信号处理课程设计课程设计|计算机课程设计|精品课程网站设计 一、??????? 设计目的 ? 熟悉filter的设计及作用 傅氏变换后的应用情况 ? 二、??????? 内容 ? 1、??????????? 产生100、200、300Hz的正弦信号并叠加,前后情况对比 2、??????????? 设计在150—250之间的带通 3、??????????? 合成信号在通过filter前后的情况对比 ? 三、??????? 设计内容 ? (一)、设计思路 1、产生100Hz、200Hz、300Hz的正弦信号 y1=A*sin(2*pi*f1*t); y2=A*sin(2*pi*f2*t); y3=A*sin(2*pi*f3*t); ? 2、将三个正弦信号叠加 ? y123=y1+y2+y3; ? 3、通过fft函数进行傅立叶变换,转换成频域函数 ? y1_fft=fft(y1); y2_fft=fft(y2); y3_fft=fft(y3); y123_fft=fft(y123); ? 4、带通滤波器的设计 ? f_N=3000; %取样频率 f_P=[150,250]; ????????????????????????????%通带截止频率 f_S=[100,300];? ???????????????????????????%通带起始频率 R_p=3;???????????????????? ??????????????%通带内波动 R_s=25;???????????????????? ?????????????%阻带内最小衰减 Ws=f_S/(f_N/2);????????? ??????????????%阻带起始角频率 Wp=f_P/(f_N/2);?????????? ?????????????%阻带截止角频率 [n,Wn]=buttord(Wp,Ws,R_p,R_s);????? %利用巴特沃斯滤波器设计函数求出滤波器阶数和截至频率 [b,a]=butter(n,Wn);?? ?????????????????%利用butter函数计算滤波器系数 freqz(b,a,1000,3000);????????????????? %计算幅频响应 ? 5、合并信号经过带通滤波器后的时域和频域 ? y_filter=filter(b,a,y123);??????? ?%计算合并信号通过滤波器后的时域信号 y_filter_fft=fft(y_filter);?? ?????%计算合并信号通过滤波器后的频域信号 ? ?(二)、程序代码分析 clear; ? %定义参数及函数方程 fs=1800; t=1/fs:1/fs:0.05 f1=100; f2=200; ???????????????????? f3=300; ???????????????????? A=1; ???????????????????? y1=A*sin(2*pi*f1*t); ???????????????????? y2=A*sin(2*pi*f2*t); ???????????????? ????y3=A*sin(2*pi*f3*t); ???????????????????? y123=y1+y2+y3; ???????????????????? ???????????????????? %时域图 ???????????????????? figure(1); ???????????????????? subplot(4,2,1); ???????????????????? plot(t,y1,b); ???????????????????? title(y1的时域波形图); ???????????????????? xlabel(时间t); ???????????????????? ylabel(幅度A); ???????????????????? subplot(4,2,2); ???????????????????? plot(t,y2,b); ???????????????????? title(y2的时域波形图); ???????????????????? xlabel(时间t); ???????????????????? ylabel(幅度A); ???????????????????? subplot(4,2,3); ???????????????????? plot(t,y3,b); ???????????????????? title(y3的时域波形图); ???????????????????? xlabel(
显示全部
相似文档