Matlab程序小波分析时间序列数据

  • q2_431020
    了解作者
  • 2MB
    文件大小
  • rar
    文件格式
  • 0
    收藏次数
  • VIP专享
    资源类型
  • 0
    下载次数
  • 2022-05-28 08:46
    上传日期
以美国某气象站1894~2010年连续的年降水量为例,试应用小波分析,完成如下任务:①小波变换系数;②绘制小波系数实部等值线图;③绘制小波系数模和模方等值线图;④绘制小波方差图;以及⑤绘制不同时间尺度的小波实部过程线。所谓年降水量时间序列的多时间尺度是指:年降水量在演化过程中,并不存在真正意义上的变化周期,而是其变化周期随着研究尺度的不同而发生相应的变化,这种变化一般表现为小时间尺度的变化周期往往嵌套在大尺度的变化周期之中。也就是说,年降水量变化在时间域中存在多层次的时间尺度结构和局部变化特征。Details: http://blog.sciencenet.cn/blog-1148346-794768.html
WavletAnalysis.rar
  • WavletAnalysis
  • main.asv
    2.5KB
  • ePrec.mat
    734B
  • 6_var.jpg
    100.6KB
  • main.m
    3.2KB
  • Prec.mat
    718B
  • entendvalue.m
    532B
  • pictures
  • 7_55a_coefficient.jpg
    112.6KB
  • 1.png
    143.9KB
  • 2.png
    200.1KB
  • 0_trend.png
    41KB
  • 8_30a_coefficient.jpg
    128.4KB
  • 4_mod_contourf.jpg
    447KB
  • 3_real_contourf.jpg
    501.2KB
  • 5_mod_square_contourf.jpg
    446.2KB
  • LI_Var.m
    369B
  • Year.mat
    363B
  • cePrec.mat
    123.4KB
  • LI_period_plot.asv
    1KB
  • LI_plot.m
    889B
  • LI_period_plot.m
    630B
  • LI_contourf.m
    1.1KB
