matlab与信号处理实践,信号处理与MATLAB实践 实验报告.docx

北京工业大学

信号处理与MATLAB实践

信号处理与MATLAB实践

实验报告

班级：140200

学生姓名：陈艺飞

学号指导教师：孙中华、张延华

完成时间：2016.12

信号处理与MATLAB实践实验报告 陈艺PAGE \* MERGEFORMAT8

目录

TOC \o "1-3" \h \z \u HYPERLINK \l "_Toc469594113" 一、任选题 PAGEREF _Toc469594113 \h 3

HYPERLINK \l "_Toc469594114" 1、作业要求 PAGEREF _Toc469594114 \h 3

HYPERLINK \l "_Toc469594115" 2、流程图 PAGEREF _Toc469594115 \h 3

HYPERLINK \l "_Toc469594116" 3、解题步骤 PAGEREF _Toc469594116 \h 4

HYPERLINK \l "_Toc469594117" 4、运行环境 matlab PAGEREF _Toc469594117 \h 8

HYPERLINK \l "_Toc469594118" 5、运行代码汇总 PAGEREF _Toc469594118 \h 8

HYPERLINK \l "_Toc469594119" 二、双号必选题： PAGEREF _Toc469594119 \h 10

HYPERLINK \l "_Toc469594120" 1、题目1 PAGEREF _Toc469594120 \h 10

HYPERLINK \l "_Toc469594121" 2、题目2 PAGEREF _Toc469594121 \h 15

HYPERLINK \l "_Toc469594122" 三、心得与体会 PAGEREF _Toc469594122 \h 17

一、任选题

1、作业要求

作业-1太阳黑子活动周期的分析

太阳黑子的活动是周期的，大约每11年达到一个爆发高峰。试证明这一点。

设计提示：

①首先下载太阳黑子的数据，可以从比利时皇家天文台(Royal Observatory of Belgium)的太阳影响数据分析中心(Solar Influences Data Analysis Center-SIDC)下载，网址是： HYPERLINK "http://sidc.oma.be/index.php3" http://sidc.oma.be/index.php3

下载数据的时间段可以从1741年1月一直到当前。

②以横坐标表示年份，纵坐标表示黑子出现的数量，绘制Wolfer图。

③、应用FFT技术分析Wolfer数，在复平面上可以直接绘制出由Y给出的傅立叶系数的分布图。

④、绘制周期图(Y的模的平方被定义为功率，功率与频率的关系曲线则被定义为周期图)。

注：要获得太阳黑子的活动周期，你可以使用月度数据。该数据可以通过点击太阳影响数据分析中心网页左边的导航条Sunspots→download of data→monthly and monthly smoothed sunspot number下载。下载的数据文件的第1列是年和月，第3列是该月太阳黑子的平均数，第4列是平滑后该月太阳黑子的平均数。为了分析太阳黑子的活动规律，即可以取任一时间段的数据，也可以取下载的全部数据。为分析方便，建议用Excel将数据文件读进表格并且存到MATLAB工作路径。

Wofler数时域图FFT

Wofler数

时域图

FFT变换

功率谱密度

(功率-频率)

功率谱密度

(功率-周期)

标出周期

标出周期

3、解题步骤

①下载自 HYPERLINK "http://sidc.oma.be/silso/datafiles" http://sidc.oma.be/silso/datafiles 页面内的Total sunspot number——Yearly mean total sunspot number [1700 - now]，CSV格式，并将其根据Matlab语法要求将其重命名为csv。

②保存太阳黑子数据(csv)至MATLAB工作路径

③装载太阳黑子的数据

load csv; %读取太阳黑子数据

④以横坐标表示年份，纵坐标表示黑子出现的数量，绘制Wolfer图。

year=:,1); %读取年份信息

wolfer=:,2);