MATLAB实验六--答案
河南财经政法大学数学与信息科学学院
实验报告 （6）A的行列式 （7）A的迹 3、求向量组 1 det(A) trace(A)
1, 1, 2, 4 , 2 0, 3, 1, 2 ， 3 3, 0, 7, 14 ,
T T T T T
4 1, 1, 2, 0 , 5 2, 1, 5, 6
m=[2 a -1; a -1 1;4 5 -5] b=[1;2;-1] [m,n]=size(A); y=[]; if norm(b)>0 %非齐次方程组% %有唯一解% if rank(A)==rank([A,b]) if rank(A)==n x=A\b; else x=A\b; y=null(A,'r'); end else x=[]; end else x=zeros(n,1); if rank(A)<n disp('原方程有无穷个解，基础解系为y'); y=null(A,'r'); end 河南财经政法大学数学与信息科学学院 6 %非0解% %齐次方程组% %0解% disp('原方程有零解x'); %方程无解% disp('此方程无解'); %方程有无穷多个解，基础解系% disp('原方程有无穷个解，其特解为x，齐次方程的基础解系为y'); disp('原方程有唯一解x');
实验报告 end
syms t; A=[2,t,-1;t,-1,1;4,5,-5]; B=[2,t,-1,1;t,-1,1,2;4,5,-5,-1]; C=solve(det(A)); fprintf('有唯一解时 取值不为：'); C(1,1) C(2,1) B(1,2)= C(1,1); A(1,2)=B(1,2); B(1,2)=C(1,1); A(2,1)=B(1,2); if rank(A)==rank(B) fprintf('有无穷多组解时 取值为：'); C(1,1) D=reff(B) else C(1,1) fprintf('时，无解'); end B(1,2)= C(2,1); A(1,2)=B(1,2); B(1,2)=C(2,1); A(2,1)=B(1,2); if rank(A)==rank(B) fprintf('有无穷多组解时 取值为：'); C(2,1) D=reff(B) else C(2,1) fprintf('时，无解'); end
disp('rank(A)= rank(B)=2<3,所以λ=1时次方程组有无穷多组解.') 河南财经政法大学数学与信息科学学院 4
实验报告 format rat C=null(A,'r'); [R1,jb]=rref(A); syms k1 k2 x1 x2 D1=k1*R1(:,1)+k2*R1(:,2); disp('齐次方程组的为通解为X1') X1=[x1*k1;x2*k2;0] R2=rref(B); D2=R2(:,4); disp('此非齐次方程组的一个特解为X2') X2=[1;-1;0] disp('此非齐次方程组的通解为X') X=[x1*k1+1;x2*k2-1;0] disp('当λ=-4/5时:') clear A=[2 -4/5 -1;-4/5 -1 1;4 5 -5]; B=[2 -4/5 -1 1;-4/5 -1 1 2;4 5 -5 -1]; rA=rank(A) rB=rank(B) %rA =2 %rB =3 %求解齐次方程的有理基
1 1 1 8、设矩阵 A= 1 1 1 ，求一正交矩阵 Q 使得 Q-1AQ 为对角矩阵 1 1 1
程序如下： clear A=[1,1,1;1,1,1;1,1,1]; [V,D]=eig(A) 河南财经政法大学数学与信息科学学院 7 %按Q'AQ算的
实验报告 V'*V Q=V Q'AQ %判断可知V是正交阵 %把V赋给Q %Q'AQ是对角阵
程序如下： clear a1=[1 -1 2 4]'; a2=[0 3 1 2]'; a3=[3 0 7 14]'; a4=[1 -1 2 0]'; a5=[2 1 5 6]'; A=[a1 a2 a3 a4 a5] b=rank(A) [R,jb]=rref(A) A(:,jb) %b=3,故向量组线性相关 %求极大无关组所在列 %向量组的极大无关组
河南财经政法大学数学与信息科学学院
2
实验报告
2 x 1 3 x 2 4 x 3 1 5.求方程 5 x 1 4 x 2 x 3 4 的解 7 x1 3 x 2 2 x 3
（1）求齐次线性方程组的基础解系、并化为阶梯型 clear all A=[2 3 4;5 4 1;1 3 2]; b=[1 4 7]'; [m,n]=size(A); R=rank(A); B=[A b]; Rr=rank(B); format rat if R==Rr&&R==n % n为未知数的个数，判断是否有唯一解 x=A\b elseif R==Rr&&R<n x=A\b C=null(A, r) %应的基础解系 % 这种情形方程组通解xx=k(p)*C(:,P)(p=1…n-R) else X= 'No solution' end %判断是否有无穷解 %求特解 %求AX=0的基础解系,所得C为n-R列矩阵，这n-R列即为对 %输入矩阵A,b
disp('rank(A)< rank(B),所以λ=-4/5时此方程组无解.')
clear syms y; A=[2 y -1;y -1 1;4 5 -5]; b=[1 2 -1]'; DA=det(A) solve(DA,'y') %下面求y取不同值时，解的情况： %1、当y不等于-4\5且不等于1时，方程有唯一解，求解为： X=A\b% %当y= -4/5时： syms y; y= -4/5; A=[2 y -1;y -1 1;4 5 -5]; b=[1 2 -1]'; if rank(A)==rank([A,b]) fprintf('y= -4/5时，方程组有无穷多解\n ' ) X=A\b else fprintf('y= -4/5时，方程组无解\n ' ) 河南财经政法大学数学与信息科学学院 5
的秩，并判断其线性相关性，求极大无关组。
clc a1=[1,-1,2,4]'; a2=[0,3,1,2]'; a3=[3,0,7,14]'; a4=[1,-1,2,0]'; a5=[2,1,5,6]'; a=[a1,a2,a3,a4,a5]; [R,jb]=rref(a) r=length(jb) %求矩阵a的秩 if (r==4) %满秩 fprintf('这五个向量线性无关\n') else fprintf('这五个向量线性相关\n') end fprintf('极大无关组为：') a(:,jb) 4、已知 A=[1 2 3;6 3 4;6 3 1] 程序如下： clear A=[1 2 3;6 3 4;6 3 1] [V,D]=eig(A) b=rank(V) %rank(V)=3,V满秩，故A可对角化 求 A 的特征值、特征向量以及能否对角化
实验报告 end %当y=1时： syms y; y=1; A=[2 y -1;y -1 1;4 5 -5]; b=[1 2 -1]'; if rank(A)==rank([A,b]) fprintf('y=1时，方程组有无穷多解') else fprintf('y=1时，方程组无解') end format rat %指定有理格式输出 B=null(A) %求解空间的有理基 C=rref(A) syms k1 k2 X=k1*B(:,1)'+k2*B(:,2)'%写出方程组的通解
实验报告
实验六：线性代数的运算
实验目的：
1、熟悉 MATLAB 的矩阵初等运算； 2、掌握求矩阵的秩、逆、特征值、特征向量、化最简阶梯形的命令； 3、会用 MATLAB 求解线性方程组及化二次型为标准型.
பைடு நூலகம்
实验内容
1. 2. 3. 4. 5. 求矩阵的秩、逆 判断向量组是否线性相关，并找出向量组中的一个最大线性无关组 求解线性方程组（用分解） 求矩阵的特征值、特征向量以及矩阵能否对角化 化二次型为标准型.
[V,D]=eig(ones(3)) a1=null(ones(3)) a2=null(ones(3)-3*eye(3)) B=[a1 a2] B'*ones(3)*B
9、判断
2 2 2 f 2 x12 3 x2 9 x3 19 x4 2 x1 x2 4 x1 x3 2 x1 x4 6 x2 x4 12 x3 x4 是否正定。
河南财经政法大学数学与信息科学学院
3
实验报告
x 0 6.求行列式 0 y
程序如下： clear syms x y
y x 0 0
0 y x 0
0 0 y x
A=[x y 0 0;0 x y 0;0 0 x y;y 0 0 x] det(A)
syms x y; A=x*eye(4); for i=1:3 A(i,i+1)=y; end A(4,4)=x; A(4,1)=y; A det(A)
7.
2 x1 x 2 x 3 1 取何值时，线性方程组 x1 x 2 x 3 2 无解、唯一解、无穷多解？并求出有无 4 x 5 x 5 x 1 2 3 1
MATLAB全部实验及答案
MATLAB全部实验及答案实验一、MATLAB基本操作实验内容及步骤4、有关向量、矩阵或数组的一些运算（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.^B5、 已知⎪⎭⎫⎝⎛+⋅=-334sin 234πt e y t 若需要计算t ∈[-1,1]，取间隔为0.01，试计算出相对应的y 值。
MATLAB实验指导书(附答案)
MATLAB基础实验指导书漳州师范学院物电系2010年10月目录实验一MATLAB环境的熟悉与基本运算 (2)实验二MATLAB数值运算 (8)实验三MATLAB语言的程序设计 (12)实验四MATLAB的图形绘制 (16)实验五采用SIMULINK的系统仿真 (20)实验六MATLAB在电路中的应用 (25)实验七MATLAB在信号与系统中的应用 (27)实验八MATLAB在控制理论中的应用 (29)实验一 MATLAB环境的熟悉与基本运算一、实验目的：1．熟悉MATLAB开发环境2．掌握矩阵、变量、表达式的各种基本运算二、实验基本知识：1.熟悉MATLAB环境:MATLAB桌面和命令窗口、命令历史窗口、帮助信息浏览器、工作空间浏览器文件和搜索路径浏览器。
2.掌握MATLAB常用命令3.MATLAB变量与运算符变量命名规则如下：（1）变量名可以由英语字母、数字和下划线组成（2）变量名应以英文字母开头（3）长度不大于31个（4）区分大小写MATLAB中设置了一些特殊的变量与常量，列于下表。
MATLAB运算符，通过下面几个表来说明MATLAB的各种常用运算符表2 MATLAB算术运算符表3 MATLAB关系运算符表4 MATLAB逻辑运算符表5 MATLAB特殊运算4.MATLAB的一维、二维数组的寻访表6 子数组访问与赋值常用的相关指令格式5.MATLAB的基本运算表7 两种运算指令形式和实质内涵的异同表6.MATLAB的常用函数表8 标准数组生成函数表9 数组操作函数三、实验内容1、学习使用help命令，例如在命令窗口输入help eye，然后根据帮助说明，学习使用指令eye（其它不会用的指令，依照此方法类推）2、学习使用clc、clear，观察command window、command history和workspace等窗口的变化结果。
3、初步程序的编写练习，新建M-file，保存（自己设定文件名，例如exerc1、exerc2、exerc3……），学习使用MATLAB的基本运算符、数组寻访指令、标准数组生成函数和数组操作函数。
MATLAB)课后实验答案-精简版.do
实验一1. 先求下列表达式的值，然后显示MATLAB 工作空间的使用情况并保存全部变量。
(1) 0122sin851z 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 --+=++=--(4)2242011122123t t z t t t t t ⎧≤<⎪=-≤<⎨⎪-+≤<⎩，其中t =0:0.5:2.52. 已知：34787,2033657327A B ⎡⎤⎡⎤⎢⎥⎢⎥==⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦⎣⎦求下列表达式的值：(1) A+6*B 和A-B+I （其中I 为单位矩阵） (2) A*B 和A.*B (3) A^3和A.^3 (4) A/B 及B\A(5) [A,B]和[A([1,3],:);B^2] 3.设有矩阵A 和 B123453166789101769,111213141502341617181920970212223242541311A B ⎡⎤⎡⎤⎢⎥⎢⎥-⎢⎥⎢⎥⎢⎥⎢⎥==-⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦(1) 求它们的乘积C 。
(2) 将矩阵C 的右下角3×2子矩阵赋给D 。
(3) 查看MATLAB 工作空间的使用情况。
(1) 求[100,999]之间能被21整除的数的个数。
(2) 建立一个字符串向量，删除其中的大写字母。
ch='ABC123d4e56Fg9';则要求结果是：实验二1. 设有分块矩阵33322322E R A O S ⨯⨯⨯⨯⎡⎤=⎢⎥⎣⎦，其中E 、R 、O 、S 分别为单位矩阵、随机矩阵、零矩阵和对角阵，试通过数值计算验证22E R RS A O S +⎡⎤=⎢⎥⎣⎦。
解: M 文件如下；E=eye(3); R=rand(3,2); O=zeros(2,3);S=diag([1 2]);A=[E R a=[E,R+R*SO S] O ,S^2] A^2=a 2. 产生5阶希尔伯特矩阵H 和5阶帕斯卡矩阵P ，且求其行列式的值Hh 和Hp 以及它们的条件数Th 和Tp ，判断哪个矩阵性能更好。
Matlab实验报告六(三次样条与分段线性插值)
实验名称插值与拟合
所属课程数学软件与实验
实验类型综合型实验
专业信息与计算科学
班级
学号
姓名
指导教师
一、实验概述
【实验目的】
学会在matlab环境下使用几种不同的插值法和拟合两种方法构造函数依据已经知道的某些特殊点来推测实际问题中需要知道但又不便于测量出来的量。
【实验原理】
1.z=interp2(x0,y0,z0,x,y,’method’): 要求x0,y0单调；x, y可取为矩阵, 或x取行向量, y取为列向量, x,y的值分别不能超出x0,y0的范围。
2.分段线性插值与计算量与n无关;n越大, 误差越小.
3.三次样条插值比分段线性插值更光滑。
4.‘linear’ : 分段线性插值；‘spline’ : 三次样条
二、实验内容
问题1 对函数, x([-5,5], 分别用分段线性插值和三次样条插值作插值（其中插值节点不少于20）, 并分别作出每种插值方法的误差曲线.
1180 1320 1450 1420 1400 1300 700 900];
mesh(x,y,z)
xi=0:20:2800;
yi=0:20:2400;
zi=interp2(x,y,z,xi',yi,'cubic');
mesh(xi,yi,zi)
3.结果
4.结论及分析
通过实验,结果正确,分析无误。
三、实验小结
1270 1500 1200 1100 1350 1450 1200 1150
1230 1390 1500 1500 1400 900 1100 1060
1180 1320 1450 1420 1400 1300 700 900
Matlab实验指导书(含答案)汇总
实验Matlab 操作环境熟悉、实验目的1．初步了解Matlab 操作环境。
2．学习使用图形函数计算器命令funtool 及其环境。
二、实验内容熟悉Matlab 操作环境，认识命令窗口、内存工作区窗口、历史命令窗口；学会使用format 命令调整命令窗口的数据显示格式；学会使用变量和矩阵的输入，并进行简单的计算；学会使用who 和whos 命令查看内存变量信息；学会使用图形函数计算器funtool ，并进行下列计算：1．单函数运算操作。
求下列函数的符号导数(1) y=sin(x);(2) y=(1+x)A3*(2-x); 求下列函数的符号积分(1) y=cos(x);(2) y=1/(1+x^2);(3) y=1/sqrt(1-x^2);(4) y=(x-1)/(x+1)/(x+2);求反函数(1) y=(x-1)/(2*x+3);(2) y=exp(x);(3) y=log(x+sqrt(1+x^2));代数式的化简(1) (x+1)*(x-1)*(x-2)/(x-3)/(x-4);(2) sin (x)A2+cos(x)A2;(3) x+sin(x)+2*x-3*cos(x)+4*x*sin(x);2．函数与参数的运算操作。
从y=xA2通过参数的选择去观察下列函数的图形变化(1) y1=(x+1)A2⑵ y2=(x+2)A2(3) y3=2*x^2⑷ y4=x^2+2(5) y5=x^4⑹ y6=x^2/23．两个函数之间的操作求和(1) sin(x)+cos(x)(2) 1+x+x^2+x^3+x^4+x^5乘积(1) exp(-x)*sin(x)(2) sin(x)*x商(1) sin(x)/cos(x);⑵ x/(1+x^2);(3) 1/(x-1)/(x-2);求复合函数(1) y=exp(u) u=sin(x)(2) y=sqrt(u) u=1+exp(xA2)(3) y=sin(u) u=asin(x)(4) y=sinh(u) u=-x三、设计提示1．初次接触Matlab 应该注意函数表达式的文本式描述。
实验六、利用MATLAB计算复变函数在孤立奇点处的留数及进行复积分
实验六、利用MATLAB 计算复变函数在孤立奇点处的留数及进行复积分计算一、本实验教学的作用：熟悉MATLAB 基本命令与操作，利用MATLAB 计算复变函数在孤立奇点处的留数；利用MATLAB 计算围线积分。
通过实验具体操作，培养学生综合实践能力。
二、本实验教学目的及学生能力标准：会利用MATLAB 计算复变函数在孤立奇点处的留数；会利用MATLAB 计算围线积分。
三、实验内容：四、相关知识1.在孤立起点处的留数----通过求极限的方法计算留数假设已知奇点α和重数m 则用下面的MATLAB 语句求出相应的留数B=limit(F*(z-α),z, α) 单奇点B=limit(F*(z-α)*m ,z,m -1)/prod(1:m-1),z,α m 重奇点例1 计算z e z z z z f 23)3cos()1(1)(-+-=π在孤立起点处的留数 解 函数)(z f 在0=z 是三重奇点，在1=z 是简单奇点>>syms z>>f=cos(z+pi/3)*exp(-2*z)/z^3*(z-1);>>limit(diff(f*z^3,z,2)/prod(1;2),z,0);>> limit((f*(z-1),z,1)ans=-1/4-1/2*3^(1/2)1/2*exp(-2)*cos(1)-1/2*3^(1/2)*exp(-2)*sin(1)学生练习1 计算z e z z z z f 23)3sin()1(1)(-+-=π在孤立起点处的留数 提示 函数)(z f 在0=z 是三重奇点，在1=z 是简单奇点>>syms z;>>f=sin(z+pi/3)*exp(-2*z)/z^3*(z-1);>>limit(diff(f*z^3,z,2)/prod(1;2),z,0);>> limit((f*(z-1),z,1)ans=-1/4*3^(1/2)+1/2-1/2*exp(-2)*sin(1)+1/2*3^(1/2)*exp(-2)*cos(1)例2 计算3542)(zz i z z f ++=在孤立起点处的留数 解 函数)(z f 在0=z 是三重奇点，在i i z 2,2-=是简单奇点>>syms z;>>f= (z+2*i3) /z^5+z^3);>>limit(diff(f*z^3,z,2)/prod(1;2),z,0);>> limit((f*(z-2*i),z,2i)>> limit((f*(z+2*i),z,-2i)ans=i/8-i/8例3 计算函数1)(2-=z e z f z在∞=z 处的留数 解 函数)(z f 在扩充复平面有三个极点：∞-=,1,1z>>syms z>>z1=exp(z)/(z^2-1);>>B1=limit(z1*(z-1),z,1)>>B2=limit(z1*(z+1),z,-1)>>B=B1+B2ansB1=1/2*exp(1)B2=-1/2*exp(-1)B=1/2*exp(1)-1/2*expp(-1)学生练习2 计算函数4sin )(zz z z f +=在0=z 处的留数 提示 函数)(z f 在0=z 是四重奇点>>syms z;>>f= (sin(z)+z)/z^4;>>limit(diff(f*z^4,z,3)/prod(1;3),z,0);ans=-1/6学生练习3 计算下列函数在奇点处的留数：（1） z z z 212-+ （2）14-z z 解 在Matlab 命令窗口键入：>> [r1,p1,k1]=residue([1,1],[1,-2,0])r1 =1.5000-0.5000p1 =2k1 =[ ]>> [r2,p2,k2]=residue([1 0],[1 0 0 0 -1])r2 =0.25000.2500-0.2500 + 0.0000i-0.2500 - 0.0000ip2 =-1.00001.00000.0000 + 1.0000i0.0000 - 1.0000ik2 =[ ]反之：>> [B,A]=residue([0.2500 0.2500 -0.2500 -0.2500],[-1 1 i -i],[])B =0 0 1 0A =1 0 0 0 -12.求积分2.1 非闭合路径的积分非闭合路径的积分，用函数int 求解，方法同微积分部分的积分。
matlab实验内容答案
m a t l a b实验内容答案(总16页) -本页仅作为预览文档封面，使用时请删除本页-实验报告说明：matlab 课程实验需撰写8个实验报告，每个实验报告内容写每次实验内容中标号呈黑体大号字显示的题目。
第一次实验内容：实验一 MATLAB 运算基础一、实验目的1.熟悉启动和退出MATLAB 的方法。
2.熟悉MATLAB 命令窗口的组成。
3.掌握建立矩阵的方法。
4.掌握MATLAB 各种表达式的书写规则以及常用函数的使用。
二、实验内容1.先求下列表达式的值，然后显示MATLAB 工作空间的使用情况并保存全部变量。
（1）22sin8511z e ︒=+ （2）12ln(2z x =，其中2120.455i +⎡⎤=⎢⎥-⎣⎦（3）0.30.33sin(0.3), 3.0, 2.9, 2.8,,2.8,2.9,3.02a ae e z a a --=+=--- 提示：利用冒号表达式生成a 向量，求各点的函数值时用点乘运算。
（4）2220141122123t t z t t t t t ⎧≤<⎪=-≤<⎨⎪--≤<⎩，其中t ＝0：：提示：用逻辑表达式求分段函数值。
2.已知12344347873657A -⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦，131203327B -⎡⎤⎢⎥=⎢⎥⎢⎥-⎣⎦求下列表达式的值：（1） A+6=B 和A-B+I(其中I 为单位矩阵)。
（2） A*B 和A.*B 。
（3） A^3和A^.3 。
（4） A/B 和B\A 。
（5）[A ，B]和[A([1,3],;);B^2] 。
3.设有矩阵A 和B12345678910111213141516171819202122232425A ⎡⎤⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎣⎦， 30161769023497041311B ⎡⎤⎢⎥-⎢⎥⎢⎥=-⎢⎥⎢⎥⎢⎥⎣⎦（1） 求它们的乘积C 。
（2） 将矩阵C 的右下角3×2子矩阵赋给D（3） 查看MATLAB 工作空间使用情况。
(完整版)MATLAB)课后实验答案[1]
1 + e2 (2) z = 1 ln( x + 1 + x 2 ) ，其中 x = ⎡⎢ 2⎣-0.45 ⎦2 2 ⎪t 2 - 2t + 1 2 ≤ t <3 ⎨实验一MATLAB 运算基础1. 先求下列表达式的值，然后显示 MATLAB 工作空间的使用情况并保存全部变量。
(1) z = 2sin 8501221 + 2i ⎤5 ⎥(3) z = e 0.3a - e -0.3asin(a + 0.3) + ln 0.3 + a ,a = -3.0, - 2.9, L , 2.9, 3.03⎧t 2 0 ≤ t < 1 (4) z = ⎪t 2 - 11 ≤ t <2 ，其中 t=0:0.5:2.5 4⎩解：M 文件:z1=2*sin(85*pi/180)/(1+exp(2))x=[2 1+2*i;-.45 5];z2=1/2*log(x+sqrt(1+x^2))a=-3.0:0.1:3.0;3=(exp(0.3.*a)-exp(-0.3.*a))./2.*sin(a+0.3)+log((0.3+a)./2)t=0:0.5:2.5;z4=(t>=0&t<1).*(t.^2)+(t>=1&t<2).*(t.^2-1)+(t>=2&t<3) .*(t.^2-2*t+1)4.完成下列操作：(1)求[100,999]之间能被21整除的数的个数。
(2)建立一个字符串向量，删除其中的大写字母。
解：(1)结果：m=100:999;n=find(mod(m,21)==0);length(n)ans=43(2).建立一个字符串向量例如：ch='ABC123d4e56Fg9';则要求结果是：ch='ABC123d4e56Fg9';k=find(ch>='A'&ch<='Z');ch(k)=[]ch=⎣O2⨯3⎥，其中E、R、O、S分别为单位矩阵、随机矩阵、零矩S⎦阵和对角阵，试通过数值计算验证A=⎢⎥。
matlab实验指导答案详解(非常详细正确)
实验一MATLAB 工作环境熟悉及简单命令的执行一、实验目的：熟悉MATLAB 的工作环境，学会使用MATLAB 进行一些简单的运算。
二、实验内容：MATLAB 的启动和退出，熟悉MATLAB 的桌面（Desktop ），包括菜单（Menu ）、工具条（Toolbar ）、命令窗口(Command Window)、历史命令窗口、工作空间(Workspace)等；完成一些基本的矩阵操作；学习使用在线帮助系统。
三、实验步骤：1、启动MATLAB ，熟悉MATLAB 的桌面。
2、在命令窗口执行命令完成以下运算，观察workspace 的变化，记录运算结果。
（1）（365-52⨯2-70）÷3 >>(365-52*2-70)/3 ans = 63.6667（2）>>area=pi*2.5^2 area = 19.6350（3）已知x=3，y=4，在MATLAB 中求z ：()232y x y x z -= >>x=3 >>y=4>>z = x ^2 * y ^3 / (x - y) ^2 z = 576（4）将下面的矩阵赋值给变量m1,在workspace 中察看m1在内存中占用的字节数。
m1=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡11514412679810115133216 执行以下命令>>m1 =[16 2 3 13 ; 5 11 10 8 ; 9 7 6 12 ; 4 14 15 1 ] >>m1( 2 , 3 ) ans = 10 >>m1( 11 ) ans = 6>>m1( : , 3 ) ans =3 10 6 15>>m1( 2 : 3 , 1 : 3 ) ans =5 11 10 9 7 6>>m1( 1 ,4 ) + m1( 2 ,3 ) + m1( 3 ,2 ) + m1( 4 ,1) ans = 34（5）执行命令>>help abs查看函数abs 的用法及用途，计算abs( 3 + 4i ) （6）执行命令>>x=0:0.1:6*pi; >>y=5*sin(x); >>plot(x,y)（6）运行MATLAB 的演示程序，>>demo ，以便对MATLAB 有一个总体了解。
MATLAB)课后实验答案-精简版
实验一 MATLAB 运算基础1. 先求下列表达式的值，然后显示MA TLAB 工作空间的使用情况并保存全部变量。
(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.022aaee a z a a --+=++=--(4) 2242011122123t t z t t t t t ⎧≤<⎪=-≤<⎨⎪-+≤<⎩，其中t =0:0.5:2.52. 已知：1234413134787,2033657327A B --⎡⎤⎡⎤⎢⎥⎢⎥==⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦⎣⎦求下列表达式的值：(1) A+6*B 和A-B+I （其中I 为单位矩阵） (2) A*B 和A.*B(3) A^3和A.^3(4) A/B及B\A(5) [A,B]和[A([1,3],:);B^2]3. 设有矩阵A 和B1234530166789101769,111213141502341617181920970212223242541311A B ⎡⎤⎡⎤⎢⎥⎢⎥-⎢⎥⎢⎥⎢⎥⎢⎥==-⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦(1) 求它们的乘积C 。
(2) 将矩阵C 的右下角3×2子矩阵赋给D 。
(3) 查看MA TLAB 工作空间的使用情况。
4. 完成下列操作：(1) 求[100,999]之间能被21整除的数的个数。
(2) 建立一个字符串向量，删除其中的大写字母。
(2). 建立一个字符串向量 例如：ch='ABC123d4e56Fg9';则要求结果是：实验二 MATLAB 矩阵分析与处理1. 设有分块矩阵33322322E R A O S ⨯⨯⨯⨯⎡⎤=⎢⎥⎣⎦，其中E 、R 、O 、S 分别为单位矩阵、随机矩阵、零矩阵和对角阵，试通过数值计算验证22E R R S A OS +⎡⎤=⎢⎥⎣⎦。