第七講 數(shù)值數(shù)組及向量化運(yùn)算

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

隨機(jī)矩陣(均勻分布、正態(tài)分布)對角陣、單位陣、魔方數(shù)組第3頁，共24頁，2023年，2月20日，星期三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)第4頁，共24頁，2023年，2月20日，星期三XidianUniversity5reshape,repmatflipud,fliplr,rot90第5頁，共24頁，2023年，2月20日，星期三XidianUniversity6數(shù)組運(yùn)算A’,A.’,A+B,A-BA*B,A.*B,A/B,A./B,A\B,A.\BA.^a,a.^A第6頁，共24頁，2023年，2月20日，星期三XidianUniversity7非數(shù)NaN,inf,isnan,isinf第7頁，共24頁，2023年，2月20日，星期三XidianUniversity8特殊操作findind2sub%單下標(biāo)換全下標(biāo)第8頁，共24頁，2023年，2月20日，星期三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)第9頁，共24頁，2023年，2月20日，星期三XidianUniversity10關(guān)系操作<,<=,>,>=,==,~=第10頁，共24頁，2023年，2月20日，星期三XidianUniversity11邏輯操作&,|,~,xor&&,||all,anyfalse,true,logicalisempty,isfinite,isprime,isreal,isletter,isspace,isinf,isnanisa,ishandle,ischar,islogical,isglobal,isnumeric第11頁，共24頁，2023年，2月20日，星期三XidianUniversity12數(shù)值計算差分diff梯度gradient求和sum:sum(A,1),sum(A,2)累計求和cumsum近似數(shù)值積分trapz累計近似數(shù)值積分cumtrapz第12頁，共24頁，2023年，2月20日，星期三XidianUniversity13數(shù)值積分trapz(x,y,tol)quad(fun,a,b,tol)dblquad(fun,xmin,ymin,xmax,ymax,tol)第13頁，共24頁，2023年，2月20日，星期三XidianUniversity14求極值一元函數(shù)極值fminbnd(fun,x1,x2,options)

fminbnd(@(x)x^2+2*x-3,-1,1)單純形法求多元函數(shù)的極值fminsearch(fun,x0,options)第14頁，共24頁，2023年，2月20日，星期三XidianUniversity15常微分方程數(shù)值略第15頁，共24頁，2023年，2月20日，星期三XidianUniversity16矩陣和代數(shù)方程高斯列主元消去法[R,ci]=rref(A)特征值、特征向量、行列式X=null(A);A矩陣零空間的全部正交基AX=0Z=orth(A);A矩陣值空間的全部正交基span(A)=span(Z)第16頁，共24頁，2023年，2月20日，星期三XidianUniversity17代數(shù)方程的求解fsolve(fun,x0)fzero(fun,x0)對分法第17頁，共24頁，2023年，2月20日，星期三XidianUniversity18概率分布圖二項式分布Binomialbino卡方分布Chisquarechi2指數(shù)分布ExponentialexpF分布Ff幾何分布Geometricgeo正態(tài)分布Normalnorm泊松分布PoissonpoissT分布Tt均勻分布Uniformunif離散均勻分布DiscreteUniformunid第18頁，共24頁，2023年，2月20日，星期三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)第19頁，共24頁，2023年，2月20日，星期三XidianUniversity20統(tǒng)計指令min(X)/max(X)mean(X)std(X)var(X)cov(X)corrcoef(X)hist/histfit第20頁，共24頁，2023年，2月20日，星期三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]);第21頁，共24頁，2023年，2月20日，星期三XidianUniversity22多項式運(yùn)算c=conv(p1,p2)[q,r]=deconv(p1,p2)poly2str(q,‘x’)y=polyval(p,x)第22頁，共24頁，2023年，2月20日，星期三XidianUniversity23多項式擬合p=polyfit(x,y,n)第23頁，共24頁，2023年，2月20日，星期三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);