资源描述
第4章Matlab的其他函数库,第1节数据分析数据库（datafun）4.1.1基本的数据分析,基本数据处理功能是按列进行，行向则表示数据的不同的样本。P5859表例data=154498367;158998175;1551006886;145637596;145637596;.141556575;155566485;147898777;1479654100;145607667std(标准差)指各列中N个元素与该列平均值的平方和开方trapz（求积分）梯形法求积分近似求元素和,4.1.2用于场论的数据分析,cradient求二维场和三维场的近似梯度del2求二维场和三维场的拉普拉斯算子cross两个相量的矢量积（叉乘）cot两个相量的数量积（点乘）,4.1.3用于随机数据分析的函数,rand（m，n）产生在0和1之间均匀分布的m行n列的随机数矩阵randn（m,n）产生正态分布的m行n列随机数矩阵（均值为0，标准差为1）x=rand(1,1000)；hist(x)y=randn(1,1000)；hist(y,50),4.1.4用于相关分析和傅立叶分析的函数,用于信号的时域分析和频域处理。corrcoef给出两个同长信号的相关系数cov（x,y）给出协方差矩阵conv（x,y）给出x,y的卷积filter（b,a,y）根据输入信号x和线性系统求输出信号的函数Xfft（x，N）求出时域信号x的离散傅立叶变换Xt=0:.001:3;u=sin(300*t)+2*cos(200*t);U=fft(u);plot(abs(U),第2节矩阵的分解与变换（matfun函数库）4.2.1线性方程组的系数矩阵,det（a）求方阵a的行列式rank（a）秩trace（a）矩阵主对角线上元素的和（迹）inv（a）逆阵pinv(a)(P62-63),4.2.2矩阵的分解p63,三角分解(lu分解)l,u=lu(a)正交分解（qr分解）q,r=qr(a)奇异值分解（svd分解）u,s,v=svd(a),4.2.3矩阵的特征值分析,e,r=eig(a)p=poly(a)roots(p),4.2.4特殊矩阵库(specmat)elmat,常用的矩阵函数,第3节多项式函数（polyfun）,一元高次代数多项式,Matlab中的表示方法,aa(1),a(2),a(n),a(n+1)幂次隐含在系数元素离向量右端的元素间隔中例：a(x)=3x3+2x+1则a3,0,2,1,多项式行向量的构造直接键入，无该次幂则取零用poly(A)命令来生成矩阵A的特征多项式。如果A是形如a0a1an-1an的向量，由命令poly(A)可生成(x-a0)(x-a1)(x-an-1)(x-an)所对应的多项式。,例：已知向量A=13480，用此向量构造一多项式并显示结果。A=1-34-80;PA=poly(A)PA=11132606-2720PAX=poly2sym(PA,X)PAX=X3+113*X2+2606*X-2720,4.3.2多项式运算,多项式的加减运算运算符进行加减运算的多项式应该具有相同的阶次，如果阶次不同，低阶的多项式必须用零添补至高阶多项式的阶次。例：求两个多项式a(x)=5x4+4x3+3x2+2x+1和b(x)=3x2+0 x+1的和。a=54321;b=301;c=a+00bc=54622,手工数再补零的方式不可取,functiony=polyadd(x1,x2)n1=length(x1);n2=length(x2);ifn1n2x2=zeros(1,n1-n2),x2;elseifn1fplot(sin,0,2*pi)fplot(sin,0,2*pi)句柄函数fplot(x)sin(x.2+y.2),0,2*pi)匿名函数,ezplot,ezplot(sin)默认区间-2*pi,2*piezplot(sin,0,2*pi)%指定区间画图ezmesh(x,y)sin(x.2+y.2)ezsurf(x,y)sin(x.2+y.2),任意函数的数值积分定积分子程序quad（函数名,初值x0,终值xf）微分方程数字解Ode23(fname,自变量初值x0，自变量终值xf,因变量初值y0)ode45(fname,自变量初值x0，自变量终值xf,因变量初值y0),任意函数的数值积分,第5节字符串函数库p77,字符串赋值size(s);abs(s)setstr(abs(s);dec2hex(abs(s)字符串语句的执行eval字符串输入输出disp（）input()sprintf(),第6,7,8节,
展开阅读全文