[指南]数学软件应用复习

《[指南]数学软件应用复习》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、数学软件应用期末复习1、R软件基础知识1)用c()函数建立向量x<-c(10.4,5.6,3.1,6.4,21.4),用重复函数rep()建立由一些基本元素重复而成的向量w=rep(c(l,2,3),c(4,3,2));w#将1,2,3分别重复4,3,2次2)用m:ri建立从m到n的白然序列(【句量)a:b表示从a开始，逐项加1(或减1),直到b为止。a=l:5;a，用seq()建立等差数列(向量)seq(from=valuel,to=value2,by=value3)即从valuel开始，到value2

2、结束,间隔为value3«seq(-5,5,by=.2)->sl;sl3)用array()>X<-array(data=l:20,dim=c(4,5));l〜20,四行五夕U或matrix()A<-matrix(l:15,nrow=3,ncoI=5,byrovv=TRUE)•建立矩阵、访问数纟ft元素：用A[i,j](A[,j]、A[i,])略写某一维的下标，则表示该维全选；a[,,]或a[]都表示整个数组。:用t()转置矩阵,用rbind()纵向连接矩阵cbind(A,B)把其自变量横向拼接；用det(

3、)计算行列式。4)求A的逆矩阵用solve(A)；求方程组Ax=b可用x=solve(A,b)或x=qr.solve(A,b)5)求对称矩阵A的特征值、特征向量：可用ev=eigen(A);ev$values;evSvectors完成6)建立数据框对象:用data.frame(),也可用read.table()、read.csv()读取外部文件数据建立数据椎。7)用R进行回归分析的儿个配套指令2、Matlab基本操作1)MatLab的向量、矩阵多种定义方式c=[4;7;l;0;5]用空格和逗号分隔牛•成行

4、向量，用分号分隔生成列向量。x=a:步长：b课件22)MatLab常用的数组操作：点运算，查找、最人最小、求和等3、MatLab的常用线性代数运算：1)行列式a=det(A),逆矩阵inv(A),秩rank(A),矩阵乘A*B,矩阵除A/B,最简阶梯型rrefA=[213;132];formatratrref(A)format2)线性方程组求解(特解、通解p5859、非负解)x=lsqnonneg(A,b)nJ求出方程组Ax=b一个非负最小二乘解(特解)；斗+壬+土今兀・屯=7求线性方程组

5、3^42^4^

6、+^-3^5=-2的通解Xj+2马+2x4+6x5=235j

7、4xj+3X,4-3x4-x5=12（1）先判斯线性方程组解的情况，输入并运行以下代码：a=[l1211-3：01226;5433-1];b=[7-22312]*；三.Matlabft用宾■解：111;3rl=rank(a),r2=rank([ab])运行结果为rl=r2=2<5»说明该线性方程组有无穷多个解.(2)输入以下代码求出线性方程组的一个特解xO=ab运行结果为：xO=(3.1667.0,0.0,3.8333)'(3)输入代码求出导

8、出组aX=0的一个基础解系xx=null(a)运行结果是：Xx=0.75300.0176-0.0000-0,4167-0.7464-0.0000-0.30430.4533-0.7071-0.30430.45330・70710.2723-0.1778-0.0000(4)线性方程组的通解为31667、r0.7530、f0.0176、r0、0-0.4167-0.74640X二0-0.3043+60.4533+G-0.70710-030430.453307071、3・8333,、0.2723丿T1778丿、0>其

9、中q,C2,C3为任意常数。3)特征值和特征向量等[V,D]=eig(A)v是特征值组成的矩阵，E=eig(A)特征向量4、MatLab得常用微积分计算(符号与数值)：课件51)导数/偏导数diff,原函数/定积分int,极限limit等2)函数的局部最优值3)MatLab微分方程求解(符号与数值)课件83、MatLab流程控制语句和自定义函数的编写课件34、MatLab基木的曲线曲面绘制1)用plot指令绘制2维曲线图，能对自己定义的函数进行作图2)用meshgrid指令生成曲而网格点的底而投影坐标。用

10、mesh或surf指令绘制3维曲面5、MatLab优化问题1)线性规划2)二次规划3)非线性规划6、MatLab插值•拟合1)用interpl作线性插值或样条插值曲线2)用lsqcurvefit作

11、11

12、线拟合模拟练习limXT8(x-2symsx;f=((x-2)/(x+2))Axlimit(f,x,inf)x2xdxsymsx;x2e~Adxint(xA2*log(x))symsx;int(xA2*exp(-x),l