第七讲 数值数组及向量化运算优秀PPT.ppt

《第七讲 数值数组及向量化运算优秀PPT.ppt》由会员分享，可在线阅读，更多相关《第七讲 数值数组及向量化运算优秀PPT.ppt（24页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、第七讲 数值数组及向量化运算Xidian University1现在学习的是第1页，共24页Xidian University2数组创建n逐个构造n等差数列型:,linspacen等比数列型 logspace,semilogx/yn工作空间浏览器nload/save命令现在学习的是第2页，共24页Xidian University3n全1矩阵、全0矩阵、空矩阵、随机矩阵(均匀分布、正态分布)n对角阵、单位阵、魔方数组现在学习的是第3页，共24页Xidian University4数组的访问nA,A(:),A(:,:)nA(r,c),A(r,:),A(:,c)nA(I)nA(L)%sizeof(

2、A)=sizeof(L)nA(2:end,3:end)现在学习的是第4页，共24页Xidian University5nreshape,repmatnflipud,fliplr,rot90现在学习的是第5页，共24页Xidian University6数组运算nA,A.,nA+B,A-BnA*B,A.*B,A/B,A./B,AB,A.BnA.a,a.A现在学习的是第6页，共24页Xidian University7非数nNaN,inf,nisnan,isinf现在学习的是第7页，共24页Xidian University8特殊操作nfindnind2sub%单下标换全下标现在学习的是第8页，共

3、24页Xidian University9空数组na=nb=rand(0)nc=rand(1,0)nsize(a)size(b)size(c)nlength(a)length(b)length(c)nisempty(a)isempty(c)现在学习的是第9页，共24页Xidian University10关系操作n,=,=,=现在学习的是第10页，共24页Xidian University11逻辑操作n&,|,xorn&,|nall,anynfalse,true,logicalnisempty,isfinite,isprime,isreal,isletter,isspace,isinf,isn

4、annisa,ishandle,ischar,islogical,isglobal,isnumeric现在学习的是第11页，共24页Xidian University12数值计算n差分 diffn梯度 gradientn求和 sum:sum(A,1),sum(A,2)n累计求和 cumsumn近似数值积分 trapzn累计近似数值积分 cumtrapz现在学习的是第12页，共24页Xidian University13数值积分ntrapz(x,y,tol)nquad(fun,a,b,tol)ndblquad(fun,xmin,ymin,xmax,ymax,tol)现在学习的是第13页，共24页

5、Xidian University14求极值n一元函数极值 fminbnd(fun,x1,x2,options)fminbnd(x)x2+2*x-3,-1,1)n单纯形法求多元函数的极值fminsearch(fun,x0,options)现在学习的是第14页，共24页Xidian University15常微分方程n数值略现在学习的是第15页，共24页Xidian University16矩阵和代数方程n高斯列主元消去法nR,ci=rref(A)n特征值、特征向量、行列式nX=null(A);nA矩阵零空间的全部正交基nAX=0nZ=orth(A);nA矩阵值空间的全部正交基nspan(A)=

6、span(Z)现在学习的是第16页，共24页Xidian University17代数方程的求解nfsolve(fun,x0)nfzero(fun,x0)n对分法现在学习的是第17页，共24页Xidian University18概率分布图二项式分布二项式分布Binomialbino卡方分布卡方分布Chisquarechi2指数分布指数分布ExponentialexpF分布分布Ff几何分布几何分布Geometricgeo正态分布正态分布Normalnorm泊松分布泊松分布PoissonpoissT分布分布Tt均匀分布均匀分布Uniformunif离散均匀分布离散均匀分布Discrete Uni

7、formunid现在学习的是第18页，共24页Xidian University19nnorm,chi2,bino,f,geo,exp,poiss,t,unif,unidnnormpdf(x,mu,sigma)normcdf(x,nu,sigma)normrnd(nu,sigma,m,n)nrand(m,n)nrandn(m,n)现在学习的是第19页，共24页Xidian University20统计指令nmin(X)/max(X)nmean(X)nstd(X)nvar(X)ncov(X)ncorrcoef(X)nhist/histfit现在学习的是第20页，共24页Xidian Univer

8、sity21P159 例4.3-5nmu=2;nsigma=0.5;nx=randn(1000,1);ny1=sigma*x+mu;ny2=sigma*(x+mu);nsubplot(2,1,1);n histfit(y1),;naxis(-5 5 0 100);nsubplot(2,1,2);n histfit(y2);n axis(-5 5 0 100);现在学习的是第21页，共24页Xidian University22多项式运算nc=conv(p1,p2)nq,r=deconv(p1,p2)npoly2str(q,x)ny=polyval(p,x)现在学习的是第22页，共24页Xidi

9、an University23多项式拟合np=polyfit(x,y,n)现在学习的是第23页，共24页Xidian University24最小二乘法x0=(-10:0.5:10);x0=(-10:0.5:10);a0=2-2 1 7;y0=polyval(a0,x0);y1=y0+(rand(size(x0)-0.5)*800;m=length(x0);n=3;X=zeros(m,n+1);for k=1:n X(:,n-k+1)=x0.k;endX(:,n+1)=ones(m,1);aT=(Xy1);y=polyval(aT,x0);plot(x0,y0,:,x0,y1,.,x0,y,-);hold on;plot(x0 x0,y0 y1)disp(aT)legend(source data,noise added,fitting data);现在学习的是第24页，共24页