第七講數(shù)值數(shù)組及向量化運算XidianUniversity1第一頁，共二十四頁，2022年，8月28日XidianUniversity2數(shù)組創(chuàng)建逐個構(gòu)造等差數(shù)列型:,linspace等比數(shù)列型logspace,semilogx/y工作空間瀏覽器load/save命令第二頁，共二十四頁，2022年，8月28日XidianUniversity3全1矩陣、全0矩陣、空矩陣、

隨機矩陣(均勻分布、正態(tài)分布)對角陣、單位陣、魔方數(shù)組第三頁，共二十四頁，2022年，8月28日XidianUniversity4數(shù)組的訪問A,A(:),A(:,:)A(r,c),A(r,:),A(:,c)A(I)A(L)%sizeof(A)==sizeof(L)A(2:end,3:end)第四頁，共二十四頁，2022年，8月28日XidianUniversity5reshape,repmatflipud,fliplr,rot90第五頁，共二十四頁，2022年，8月28日XidianUniversity6數(shù)組運算A’,A.’,A+B,A-BA*B,A.*B,A/B,A./B,A\B,A.\BA.^a,a.^A第六頁，共二十四頁，2022年，8月28日XidianUniversity7非數(shù)NaN,inf,isnan,isinf第七頁，共二十四頁，2022年，8月28日XidianUniversity8特殊操作findind2sub%單下標(biāo)換全下標(biāo)第八頁，共二十四頁，2022年，8月28日XidianUniversity9空數(shù)組a=[]b=rand(0)c=rand(1,0)size(a)size(b)size(c)length(a)length(b)length(c)isempty(a)isempty(c)第九頁，共二十四頁，2022年，8月28日XidianUniversity10關(guān)系操作<,<=,>,>=,==,~=第十頁，共二十四頁，2022年，8月28日XidianUniversity11邏輯操作&,|,~,xor&&,||all,anyfalse,true,logicalisempty,isfinite,isprime,isreal,isletter,isspace,isinf,isnanisa,ishandle,ischar,islogical,isglobal,isnumeric第十一頁，共二十四頁，2022年，8月28日XidianUniversity12數(shù)值計算差分diff梯度gradient求和sum:sum(A,1),sum(A,2)累計求和cumsum近似數(shù)值積分trapz累計近似數(shù)值積分cumtrapz第十二頁，共二十四頁，2022年，8月28日XidianUniversity13數(shù)值積分trapz(x,y,tol)quad(fun,a,b,tol)dblquad(fun,xmin,ymin,xmax,ymax,tol)第十三頁，共二十四頁，2022年，8月28日XidianUniversity14求極值一元函數(shù)極值fminbnd(fun,x1,x2,options)

fminbnd(@(x)x^2+2*x-3,-1,1)單純形法求多元函數(shù)的極值fminsearch(fun,x0,options)第十四頁，共二十四頁，2022年，8月28日XidianUniversity15常微分方程數(shù)值略第十五頁，共二十四頁，2022年，8月28日XidianUniversity16矩陣和代數(shù)方程高斯列主元消去法[R,ci]=rref(A)特征值、特征向量、行列式X=null(A);A矩陣零空間的全部正交基AX=0Z=orth(A);A矩陣值空間的全部正交基span(A)=span(Z)第十六頁，共二十四頁，2022年，8月28日XidianUniversity17代數(shù)方程的求解fsolve(fun,x0)fzero(fun,x0)對分法第十七頁，共二十四頁，2022年，8月28日XidianUniversity18概率分布圖二項式分布Binomialbino卡方分布Chisquarechi2指數(shù)分布ExponentialexpF分布Ff幾何分布Geometricgeo正態(tài)分布Normalnorm泊松分布PoissonpoissT分布Tt均勻分布Uniformunif離散均勻分布DiscreteUniformunid第十八頁，共二十四頁，2022年，8月28日XidianUniversity19norm,chi2,bino,f,geo,exp,poiss,t,unif,unidnormpdf(x,mu,sigma)

normcdf(x,nu,sigma)

normrnd(nu,sigma,m,n)rand(m,n)randn(m,n)第十九頁，共二十四頁，2022年，8月28日XidianUniversity20統(tǒng)計指令min(X)/max(X)mean(X)std(X)var(X)cov(X)corrcoef(X)hist/histfit第二十頁，共二十四頁，2022年，8月28日XidianUniversity21P159例4.3-5mu=2;sigma=0.5;x=randn(1000,1);y1=sigma*x+mu;y2=sigma*(x+mu);subplot(2,1,1);histfit(y1),;axis([-550100]);subplot(2,1,2);histfit(y2);axis([-550100]);第二十一頁，共二十四頁，2022年，8月28日XidianUniversity22多項式運算c=conv(p1,p2)[q,r]=deconv(p1,p2)poly2str(q,‘x’)y=polyval(p,x)第二十二頁，共二十四頁，2022年，8月28日XidianUniversity23多項式擬合p=polyfit(x,y,n)第二十三頁，共二十四頁，2022年，8月28日XidianUniversity24最小二乘法x0=(-10:0.5:10)';x0=(-10:0.5:10)';a0=[2-217];y0=polyval(a0,x0);y1=y0+(rand(size(x0))-0.5)*800;m=length(x0);n=3;X=zeros(m,n+1);fork=1:nX(:,n-k+1)=x0.^k;endX(:,n+1)=ones(m,1);

aT=(X\y1)';y=polyval(aT,x0);