介绍威斯康星大学气候学副教授Dan Vimont自编的MATLAB函数库

hillside

Climatological Analysis(《气候学分析》课程)

主讲：Dan Vimont

email:  dvimont@wisc.edu

Course Website（课程网址）:
http://www.aos.wisc.edu/~dvimont/aos575

Meeting Time and Location:
Class meets MWF at 9:55-10:45 a.m.,in AO&SS 811
Office Hours:
By appointment

Course Syllabus:  (PDF)
Course Schedule:  (PDF)
Paper Discussions:  (PDF)
Final Project Description:  (PDF)
Lecture Notes
We will be using lecture notes provided courtesy of Profs. Dennis Hartmannand John M. Wallace (Department of Atmospheric Sciences, Unversity of Washington).  Allchapters are in PDF format.

Chapter1:  Review of Basic Statistics
BivariateNormal Distribution Notes
T-testnotes
Chapter2:  Composite Analysis
Chapter3:  Linear Regression
Chapter4:  Multivariate Analysis
Chapter6:  Time Series Analysis
Chapter6a:  Time series analysis, Part 1
Chapter6b:  Time series analysis, Part 2
Chapter6c:  Cross-spectral analysis
Chapter6d:  Other spectral techniques (SST, Multi-Taper, Maximum Entropy,etc.)

作者称：I also have a lot ofmatlab functions, including my matlab libraries, here:
http://www.aos.wisc.edu/~dvimont/matlab
EOF Analysis
eof_routine.m
cosweight.m
Suppose you want to perform EOF analysis on tropical Pacific SST from the NCEPreanalysis.  Here's a quick way to do that, assuming you've downloaded mymatlab libraries, above:

lims = [120 285 -30 30];
tim = get_time(1948,2002,1948);
[sst, lat, lon] = getnc2('skt.mon.mean.nc', 'skt', lims, 1, tim);
land = getnc2('land.sfc.gauss.nc', 'land', lims, 1, 1);
sst = remove_anncyc(sst);
sst = cosweight(sst, lat);
[lam, lds, pcs, per] = eof_routine(sst(:,land==0));

%  Get a scaled map of the leading EOF
eof1 = repmat(NaN, [length(lat) length(lon)]);
eof1(land == 0) = lds(:,1)*sqrt(lam(1));
eof1 = uncosweight(eof1, lat);
chi2example.m:  
1.  A demo for using the Chi2 distribution to determine confidence limitson variance
2.  Shows how to construct an AR1 time series, compute its power spectrum,and compare the spectrum to a theoretical AR1 power spectrum with confidencelimits derived from the chi2 distribution.
Example1_16_Feb.m:
This example demonstrates univariate and multivariate linear regression inMATLAB.  Regression is used to predict summer ozone concentrations alongthe Gulf Coast from meteorologicalvariables.  More information about this data set can be found from "Case Studies in Environmental Statistics",from the Geophysical Statistics Project, through NCAR.  This exampleshould provide enough information to solve problem 1 on HW2.
This routine uses:  Regression_example.mat

Example2_16_Feb.m:
This example demonstrates how to generate regression maps of SST or 500mbheight onto the CTI.  This requires accounts on the 14th floor computer toaccess the data.  This example demonstrates how you might solve Problem 2from HW2 (2005)

HW2_prob3_code.m:
This example demonstrates how you would go about solving problem 3 on HW2(2005).

eof_example.m:
This example demonstrates EOF/PC analysis on a simple matrix.  Thiscalculates the EOF's in two ways:  by eigenanalysis of the covariancematrix, and by singular value decomposition of the data matrix.
define_pdo.m:
This is the code used to calculate the PDO time series in Homework 3. Note that it uses a bunch of routines from my MATLAB libraries, at:
     http://www.aos.wisc.edu/~dvimont/matlab

define_aleutianlow.m:
This is the code used to calculate the Aleutian Low time series in Homework3.  Note that it uses a bunch of routines from my MATLAB libraries, at:
     http://www.aos.wisc.edu/~dvimont/matlab

weekly_extremes.m:This routine generates anAR(1) time series, then calculates extreme values in that time series relativeto the rest of the week.

These two routines can be used in Homework 4:autocorr.m
partial_autocorr.m

spectral_analysis.m:Fun with spectralanalysis, including generating AR1 and AR2 null hypotheses.  There's alsoa routine for calculating alph1 and alph2 for the AR2 null hypothesis. Hey, check out the similarity between calculating alph1 and alph2 herein, and"partial_autocorr.m", above.  You can do this using either.

filter_nr.m:Non-recursive filteringroutine.
freq_stuff.m:
This routine calculates the frequency response function for a couple offilters.  Also, the IRF, etc...
tutorial.m:This is the companion tothe tutorial handed out in class.

EOF analysis example scripts:  EOF analysis applied to the followingimage:


First, download the image to your ~/matlab/Data directory.  Then, downloadthe following scripts to an appropriate directory:
eof_pics.m
plot_island_demo.m
subplot2.m
global_axes.mfigure_tall.m
The EOF analysis example is in 'eof_pics.m', and the plotting routine is in'plot_island_demo.m'.  The other routines are used in the plottingroutines.

Here are some other routines that might be useful:
wysiwyg.m
worldlo.m
worldlo.mat
corr.m

注：Dan Vimont副教授在页面要求：If you use these files, please acknowledge their source at http://www.aos.wisc.edu/~dvimont/matlab。附件matlab.tar.gz为作者提供下载的全部自编函数库。

liutaoz21

顶楼主！感谢分享！Mark~

kongfeng0824

真是好贴！！！！！！！！！！！！！！

ssl

qxtlyf

真是好贴！！！！！！！！！！！！！！

郭衡

价廉物美，多谢多谢！

白云从空

这个好厉害，真的是不错

zihan6688

极其好的帖子。。。。

飞天猫

这个好厉害，真的是不错

XIN如止水

感谢lz分享~~