内容介绍
% Created by LI Xu % 15 May, 2014 % Version 1.0 % Main program for Wavlet Analysis clear; clc; load Prec load ePrec load cePrec load Year clear wname % Location of the Original Arrary in % the Extended Arrary [st, ed]=entendvalue(Prec, ePrec); % Retrieve the Coefs of Original Arrary Orig_coefs=coefs(:, st:ed); clear coefs ePrec st ed % Retrieve the real part of Orig_coefs realpart=real(Orig_coefs); % Drawing the trend of the Precipitation % LI_plot(Year, Prec, 600); % xlabel('Year/a', 'fontangle', 'italic', 'fontsize', 20, 'fontweight', 'b'); % ylabel('Annual Precipitation/mm', 'fontangle', 'italic', 'fontsize', 20, 'fontweight', 'b'); % legend('Observation', 'Trend Line'); % export_fig 0_trend.png % % Drawing the contourf map of the real part % text_label=-500:100:600; % text_x=1904:10:2004; % % LI_contourf(Year, realpart, scales, text_label, text_x, 14); % xlabel('Year/a', 'fontangle', 'italic', 'fontsize', 20, 'fontweight', 'b'); % ylabel('Period Scale/a', 'fontangle', 'italic', 'fontsize', 20, 'fontweight', 'b'); % export_fig 3_real_contourf.png % Drawing the modulus map of the Wavelet Coefficients % mod_coef=abs(Orig_coefs); % text_label=0:50:620; % text_x=1904:10:2004; % % LI_contourf(Year, mod_coef, scales, text_label, text_x, 14); % xlabel('Year/a', 'fontangle', 'italic', 'fontsize', 20, 'fontweight', 'b'); % ylabel('Modulus/dimensionless', 'fontangle', 'italic', 'fontsize', 20, 'fontweight', 'b'); % export_fig 4_mod_contourf.png % Drawing the modulus square of the Wavelet Coefficients % mod_square=(abs(Orig_coefs)).^2; % text_label=0:15000:380000; % text_x=1904:10:2004; % % LI_contourf(Year, mod_square, scales, text_label, text_x, 14); % xlabel('Year/a', 'fontangle', 'italic', 'fontsize', 20, 'fontweight', 'b'); % ylabel('Modulus Square/dimensionless', 'fontangle', 'italic', 'fontsize', 20, 'fontweight', 'b'); % export_fig 5_mod_square_contourf.png % Plot Wavelet Variance % http://wenku.baidu.com/view/bec1ebc90c22590102029dce.html Wavlet_var=sum(abs(Orig_coefs).^2, 2); Wavlet_var=Wavlet_var./length(Wavlet_var); LI_Var(scales, Wavlet_var); set(gca, 'xgrid', 'on'); xlabel('Period Scale/a', 'fontangle', 'italic', 'fontsize', 20, 'fontweight', 'b'); ylabel('Wavelet Variance/dimensionless', 'fontangle', 'italic', 'fontsize', 20, 'fontweight', 'b'); export_fig 6_var.png % Plot the wavelet coefficient for 55a-period scale % a55_coef=realpart(55, :); % LI_period_plot(Year, a55_coef); % legend('55a-Period Scale', 'location', 'north'); % % xlabel('Year/a', 'fontangle', 'italic', 'fontsize', 20, 'fontweight', 'b'); % ylabel('Wavelet Coefficient/dimensionless', 'fontangle', 'italic', 'fontsize', 20, 'fontweight', 'b'); % export_fig 7_55a_coefficient.png % Plot the wavelet coefficient for 30a-period scale % a30_coef=realpart(30, :); % LI_period_plot(Year, a30_coef); % legend('30a-Period Scale', 'location', 'north'); % % xlabel('Year/a', 'fontangle', 'italic', 'fontsize', 20, 'fontweight', 'b'); % ylabel('Wavelet Coefficient/dimensionless', 'fontangle', 'italic', 'fontsize', 20, 'fontweight', 'b'); % export_fig 8_30a_coefficient.png % Close All Figures close all;
评论
    相关推荐
    • swim.html:http
      swim.html 来源
    • Html2Xaml:http的分叉
      Html2Xaml 这不是将 html 转换为 xaml 的好方法,我认为最好将其转换为列表,并通过 ListView 显示内容
    • html-redirect:生成 HTML 重定向页面
      htmlRedirect ( 'http://example.com' ) . pipe ( fs . createWriteStream ( 'example1.html' ) ) ; 示例1.html: <!DOCTYPE html > < html > < head > < title > (redirect) </ ...
    • Lingulo-Responsive-Tutorial:httpHTML5自适应网站教程
      语言响应式教程 HTML5自适应网站教程, 此HTML5网站是Christoph Anastasiades创建HTML5 / CSS3教程的一部分( )随意使用用于您的个人或商业项目...为了最大程度地减少HTTP请求的数量,请将所有CSS文件合并为一个文件
    • Cryptoquips:这个简单的应用程序从http解析html
      密码学 这个简单的应用程序解析html,并在卡视图中显示列出的Cryptoquip图像。
    • Lightweight C-HTTP & HTML Wrapper:轻量级的C-HTTPHTML包装器-开源
      它是一个C-HTTPHTML包装器,用于自动处理Web数据。 一些示例是:-接收网站以检查其是否正常运行-从例如下载您的个人数据。 网上银行,帐单...-处理网站中包含的数据,例如天气数据-接收大量数据,尤其是对网站...
    • slackify-html
      slackify-htmlhtml转换为 用法示例: var slackify = require('slackify-html'); var text = slackify('this <a>link</a> is <b>important</b>');...// text variable contains 'this <http> is *important*'
    • HTML5-Examples:httpHTML 5示例-html source website
      HTML5-示例 HTML 5示例。 版权所有2012:SKJ Technologies Private Ltd 这些示例已根据“知识共享署名-相同的ShareAlike 3.0未移植”(CC BY-SA 3.0)许可证进行共享。 有关CC BY-SA 3.0许可证的更多详细信息,请...
    • html定时抓取程序
      利用http协议,定时抓取指定网页内容,保存为.htm文件的程序;利用http协议,定时抓取指定网页内容,保存为.htm文件的程序;
    • html日期控件
      http://hiyu2218.bokele.com/?CH=6441