小波分析—时间序列的多时间尺度分析

sam_doggy

最近刚开始学习小波分析,看到很好的一个博客,博主申明可以随意转载了,就转来分享.
原链接如下 http://blog.sciencenet.cn/blog-1148346-794768.html 不全篇转载了

sam_doggy

下面是matlab新手对作者源程序的一些注解一边看一边更新:
在LI_plot中,作者使用了bar(条状图) 函数和regress(线性回归)函数以及自编函数PlotAddText,
因为作者作图的需要,对bar函数的参数进行了调控,set 命令是调节坐标轴等参数,[0, 255, 255]./255这个颜色是天青色

nunu18

楼主这个程序是matlab的是吧？

sam_doggy

nunu18 发表于 2018-3-3 12:19
楼主这个程序是matlab的是吧？

是matlab的，我还在慢慢看，很多不懂的，很小白

sam_doggy

%z作者的Li_plot函数代码
% Created by LI Xu
% Version 1.0
% Feburary 27, 2016

% Description:
% Plot precipitation

% If you have any question about this code,
% please do not hesitate to contact me via E-mail:
% jeremy456@163.com

% Blog:
% http://blog.sciencenet.cn/u/lixujeremy
% http://lixuworld.blogspot.com/

function LI_plot(years, prec)
    clf;
    handle=bar(years, prec);
    set(gca, 'xlim', [min(years)-1, max(years)+1]);
    set(gca, 'fontname', 'Arial', 'fontsize', 15);
    set(handle, 'facecolor', [0,1,1]);
    set(gcf, 'color', 'w');
    set(gca, 'ylim', [600, 1800]);
    hold on;
    % Trend Line
    X=[ones(size(years)), years];
    [b, ~, ~, ~, stats]=regress(prec, X);
   
    b1=sprintf('%.3f', b(1));
    b2=sprintf('%.3f', b(2));
    str1=['{\ity}=', b2, '{\itx}+', b1];
    pval=stats(3);
    pval=sprintf('%.3f', pval);
    str2=['{\itp}=', pval];
    instr={str1, str2};
    instratts.font='Arial';
    instratts.xscale=0.1;
    instratts.yscale=0.85;
    instratts.fontsize=30;
    instratts.interval=0.1;
    PlotAddText(gca, instr, instratts);
    xlabel('{\itYear/a}', 'fontname', 'Times New Roman', 'fontsize', 20);
    ylabel('{\itAnnual Precipitation/mm}', 'fontname', 'Times New Roman', 'fontsize', 20);
    xlim=get(gca, 'xlim');
    yy=b(1)+b(2)*xlim;
    linehand=plot(xlim, yy, 'k-');
    set(linehand, 'linewidth', 2);
    hold off;
end

sam_doggy

读这个代码的%trend line行以前，我对于图像的细微调整大多采用xlim,axis,xlabel这几个简单的函数，以及figure界面的鼠标单击，调整对于图像句柄的调用几乎不会。事实上，简单的界面和只需要调整一次以及自娱自乐的那些东西直接在figure界面里调整倒是也够了。
今天又看了一些相关的课件，认为主要需要掌握的是get，set两个函数就行。get是查看参数，set是设定。
关于这方面的理解我把常见的一些图和我之前做的图采用get多看了一些参数。
如果有和我一样的小白建议多get一些参数看看。
下面是我认为比较简单的课件。

scirencc

这个代码能跑吗

最初的向往

很有聊，谢谢楼主分析

青苹果派派

谢谢分享！

sam_doggy

scirencc 发表于 2018-3-9 19:18
这个代码能跑吗

这是函数，main我放网盘了啊