统计信号分析与处理报告讲义.docxVIP

中国地质大学（武汉）统计信号分析处理实验报告小组成员: 魏彦斌 马全林 陈飞班号: 075132 _院系：__机电学院专业：_通信工程指导教师：_侯强老师一实验内容实验一、地震时间间隔的密度估计；实验二、地震震级与频度回归分析；实验三、地震空间分布聚类分析；实验四、地震优势深度聚类分析；二.实验要求及结果。 实验一、地震时间间隔的密度估计；读入数据后，把时间列（第一二列）转换成数据格式，然后分：3级以下，3到4级，4到5级，5到6级，6级以上等6个部分分别计算地震之间的时间差t，这个t就是地震时间间隔，而且是一个随机变量，对t这个随机变量进行密度估计。代码：%% clear all;close all; clc;filename = 中国地震台网(CSN)地震目录(1970-01-01至2015-09-31).xls;sheet = 1;xRange = A5:A8462;% xRange = A3:A8462;x2Range = B5:B8462;% x2Range = B3:B8462;yRange = H3:H8462;% [~,x]= xlsread(filename, sheet, xRange);% [~,x2] = xlsread(filename, sheet, x2Range);[~,x]= xlsread(filename, sheet, xRange);[ttt,x2] = xlsread(filename, sheet, x2Range);ml= xlsread(filename, sheet, yRange); %读取数据%% X = x(~isnan(ml));X2 = x2(~isnan(ml));n = length(X2); %去掉无数据的日期和时间 for i = 1:n %将日期时间转化为数值形式 Xyy(i) = str2double(X{i}(1:4)); %年 Xmm(i) = str2double(X{i}(6:7)); %月 Xdd(i) = str2double(X{i}(9:10)); %日 XHH(i) = str2double(X2{i}(1:2)); %时 XMM(i) = str2double(X2{i}(4:5)); %分 XSS(i) = str2double(X2{i}(7:8)); %秒 endxx = datenum(Xyy,Xmm,Xdd,XHH,XMM,XSS); %将时间转化为数值形式ML = ml(~isnan(ml) ); %去掉无数据项a=1; b=1;c=1;d=1; e=1;for i=1:n if ML(i)=3.0 t_3(a)=xx(i); a=a+1; elseif ML(i)3.0ML(i)=4.0 t_34(b)=xx(i); b=b+1; elseif ML(i)4.0ML(i)=5.0 t_45(c)=xx(i); c=c+1; elseif ML(i)5.0ML(i)=6.0 t_56(d)=xx(i); d=d+1; else t_6(e)=xx(i); e=e+1; end;end;%求个部分时间差for i=1:(length(t_3)-1) tt_3(i)=t_3(i)-t_3(i+1);endfor i=1:(length(t_34)-1) tt_34(i)=t_34(i)-t_34(i+1);endfor i=1:(length(t_45)-1) tt_45(i)=t_45(i)-t_45(i+1); endfor i=1:(length(t_56)-1) tt_56(i)=t_56(i)-t_56(i+1);endfor i=1:(length(t_6)-1) tt_6(i)=t_6(i)-t_6(i+1);endx=linspace(min(tt_3)-1,max(tt_3),1200);p=Parzen(tt_3,x,15,[]);plot(x,p);grid on;figure(2);x=linspace(min(tt_34),max(tt_34),50)