信号处理课程设计实验报告.docVIP

信号处理课程设计实验报告 专业：电子信息工程 班级： XXX 姓名： XXX 学号：XXX XX年XX月XX日 一、课程设计目的 信号处理课程设计是学习完信号与系统、数字信号处理及DSP课程后全面复习课程所学知识分析和解决问题的能力 二、课程设计的内容和要求 1 要求： 课程设计 ⑴ 课程设计题目和题目设计要求； ⑵ 设计思想和系统功能结构及功能说明； ⑶ 设计中关键部分的详细描述和介绍，采用流程图描述关键模块的设计思路； ⑷ 总结，包括设计过程中遇到的问题和解决方法，心得体会等； ⑸ 参考文献； ⑹ 程序源代码清单。 2 课程设计题目和内容： 实验球磨机描述 图2 磨机滚筒尺寸 实验用球磨机内壁上安装有锰钢衬板。实验过程中使用Φ30mm、Φ20mm和Φ15mm三种锰钢钢球。磨机内所装钢球的最大重量为80kg，矿物处理量为10kg/h。实验中处理的矿石为经过破碎的铜矿石，颗粒尺寸小于6mm。 实验数据来源于XMQL 420×450格子型球磨机数据采集系统。其中振动加速度传感器安装在球磨机中部筒体外表面上，采样频率为51,200 Hz，声音传感器安装在距球磨机入口1/3处。该实验球磨机最大钢球装载量为80Kg，磨粉能力为10Kg/小时，转速为53转/分钟。研磨实验方案是在球磨机筒体内添加不同配比(直径30、20和15mm)的钢球、铜矿石和水，混合均匀后研磨运行1分钟后清洗球磨机。 图2.4 实验球磨机信号采集 （1）振动加速度信号读取：txt文件读取并保存为mat文件 实验程序：file_path=cd(E:\yzm );%设置为当前目录 D=dir(*.txt); m=length(D); for i=1:m i if D(i).isdir==0 filename=D(i).name; fid = fopen([filename]); while ~feof(fid) for i=1:17 tline = fgetl(fid); end a = fscanf(fid,%f %f,[2 inf]); end data=a; eval([save ZD filename(1:8) data]) status = fclose(fid); end end plot(data(1:end,1),data(1:end,2)) 实验结果： Y=wavread(E:\yzm\NVOICE16.wav); [i,j]=size(Y); m=length(Y); FS=8000; t=(0:(m-1))*(1/FS); t=t; [Y,FS,NBITS]=wavread(E:\yzm\NVOICE16.wav); plot(t,Y); 实验结果： 振 声 数 据 图 像 ⑵ 任务2：振动和振声数据预处理：降低采样频率，提出粗大误差和噪声 实验程序：%时域信号显示与误差剔除 function tsignal(file) %clear all; %clc; load(file) m=length(file); filename=file(1:m-4); test=eval(filename); %y=test(1024001:2:1536000,2);%20-30秒的采样数据 y=test(1:2:424372,1);%20-30秒的采样数据 ymean=mean(y) ystd=std(y); outliers=abs(y-ymean)3*ystd; %剔除离群点数据 nout=sum(outliers) y(outliers==1)=[]; yn=length(y) ymin=min(y) ymax=max(y) figure plot(y) % axis([0 60 -0.2 0.2]) grid title(NVOICE16.mat T signal) xlabel(采样点数) ylabel(振动幅值) 震动剔除实验结果： 振声剔除实验结果： （3）任务3：振动和振声数据频谱，请画出其傅立叶变换图； 实验程序：振动信号：file_path=D:\work\ballmill\MJ_DATA\ZD2_data\; D=dir(file_path); mD=length(D); PP=[];FF=[];SS=[]; for i=1:mD i if D(i).