《常用统计分析软件matlab课件》由会员分享,可在线阅读,更多相关《常用统计分析软件matlab课件(50页珍藏版)》请在金锄头文库上搜索。
1、Introduction顾世梁扬州大学农学院 87979358(o), 87328797(h) 1.1 Matlab 的特点: 1. 矩阵处理无所不能; 2. 计算功能包罗万象; 3. 画图功能丰富多彩; 4. 操作应用方便直观.1 简介1.由Matrix Laboratory 1967年推出。 Mathworks于1984年 发布Matlab。 2. Matlab 不断改进,现 在为7.4版本。 3. 有些网站有下载。1.2 发生发展1.3 安装硬件要求不太高, 一般的台式机、笔 记本均可安装。2 基本操作 2.1 视窗环境2.1.1 命令窗口: 命令和数据的输入输出2.1.2 M文件窗口
2、:编辑源程序文件和调试程序2.1.3 工作空间窗口:存放变量相关信息2.1.4 当前目录窗口:存放M或函数文件的工作目录 2.1.5 命令历史窗口:命令的历史纪录 2.1.6 起始面板窗口:工具箱组合 2.1.7 图形窗口:画图 2.1.8 帮助窗口:功能强大的帮助 2.2 赋值1) x= 1 2 3 4(直接键盘输入);2) x=2:.5:20, x=a:i:b, x=a:b, x=a:i:b;3) x=linspace(1,10,20), x=linspace(1,20);4) x=rand(3,4), x=randn(m,n),x=random(name,a,b,c,m,n),x=Ran
3、dom(unid,k,m,n)5)x,y= meshgrid(-2:.1:2)6) import (通过array editor)7) Load data.dat, data.txt, data.csv, load filename or load filename8) xlsread exam.xls, xlsread(filename)9) 复制与粘贴除此以外:Matlab可读写几十种数据、 文字、图形、音像、压缩文件。imread -Read image from graphics file. p=imread(sxh.jpg), image(p)wavread -Read (WAV)
4、sound file.Sd=wavread(music.wav), sound(Sd, 22500)aviread -Read movie (AVI) file. mv=Aviread(vip.avi), movie(mv, 2) 2.3 矩阵基本数据操作函数 极大、极小:max, min 总和、平均:sum(x), mean(x) 中位数:median, iqr, prctile(x, p) 偏度系数skewness, 峰度系数 kurtosis 排序:sort 乘积:prod 标准差:std, std(a,1) 方差:var, var(x,1); 协方差 cov 相关系数:corrcoef
5、(x, y), corr(x) 2.4 矩阵运算2.4.1 矩阵运算:加减: a+b, a-b乘法: a*b, a2除法:左除 ab, 右除 a/b2.4.2 数组运算乘法: a.*b除法: a./b, a.b幂函数:a.2, a.b2.4.3 矩阵函数运算 C=sqrtm(A), B=expm(A), logm(A), C=funm(A, fun), fun=sqrtm,sin 2.4.4 逻辑运算 逻辑关系:A=b, Ab, ab, a=b, a=b 逻辑与:A stem(x, y , filled)画竖线 Text(x, y, str), (*str必须为列向量)注 Text(x, y,
6、 str, fontsize,15) Hold on; hold off; Subplot(2,3,2), subplot(2,3,6), subplot(1,1) Plot 绘图函数的参数参数意义参数意义 rRed-红色-实线实线 ggreen-绿色-虚线线 bBlue-蓝色:点线线 yYellow-黄色-.点划线线 mMagenta-深红o圆圆圈 cCyan-青蓝x叉号 wWhite-白色+加号 kBlack-黑色s方形 *-星号d菱形.点号 3.2 三维画图mesh(x, y, z)surf(x, y, z)scatter3(x, y, z)WaterfallStem3(x, y, z)
7、, stem3(x, y, z, filled)Text(x, y, z, str), text(x, y, z, str, fontsize, 12) 3.3 统计作图 bar(x), bar3(x), area(x), pie(x), pie3(x) hist(y), hist(y, x), stem(x, y), stem3(x, y, z) scatter(x, y), scatter(x, y, s, c), scatter3(x, y, z) text(x, y, str), text(x, y, z, str,fontsize,15) fplot(fun, a b); fplot
8、(100/(1+20*exp(-.4*x),0 20 ezplot(sin(3*t)*cos(t),sin(3*t)*sin(t),0,pi) ezplot3(cos(t),t.*sin(t),sqrt(t),0,4*pi)fplot(tpdf(x,2),-4 4), hold on fplot(normpdf(x),-4 4) 4 符号运算 4.1 多项式运算 多项式赋值 p=2 1 3 5 由根生成多项式 poly(p) 多项式转换 poly2sym(p) 多项式加、减 p+q, 多项式求根 roots(p) 多项式乘 、除conv(p, q),deconv(p, q) 多项式取值 pol
9、yval(p, x) 多项式拟合 polyfit(x, y, k) 4.2 符号运算 泰勒级数:taylor(sin(x), taylor(exp(x),n); 级数和: symsum(k), symsum(k2,0,n); 求极限:limit(xx), limit(sin(x)/x,0); limit(1+1/x)x), limit(1+1/x)x,inf); 方程精简: pretty, simplify, simple; 微、积分: diff, int, quad解方程(solve), solve(a*x2+b*x+c=0) x y=solve(3*x+4*y-10,4*x-7*y+3)
10、x y z=solve(4*x-6*y+15*z-3, x+3*y-7*z+16,3*x-17*y-20*z-5) 解微分方程(dsolve), dsolve(Dy+y+c,) y=dsolve(Dy-b*y+c*y2,y(0)=k/(1+a) syms c b k; y=subs(y,c b/k); pretty(y)解非线性方程组(fsolve)function F=myfun(x) F = 2*x(1) - x(2) - exp(-x(1);-x(1) + 3*x(2) - exp(-x(2); %将其存入myfun.m x0=1 1; fsolve(myfun,x0)5 函数和编程 5
11、.1自定义函数1) 编辑m文件,2) function开头,3) 定义 依变量(输出变量)名=函数名(自变量,参数) ,4) 按 (数组而非矩阵运算)格式输入函数 ,5) 存贮。Peaks函数:二维正态分布概率密度函数: 5.2 编程程序文件的几种结构 顺序结构(sample1) 循环结构(sample2) 分支结构(sample3) 在m文件中Run运行在命令窗口中运行http:/bbs.bc- 统计分析 6.1 简单数据计算 max, min, mean, median, sum, prod; geomean, harmmean; nanmax, nanmean, nanmedian; n
12、anmin, nanstd, nansum; factorial, primes, perms(x1 x2 x3 x4; mod(x,n), round, rem; fix, floor, ceil, sign 6.2 概率分布 6.2.1 计算概率密度函数值(pdf): f=binopdf(x, n, p), normpdf(x, mu,), tpdf(x, v), chi2pdf(x, v), fpdf(x, v1, v2), 6.2.2 计算累积分布函数(cdf, 概率值): P=binocdf(x, n, p), normcdf(x, u, ), tcdf(x, v), fcdf(F,
13、 v1,v2), chi2cdf(x, v), 6.2.3 产生服从一定分布的随机数 rand(3, 4), randn(5, 6), normrnd(mu, ), y=random(name,a1,a2,a3,m,n), y=random(Nomal, 0,1,3,5), random(poiss,np,m,n), y=random(unif,a,b,m,n), y=random(unid, k, m,n), 用help random查看name 6.3 常用统计计算 mean, median, geomean, harmmean; max, min, sum, sort; var(x),v
14、ar(x,1), std(x), std(x,1), cov; corrcoef, corr; kurtosis, skewness , size; tabulate; nanmax, nanmean, nanmedian; nanmin, nanstd, nansum; range, iqr, prctile(x, p), quantile(x, p) 6.4 假设测验ztest(利用正态分布测验) h, p, ci, z=ztest(x, m, , , tail)ttest (单个样本的 t 测验) h, p, ci, stat=ttest(x, m, , tail)ttest2 (两个样
15、本的 t 测验) h, p, ci, stat=ttest2(x1, x2, , tail) 6.5 方差分析 anova1 (单向分组资料方差分析) p, table, stats=anova1(x) multcompare(stats)anova2 (两向分组资料方差分析) p, table, stat=anova2(x) p, table, stat=anova2(x, n) multcompare(stat,display,ctype, est, dim)Anovan (多因素资料的方差分析) Anovan(x, group, model) 三因素model=1 2 3 4 5 6 7 (三因素方差分析编码表)数值值含义义 1A(主效) 2B(主效) 3A*B(互作) 4C(主效) 5A*C(互作) 6B*C(互作) 7A*B*C(互作)四因素方差分析编码表(model)数 值含 义数 值含 义 1A(主效)9A*D 2B(主效)10B*D 3A*B(互