MTLB实验练习题计算机南邮MATLAB数学实验大作业答案
MATLAB实验：运算基础,并且附有答案
实验二、MATLAB运算基础一、实验目的掌握MATLAB各种表达式的书写规则及常用函数的使用。
掌握MATLAB中字符串、元胞数组和结构的常用函数的使用。
二、实验内容及步骤1、设有矩阵A和B，A=[1 2 3 4 5;6 7 8 9 10;11 12 13 14 15;16 1718 19 20;21 22 23 24 25]，B=[3 0 16;17 -6 9;0 23 -4;9 7 0;4 1311]1）求它们的乘积C >>C=A*B2）将矩阵C的右下角3x2子矩阵赋给D >>I=[3 4 5];J=[2 3];D=C(I,J)也可以用>>D=C([3 4 5],[2 3])D =520 397705 557890 7172、完成下列操作1）求[100，999]之间能被61整除的数及其个数（提示：先利用冒号表达式，再利用find和length函数。
）>> a=100:999;find(rem(a,61)==0)ans =23 84 145 206 267 328 389 450 511572 633 694 755 816 877>> b=a(ans)b =122 183 244 305 366 427 488 549 610 671 732 793 854 915 976>> length(b)ans =152）建立一个字符串向量，删除其中的大写字母（提示：利用find函数和空矩阵。
）a=’I am maying’;a( find(a>’A’&a<’Z’))=[]3、已知A=[23 10 -78 0；41 -45 65 5；32 5 0 32；6 -54 92 14]，取出其前3行构成矩阵B，其前两列构成矩阵C，其左下角3x2子矩阵构成矩阵D，B与C的乘积构成矩阵E，分别求E<D、E&D、E|D、~E|~D。
MATLAB全部实验及答案
MATLAB全部实验及答案实验一、MATLAB基本操作实验内容及步骤1、命令窗口的简单使用（1）简单矩阵的输入（2）求[12+2×(7-4)]÷32的算术运算结果2、有关向量、矩阵或数组的一些运算（1）设A=15；B=20；求C=A+B与c=a+b？（2）设A=[1 2 3;4 5 6;7 8 9]，B=[9 8 7;6 5 4;3 2 1]；求A*B与A.*B？A*B就是线代里面的矩阵相乘 A.*B是对应位置的元素相乘（3）设a=10，b=20；求i=a/b=0.5与j=a\b=2？（4）设a=[1 -2 3;4 5 -4;5 -6 7]请设计出程序，分别找出小于0的矩阵元素及其位置（单下标、全下标的形式），并将其单下标转换成全下标。
clear,clca=[1 -2 3;4 5 -4;5 -6 7];[x,y]=find(a<0);c=[];for i=1:length(x)c(i,1)=a(x(i),y(i));c(i,2)=x(i);c(i,3)=y(i);c(i,4)=(y(i)-1)*size(a,2)+x(i);endc（5）在MATLAB命令行窗口运行A=[1,2;3,4]+i*[5,6;7,8]；看结果如何？如果改成运行A=[1,2;3,4]+i[5,6;7,8]，结果又如何？前面那个是虚数矩阵，后面那个出错（6）请写出完成下列计算的指令：a=[1 2 3;3 4 2;5 2 3]，求a^2=？,a.^2=？a^2= 22 16 1625 26 2326 24 28a.^2=1 4 99 16 425 4 9（7）有一段指令如下，请思考并说明运行结果及其原因clearX=[1 2;8 9;3 6];X( : ) 转化为列向量（8）使用三元组方法，创建下列稀疏矩阵2 0 8 00 0 0 10 4 0 06 0 0 0方法一：clear,clcdata=[2 8 1 4 6];ir=[1 1 2 3 4 ];jc=[1 3 4 2 1];s=sparse(ir,jc,data,4,4);full(s)方法二：不用三元组法clear,clca=zeros(4,4);a(1,[1,3])=[2,8];a(2,4)=1;a(3,2)=4;a(4,1)=6;a（9） 写出下列指令的运行结果>> A = [ 1 2 3 ]; B = [ 4 5 6 ];>> C = 3.^A>> D = A.^B3、 已知⎪⎭⎫ ⎝⎛+⋅=-334sin 234πt e y t 若需要计算t ∈[-1,1]，取间隔为0.01，试计算出相对应的y 值。
南邮MATLAB数学实验答案(全)[精品文档]
![南邮MATLAB数学实验答案(全)[精品文档]](https://img.taocdn.com/s3/m/3d22b17eb307e87101f6963b.png)
第一次练习教学要求：熟练掌握Matlab 软件的基本命令和操作，会作二维、三维几何图形，能够用Matlab 软件解决微积分、线性代数与解析几何中的计算问题。
补充命令vpa(x,n) 显示x 的n 位有效数字，教材102页fplot(‘f(x)’,[a,b]) 函数作图命令，画出f(x)在区间[a,b]上的图形 在下面的题目中m 为你的学号的后3位（1-9班）或4位（10班以上） 1.1 计算30sin limx mx mx x →-与3sin limx mx mxx →∞- syms xlimit((902*x-sin(902*x))/x^3) ans =366935404/3limit((902*x-sin(902*x))/x^3,inf) ans = 0 1.2 cos1000xmxy e =，求''y syms xdiff(exp(x)*cos(902*x/1000),2) ans =(46599*cos((451*x)/500)*exp(x))/250000 - (451*sin((451*x)/500)*exp(x))/250 1.3 计算221100xy e dxdy +⎰⎰dblquad(@(x,y) exp(x.^2+y.^2),0,1,0,1) ans = 2.13941.4 计算4224x dx m x +⎰ syms xint(x^4/(902^2+4*x^2)) ans =(91733851*atan(x/451))/4 - (203401*x)/4 + x^3/12 1.5 (10)cos ,x y e mx y =求syms xdiff(exp(x)*cos(902*x),10) ans =-356485076957717053044344387763*cos(902*x)*exp(x)-3952323024277642494822005884*sin(902*x)*exp(x)1.6 0x =的泰勒展式(最高次幂为4).syms xtaylor(sqrt(902/1000+x),5,x) ans =-(9765625*451^(1/2)*500^(1/2)*x^4)/82743933602 +(15625*451^(1/2)*500^(1/2)*x^3)/91733851-(125*451^(1/2)*500^(1/2)*x^2)/406802 + (451^(1/2)*500^(1/2)*x)/902 +(451^(1/2)*500^(1/2))/500 1.7 Fibonacci 数列{}n x 的定义是121,1x x ==12,(3,4,)n n n x x x n --=+=用循环语句编程给出该数列的前20项（要求将结果用向量的形式给出）。
（完整版）MATLAB）课后实验答案[1]
![（完整版）MATLAB）课后实验答案[1]](https://img.taocdn.com/s3/m/306ff29b690203d8ce2f0066f5335a8102d266cb.png)
（完整版）MATLAB）课后实验答案[1]实验⼀ MATLAB 运算基础1. 先求下列表达式的值，然后显⽰MATLAB ⼯作空间的使⽤情况并保存全部变量。
(1) 0122sin 851z e =+(2) 21ln(2z x =，其中2120.455i x +??=?- (3) 0.30.330.3sin(0.3)ln , 3.0, 2.9,,2.9,3.022a a e e az a a --+=++=--L (4) 2242011122123t t z t t t t t ?≤=-≤，其中t =0:0.5:2.5 解：4. 完成下列操作：(1) 求[100,999]之间能被21整除的数的个数。
(2) 建⽴⼀个字符串向量，删除其中的⼤写字母。
解：(1) 结果：(2). 建⽴⼀个字符串向量例如：ch='ABC123d4e56Fg9';则要求结果是：实验⼆ MATLAB 矩阵分析与处理1. 设有分块矩阵33322322E R A O S=?，其中E 、R 、O 、S 分别为单位矩阵、随机矩阵、零矩阵和对⾓阵，试通过数值计算验证2 2E R RS A O S +??=。
解: M ⽂件如下；5. 下⾯是⼀个线性⽅程组：1231112340.951110.673450.52111456x x x ??=???(1) 求⽅程的解。
(2) 将⽅程右边向量元素b 3改为0.53再求解，并⽐较b 3的变化和解的相对变化。
(3) 计算系数矩阵A 的条件数并分析结论。
解： M ⽂件如下：123d4e56g9实验三选择结构程序设计1. 求分段函数的值。
2226035605231x x x x y x x x x x x x ?+-<≠-?=-+≤<≠≠??--?且且及其他⽤if 语句实现，分别输出x=-5.0,-3.0,1.0,2.0,2.5,3.0,5.0时的y 值。
解：M ⽂件如下：2. 输⼊⼀个百分制成绩，要求输出成绩等级A、B、C、D、E。
MATLAB实验二运算基础答案
实验二、MA TLAB运算基础一、实验目的掌握MA TLAB各种表达式的书写规则及常用函数的使用。
掌握MA TLAB中字符串、元胞数组和结构的常用函数的使用。
二、实验内容及步骤1、设有矩阵A和B，A=[1 2 3 4 5;6 7 8 9 10;11 12 13 14 15;16 1718 19 20;21 22 23 24 25]，B=[3 0 16;17 -6 9;0 23 -4;9 7 0;4 1311]1）求它们的乘积C2）将矩阵C的右下角3x2子矩阵赋给D>> A=[1:1:5;6:1:10;11:1:15;16:1:20;21:1:25];B=[3 0 16;17 -6 9;0 23 -4;9 7 0;4 13 11];...C=A*B,D=C(3:5,2:3)C =93 150 77258 335 237423 520 397588 705 557753 890 717D =520 397705 557890 7172、完成下列操作1）求[100，999]之间能被61整除的数及其个数（提示：先利用冒号表达式，再利用find和length函数。
）>> A=[100:999];B=find(rem(A,61)==0),B_total=length(B)B =23 84 145 206 267 328 389 450 511 572633 694 755 816 877B_total =152）建立一个字符串向量，删除其中的大写字母（提示：利用find函数和空矩阵。
）>> a=['MA TLAB is important'], b=abs(a); c=find(b<=90 & b>=65) , a(c)=[],a =MA TLAB is importantc = 1 2 3 4 5 6a =is important⑶已知A=[23 10 -78 0；41 -45 65 5；32 5 0 32；6 -5492 14]，取出其前3行构成矩阵B，其前两列构成矩阵C，其左下角3x2子矩阵构成矩阵D，B与C的乘积构成矩阵E，分别求E<D、E&D、E|D、~E|~D。
(完整版)MATLAB)课后实验答案[1]
![(完整版)MATLAB)课后实验答案[1]](https://img.taocdn.com/s3/m/a428921bf524ccbff12184eb.png)
实验一 MATLAB 运算基础1.先求下列表达式的值，然后显示 MATLAB 工作空间的使用情 况并保存全部变量解:4.完成下列操作：(1) 求[100,999] 之间能被21整除的数的个数。
(1) z i 2sin 85° 1 e 2Z 2 卯(x L)，其中x2 0.45 2i 5Z 3 0.3a 0.3a e e 2sin (a 0.3)3.0, 2.9,L ,2.9, 3.0 Z 4t 2t 2 t 2 1 2t 其中 t=0:0.5:2.5(2) 建立一个字符串向量，删除其中的大写字母。
解: (1)结果：m=100:999;n=fin d(mod(m,21)==0);len gth( n)ans =43(2).建立一个字符串向量例如:ch二'ABC123d4e56Fg9：则要求结果是: ch二'ABC123d4e56Fg9：k=fi nd(ch>='A'&ch<='Z'); ch(k)=[] ch =123d4e56g9 实验二MATLAB 矩阵分析与处理1.设有分块矩阵A E 3 3 °2 3 R 3 2S 2 2 ，其中E 、R 、0、S 分别为单位矩阵、随机矩阵、零矩阵和对角阵，试通过数值计算验证A 2 E R RS 0 S 22 3 4 1 1 1 3 4 5 1 1 1 4 5 6 x , 0.95 x 2 0.67 x 3 0.52(1) 求方程的解。
(2) 将方程右边向量元素 b 3改为0.53再求解，并比较 b 3的变化和解的相对变化。
(3)计算系数矩阵A 的条件数并分析结论。
解：M 文件如下： 解：M 文件如下; 5.下面是一个线性方程组:实验三选择结构程序设计1. 求分段函数的值。
x2x 6 x 0且x 3y x2 5x 6 0 x 5且x 2 及x 3x2x 1 其他用if语句实现，分别输出x=-5.0,-3.0,1.0,2.0,2.5,3.0,5.0 时的y值。
MATLAB数学实验第二版课后练习题含答案
MATLAB数学实验第二版课后练习题含答案课后练习题MATLAB数学实验第二版的课后练习题如下：第一章课后练习题1.编写MATLAB程序，计算并输出下列公式的结果：y = \\frac{1}{\\sqrt{2\\pi\\sigma^2}} e^{-\\frac{(x-\\mu)^2}{2\\sigma^2}}其中，x, $\\mu$, $\\sigma$ 分别由用户输入。
要求输出结果精确至小数点后两位。
答案如下：x=input('请输入 x 的值:');mu=input('请输入 mu 的值:');sigma=input('请输入 sigma 的值:');y=1/sqrt(2*pi*sigma^2) *exp(-(x-mu)^2/ (2*sigma^2));fprintf('y = %.2f\', y);2.编写MATLAB程序，求解下列方程的解：4x + y = 11\\\\x + 2y = 7答案如下：A= [4,1;1,2];B= [11;7];X=inv(A) *B;fprintf('x = %.2f, y = %.2f\', X(1), X(2));第二章课后练习题1.编写MATLAB程序，计算下列多项式的值：P(x) = x^4 - 2x^3 + 3x^2 - x + 1其中，x 由用户输入。
要求输出结果精确至小数点后两位。
答案如下：x=input('请输入 x 的值:');y=x^4-2*x^3+3*x^2-x+1;fprintf('P(%.2f) = %.2f\', x, y);2.编写MATLAB程序，绘制下列函数的图像：f(x) = \\begin{cases} x + 1, & x < 0 \\\\ x^2, & 0 \\leq x < 1 \\\\ 2x - 1, & x \\geq 1 \\end{cases}答案如下：x=-2:0.01:2;y1=x+1;y2=x.^2.* ((x>=0) & (x<1));y3=2*x-1;plot(x,y1,x,y2,x,y3);legend('y1 = x + 1','y2 = x^2','y3 = 2x - 1');总结本文提供了《MATLAB数学实验第二版》的部分课后练习题及其答案。
matlab数学实验考试题及答案
matlab数学实验考试题及答案一、选择题（每题2分，共10分）1. MATLAB中用于生成0到1之间均匀分布的随机数的函数是？A. randB. randiC. randnD. randperm答案：A2. 下列哪个命令可以计算矩阵的行列式？A. detB. rankC. eigD. inv答案：A3. MATLAB中用于求解线性方程组的命令是？A. solveB. linsolveC. fsolveD. ode45答案：A4. 在MATLAB中，如何创建一个3x3的单位矩阵？A. eye(3)B. ones(3)C. zeros(3)D. identity(3)答案：A5. MATLAB中用于绘制二维图形的函数是？A. plotB. surfC. meshD. contour答案：A二、填空题（每题3分，共15分）1. MATLAB中，使用________函数可以计算矩阵的迹。
答案：trace2. 若要在MATLAB中创建一个从1到10的向量，可以使用________函数。
答案：1:103. MATLAB中，使用________函数可以计算矩阵的特征值。
答案：eig4. 若要在MATLAB中绘制一个正弦波，可以使用________函数。
答案：sin5. MATLAB中，使用________函数可以计算矩阵的逆。
答案：inv三、简答题（每题10分，共20分）1. 描述MATLAB中如何使用循环结构来计算并打印1到100之间所有奇数的和。
答案：可以使用for循环结构，初始化一个变量sum为0，然后遍历1到100之间的每个数，使用模运算符判断是否为奇数，如果是，则将其加到sum上，最后打印sum的值。
2. 简述MATLAB中如何使用条件语句来检查一个数是否为素数，并打印出所有小于100的素数。
答案：可以使用for循环遍历2到99之间的每个数，对于每个数，使用一个while循环检查它是否有除1和它本身之外的因数，如果没有，则使用if语句判断该数是否为素数，如果是，则打印该数。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
“M A T L A B ”练习题 要求：抄题、写出操作命令、运行结果，并根据要求，贴上运行图。
1、求230x e x -=的所有根。
（先画图后求解）（要求贴图）>> solve('exp(x)-3*x^2',0)ans =-2*lambertw(-1/6*3^(1/2))-2*lambertw(-1,-1/6*3^(1/2))-2*lambertw(1/6*3^(1/2))2、求下列方程的根。
1) 5510x x ++=a=solve('x^5+5*x+1',0);a=vpa(a,6)a =1.10447+1.05983*i-1.00450+1.06095*i-.-1.00450-1.06095*i1.10447-1.05983*i2）1sin02x x-=至少三个根>> fzero('x*sin(x)-1/2', 3)ans =2.9726>> fzero('x*sin(x)-1/2',-3)ans =-2.9726>> fzero('x*sin(x)-1/2',0)ans =-0.74083）2sin cos0x x x-=所有根>> fzero('sin(x)*cos(x)-x^2',0) ans =>> fzero('sin(x)*cos(x)-x^2',0.6)0.70223、求解下列各题：1）30sin limx x x x ->- >> sym x;>> limit((x-sin(x))/x^3)ans =1/62） (10)cos ,x y e x y =求>> sym x;>> diff(exp(x)*cos(x),10)ans =(-32)*exp(x)*sin(x)3）21/20(17x e dx⎰精确到位有效数字）>> sym x; >> vpa((int(exp(x^2),x,0,1/2)),17)0.4）42254x dx x+⎰ >> sym x;>> int(x^4/(25+x^2),x)ans =125*atan(x/5) - 25*x + x^3/35）求由参数方程arctan x y t⎧⎪=⎨=⎪⎩dy dx 与二阶导数22d y dx 。
>> sym t;>> x=log(sqrt(1+t^2));y=atan(t);>> diff(y,t)/diff(x,t)ans =1/t6）设函数y =f (x )由方程xy +e y = e 所确定，求y ′(x )。
>> syms x y;f=x*y+exp(y)-exp(1); >> -diff(f,x)/diff(f,y) ans =-y/(x + exp(y))7）sin2xe xdx+∞-⎰>> syms x;>> y=exp(-x)*sin(2*x);>> int(y,0,inf)ans =2/58）08x=展开（最高次幂为）>> syms xf=sqrt(1+x);taylor(f,0,9)ans =- (429*x^8)/32768 + (33*x^7)/2048 - (21*x^6)/1024 + (7*x^5)/256 - (5*x^4)/128 +x^3/16 - x^2/8 + x/2 + 19） 1sin (3)(2)x y e y =求>> syms x y;>> y=exp(sin(1/x));>> dy=subs(diff(y,3),x,2)dy =-0.582610）求变上限函数2x x⎰对变量x 的导数。
>> syms a t; >> diff(int(sqrt(a+t),t,x,x^2))Warning: Explicit integral could not be found.ans =2*x*(x^2 + a)^(1/2) - (a + x)^(1/2)4、求点(1,1,4)到直线L :31102x y z --==- 的距离 >> M0=[1,1,4];M1=[3,0,1];M0M1=M1-M0;v=[-1,0,2];d=norm(cross(M0M1,v))/norm(v)d =1.09545、已知22()2(),x f x μσ--=分别在下列条件下画出()f x 的图形：（要求贴图）(1)1,011σμ=时＝，－，，在同一坐标系里作图>> syms x;>> fplot('(1/sqrt(2*pi))*exp(-((x)^2)/2)',[-3,3],'r')>> hold on>> fplot('(1/sqrt(2*pi))*exp(-((x-1)^2)/2)',[-3,3],'y')>> hold on>> fplot('(1/sqrt(2*pi))*exp(-((x+1)^2)/2)',[-3,3],'g')>> hold off(2)0,124μσ＝时＝，，，在同一坐标系里作图。
>> syms x;fplot('(1/sqrt(2*pi))*exp(-((x)^2)/2)',[-3,3],'r')hold onfplot('(1/(sqrt(2*pi)*2))*exp(-((x)^2)/(2*2^2))',[-3,3],'y') hold onfplot('(1/(sqrt(2*pi)*4))*exp(-((x)^2)/(2*4^2))',[-3,3],'g') hold off6、画下列函数的图形：（要求贴图）（1）sin020cos02 4x u tty u tutz⎧⎪=≤≤⎪=⎨≤≤⎪⎪=⎩>> ezmesh('u*sin(t)','u*cos(t)','t/4',[0,20,0,2]) (2) sin()03,03z xy x y=≤≤≤≤>> x=0:0.1:3;y=x;[X Y]=meshgrid(x,y);Z=sin(X*Y);>> mesh(X,Y,Z)（3）sin(3cos)02 cos(3cos)02 sinx t uty t uuz uππ=+⎧≤≤⎪=+⎨≤≤⎪=⎩ezmesh('sin(t)*(3+cos(u))','cos(t)*(3+cos(u))','sin(u)',[0,2*pi,0,2*pi])7、 已知422134305,203153211A B -⎛⎫⎛⎫ ⎪ ⎪=-=-- ⎪ ⎪ ⎪ ⎪-⎝⎭⎝⎭，在MATLAB 命令窗口中建立A 、B 矩阵并对其进行以下操作：(1) 计算矩阵A 的行列式的值det()A>> A=[4,-2,2;-3,0,5;1,5,3];>> det(A)ans =-158(2) 分别计算下列各式：1122,*,.*,,,,T A B A B A B AB A B A A --- >> A=[4,-2,2;-3,0,5;1,5,3];B=[1,3,4;-2,0,-3;2,-1,1]; >> 2*A-Bans =7 -7 0-4 0 130 11 5>> A*Bans =12 10 247 -14 -7-3 0 -8>> A.*Bans =4 -6 86 0 -152 -5 3>> A*inv(B)ans =-0.0000 -0.0000 2.0000 -2.7143 -8.0000 -8.1429 2.4286 3.0000 2.2857 >> inv(A)*Bans =0.4873 0.4114 1.00000.3671 -0.4304 0.0000-0.1076 0.2468 0.0000>> A*Aans =24 2 4-7 31 9-8 13 36>> A'ans =4 -3 1-2 0 52 5 3>>8、在MATLAB中分别利用矩阵的初等变换及函数rank、函数inv求下列矩阵的秩：(1)16323540,11124A-⎛⎫⎪=-⎪⎪--⎝⎭求 rank(A)=?>> A=[1,-6,3,2;3,-5,4,0;-1,-11,2,4];>> rank(A)ans =3(2) 35011200,10201202B ⎛⎫ ⎪⎪= ⎪ ⎪⎝⎭求1B -。
>> B=[3,5,0,1;1,2,0,0;1,0,2,0;1,2,0,2]>> inv(B)ans =2.0000 -4.0000 -0.0000 -1.0000-1.0000 2.5000 0.0000 0.5000-1.0000 2.0000 0.5000 0.50000 -0.5000 0 0.50009、在MATLAB 中判断下列向量组是否线性相关，并找出向量组1(1132),T α=234(1113),(5289),(1317)T T T ααα=--=-=-中的一个最大线性无关组。
>> a1=[1 1 3 2]'a2=[-1 1 -1 3]'a3=[5 -2 8 9]'a4=[-1 3 1 7]'A= [a1, a2 ,a3 ,a4] ;[R jb]=rref(A) a1 =1132a2 =-11-13a3 =589a4 =-1317R =1.0000 0 0 1.0909 0 1.0000 0 1.7879 0 0 1.0000 -0.0606 0 0 0 0 jb =1 2 3>> A(:,jb)1 -1 51 1 -23 -1 82 3 910、在MATLAB中判断下列方程组解的情况，若有多个解，写出通解。
（1）123412341234123442020 3720 31260 x x x xx x x xx x x xx x x x-+-=⎧⎪--+=⎪⎨++-=⎪⎪--+=⎩一：>> A=[1,-1,4,2;1,-1,-1,2;3,1,7,-2;1,-3,-12,6]; >> rank(A)ans =3>> rref(A)ans =1 0 0 00 1 0 -20 0 1 00 0 0 0二：>> A=[1,-1,4,2;1,-1,-1,2;3,1,7,-2;1,-3,-12,6]; >> format ratn=4;RA=rank(A)RA =3>> if(RA==n)fprintf('%方程只有零解')elseb=null(A,'r')endb =0 21 >> syms kX=k*bX =2*kk(2)123123123123234245 38213 496x x xx x xx x xx x x++=⎧⎪-+=-⎪⎨+-=⎪⎪-+=-⎩>> A=[2 3 1;1 -2 4;3 8 -2;4 -1 9]; b=[4 -5 13 -6]';B=[A b];>> n=3;>> RA=rank(A)RA =2>> RB=rank(B)RB =2rref(B)ans =1 02 -1 0 1 -1 2 0 0 0 0 0 0 0 0 >> format ratif RA==RB&RA==n %判断有唯一解X=A\belseif RA==RB&RA<n %判断有无穷解X=A\b %求特解C=null(A,'r') %求AX=0的基础解系else X='equition no solve' %判断无解endWarning: Rank deficient, rank = 2, tol = 8.9702e-015. X =3/2-1/2C =-21111、求矩阵211020413A-⎛⎫⎪= ⎪⎪-⎝⎭的逆矩阵1A-及特征值和特征向量。