Matlab练习题
MATLAB习题及参考答案经典.doc
习题：1， 计算⎥⎦⎤⎢⎣⎡=572396a 与⎥⎦⎤⎢⎣⎡=864142b 的数组乘积。
2， 对于B AX =，如果⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=753467294A ，⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=282637B ，求解X 。
3， 已知：⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=987654321a ，分别计算a 的数组平方和矩阵平方，并观察其结果。
4， 角度[]604530=x ，求x 的正弦、余弦、正切和余切。
(应用sin,cos,tan.cot)5， 将矩阵⎥⎦⎤⎢⎣⎡=7524a 、⎥⎦⎤⎢⎣⎡=3817b 和⎥⎦⎤⎢⎣⎡=2695c 组合成两个新矩阵： （1）组合成一个4⨯3的矩阵，第一列为按列顺序排列的a 矩阵元素，第二列为按列顺序排列的b 矩阵元素，第三列为按列顺序排列的c 矩阵元素，即 ⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡237912685574（2）按照a 、b 、c 的列顺序组合成一个行矢量，即 []2965318772546， 将(x -6)(x -3)(x -8)展开为系数多项式的形式。
(应用poly,polyvalm)7， 求解多项式x 3-7x 2+2x +40的根。
(应用roots)8， 求解在x =8时多项式(x -1)(x -2) (x -3)(x -4)的值。
(应用poly,polyvalm)9， 计算多项式9514124234++--x x x x 的微分和积分。
(应用polyder,polyint ，poly2sym)10， 解方程组⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡66136221143092x 。
(应用x=a\b)11， 求欠定方程组⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡5865394742x 的最小范数解。
(应用pinv) 12， 矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=943457624a ，计算a 的行列式和逆矩阵。
(应用det,inv)13， y =sin(x )，x 从0到2π，∆x =0.02π，求y 的最大值、最小值、均值和标准差。
matlab考试题目及答案
matlab考试题目及答案1. 题目：编写一个MATLAB函数，实现计算并返回一个向量中所有元素的平方和。
答案：以下是一个简单的MATLAB函数，用于计算并返回一个向量中所有元素的平方和。
```matlabfunction sumOfSquares = calculateSumOfSquares(vector)sumOfSquares = sum(vector.^2);end```2. 题目：给定一个3x3的矩阵A，使用MATLAB编写代码，求出矩阵A 的转置。
答案：可以通过简单的转置操作来求得矩阵A的转置。
```matlabA = [1 2 3; 4 5 6; 7 8 9];A_transpose = A';```3. 题目：编写一个MATLAB脚本，实现对一个二维数组进行排序，并返回排序后的数组。
答案：以下是一个MATLAB脚本，用于对一个二维数组进行排序，并返回排序后的数组。
```matlabfunction sortedArray = sort2DArray(array)sortedArray = sort(array(:));end```4. 题目：给定一个向量x，使用MATLAB编写代码，计算并返回向量x的元素个数。
答案：可以通过内置函数`numel`来计算向量x的元素个数。
```matlabx = [1, 2, 3, 4, 5];numElements = numel(x);```5. 题目：编写一个MATLAB函数，实现计算并返回两个向量元素的点积。
答案：以下是一个简单的MATLAB函数，用于计算两个向量的点积。
```matlabfunction dotProduct = calculateDotProduct(vector1, vector2)dotProduct = dot(vector1, vector2);end```6. 题目：给定一个矩阵B，使用MATLAB编写代码，求出矩阵B的行列式。
matlab试题及答案
matlab试题及答案# MATLAB试题及答案一、选择题1. MATLAB的基本数据单位是：A. 矩阵B. 向量C. 标量D. 数组答案：A2. 下列哪个命令可以用来绘制函数图形？A. `plot`B. `graph`C. `draw`D. `chart`答案：A3. MATLAB中，以下哪个是正确的矩阵转置操作？A. `transpose(A)`B. `A'`C. `A^T`D. `flip(A)`答案：B二、简答题1. 简述MATLAB中矩阵的基本操作。
答案：在MATLAB中，矩阵是最基本的数据结构，可以进行加、减、乘、除等基本运算。
矩阵的创建可以使用方括号`[]`，例如`A = [1 2;3 4]`。
矩阵的转置使用单引号`'`，例如`A'`。
矩阵的求逆使用`inv`函数，例如`inv(A)`。
2. MATLAB中如何实现循环结构？答案：MATLAB中实现循环结构主要有两种方式：`for`循环和`while`循环。
`for`循环用于已知迭代次数的情况，例如：```matlabfor i = 1:5disp(i);end````while`循环用于迭代次数未知的情况，例如：```matlabi = 1;while i <= 5disp(i);i = i + 1;end```三、计算题1. 给定矩阵A和B，请计算它们的乘积C，并求C的行列式。
A = [1 2; 3 4]B = [5 6; 7 8]答案：首先计算矩阵乘积C：```matlabC = A * B;```然后计算C的行列式：```matlabdetC = det(C);```结果为：```matlabC = [19 22; 43 50]detC = -16```2. 编写一个MATLAB函数，计算并返回一个向量的范数。
答案：```matlabfunction norm_value = vector_norm(v)norm_value = norm(v);end```四、编程题1. 编写一个MATLAB脚本，实现以下功能：- 随机生成一个3x3的矩阵。
Matlab 考题题整理 带答案
MATLAB 考试试题(1)产生一个1x10的随机矩阵，大小位于（-5 5），并且按照从大到小的顺序排列好！（注：要程序和运行结果的截屏）答案：a=10*rand(1,10)-5;b=sort(a,'descend')1.请产生一个100*5的矩阵，矩阵的每一行都是[1 2 3 4 5]2. 已知变量：A=’ilovematlab’；B=’matlab’, 请找出：（A）B在A中的位置。
（B）把B放在A后面，形成C=‘ilovematlabmatlab’3. 请修改下面的程序，让他们没有for循环语句！A=[1 2 3; 4 5 6; 7 8 9];[r c]=size(A);for i=1:1:rfor j=1:1:cif (A(i,j)>8 | A(i,j)<2)A(i,j)=0;endendend4. 请把变量A=[1 2 3; 4 5 6; 7 8 9]写到文件里(output.xls)，写完后文件看起来是这样的1 2 3 4 5 6 7 8 95.试从Yahoo网站上获得微软公司股票的2008年9月的每日收盘价。
6.编写M文件，从Yahoo网站批量读取60000.SH至600005.SH在2008年9月份的每日收盘价（提示：使用字符串函数）。
7. 将金牛股份（000937）2005年12月14日至2006年1月10日的交易记录保存到Excel中，编写程序将数据读入MATLAB中，进一步将数据读入Access数据库文件。
8.已知资产每日回报率为0.0025，标准差为0.0208，资产现在价值为0.8亿，求5％水平下资产的10天在险价值(Var)。
9.a=[1 2 3 4 5],b=a(1)*a(5)+a(2)*a(4)+a(3)*a(3)+a(4)*a(2)+a(5)*a(1).试用MATLAB中最简单的方法计算b，注意最简单哦。
1、求下列联立方程的解3x+4y-7z-12w=45x-7y+4z+ 2w=-3x +8z- 5w=9-6x+5y-2z+10w=-8求系数矩阵的秩；求出方程组的解。
matlab开卷考试题及答案
matlab开卷考试题及答案MATLAB开卷考试题及答案一、选择题（每题2分，共20分）1. MATLAB中用于创建一个3x3的单位矩阵的命令是：A. `eye(3)`B. `ones(3)`C. `zeros(3)`D. `identity(3)`答案：A2. 下列哪个函数用于计算矩阵的逆？A. `inv()`B. `det()`C. `eig()`D. `trace()`答案：A3. MATLAB中，用于绘制二维图形的函数是：A. `plot()`B. `surf()`C. `mesh()`D. `bar()`答案：A4. 如何在MATLAB中计算向量的范数？A. `norm()`B. `abs()`C. `length()`D. `size()`答案：A5. MATLAB中，用于创建一个随机数矩阵的函数是：A. `rand()`B. `randi()`C. `randn()`D. `randperm()`答案：A6. 下列哪个命令用于在MATLAB中求解线性方程组？A. `solve()`B. `fsolve()`C. `fzero()`D. `linsolve()`答案：A7. MATLAB中，用于计算矩阵的行列式的函数是：A. `det()`B. `trace()`C. `eig()`D. `inv()`答案：A8. 如何在MATLAB中创建一个全1的3x3矩阵？A. `ones(3,3)`B. `eye(3,3)`C. `zeros(3,3)`D. `rand(3,3)`答案：A9. MATLAB中，用于计算矩阵的特征值的函数是：A. `eig()`B. `eigs()`C. `svd()`D. `qr()`答案：A10. 下列哪个函数用于计算矩阵的奇异值分解？A. `eig()`B. `svd()`C. `qr()`D. `lu()`答案：B二、填空题（每题3分，共30分）1. MATLAB中，用于创建一个3x3的零矩阵的命令是 `____`。
matlab考试题
matlab考试题一、选择题1. 在MATLAB中，下列哪个命令可以用于绘制三维曲面图？A. plot3B. surfC. meshD. contour32. 如何在MATLAB中创建一个名为“myArray”的数组，其元素从1到100？A. myArray = 1:100;B. myArray = array(1:100);C. array = {1, 2, 3, ..., 100};D. array = 1:100, name = 'myArray';3. 在MATLAB中，如何计算矩阵A和矩阵B的点积？A. A * BB. A.* BC. dot(A, B)D. A + B4. 若要在MATLAB中定义一个函数，使得输入参数x的平方和x的和相加，应该如何编写该函数？A. function y = x^2 + xB. function y = sum(x^2, x)C. function y = x^2 + sum(x)D. function y = x.^2 + x5. 在MATLAB中，如何清空当前工作空间的所有变量？A. clear allB. clearC. delete allD. clear workspace二、填空题1. 在MATLAB中，使用 ________ 命令可以导入外部数据文件。
2. 要创建一个从10到20的行向量，可以使用 MATLAB 代码________。
3. 在MATLAB中，使用 ________ 函数可以计算矩阵的行列式。
4. 若要在MATLAB中绘制一个单位圆，可以使用 ________ 命令。
5. 为了在MATLAB中保存工作，可以使用 ________ 命令。
三、简答题1. 请简述MATLAB的主要功能和应用领域。
2. 描述在MATLAB中如何进行矩阵运算。
3. 说明在MATLAB中创建和使用函数的基本步骤。
4. 讨论在MATLAB中进行数据分析和可视化的常用方法。
Matlab练习
Matlab 练习习题1． 设a=1.2,b=-4.6,c=8.0,e=-4.0,计算⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡+=d bc e a t ππ22arctan2． 设a=5.67,b=7.811,计算)lg(b a e ba ++3． 已知园半径为15，求其直径、周长和面积。
4． 已知三角形三边a=8.5,b=14.6,c=18.5,求三角形面积2/)(,))()(((c b a s c s b s a s s area ++=---=5． 下列命令执行后，L1、L2、L3、L4的值分别是多少？A=1:9；B=10-A ； L1=A==B ； L2=A<=5； L3=A>3&A<7;L4=find(A>3&A<7)； 习题1． 设矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=321212113A⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---=111012111B求：（1）2A+B ；（2）4A 2-3B 2；（3）AB ；（4）BA ；（5）AB-BA2.设三阶矩阵A 、B ，满足A -1BA=6A+BA⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡=71000410031A 求矩阵B2． 设（2E-C -1B ）A T =C -1,其中E 是4阶单位矩阵，A T 是4阶矩阵A 的转置，⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡---=1021000032231021B ⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=1021000002101021C求矩阵A3． 有一4阶魔方矩阵a ，找出矩阵中大于7的元素，并将它们重新排列成列向量b 。
4． 给定矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=054000031A ,删去整行和整列的0。
5． 设2阶矩阵A 、B 、X ，满足X-2A=B-X,⎥⎦⎤⎢⎣⎡--=2112A⎥⎦⎤⎢⎣⎡--=0220B 求矩阵X6． 求矩阵的主对角元素、逆矩阵、行列式的值、秩、特征值和特征向量⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----=163053064A ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--=1124111221B7． 分别用矩阵求逆、矩阵除法求方程组的解⎪⎩⎪⎨⎧=+--=+-=+-1074453932z y x z y x z y x 8． 已知多项式P 1(x)=3x+2,P 2(x)=5x 2-x ＋2，P 3(x)=x 2-0.5,求：（1） P(x)=P 1(x)+P 2(x)+P 3(x) （2） P(x)=P 1(x)*P 2(x)*P 3(x) （3） P(x)=0的全部根计算x i =0.2*i,i=0,1,2各点上的P(x i )。
(完整word版)Matlab考试题库+答案
填空题1、标点符号; %—用来表示该行为注释行。
可以使命令行不显示运算结果，2、x为0～4pi，步长为0.1pi的向量，使用命令 x=0:0.1*pi:4*pi创建。
3、输入矩阵A=，使用全下标方式用A(2,2)取出元素“-5”，使用单下标方式用A(5)取出元素“-5”。
4、符号表达式sin(2*a+t)+m中独立的符号变量为t。
5、M脚本文件和M函数文件的主要区别是M脚本文件没有函数定义和M函数文件有函数定义_______。
6. 设x是一维数组，x的倒数第3个元素表示为x(_end-2_)设y为二维数组，要删除y的第34行和48列，可使用命令y(34,:)=[] ；y(:,48)=[]；7. 将变量x以Ascii文本格式存储到文件fname.txt，应使用命令save _x ；8. 在while 表达式, 语句体, End 循环语句中，表达式的值非零时表示循环条件为真，语句体将被执行，否则跳出该循环语句；9.要从键盘读入一个字符串并赋值给变量x，且给出提示“Who is she?”，应使用命令x=input(‘Who is she?’,’s’)；10．设A=和B=和C=均为m*n矩阵，且存在于WorkSpace中，要产生矩阵D=，可用命令D=(A-C)/B.^C，计算可用命令det(inv(A’*B)11. 在MATLAB命令窗口中的“>>”标志为MATLAB的命令行提示符，“│”标志为输入提示符。
12.已知A=[1 2 3；4 5 0；7 8 9]；B=[1 0 3；1 5 0；0 1 2]；写出下列各指令运行的结果。
A+B；A.*B；A==B ；ans= [2,2,6;5,10,0;7,9,11]ans= [1,0,9;4,25,0;0,8,18]ans= [1,0,1;0,1,1;0,0,0]13.已知A是矩阵，求A的对角矩阵函数是 diag ；求A的下三角矩阵函数是 tril。
(完整版)matlab基础练习题(带答案)
Matlab 基础练习题常量、变量、表达式1、 MATLAB 中，下面哪些变量名是合法的？（ ）（A ）_num （B ）num_ （C ）num- （D ）-num 2、 在MA TLAB 中，要给出一个复数z 的模，应该使用（ ）函数。
（A ）mod(z) （B ）abs(z) （C ）double(z) （D ）angle(z) 3、 下面属于MATLAB 的预定义特殊变量的是？（ ）（A ）eps （B ）none （C ）zero （D ）exp4、 判断：在MA TLAB 的内存工作区中，存放一个英文字符 'a' 需要占用1个字节，存放一个中文字符‘啊’需要占用2个字节。
（ 错，都是2个字节 ）5、 判断：MA TLAB 中，i 和j （ 对 ）6、 判断：MA TLAB 中，pi 代表圆周率，它等于3.14。
（ 错，后面还有很多位小数 ）7、 在MA TLAB 中，若想计算的51)3.0sin(21+=πy 值，那么应该在MA TLAB 的指令窗中输入的MA TLAB 指令是__y1=2*sin(0.3*pi)/(1+sqrt(5))_。
8、 在MA TLAB 中，a = 1，b = i ，则a 占_8__个字节，b 占_16_个字节，c 占________字节。
9、 在MA TLAB 中，inf 的含义是__无穷大__，nan 的含义是__非数（结果不定）___。
数组1、 在MA TLAB 中，X 是一个一维数值数组，现在要把数组X 中的所有元素按原来次序的逆序排列输出，应该使用下面的（ ）指令。
（A ）X[end:1] （B ）X[end:-1:1] （C ）X (end:-1:1) （D ）X(end:1) 2、 在MA TLAB 中，A 是一个字二维数组，要获取A 的行数和列数，应该使用的MATLAB的命令是（ ）。
（A ）class(A) （B ）sizeof(A) （C ）size(A) （D ）isa(A)3、 在MATLAB 中，用指令x=1:9生成数组x 。
matlab练习
一、填空1、命令clear 、clf 、clc 、who 、whos 的含义分别是 。
2、若矩阵2310,5112A B ⎛⎫⎛⎫== ⎪ ⎪⎝⎭⎝⎭,则A.*B= ,A*B= 。
3、生成4行3列的元素全为1的矩阵的命令是 ；生成4阶全零矩阵的命令是 。
生成对角线元素为[2 4 6]矩阵的命令是 。
4、假定A 是一个6阶方阵，选取矩阵A 第3行的指令是 ，选取第4行第2列元素的指令是 。
5、生成一个从1到50，步长为10的等差数列构成的数组，可以使用的命令是 。
6、求x 的平方根使用的命令是 ，round(pi)的结果是7、命令subplot(2,2,2)的功能是把图形窗口分为 个子图，并把第 个子图作为当前图形窗口。
8、matlab 中，求解线性方程组时，矩阵行变换化简的命令是 。
9、在用命令 p=polyfit （x ，y ，n ）对数据进行多项式拟合时，参数n 的含义为 。
10、求函数4y x cos7x =的5阶导数的命令是 。
11、用符号法求解微分方程2xy -3y =x ,y(1)=0,y(5)=0'''的解的指令是 。
12、Matlab 可以输入字母、汉字，但是M 文件中标点符号必须（） 状态下输入。
14、求xe 的命令是（），求x 的自然对数lnx 的命令（）。
15、画图时，在x 轴旁边加注文字说明的命令是( )，图名的标注命令是（），图例标注命令是（），在鼠标指定位置上标注的命令是（），将一个图形窗口分成多个子图的命令是（），画空间曲线的命令是（）。
16、绘制三维空间曲面图形的命令是（）和（），生成格点矩阵的命令是( )17、用矩形法、复合梯形公式、复合辛普生公式求1204dx 1x +⎰的定积分。
18、下列叙述正确的是：（ ）(A) matlab 中的变量是不区分大小写的。
(B) 在matlab 所输入的命令后如果输入逗号，则不显示执行结果。
(C) matlab 的变量必须先定义其类型才可使用。
matlab试题及答案开卷
matlab试题及答案开卷1. MATLAB基础操作在MATLAB中，如何创建一个名为“myMatrix”的3x3矩阵，其元素从1到9按顺序排列？答案：在MATLAB中，可以使用以下命令创建名为“myMatrix”的3x3矩阵：```matlabmyMatrix = [1 2 3; 4 5 6; 7 8 9];```或者使用内置函数`magic`：```matlabmyMatrix = magic(3);```2. 矩阵运算给定两个矩阵A和B，其中A = [2 3; 4 5]，B = [6 7; 8 9]，计算A+B的结果。
答案：在MATLAB中，可以使用加号`+`来计算两个矩阵的和：```matlabA = [2 3; 4 5];B = [6 7; 8 9];C = A + B;```计算结果C为：```matlabC = [8 10; 12 14];```3. 向量操作创建一个从0到1的等差数列，步长为0.1。
答案：在MATLAB中，可以使用`linspace`函数或`:`操作符来创建等差数列：```matlabv = 0:0.1:1;```或者使用`linspace`函数：```matlabv = linspace(0, 1, 11);```两种方法都将得到一个包含11个元素的向量，从0开始，到1结束，步长为0.1。
4. 条件语句编写一个MATLAB脚本，判断一个给定的数n是否为素数。
答案：可以通过以下MATLAB脚本来判断一个数n是否为素数：```matlabfunction isPrime = isPrimeNumber(n)if n <= 1isPrime = false;elsefor i = 2:sqrt(n)if mod(n, i) == 0isPrime = false;return;endendisPrime = true;endend```调用此函数时，传入一个数值n，函数将返回一个布尔值，表示n 是否为素数。
matlab试题及答案
matlab试题及答案一、单项选择题（每题2分，共10分）1. MATLAB中用于创建向量的命令是：A. vectorB. arrayC. linspaceD. colon答案：D2. 在MATLAB中，以下哪个函数用于计算矩阵的行列式？A. detB. rankC. invD. eig答案：A3. MATLAB中用于进行矩阵转置的运算符是：A. 'B. .C. ^D. !答案：A4. 若A是一个3x3的矩阵，执行命令A(2,:)=0;后，矩阵A的第二行将变为：A. [0 0 0]B. [1 0 0]C. [0 1 0]D. [0 0 1]答案：A5. MATLAB中，以下哪个函数用于绘制三维曲面图？A. plotB. surfC. barD. hist答案：B二、填空题（每题3分，共15分）1. MATLAB中，使用________命令可以清除所有变量。
答案：clear2. 若要在MATLAB中创建一个从0到1的100个元素的向量，可以使用命令________。
答案：linspace(0,1,100)3. MATLAB中，使用________函数可以计算矩阵的特征值。
答案：eig4. 在MATLAB中，________函数用于计算两个矩阵的点乘。
答案：dot5. 若要在MATLAB中绘制一个圆，可以使用________函数。
答案：plot三、简答题（每题5分，共20分）1. 请解释MATLAB中矩阵索引的概念。
答案：在MATLAB中，矩阵索引指的是通过行号和列号来访问矩阵中特定元素的过程。
例如，A(2,3)表示访问矩阵A的第二行第三列的元素。
2. MATLAB中如何实现矩阵的元素乘法？答案：在MATLAB中，矩阵的元素乘法可以通过使用点乘运算符（.*）来实现。
例如，C = A .* B，其中A和B是相同大小的矩阵。
3. 请说明MATLAB中如何使用循环结构。
答案：MATLAB中可以使用for循环和while循环两种循环结构。
MATLAB平时练习题(附答案)
一、填空题1、MATLAB常用操作界面包括、工作空间窗口、、、内存数组编辑器、M文件编辑/调试器、帮助导航/浏览器、图形窗口等。
2、MATLAB中Inf或inf表示、NaN或nan表示、nargout表示。
3、MATLAB中逗号主要用作；用作输入量与输入量之间的分隔符；。
4、工作空间浏览器主要用于内存变量的、和。
5、MATLAB实现将全下标转换为单下标的指令为、据单下标换算出全下标的指令为。
6、二维数组的标识有、、“逻辑1”标识。
二、简述下列命令在MATLAB中的作用1.clf %清除图对象clear %清除工作空间内的所有变量clc %清除当前屏幕上显示的所有内容，但不清除工作空间中的数据2.ceil %沿+∞方向取整factor %符号计算的因式分解3.sph2cart %球坐标变为直角坐标pow2 %2的幂4.digits(d) %设置当前的可变算术精度的位数为整数 d 位vpa(pi,100) %设置当前pi的小数位数为1005.box on %打开框状坐标轴开grid off %网格关6.logspace %对数分度向量cat %串接成高维数组7.sym2poly %符号多项式转变为双精度多项式系数向量poly2sym %双精度多项式系数转变为向量符号多项式8.plot3 %三维线图poly2str %以习惯方式显示多项式9.bar %二维直方图pie %二维饼图10.zoom on %打开图形缩放模式 edit %M文件编辑11.whos %对当前工作空间变量的信息进行列表figure %生成图形窗口12.cart2sph %直角坐标变为球坐标pol2cart % 极或柱坐标变为直角坐标13.diff %数值差分、符号微分14.ezplot3 %画三维曲线的简捷指令fix %向零取整factor %符号计算的因式分解三简答题1.叙述MATLAB的主要窗口包括哪些？其作用是什么？命令窗口（Commend Window）——供输入命令，实现计算或绘图功能。
Matlab考试题库及答案(教师出卷参考专用)
Matlab考试题库及答案(教师出卷参考专用)一、选择题1.以下哪个函数用于在Matlab中创建一个图形窗口？A. figureB. plotC. graphD. window答案：A2.在Matlab中，以下哪个选项可以用来定义一个矩阵？A. A = [1 2 3; 4 5 6]B. A = (1, 2, 3, 4, 5, 6)C. A = {1, 2, 3, 4, 5, 6}D. A = 1 2 3; 4 5 6答案：A3.以下哪个函数用于求解线性方程组Ax=b？A. solveB. linsolveC. solve(A, b)D. linsolve(A, b)答案：D4.在Matlab中，如何计算矩阵A和矩阵B的乘积？A. A BB. A \ BC. A . BD. A .\ B答案：A5.以下哪个函数用于在Matlab中绘制三维散点图？A. scatterB. scatter3C. plot3D. bar3答案：B二、填空题1.在Matlab中，要创建一个名为"myfig"的图形窗口，可以使用______函数。
答案：figure('Name', 'myfig')2.在Matlab中，要计算矩阵A的行列式，可以使用______函数。
答案：det(A)3.在Matlab中，若要计算变量x的平方，可以使用______运算符。
答案：.^24.在Matlab中，若要计算矩阵A的逆矩阵，可以使用______函数。
答案：inv(A)5.在Matlab中，要绘制一个正弦波形，可以使用______函数。
答案：plot(sin(x))三、判断题1.在Matlab中，矩阵的索引从1开始计数。
（）答案：正确2.在Matlab中，可以使用逻辑运算符"&&"和"||"。
（）答案：错误3.在Matlab中，矩阵乘法满足交换律。
MATLAB习题及答案
填空题1. MATLAB于1984年由美国Mathworks公司推出，其后每年更新（两次。
2. MATLAB是一种以（矩阵）运算为基础的交互式程序设计语言。
3. MATLAB具有卓越的数值计算能力和符号计算、文字处理、可视化建模仿真和实时控制等众多功能，其每个变量代表一个（矩阵），每个元素都看作（复数）。
4.通过命令（help）、（lookfor），可以查找所有命令或函数的使用方法。
5.执行语句a=1:2:10，得到的一维数组是（1 3 5 7 9）.6.执行语句b=linspace（1,10,10）后，一维数组b包含（10）个元素，最大值是10）7.函数rem（）的功能是取（余）数。
8.若p=[1 0 0;1 1 0],则p|〜p=（[1 1 1;1 1 1]）.（注：填空时请用本题的p的方式表示结果）9.若p=[1 0 0;1 1 0],则all（p）=（[1 0 0]）.10.矩阵的加减运算，要求相加减的矩阵阶数相同。
若A=[1 2 3 4；2 3 1 8]，则执行语句：[n,m]=size（A），则n=（2 ）,m=（4 ）.11.对于一维矩阵，求其长度的函数是（length（））.12.数组和数组之间的运算，尤其是对于乘除运算和乘方运算，如果采用点方式进行计算，表明是数组的（元素）之间的运算关系。
13.求矩阵运算A*B时，要求在维度上，A的（列）数与B的（行）数相等。
二、判断题1.MATLAB只有一种数据类型，一种标准的输入输出语句，不需编译，可直接运行。
（对2.MATLAB的特殊常量是一些预选定义好的数值变量。
（对3.MATLAB变量名不区分大小写。
（错4.i是特殊常量。
（对5.NAN是非数。
（对6.MATLAB中所有的变量都表示一个矩阵或一个向量。
（对7.MATLAB中变量不需要先定义后使用，会自动根据实际赋值的类型对变量类型进行定义。
（对8.clc命令可以从内存中删除一个、多个和所有变量。
MATLAB考试试题及答案
MATLAB考试试题及答案一、选择题（每题5分，共25分）1. 在MATLAB中，下列哪个命令用于创建一个行向量？A. v = [1; 2; 3]B. v = [1 2 3]C. v = [1, 2, 3]D. v = (1, 2, 3)答案：B2. 在MATLAB中，下列哪个命令用于计算矩阵A的行列式？A. det(A)B. det(A')C. det(inv(A))D. det(A^2)答案：A3. 在MATLAB中，下列哪个命令用于计算矩阵A的逆？A. inv(A)B. A^(-1)C. pinv(A)D. A\B答案：A4. 在MATLAB中，下列哪个命令用于求解线性方程组Ax= b？A. A\bB. A/BC. B/AD. A^-1b答案：A5. 在MATLAB中，下列哪个命令用于绘制二维图形？A. plot(x, y)B. scatter(x, y)C. bar(x, y)D. pie(x, y)答案：A二、填空题（每题5分，共25分）6. 在MATLAB中，可以使用______命令创建一个等差数列。
答案：linspace7. 在MATLAB中，可以使用______命令创建一个等比数列。
答案：logspace8. 在MATLAB中，可以使用______命令计算矩阵A的特征值。
答案：eig(A)9. 在MATLAB中，可以使用______命令计算矩阵A的特征向量。
答案：eigenvector(A)10. 在MATLAB中，可以使用______命令计算矩阵A的奇异值。
答案：svd(A)三、解答题（每题25分，共75分）11. 编写MATLAB程序，求解以下线性方程组：2x + 3y - z = 1x - y + 2z = 03x + 2y - 4z = -3答案：```A = [2 3 -1; 1 -1 2; 3 2 -4];b = [1; 0; -3];x = A\b;disp('解为：');disp(x);```12. 编写MATLAB程序，绘制以下函数的图形：y = sin(x) + cos(x)，x ∈ [0, 2π]答案：```x = linspace(0, 2pi, 100);y = sin(x) + cos(x);plot(x, y);title('y = sin(x) + cos(x)');xlabel('x');ylabel('y');grid on;```13. 编写MATLAB程序，计算以下矩阵的特征值和特征向量：A = [1 2 3; 4 5 6; 7 8 9]答案：```A = [1 2 3; 4 5 6; 7 8 9];[V, D] = eig(A);disp('特征值：');disp(diag(D));disp('特征向量：');disp(V);```14. 编写MATLAB程序，使用牛顿迭代法求解方程f(x) = x^3 - 4x + 2 = 0在x = 1附近的根。
matlab考试试题
matlab考试试题一、选择题（每题2分，共20分）1. MATLAB的全称是：A. Matrix LaboratoryB. Microprocessor Application Tool SetC. Micro Advanced Tool SetD. None of the above2. 在MATLAB中，用于生成1到10的行向量的命令是：A. `1:10`B. `linspace(1,10,1)`C. `logspace(1,10,1)`D. `1.10`3. MATLAB中的基本数据类型是：A. 向量B. 矩阵C. 字符串D. 单元数组4. 下列哪个命令可以计算矩阵的行列式？A. `det`B. `diag`C. `dot`D. `diff`5. 在MATLAB中，用于绘制三维曲线的函数是：A. `plot`C. `surf`D. `plot3`6. MATLAB中，创建一个空矩阵的命令是：A. `zeros(3)`B. `ones(3)`C. `eye(3)`D. `rand(3)`7. 下列哪项不是MATLAB的编程范式？A. 函数式编程B. 面向对象编程C. 过程式编程D. 事件驱动编程8. MATLAB中，用于求解线性方程组的函数是：A. `solve`B. `fsolve`C. `linsolve`D. `linspace`9. 在MATLAB中，用于计算向量范数的函数是：A. `norm`B. `abs`C. `sum`D. `mean`10. MATLAB中，用于生成随机数的函数是：A. `rand`B. `randn`D. All of the above二、简答题（每题5分，共30分）11. 简述MATLAB中矩阵的基本操作，至少列出三种。
12. 解释MATLAB中的M文件和脚本文件的区别。
13. 在MATLAB中，如何实现循环结构？请举例说明。
14. 描述MATLAB中函数的作用域，包括局部变量和全局变量。
matlab笔试试题及答案
matlab笔试试题及答案MATLAB笔试试题及答案1. MATLAB中，如何创建一个名为`myMatrix`的3x3矩阵，其元素从1到9依次排列？答案：可以使用以下命令创建矩阵：```matlabmyMatrix = [1 2 3; 4 5 6; 7 8 9];```2. 在MATLAB中，如何计算向量`v`的范数？答案：可以使用`norm`函数来计算向量的范数。
例如，如果`v`是一个向量，则：```matlabv = [1, 2, 3];norm_v = norm(v);```3. MATLAB中，如何将一个矩阵`A`转置？答案：可以使用单引号(`'`)来转置矩阵`A`：```matlabA = [1 2 3; 4 5 6; 7 8 9];A_transpose = A';```4. 如果要找到矩阵`B`中的最大元素，应该如何操作？答案：可以使用`max`函数来找到矩阵中的最大元素：```matlabB = [10 20 30; 40 50 60; 70 80 90];max_element = max(B(:));```5. 在MATLAB中，如何实现一个简单的for循环，计算从1到10的和？答案：可以使用以下代码实现：```matlabsum = 0;for i = 1:10sum = sum + i;end```6. MATLAB中，如何创建一个二维数组，其元素为当前索引的乘积？答案：可以使用以下代码创建：```matlabC = zeros(3,3);for i = 1:3for j = 1:3C(i,j) = i * j;endend```7. 如何在MATLAB中使用条件语句？答案：可以使用`if`语句来实现条件判断。
例如：```matlabx = 5;if x > 0disp('x is positive');elsedisp('x is not positive');end```8. 在MATLAB中，如何使用`plot`函数绘制一个正弦波？答案：可以使用以下代码绘制正弦波：```matlabx = 0:0.1:2*pi;y = sin(x);plot(x, y);title('Sine Wave');xlabel('x');ylabel('sin(x)');```9. MATLAB中如何使用`fft`函数计算一个信号的快速傅里叶变换？答案：可以使用以下代码计算信号的FFT：```matlabt = 0:0.001:1;x = cos(2*pi*50*t) + randn(size(t));X = fft(x);f = linspace(0, 1000, length(x));plot(f, abs(X));title('Frequency Domain');xlabel('Frequency (Hz)');ylabel('Amplitude');```10. 如何在MATLAB中使用`subplot`函数创建多个子图？答案：可以使用以下代码创建多个子图：```matlabsubplot(2,1,1); % 第一个子图plot(1:10, (1:10).^2);title('Square of Numbers');subplot(2,1,2); % 第二个子图plot(1:10, (1:10).^3);title('Cube of Numbers');```。
matlab考试题B卷及答案
matlab考试题B卷及答案1. MATLAB基础操作题：请在MATLAB中创建一个名为`myVector`的向量，包含元素1, 2, 3, 4, 5，并计算其元素之和。
答案：首先，创建向量`myVector`可以使用以下代码：```matlabmyVector = [1 2 3 4 5];```然后，计算元素之和可以使用`sum`函数：```matlabsumOfElements = sum(myVector);```最终，`sumOfElements`的值将为15。
2. MATLAB矩阵运算题：给定一个3x3的矩阵A，其元素为A=[1 2 3;4 5 6; 7 8 9]，请计算矩阵A的转置。
答案：计算矩阵A的转置可以使用`transpose`函数或者简单的`.'`操作符。
以下是使用`transpose`函数的代码：```matlabA = [1 2 3; 4 5 6; 7 8 9];A_transpose = transpose(A);```或者使用`.'`操作符：```matlabA_transpose = A.';```得到的转置矩阵`A_transpose`将是：```1 4 72 5 83 6 9```3. MATLAB编程题：编写一个MATLAB函数，该函数接受一个整数n作为输入，并返回一个n阶的魔方矩阵。
答案：魔方矩阵是一个方阵，其中每一行、每一列以及两条对角线上的元素之和都相等。
以下是一个生成n阶魔方矩阵的MATLAB函数：```matlabfunction magicMatrix = createMagicSquare(n)if mod(n, 2) == 0error('n must be an odd number.');endmagicMatrix = zeros(n);num = 1;for i = 1:nfor j = 1:nif i == 1j = j;elsej = (j + 2 * (i - 1) - n) % n + 1;endmagicMatrix(i, j) = num;num = num + 1;endendend```调用此函数并传入一个奇数n，例如`createMagicSquare(3)`，将返回一个3阶魔方矩阵。
matlab基础考试题及答案
matlab基础考试题及答案1. MATLAB中，如何创建一个名为"myMatrix"的3x3矩阵，其元素为1到9的连续整数？答案：可以使用以下命令创建矩阵：```matlabmyMatrix = [1 2 3; 4 5 6; 7 8 9];```或者使用`magic`函数：```matlabmyMatrix = magic(3);```2. 在MATLAB中，如何计算一个向量的元素平方？答案：可以使用`.^`运算符来计算向量元素的平方：```matlabv = [1 2 3 4];v_squared = v.^2;```3. MATLAB中，如何将一个向量逆序？答案：可以使用`flip`函数或者`end`索引来逆序一个向量：```matlabv = [1 2 3 4];v_reversed = flip(v);```或者```matlabv_reversed = v(end:-1:1);```4. 在MATLAB中，如何计算矩阵的行列式？答案：可以使用`det`函数来计算矩阵的行列式：```matlabA = [1 2; 3 4];detA = det(A);```5. MATLAB中，如何使用循环结构来计算1到100的自然数之和？答案：可以使用`for`循环来计算总和：```matlabsum = 0;for i = 1:100sum = sum + i;end```或者使用`sum`函数直接计算：```matlabsum = sum(1:100);```6. 在MATLAB中，如何创建一个复数？答案：可以通过指定实部和虚部来创建一个复数：```matlabcomplexNumber = 3 + 4i;```7. MATLAB中，如何绘制一个正弦波？答案：可以使用`plot`函数和`sin`函数来绘制正弦波：```matlabx = linspace(0, 2*pi, 100);y = sin(x);plot(x, y);```8. 如何在MATLAB中计算两个矩阵的点积？答案：可以使用`dot`函数或`*`运算符来计算两个矩阵的点积：```matlabA = [1 2; 3 4];B = [5 6; 7 8];dotProduct = dot(A, B);```或者```matlabdotProduct = A * B;```9. MATLAB中，如何使用条件语句来检查一个数是否为偶数？答案：可以使用`if`语句和模运算符来检查一个数是否为偶数：```matlabnum = 10;if mod(num, 2) == 0disp('The number is even.');elsedisp('The number is odd.');end```10. 在MATLAB中，如何将一个字符串转换为小写？答案：可以使用`lower`函数来将字符串转换为小写：```matlabstr = 'Hello World';lowerStr = lower(str);```结束语：以上是MATLAB基础考试题及答案，希望这些题目能帮助你更好地理解和掌握MATLAB的基础知识。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
Matlab练习题１已知矩阵11 12 13 1421 22 23 24A= 31 32 33 3441 42 43 44(1)A(:,1) (2) A(2,:) (3) A(:,2:3)(4) A(2:3,2:3) (5) A(:,1:2:3) (6) A(2:3)(7) A(:) (8) A(:,:) (9) ones(2,2)(10) eye(2) (11) [A,[ones(2,2);eye(2)]](12) diag(A) (13) diag(A,1) (14) diag(A,-1)(15) diag(A,2)2(1)输入如下矩阵Ａ０π/3Ａ＝π/6 π/2(2) 求矩阵B1,B1中每一元素为对应矩阵Ａ中每一元素的正弦函数(3) 求矩阵B2, B2中每一元素为对应矩阵Ａ中每一元素的余弦函数(4) 求B12+B22(5) 求矩阵Ａ的特征值与特征矢量：称特征矢量为Ｍ，而特征值矩阵为Ｌ(6) 求Msin(L)M-1(7)使用funm命令求矩阵A的正弦函数（结果应该与（６）同）(8)求cosA(9) 证明sin2A+cos2A=I3 按题目要求用MATLAB命令完成下列矩阵运算(1) 使用rand命令产生5个2x2随机矩阵A,B,C,D,E(2) 求矩阵F(使用和不使用inv命令两种情况) F=A-1[B+C-1(D-1E)]4 手算和上机分别求A.*B‘ A.\B‘ A‘.\B其中A=[1;1;1] B=[2,3,4]5 已知A=[2 7 6;9 0 –10;3 0.5 6]; B=[8 0.2 0;3 2 5;4 0 7];求(1)A|B, A&B, A~B, A>B, A>=B, A<B, A<=B, A==B, A~=B(2)元素值为零的元素标号(3)元素值大于6的元素标号6 某专业有三名研究生,本学期选修了四门课程,若这些研究生的姓名,学号,性别,出生年月,课程名称,考试成绩可任意假定,(1)分别用结构型变量和细胞型变量表示以上信息;(2)举例说明查阅以上任何一条信息的方法;(3)求每一个研究生的平均成绩.？7 已知矩阵A=[1.2 3 5 0.9;5 1.7 5 6;3 9 0 1;1 2 3 4],求(1)A的特征多项式(2)特征多项式中未知数为20 时的值？(3)特征多项式的根？(4)特征多项式的导数8 已知五个数据点: (1,5.5), (2,43.1), (3,128), (4,290.7), (5,498.4)(1)用三次曲线拟合上述数据点(2)在同一图形中绘出数据点和拟合曲线(3)用适当图形表示拟合精度9 在实验中测得如下6组数据:(0,1.5), (1.5,3.4), (2.8,13), (3.8,36), (4.5,63), (4.9,78)请用三次曲线拟合以上数据并给出以下结果：(1)三次多项式的各项系数；(2)将数据点和拟合曲线以最佳效果在同一图中绘出；10 求下列函数的极限（1）lim(x2/sin2(x/3)) x 0 （2）lim((tanx-sinx)/sin3x) x 0 （3）lim(sin(a+x)-sin(a-x))/x x 0 （4）limxcos(1/x) x 0（5）lim((1+mx)n-(1+nx)m)/x2 x 0 （6）lim(1+1/n)(n+5) n ∞（7）lim(1-2/x)x x ∞（8）lim(1+cosx) 3secx x л/2 11 求下列函数的积分（1）x2/sin2(x/3) （2）(tanx-sinx)/sin3x（3）(sin(a+x)-sin(a-x))/x （4）xcos(1/x)（5）((1+mx)n-(1+nx)m)/x2 (6) cos2x(7) sinaxcosbx (8)cosaxsinbx？(9) arcsin(x/a) (10)1/(a+bsinx)(11) xarcsin(x/a) (12) x2arcsin(x/a)12 求下列函数的反函数（1）y=2sin3x （2）y=1+ln(x+2)（3）y=2x/(2x+1) （4）y=(ex+e(-x))/2（5）y=1/2(arccos(x/2)) （6）y=x+1/x13 求下列函数的定积分（1）(x+sinx)/(1+cosx) [0,pi/2] （2）ln(1+tanx ) [0,pi/4] （3）1/(1+cos2x ) [0,pi/2] （4）cos5xsinx [0,pi/2] （5）(3x4+3x2+1)/(x2+1) [-1,0] (6) x2+1/x4 [1,2] (7) tan2 x [0,pi/4] (8) 4cos4x [-pi/2,pi/2] (9) 1-sin3x [0,pi] (10)1/(11+5x) 3 [-2,1] (11) cosxcos2x [-pi/2 , pi/2] (12) (x 3 sin 2 x )/(x 4 + 2x 2 + 1) [-5,5] 14 求解代数方程？(1) ax2+bx+c=0(2) cos(2x)+sin(x)=115 解线性方程组a1x+b1y+c1z=d1a2x+b2y+c2z=d2 其中ai, bi, ci, di为常数, x, y, z为变量a3x+b3y+d3z=d316 求方程tan(x)+sin(x)=2在区间[-2π, 2π]上的全部实数解并用图示。
17 求解以下微分方程.(1) dy/dx=y+sinx y |x= π/2=1(2) d2y/dx2+bdy/dx+a2y=0 dy/dx|x= π/a=0 y|x=0=1(3) dx/dt=2x+y dy/dt=-x+3y x(0)=1 y(0)=2(4) dy/dx=y+cosx y |x= π/2=1(5) d2y/dx2+bdy/dx+a2y=0 dy/dx|x= π/a=1 y|x=0=2(6) dx/dt=3x+y dy/dt=-2x+5y x(0)=1 y(0)=218 求解以下微分方程(用符号解法）.(1) dy/dx=y/x+tany/x(2) (3x2+6xy2) dx +(6x2y+4y3)dy=0(3) dx/dt=-xcost+(1/2 )sin2t(4)y= (dy/dx) 2 -x (dy/dx)+x 2 /2(5) dy/dx=(x-y+1)/(x+y-3)19 绘制下列各种函数图形(1)绘制下列极坐标图形r=3(1-cosθ)r=2(1+cos θ)r=2(1+sin θ)r=cos3 θr=exp(4πθ)(2)求函数z的三维图形.定义区间与z的函数表达式如下:-5 ≤ x ≤ 5 , -5 ≤ y ≤51 1z= —————————————————(x+1)2+(y+1)2+1 (x-1)2+(y-1)2+120 用MATLAB命令绘制下列数学函数的图形(1) 绘制如下函数图形y(t)=1-2exp(-t)sin(t) (0 ≤ t ≤ 8)且在x轴写上―Time‖,y轴写上―Amplitude‖标号,图形标题为―Decaying-oscillating Exponential‖(2) 绘制如下图形y(t)=5exp(-0.2t)cos(0.9t-30O)+0.8exp(-2t) (0≤ t ≤30)(3)在0 ≤ t≤10区间内绘制下面图形y(t)=1.23cos(2.83t+240O)+0.625x(t)=0.625t(4)在0 ≤ t≤20区间内,且在同一图中绘制下述函数图形y1(t)=2.62exp(-0.25t)cos(2.22t+174O)+0.6y2(t)=2.62exp(-0.25t)+0.6y3(t)=0.6在y值对应的-2到3的区域,首先求y1的最小值与最大值,然后再求次大值与次小值21 用MATLAB命令构造一个101个元素的矢量,它的元素交替为1和-1,使用plot, comb 和stairs 命令分别绘制该矢量的图形22(1)利用prod命令生成与列写N的阶乘(N!),并对N=1~10进行验证(2)使用for命令代替prod 命令,完成(1)项中的功能(3)使用sum命令求N个数的连续和.并对N=1~10进行验证(4)使用for命令代替sum命令,完成(3)项中的功能23 绘制分段连续函数图x2 |x|≦5-5x+50 5﹤x﹤105x+50 -10﹤x﹤-50 |x|≧1024 求以下定积分：25 请绘出以下函数的三维面图：f(x,y)=3(1-x)2exp(-x2-(y+1)2)-10(x/5-x3-y5)exp(-x2-y2)-(1/3)exp(-(x+1)2-y2)26 请练习以下命令plot semilogx semilogy loglogpolar plotyy ezplot ezpolarfplot feather pie barcontour plot3 mesh surf27 已知微分方程dy/dx=(x+y)2 ,初始值为y(0)=1, 在x的取值范围[0, 1]内求方程的数值解并绘出函数图形。
28 已知sin0.32=0.314567,sin0.34=0.333487,sin0.36=0.352274, 用线性插值、抛物线插值和三次样条插值计算sin0.3367的值并作出误差。
29给定函数f(x)= , x [-5,5]，节点xk=-5+k,用三次样条插值求s(xk).30 给出f(x)=lnx数值表用线性插值及二次插值计算ln0.54的近似值。
X 0.4 0.5 0.6 0.7 0.8lnx -0.916291 -0.693147 -0.510826 -0.357765 -0.22314431 给出数据表如下X 0.25 0.30 0.39 0.45 0.53Y 0.5000 0.5477 0.6245 0.6708 0.7280试求三次样条插值s(x)，并满足条件：S’(0.25)=1.000,s’(0.53)=0.6868;s’’(0.25)=s’’(0.53)=01 求矩阵a =2 1 –3 –13 1 0 7–1 2 4 –21 0 –1 5的逆阵、对角线及其行列式的值。
(50)a=[2 1 -3 -1;3 1 0 7;-1 2 4 -2;1 0 -1 5]；inv(a)det(a)2求下列线性方程的解。
A=[6 3 4；-2 5 7；8 –4 –3]；B=[3；-4；-7]；X=A\B3 求垂直于向量A=（1，2，3）和B=（3，4，5）的向量，并计算三个向量组成图形的面积。
A=[1 2 3]；B=[3，4，5]；C=cross(A,B)S=dot(A,cross(B,C)) 123123123634325748437x x x x x x x x x ++=-++=---=-4 求矩阵a =1 1 10 4 -12 -2 1 的LU、QR分解。