实验2 数组矩阵及数值计算 matlab
Matlab教学第二章 MATLAB矩阵(数值)运算
为它专门建立一个M文件。其步骤为： • 第一步：使用编辑程序输入文件内容。 • 第二步：把输入的内容以纯文本方式存盘(设文件名为 mymatrix.m)。 • 第三步：在MATLAB命令窗口中输入mymatrix，就会自动建立一 个名为AM的矩阵，可供以后显示和调用。
１．建立矩阵
• 利用M文件建立矩阵：
• 3．矩阵的函数
• • • • • • • • • • • • • • • Matrix functions. expm - Matrix exponential. logm - Matrix logarithm. sqrtm - Matrix square root. funm - Evaluate general matrix function. Factorization utilities qrdelete - Delete a column or row from QR factorization. qrinsert - Insert a column or row into QR factorization. rsf2csf - Real block diagonal form to complex diagonal form. cdf2rdf - Complex diagonal form to real block diagonal form. balance - Diagonal scaling to improve eigenvalue accuracy. planerot - Givens plane rotation. cholupdate - rank 1 update to Cholesky factorization. qrupdate - rank 1 update to QR factorization.
实验二MATLAB矩阵及其运算
实验二 MATLAB矩阵及其运算学号：3111004117 姓名：颜琼专业：11电本日期：2013.10.15一、实验目的1.掌握MATLAB数据对象的特点以及数据的运算规则。
2.熟悉MATLAB中建立矩阵的方法以及矩阵处理的方法。
3.掌握MATLAB的分析方法。
二、实验环境计算机、MATLAB软件三、实验内容1.（1）新建一个.m文件，验证书本第15页例2-1；（2）用命令方式查看和保存代码中的所有变量；（3）用命令方式删除所有变量；（4）用命令方式载入变量z。
2. 将x=[4/3 1.2345e-6]在以下格式符下输出：短格式、短格式e方式、长格式、长格式e方式、银行格式、十六进制格式、+格式，并在实验报告中说明各种格式符下输出的结果有何区别。
3.完成指导书1-3题的内容。
4.完成思考练习第5题的内容。
四、实验过程及结果1、例1：（1）新建立es1.m文件，在MATLAB命令窗口中输入：x=1+2i;y=3-sqrt(17);z=(cos(abs(x+y))-sin(78*pi/180))/(x+abs(y))运行得到输出结果如下：>> es1z =-0.3488 + 0.3286i>>与例题结果一样，故，例1得到验证。
(2）用命令方式查看代码中的所有变量：在命令窗口输入：who whos用命令方式保存代码中的所有变量：在命令窗口输入：save es1(3) 用命令方式删除所有变量：在命令窗口输入：clear x y z(4) 用命令方式载入变量z：在命令窗口输入：load es1 z2、例2 ：新建立es2.m文件，在MATLAB命令窗口中输入：x = [4/3 1.2345e-6]format shortsingle(x)format short esingle(x)format longsingle(x)format long esingle(x)format banksingle(x)format hexsingle(x)format +single(x)运行得到输出结果为：>> es2x =1.3333 0.0000ans =1.3333 0.0000ans =1.3333e+000 1.2345e-006ans =1.3333334 0.0000012ans =1.3333334e+000 1.2345000e-006ans =1.33 0.00ans =3faaaaab 35a5b119ans =++>>1、求表达式的值：1)w = √2*(1+(0.34245*10^(-6)))；2)x = ((2*π*a)+((b+c)/(pi+a*b*c))-e^2)/(tan(b+c)+a),其中a=3.5,b=5， c=9.8；3)y=2*π*α^2*((1-π/4)*β-(0.8333-π/4)*α)，其中α=3.32，β=-7.9；4)z=1/2*e^(2*t)*log(t+√(1+t*t))；1）新建立了l11.m文件，在MATLAB命令窗口中输入：w = sqrt(2)*(1+(0.34245*10^(-6)))运行得到输出结果为：>> l11w =1.4142>>2）新建立l12.m文件，在MATLAB命令窗口中输入：a = 3.5;b = 5;c = -9.8;x = ((2*pi*a)+((b+c)/(pi+a*b*c))-exp(2))/(tan(b+c)+a)运行得到输出结果为：>> l12x =0.9829>>3）新建立l13.m文件，在MATLAB命令窗口中输入：a=3.32;b=-7.9;y=2*pi*a^2*((1-pi/4)*b-(0.8333-pi/4)*a)运行得到输出结果为：>> l13y =-128.4271>>4）新建立l14.m文件，在MATLAB命令窗口中输入：t=[2,1-3i;5,-0.65];z=1/2*exp(2*t)*log(t+sqrt(1+t*t))运行得到输出结果为：>> l14z =1.0e+004 *0.0057 - 0.0007i 0.0049 - 0.0027i1.9884 - 0.3696i 1.7706 - 1.0539i>>2、已知：A=[-1,5,4;0,7,8;3,61,7];B=[8,3,-1;2,5,3;-3,2,0];求表达式的值：1）A+6B和A^2-B+I(其中I为单位矩阵)；2）A*B、A.*B、和B.*A；3）A/B及B/A；4）[A,B]和[A([1，3]，：)；B^2]。
matlab-第2讲-矩阵及数值运算
2. roots —— 求多项式的根
a=[1 2 3; 4 5 6; 7 8 0]; p=poly(a) p= 1.00 -6.00 -72.00 -27.00 r=roots(p) r = 12.12 -5.73 ——显然 r是矩阵a的特征值 -0.39
• matlab规定多项式系数向量用行向量表示，一
MATLAB学习内容
第一讲：概述及基本操作
第二讲：矩阵及数值运算
第三讲：二维、三维绘图
第四讲： MATLAB的程序设计
第一天
第五讲：统计不回归分析
第六讲：揑值不拟合 第七讲：微分方程 第八讲：优化问题
第二天
第二讲 矩阵与数值运算
matlab 具有出色的数值计算能力，占 据世界上数值计算软件的主导地位
a./b= [a1/b1,a2/b2,…,an/bn] a.\b=[b1/a1,b2/a2,…,bn/an] a.^b=[a1^b1,a2^b2,…,an^bn] MATLAB (shuzu4)
返回
三、 矩阵的数组运算 数组运算指元素对元素的算术运算，
1、数组加减( .+ , .- ) a.+b a.- b
本讲内容
一、矩阵的创建不修改 二、矩阵运算 三、数组运算 四、利用矩阵来进行数值计算 – 多项式运算 – 线性方程组 – 数值统计 – 线性揑值 – 函数优化 – 微分方程的数值解
二、矩阵运算
1. 矩阵加、减（＋,－）运算
规则：
相加、减的两矩阵必须有相同的行和列，两矩 阵对应元素相加减。 允许参不运算的两矩阵之一是标量。标量不矩 阵的所有元素分别进行加减操作。
（2）数组-数组运算 当两个数组有相同维数时，加、减、乘、除、 幂运算可按元素对元素方式进行的，丌同大小戒维 数的数组是丌能进行运算的.
matlab实验二
北京工业大学Matlab实验报告**: ***学号: ************: **实验二、Matlab 的基本计算（一）实验目的1.掌握建立矩阵的方法。
2.掌握Matlab 各种表达式的书写规则以及常用函数的使用。
3.能用Matlab 进行基本的数组、矩阵运算。
4.掌握矩阵分析的方法以及能用矩阵运算或求逆法解线性方程组。
5.掌握Matlab 中的关系运算与逻辑运算。
(二）实验环境1．计算机2．MATLAB7.0集成环境（三）实验内容及要求1、熟练操作MATLAB7.0运行环境；2、自主编写程序，必要时参考相关资料；3、实验前应写出程序大致框架或完整的程序代码；4、完成实验报告。
（四）实验程序设计1.利用diag 等函数产生下列矩阵。
⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=032570800a ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=804050702b2.利用reshape 函数将1题中的a 和b 变换成行向量。
3.产生一个均匀分布在（-5，5）之间的随机矩阵（10×2），要求精确到小数点后一位。
4.已知：⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=76538773443412A ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--=731203321B求下列表达式的值：(1) B A K *611+=和I B A K +-=12（其中I 为单位矩阵）(2) B A K *21=和B A K *.22=(3) 331^A K =和3.32^A K =(4) B A K /41=和A B K \42=(5) ],[51B A K =和]2:);],3,1([[52^B A K = 5.下面是一个线性方程组：⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡52.067.095.03216/15/14/15/14/13/14/13/12/1x x x(1)求方程的解（矩阵除法和求逆法）(2)将方程右边向量元素3b 改为0.53，再求解，并比较3b 的变化和解的相对变化。
matlab实验矩阵和数组的操作
实验二矩阵和数组的操作一、实验环境计算机MATLAB软件二、实验目的1.掌握矩阵和数组的一般操作，包括创建、保存、修改和调用等。
2.学习矩阵和数组的加减运算与乘法、3.掌握对数组中元素的寻访与赋值，会对数组进行一般的操作。
三、实验内容与步骤1．用三种方法创建一个3X3矩阵，然后利用矩阵编辑器，将其扩充为4X5矩阵，并保存，试着调用它。
程序如下（1）直接输入法>> A=[3,2,1;4,5,6;7,8,9]A =3 2 14 5 67 8 9（2）利用MA TLAB提供的函数创建一个3X3的矩阵>> A=rand(3,3)A =0.4103 0.3529 0.13890.8936 0.8132 0.20280.0579 0.0099 0.1987（3）利用MA TLAB提供的“Matrix Editor”完成输入>> A=1A =1在矩阵编辑器中得如此矩阵1 0 00 0 00 0 0该成4X5 矩阵后如下1 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 0修改元素后为1 2 3 4 50 9 8 7 611 12 13 14 1510 19 18 17 16对文件进行保存使用save data 命令，用load data 命令刻把保存在文件的矩阵读到MATLAB的工作区的内存来2、建立一个等差数列，然后由它产生一个对角阵。
操作如下>> a=linspace(0,1.5,5)a =0 0.3750 0.7500 1.1250 1.5000>> B=diag(a)B =0 0 0 0 00 0.3750 0 0 00 0 0.7500 0 00 0 0 1.1250 00 0 0 0 1.50003、利用MATLAB的函数inv(A)求方阵A的逆矩阵。
操作如下>> A=[1,2;5,6]A =1 25 6>> B=inv(A)B =-1.5000 0.50001.2500 -0.2500四、练习1、创建一个5X5矩阵，提取主对角线以上的部分>> B=rand(5,5)B =0.0971 0.0344 0.1869 0.7547 0.11900.8235 0.4387 0.4898 0.2760 0.49840.6948 0.3816 0.4456 0.6797 0.95970.3171 0.7655 0.6463 0.6551 0.34040.9502 0.7952 0.7094 0.1626 0.5853 >> U=triu(B)U =0.0971 0.0344 0.1869 0.7547 0.11900 0.4387 0.4898 0.2760 0.49840 0 0.4456 0.6797 0.95970 0 0 0.6551 0.34040 0 0 0 0.5853 2、A=rand(3,3),B=magic(3,3),C=rand(3,4),计算AXBXC >> A=rand(3,3)A =0.1493 0.2543 0.92930.2575 0.8143 0.35000.8407 0.2435 0.1966>> B=magic(3)B =8 1 63 5 74 9 2>> C=rand(3,4)C =0.2511 0.3517 0.5497 0.75720.6160 0.8308 0.9172 0.75370.4733 0.5853 0.2858 0.3804>> A.*B.*C??? Error using ==> timesMatrix dimensions must agree.>> A*B*Cans =9.5982 12.7780 13.3892 13.39619.8497 12.9393 12.3754 13.12937.8080 10.2588 10.0835 11.84363.创建一个3X3的矩阵，并求其转置，逆矩阵>> D=rand(3,3)D = 0.5678 0.5308 0.12990.0759 0.7792 0.56880.0540 0.9340 0.4694>> E=conj(D)'E =0.5678 0.0759 0.05400.5308 0.7792 0.93400.1299 0.5688 0.4694>> F=inv(D)F =1.7826 1.3763 -2.16120.0529 -2.7944 3.3717-0.3102 5.4022 -4.33034、用两种方法求Ax=b的解（A为随机矩阵，b为四阶列向量）。
Matlab实验报告(二)矩阵和数组操作
Matlab实验报告（二）矩阵和数组操作一、实验目的1．掌握矩阵和数组的一般操作，包括创建、保存、修改和调用等。
2．学习矩阵和数组的加减运算与乘法。
3．掌握对数组中元素的寻访与赋值，会对数组进行一般的操作。
二、预备知识1．常用的产生特殊矩阵的函数?eye(m,n) 单位阵?rand（m,n) 随机矩阵?randn(m,n) 正态分布的随机矩阵?zeros(m,n) 零矩阵?ones(m,n) 全部元素都为1的矩阵?compan(A) 矩阵A的伴随矩阵?bankel(m,n) n维Hankel矩阵?invhilb(n) n维逆Hilbert矩阵?magic(n) n维Magic矩阵?toeplitz(m,n) Toeplitz矩阵?wilkinson(n) n维Wilkinson特征值测试矩阵?handamard(n) n 维Handamard矩阵?hilb(n) n维Hilbert矩阵?kron(A,B) Kronecker 张量积?pascal(n) n维Pascal矩阵?vander(A) 由矩阵A产生Vandermonde矩阵2．通过矩阵的结构变换，获得新矩阵表2 矩阵结构变化产生新矩阵L=tril(A) L=tril(A,k) 0 U=triu(A) U主对角线及以上的元素取矩阵A的元素，其余为0 L主对角线及以下元素取矩阵A 的元素，其余为0 L及第k条对角线及以下元素取矩阵A的元素，其余为U=triu(A,k) 0 B=rot90(A) B=rot90(A,k) B=fliplr(A) B=flipud(A) B=reshape(A,m,n) U 第k条对角线及以上的元素取矩阵A的元素，其余为矩阵A逆时针旋转90°得到B 矩阵A逆时针旋转k*90°得到B 矩阵A左右翻转得到B 矩阵A上下翻转得到B 将矩阵A的元素重新排列，得到m*n的新矩阵（m*n就等于A的行列式之积。
MATLAB基础教程 第2章 数组、矩阵及其运算
写出MATLAB表达式。 解：根据MATLAB的书写规则，以上MATLAB表达式为： （1）y=1/(a*log(1-x-1)+C1) （2）f=2*log(t)*exp(t)*sqrt(pi) （3）z=sin(abs(x)+abs(y))/sqrt(cos(abs(x+y))) （4）F=z/(z-exp(T*log(8)))
命令：X(3:-1:1)
命令：X(find(X>0.5)) 命令：X([1 2 3 4 4 3 2 1])
第二章 数组、矩阵及其运算
2.1 数组（矩阵）的创建和寻访
2. 二维数组的创建和寻访
例2-3 综合练习。将教材P.31~P.44的实例按顺序在MATLAB的 command窗口中练习一遍，观察并体会其输出结果。 （注意变量的大小写要和教材上的严格一致。）
A./B
B.\A
A的元素被B的对应元素相除
（与上相同）
第二章 数组、矩阵及其运算
2.3 数组、矩阵的其他运算
1. 乘方开方运算
数组的乘方运算与power函数 格式：c=a.^k或c=power(a,k) 例如： >> g=[1 2 3;4 5 6] >>g.^2 矩阵的乘方运算与mpower函数 格式：C=A^P或C=mpower(A,P) 注意：A必须为方阵
第二章 数组、矩阵及其运算
2.2 数组、矩阵的运算
3. 矩阵的加法、减法
运算规则是：若A和B矩阵的维数相同，则可以执行矩阵的加减运算， A和B矩阵的相应元素相加减。如果维数不相同，则MATLAB将给出
出错信息。
第二章 数组、矩阵及其运算
2.2 数组、矩阵的运算
3. 矩阵的乘法
matlab第二章实训报告
佛山科学技术学院《MATLAB教程第二章实训》报告专业姓名成绩班级学号日期一、目的1.学习matlab的数据类型2.矩阵和数组的算术运算3.字符串4.时间和日期5.结构体和元胞数组6.多维数组7.逻辑运算和关系运算8.数组的信息获取9.多项式二、步骤1.学习matlab的数据类型Matlab R2010a定义了15种基本的数据类型，包括整型、浮点型、字符型和逻辑型等。
用户甚至可以定义自己的数据类型。
Matlab内部的任何数据类型，都是按照数组的形式进行储存和运算的。
数值型包括整数和浮点数，其中整数包括有符号数和无符号数，浮点数包括单精度型和双精度型。
在默认情况下，matlab默认将所有数值都按照双精度浮点数类型来存储和操作。
（1）常数和变量Matlab的常数采用十进制表示，可以用带小数点的形式直接表示，也可以用科学记数法。
数值的表示范围是10^-309-10^309。
变量是数值计算的基本单元。
Matlab与其他的高级语言不同，变量使用是无需先定义，其名称就是第一次合法出现时的名称，因此用起来很便捷。
Matlab的变量命名有一定的规则：a.变量区分字母的大小写。
例如，“a”和“A”是不同的变量。
b.变量名不能超过63个字符，第63个字符后的字符会被忽略。
c.变量名必须以字母开头，变量名的组成可以是任意字母、数字或者下划线，但不能有空格和标点符号。
d.关键字（如if\while等）不能作为变量名。
在matlab中的所有表示符号包括函数名、文件名都是遵循变量名的命名规则。
Matlab中有一些自己的特殊变量，是由系统预先自动定义的，例如：ans——运算结果的默认变量名Pi——圆周率πEps——浮点数的相对误差Inf或inf——无穷大Nan或nan——不定值i或j——i=j=-1^1/2，虚数单位Nargin——函数的输入变量数目Nargout——函数的输出变量数目Realmin——最小的可用正实数Realmax——最大的可用正实数（2）整数和浮点数Matlab提供了8种内置的整数类型，为了在使用时提高运行速度和存储空间，应该尽量使用字节少的数据类型，可以使用类型转换函数将各种整数类型强制相互转换。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
2013年上学期MATLAB语言上机实验内容姓名：班级：学号：08实验2: 数组（矩阵）及数值计算1、指令A=ones(4,4);B=A(3,[2 3])C=A(3,:)<1写出输出结果：B = 1 1C = 0 0 0 02、指令>>B=magic(3)；写出求矩阵B的特征值D和特征向量V的语句：B=magic(3);[V,D]=eig(B)V =-0.5774 -0.8131 -0.3416-0.5774 0.4714 -0.4714-0.5774 0.3416 0.8131D =15.0000 0 00 4.8990 00 0 -4.89903、>>A=rand(100,100)；用向量化运算的方法，用1条语句实现求矩阵A中大于0.4的所有元素个数：A=rand(100,100);S=sum(sum(A>0.4))S = 59774要求在闭区间]2,0[ 上产生具有10个等距采样点的一维数组。
试用两种不同的指令实现。
A=[0:2*pi/9:2*pi]A =0 0.6981 1.3963 2.0944 2.7925 3.4907 4.1888 4.8869 5.5851 6.2832A=linspace(0,2*pi,10)A =0 0.6981 1.3963 2.0944 2.7925 3.4907 4.1888 4.8869 5.5851 6.28325、由指令rand('state',0),A=rand(3,5)生成二维数组A，试求该数组中所有大于0.5的元素的位置，分别求出它们的“全下标”和“单下标”。
A=rand(3,5)[ri,cj]=find(A>0.5);id=sub2ind(size(A),ri,cj);ri=ri';cj=cj'; disp(' ')disp('大于0.5的元素的全下标')disp(['行号 ',int2str(ri)])disp(['列号 ',int2str(cj)])disp(' ')disp('大于0.5的元素的单下标')disp(id')A =0.95012928514718 0.48598246870930 0.456467665168340.44470336435319 0.921812970744800.23113851357429 0.89129896614890 0.018503643248220.61543234810009 0.738207245810670.60684258354179 0.76209683302739 0.821407164295250.79193703742704 0.17626614449462大于0.5的元素的全下标行号 1 3 2 3 3 2 3 1 2列号 1 1 2 2 3 4 4 5 5大于0.5的元素的单下标1 3 5 6 9 11 12 13 146、在时间区间 [0,10]中，绘制t-=曲线。
要求分别采取“标量循环运算法”和15.0-cosy t2e“数组运算法”编写两段程序绘图。
t=linspace(0,10,200);N=length(t);y1=zeros(size(t));for k=1:Ny1(k)=1-exp(-0.5*t(k))*cos(2*t(k));endsubplot(1,2,1),plot(t,y1),xlabel('t'),ylabel('y1'),grid ony2=1-exp(-0.5*t).*cos(2*t);subplot(1,2,2),plot(t,y2),xlabel('t'),ylabel('y2'),grid on7、先运行指令x=-3*pi:pi/15:3*pi; y=x; [X,Y]=meshgrid(x,y); warning off; Z=sin(X).*sin(Y)./X./Y；产生矩阵Z。
（1）请问矩阵Z中有多少个“非数”数据？（2）用指令surf(X,Y,Z); shading interp观察所绘的图形。
（3）请写出绘制相应的“无裂缝”图形的全部指令。
x=-3*pi:pi/15:3*pi;y=x;[X,Y]=meshgrid(x,y);Z=sin(X).*sin(Y)./X./Y;NumOfNaN=sum(sum(isnan(Z)))subplot(1,2,1),surf(X,Y,Z),shading interp,title('有缝图')XX=X+(X==0)*eps;YY=Y+(Y==0)*eps;ZZ=sin(XX).*sin(YY)./XX./YY;subplot(1,2,2),surf(XX,YY,ZZ),shading interp,title('无缝图')y =sinint(x)y1 =1.65418、指令>A=3*ones(5,5);B=A(3,[1,2])C=A(:,3)>1写出输出结果： B =3 3C =11111 9、 指令>>B=magic(4), 按要求分别写出指令：求矩阵B 的行列式：__[R,ci]=rref(B)_____________________。
求矩阵B 的逆矩阵：__X=inv(B)_____________________。
求矩阵B 的迹： Z=orth(B) 。
求矩阵B 的特征值D 和特征向量V ： [V,D]=eig(B) 。
10、采用数值计算方法，画出dt tt x y x ⎰=0sin )(在]10 ,0[区间曲线，并计算)5.4(y 。
dt=0.0001;t=0:dt:10;t=t+(t==0)*eps;f=sin(t)./t;s=cumtrapz(f)*dt;plot(t,s,'LineWidth',3)ii=find(t==4.5);s45=s(ii)s45 =11、用quad 求取dx x e x sin 7.15⎰--ππ的数值积分，并保证积分的绝对精度为910-。
fx=@(x)exp(-abs(x)).*abs(sin(x));format longsq=quadl(fx,-10*pi,1.7*pi,1e-7)sq = 1.0878499381549812、求函数5.08.12cos 5.1)5(sin )(206.02++-=t t t e t t f t 在区间]5,5[-中的最小值点。
syms tfts=sin(5*t)^2*exp(0.06*t*t)-1.5*t*cos(2*t)+1.8*abs(t+0.5);dfdt=diff(fts,t);tmin=solve(dfdt,t)fobj3=subs(fts,t,tmin)tmin =-.60100931947716486053884417850955e-2fobj3 =.8990990814468455167020879772312413、设0)0(,1)0(,1)(2)(3)(22===+-dt dy y t y dt t dy dt t y d ，用数值法和符号法求5.0)(=t t y 。
数值法format longts=[0,1];y0=[1;0];dydt=@(t,y)[y(2);-2*y(1)+3*y(2)+1];[tt,yy]=ode45(dydt,ts,y0);y_05=interp1(tt,yy(:,1),0.5,'spline')y_05 =0.78958020790127符号法syms t;ys=dsolve('D2y-3*Dy+2*y=1','y(0)=1,Dy(0)=0','t')ys_05=subs(ys,t,sym('0.5'))ys =1/2-1/2*exp(2*t)+exp(t)ys_05 =.7895803564706055291685070521378014、求矩阵b Ax =的解，A 为4阶魔方阵，b 是)14(⨯的全1列向量。
A=magic(4);b=ones(4,1);[R,C]=rref([A,b]);x=A\b;xx=inv(A)*bxx =0.031250000000000.12500000000000-0.06250000000000-0.0156250000000015、求矩阵b Ax =的解，A 为4阶魔方阵，⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=4321b 。
A=magic(4);b=(1:4)';x=pinv(A)*bb_pinv=A*xx =0.023529411764710.123529411764710.123529411764710.02352941176471b_pinv =1.300000000000002.900000000000002.100000000000003.7000000000000016、求解二元函数方程组⎩⎨⎧=+=-0)cos(0)sin(y x y x 的解。
（用符号法解） S=solve('sin(x-y)','cos(x+y)','x','y');X=S.x,Y=S.yX =-1/4*pi1/4*piY =-1/4*pi1/4*pi17、已知有理分式)()()(x D x N x R =，其中)5.0)(3()(33++=x x x x N ，)125)(22()(232++-+=x x x x x D 。
（1）求该分式的商多项式)(x Q 和余多项式)(x r 。
clearN=conv([3 0 1 0],[1 0 0 0.5])D=conv([1 2 -2],[5 2 0 1])[Q,r]=deconv(N,D)N =3.00000000000000 0 1.000000000000001.50000000000000 0 0.50000000000000 0D =5 12 -6 -3 2 -2Q =0.60000000000000 -1.44000000000000r =-0.00000000000000 0.00000000000000 21.88000000000001 -5.34000000000000 -5.52000000000000 4.58000000000000 -2.88000000000000（2）用程序验算)()()()(x N x r x Q x D =+是否成立。