《气象统计应用11逐步回归实例课件.ppt》由会员分享，可在线阅读，更多相关《气象统计应用11逐步回归实例课件.ppt（34页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、气象统计应用气象统计应用海洋环境学院海洋环境学院1魏晓敏魏晓敏5江会常江会常9谭云霞谭云霞13杨吕玉慈杨吕玉慈安子敬安子敬井苗苗井苗苗田云菲田云菲杨挺杨挺常海洁常海洁李丹李丹田志平田志平杨雨晗杨雨晗陈昌硕陈昌硕李松乘李松乘王海蛟王海蛟姚佳伟姚佳伟陈青陈青李晓东李晓东王浩王浩叶俊安叶俊安2陈顺鑫陈顺鑫6李晓霞李晓霞10王旌旭王旌旭14尹相洲尹相洲陈晓丹陈晓丹刘鹏刘鹏王凯悦王凯悦尤樱杰尤樱杰丁芊丁芊刘琪瑶刘琪瑶王可王可岳欣健岳欣健范志伟范志伟刘珊刘珊王昕王昕张丹丹张丹丹冯嘉烨冯嘉烨潘锡飞潘锡飞王砚硕王砚硕张凯张凯3高荣璐高荣璐7彭名杨彭名杨11魏海伦魏海伦15张露张露耿雨耿雨邱静怡邱静怡魏亮魏亮周
2、茂铮周茂铮顾瑜顾瑜舒帅舒帅魏胜玮魏胜玮周曦林周曦林郭林郭林宋海凡宋海凡吴一乔吴一乔朱洪洲朱洪洲郭鹏嘉郭鹏嘉宋丽潞宋丽潞吴玥吴玥朱金国朱金国4郭润霞郭润霞8孙华庆孙华庆12熊超伟熊超伟16朱金双朱金双贺冰蕊贺冰蕊孙健翔孙健翔徐杉徐杉朱良帅朱良帅黄晶伟黄晶伟孙宁孙宁许杰许杰李时越李时越黄菊黄菊孙晓梅孙晓梅薛琦薛琦孙艺孙艺黄山黄山孙宇孙宇杨超杨超吴漾、付迪吴漾、付迪逐步回归实例气候与孢粉气候与孢粉气候植被孢粉气候变化孢粉百分比含量孢粉是什么孢粉外力搬运风,水,动物等降落(表土样本)化石孢粉(地层样本)地质作用土壤成岩孢粉特点：体积小，数量多：“花粉雨”；孢粉壁(孢粉素:C96H22O24)易于保存，
3、极难氧化，高温下也难溶解，放在浓酸里也不溶解。针刺阿拉曲孢针刺阿拉曲孢Aratrisporites tenuispinosusPlayford,1965 挺直斑马纹孢挺直斑马纹孢Zebrasporites rectusZhang，1984 华丽虎纹孢华丽虎纹孢Tigrinispora elegansChang,1974et1984 小圈环纹孢小圈环纹孢 Annulispora microannulatadeJeersey,1962 轮状花冠孢轮状花冠孢 Coronasporites coronaris(Chang)Zhang,1984 瘤状厚瘤孢瘤状厚瘤孢Ricciisporites tube
4、rculatus Lundblad1954孢粉是什么数据和方法数据和方法两类数据两类数据:surfacesamplesandstratumsamples.气候变量气候变量:TjulandPannareobservedat25to2910minShennongjia.方法：转换函数法方法：转换函数法孢粉（木本）百分比含量随海拔高度的变化孢粉（木本）百分比含量随海拔高度的变化孢粉（草本蕨类）百分比含量随海拔高度的变化孢粉（草本蕨类）百分比含量随海拔高度的变化地层样本步骤步骤1.利用地表样本，选择合适的因子（孢粉利用地表样本，选择合适的因子（孢粉类型），建立气候参数与孢粉百分比的类型），建立气候参数
5、与孢粉百分比的线性回归方程。线性回归方程。2.将地层样本的孢粉百分比含量带入回归将地层样本的孢粉百分比含量带入回归方程，建立气候参数变化的时间序列。方程，建立气候参数变化的时间序列。选择标准选择标准表土样本中含量较高表土样本中含量较高地层样本中连续性好地层样本中连续性好（1）读取文本文件）读取文本文件1.打开文件打开文件fopen2.读取文件读取文件fscanf3.关闭文件关闭文件fclose1.fopen格式：格式：fid=fopen(数据文件数据文件,r)r读读w写写a附加附加文件先要打开，才可进行读写操作。文件先要打开，才可进行读写操作。2.fscanf格式：格式：A=fscanf(fi
6、d,format,size)formatformat用以控制读取的数据格式，由用以控制读取的数据格式，由%加上格加上格式符组成，常见格式符：式符组成，常见格式符：%s%s：字符串：字符串s=fscanf(fid,%s)%d%d：十进制数：十进制数a=fscanf(fid,%5d)读读5位整数位整数%f%f：浮点数：浮点数b=fscanf(fid,%6.2f)读浮点数读浮点数size为为可可选选项项，用用于于控控制制读读取取的的数数据据，若若不不选选用用则则读读取取整整个个文文件件内内容容，若若选选用用它它的的值值可可以以是是下列值：下列值：N N 读取前读取前N N个元素到一个列向量。个元素到
7、一个列向量。inf inf 读取整个文件。读取整个文件。M,N M,N 读数据到读数据到MNMN的矩阵中，按列存放的矩阵中，按列存放 N N可以是可以是infinf，M M不可以。不可以。3.fclosen文文件件在在进进行行完完读读、写写等等操操作作后后，应应及及时时关关闭闭，以保证文件的安全。以保证文件的安全。n关闭文件命令格式为：关闭文件命令格式为：sta=fclose(fid)关闭打开的文件关闭打开的文件sta=fclose(all)关闭所有打开的文件关闭所有打开的文件sta表示关闭文件操作的返回代码，若关闭成表示关闭文件操作的返回代码，若关闭成功，返回功，返回0，否则返回，否则返回1
8、。fid1=fopen(table.txt,r);data=fscanf(fid1,%f);datadata=1491215262034234676fid=fopen(table.txt,r);title=fscanf(fid,%f,3,4);status=fclose(fid);titletitle=11220464153476926230注意：注意：fscanf在读取数据时，在读取数据时，是逐行进行的，在返回矩阵是逐行进行的，在返回矩阵时，是将数据逐列写入的时，是将数据逐列写入的写文本文件写文本文件1.打开文件打开文件fopen2.写文件写文件fprintf3.关闭文件关闭文件fclose
9、例例a=1:10;b=a;sqrt(a);fid3=fopen(s_table.dat,w);fprintf(fid3,平方根表平方根表n);fprintf(fid3,%2.0f%6.3fn,b);fclose(fid3);平方根表平方根表11.00021.41431.73242.00052.23662.44972.64682.82893.000103.162（2）读取二进制文件）读取二进制文件1.打开文件打开文件fopen2.读取文件读取文件fread3.关闭文件关闭文件fclose写二进制文件写二进制文件1.打开文件打开文件fopen2.写文件写文件fwrite3.关闭文件关闭文件fclo
10、sefwritecount=fwrite(fid,a,precision);a是待写入的矩阵，是待写入的矩阵，precision设定精度设定精度例，写读二进制文件例，写读二进制文件fid4=fopen(tob.bin,w);count=fwrite(fid4,magic(5),int32);status=fclose(fid4);生成一个文件名为生成一个文件名为tob.bin的二进制文件，的二进制文件，数据位数据位5阶方阵，类型为整型，输出阶方阵，类型为整型，输出count值为值为25。幻方：幻方：magic(n)生成一个生成一个n*n的矩阵，矩阵元的矩阵，矩阵元素是由整数素是由整数1到到n2组成的并且任何行任何列的组成的并且任何行任何列的和都相等，阶数和都相等，阶数n必须是大于等于必须是大于等于3的标量。的标量。读文件读文件fid5=fopen(tob.bin,r);data=(fread(fid5,25,int32);datadata=Columns1through2117234101124561218171319258142021215Columns22through25162239