东北大学MATLAB实验参考答案

MATLAB实验一：参考--答案

save wenjian.mat A a x load wenjian.mat A a x 小结及思考在 matlab 中，A*B 与 A .*B 以及 B./A 与 B.\A 之间有什么区别？
A*B 是 A 与 B 进行通常意义上的矩阵乘法得到的矩阵，而 A.*B 是矩阵 A 的每个元素与 B 的对应位置的元素相乘得到的矩阵。 B./A 是矩阵 B 的每一个元素除以 A 对应位置的元素得到的矩阵， B.\A 是矩阵 A 的每一个元素除以 B 对应位置的元素得到的矩阵。
（1）A=2015*eye(10) （2）A=8*(ones(8)-eye(8))
0 8 （2） 8
8 0 8

8 8 0 88
9. 设有分块矩阵 A
E33 O23
R32 ，其中 E,R,O,S 分别为单位阵、随机阵、零阵和对角阵， S 22
7. 写出下列函数的作用： clc Clear ones det triu fliplr zeros cat
清除命令窗口的记录清除已声明的变量和函数的记录生成全 1 矩阵求一个方阵的行列式取一个矩阵的上三角对矩阵进行左右旋转产生一个零矩阵拼接两个数组或者矩阵
8. 建立如下矩阵
0 0 2015 2015 0 0 （1） 0 0 2015 1010
3．
计算 1.369 sin
2
7 26.48 2.9 的值（输入程序） 10
1.369^2+sin(7/10*pi)*26.48^(1/2)/2.9
4．
1 3 3 1 1 2 输入矩阵 A 4 4 6 , B 2 2 2 ，指出下列命令的含义 6 8 9 3 3 3

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 运算基础1. 先求下列表达式的值，然后显示MATLAB 工作空间的使用情况并保存全部变量。

(1) 0122sin851z e =+ (2) 221ln(1)2z x x =++，其中2120.455i x +⎡⎤=⎢⎥-⎣⎦(3) 0.30.330.3sin(0.3)ln , 3.0, 2.9,,2.9,3.022a a e e a z a a --+=++=--(4) 2242011122123t t z t t t t t ⎧≤<⎪=-≤<⎨⎪-+≤<⎩，其中t =0:0.5:2.5 解： 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;z3=(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) 运算结果：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;z3=(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)z1 =0.2375z2 =0.7114 - 0.0253i 0.8968 + 0.3658i0.2139 + 0.9343i 1.1541 - 0.0044iz3 =Columns 1 through 40.7388 + 3.1416i 0.7696 + 3.1416i 0.7871 + 3.1416i 0.7913 + 3.1416iColumns 5 through 80.7822 + 3.1416i 0.7602 + 3.1416i 0.7254 + 3.1416i 0.6784 + 3.1416iColumns 9 through 120.6196 + 3.1416i 0.5496 + 3.1416i 0.4688 + 3.1416i 0.3780 + 3.1416iColumns 13 through 160.2775 + 3.1416i 0.1680 + 3.1416i 0.0497 + 3.1416i -0.0771 + 3.1416iColumns 17 through 20-0.2124 + 3.1416i -0.3566 + 3.1416i -0.5104 + 3.1416i -0.6752 + 3.1416iColumns 21 through 24-0.8536 + 3.1416i -1.0497 + 3.1416i -1.2701 + 3.1416i -1.5271 + 3.1416iColumns 25 through 28-1.8436 + 3.1416i -2.2727 + 3.1416i -2.9837 + 3.1416i -37.0245 Columns 29 through 32-3.0017 -2.3085 -1.8971 -1.5978Columns 33 through 36-1.3575 -1.1531 -0.9723 -0.8083Columns 37 through 40-0.6567 -0.5151 -0.3819 -0.2561Columns 41 through 44-0.1374 -0.0255 0.0792 0.1766Columns 45 through 480.2663 0.3478 0.4206 0.4841Columns 49 through 520.5379 0.5815 0.6145 0.6366Columns 53 through 560.6474 0.6470 0.6351 0.6119Columns 57 through 600.5777 0.5327 0.47740.4126Column 610.3388z4 =0 0.2500 0 1.2500 1.0000 2.2500 2. 已知：求下列表达式的值：(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]解：M 文件：A=[12 34 -4;34 7 87;3 65 7];B=[1 3 -1;2 0 3;3 -2 7]; A+6.*BA-B+eye(3)A*BA.*BA^3A.^3A/BB\A[A,B][A([1,3],:);B^2]运算结果：A=[12 34 -4;34 7 87;3 65 7];B=[1 3 -1;2 0 3;3 -2 7]; A+6.*BA-B+eye(3)A*BA.*BA^3A.^3A/BB\A[A,B][A([1,3],:);B^2]ans =18 52 -1046 7 10521 53 49ans =12 31 -332 8 840 67 1ans =68 44 62309 -72 596154 -5 241ans =12 102 468 0 2619 -130 49ans =37226 233824 48604247370 149188 60076678688 454142 118820 ans =1728 39304 -6439304 343 65850327 274625 343 ans =16.4000 -13.6000 7.600035.8000 -76.2000 50.200067.0000 -134.0000 68.0000109.4000 -131.2000 322.8000-53.0000 85.0000 -171.0000-61.6000 89.8000 -186.2000ans =12 34 -4 1 3 -134 7 87 2 0 33 65 7 3 -2 7ans =12 34 -43 65 74 5 111 0 1920 -5 403. 设有矩阵A和B(1) 求它们的乘积C。

东北大学 matlab实验

《MATLAB 语言与应用》上机实验作业第一部分2、用MATLAB 语句输入矩阵A 和B⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=1423143212344321A ， ⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡++++++++++++++++=4j 11j43j 22j34j 11j 42j 33j 24j 13j 22j 31j 41j 42j 33j 24j 1B前面给出的是44⨯矩阵，如果给出5)6,5(=A 命令将得出什么结果？ MATLAB 结果：>> A=[1 2 3 4;4 3 2 1;2 3 4 1;3 2 4 1];>> B=[1+4j 2+3j 3+2j 4+1j;4+1j 3+2j 2+3j 1+4j;2+3j 3+2j 4+1j 1+4j;3+2j 2+3j 4+1j 1+4j]; >> A A =1 2 3 4 4 3 2 1 2 3 4 1 3 2 4 1>> B B =1.0000 + 4.0000i2.0000 +3.0000i 3.0000 + 2.0000i4.0000 + 1.0000i 4.0000 + 1.0000i 3.0000 + 2.0000i 2.0000 + 3.0000i 1.0000 + 4.0000i 2.0000 + 3.0000i 3.0000 + 2.0000i 4.0000 + 1.0000i 1.0000 + 4.0000i 3.0000 + 2.0000i 2.0000 + 3.0000i 4.0000 + 1.0000i 1.0000 + 4.0000i>> A(5,6)=5 A =1 2 3 4 0 0 4 3 2 1 0 0 2 3 4 1 0 0 3 2 4 1 0 00 0 0 0 0 53、假设已知矩阵A，试给出相应的MATLAB命令，将其全部偶数行提取出来，赋给B矩阵，用magic(8)A=命令生成A矩阵，用上述命令检验一下结果是不是正确。

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 运算基础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 a z 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 分别为单位矩阵、随机矩阵、零矩阵与对角阵，试通过数值计算验证22E R RS A OS +⎡⎤=⎢⎥⎣⎦。

解: M 文件如下；5、下面就是一个线性方程组：1231112340.951110.673450.52111456x x x ⎡⎤⎢⎥⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎢⎥⎢⎥⎣⎦(1) 求方程得解。

(2) 将方程右边向量元素b 3改为0、53再求解，并比较b 3得变化与解得相对变化。

(3) 计算系数矩阵A 得条件数并分析结论。

解： M 文件如下：实验三选择结构程序设计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 值。

matlab实验二答案

实验二 MATLAB环境的熟悉与基本运算（二）一、实验目的1．熟悉MATLAB 开发环境2．掌握矩阵、变量、表达式的各种基本运算二、实验内容1．设两个复数a=1+2i，b=3-4i，计算a+b，a-b，a×b，a/b。

程序：clca=1+2i,b=3-4i，x1=a+b，x2=a-b，x3=a*b，x4=a/b2．已知矩阵112A=134245⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦，1518B=20362545⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦，求A\B。

程序clca=[1 1 2;1 3 4;2 4 5] b=[15 18;20 36;25 45] x=a\b3．已知矩阵A为四阶魔方矩阵，矩阵13572468B=35794688⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦，求A+B，A-B和A*B。

程序：clca=magic(4)b=[1 3 5 7;2 4 6 8;3 5 7 9;4 6 8 8] x1=a+bx2=a-bx3=a*b4.试用简单的语句输入下面的矩阵（至少用两种方法）：112233A 00A=0A 00A ⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦，其中11-210A =0-2100-2⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦，22-310A =0-3100-3⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦，33-410A =0-4100-4⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦注：函数diag(V ，K )用来建立一个对角矩阵，其中V 为某个向量，K 为向量V 偏离主对角线的列数，K 为零时表示V 为主对角线，K 大于零的数时表示V 在主对角线以上，K 为小于零的数时表示V 在主对角线以下。

方法一：程序：clcA11=[-2 1 0;0 -2 1;0 0 -2] A22=[-3 1 0;0 -3 1;0 0 -3]A33=[-4 1 0;0 -4 1;0 0 -4]b=zeros(3)V=[A11 b b;b A22 b;b b A33]方法二：程序：clcA11=diag([-2 -2 -2],0)+diag([1 1],1)A22=diag([-3 -3 -3],0)+diag([1 1],1)A33=diag([-4 -4 -4],0)+diag([1 1],1)b=zeros(3)A=[A11 b b;b A22 b;b b A33]5. 用rand函数生成一个6阶的方阵A，(1)查询方阵A 第2 行、第3 列的元素，(2)查询方阵A 第2 行的所有元素，(3)查询方阵A 第6 列的所有元素，(4)给方阵A 第3 行、第2 列的元素赋值为5，(5)给方阵A 第3 行的所有元素赋值为4，(6)给方阵A 第1列的所有元素赋值为1。

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]

MATLAB)课后实验答案[1]实验一 MATLAB 运算基础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.5 解：4. 完成下列操作：(1) 求[100,999]之间能被21整除的数的个数。

(2) 建立一个字符串向量，删除其中的大写字母。

解：(1) 结果：(2). 建立一个字符串向量例如：ch='ABC123d4e56Fg9';则要求结果是：实验二 MATLAB 矩阵分析与处理1. 设有分块矩阵33322322E R A O S ⨯⨯⨯⨯⎡⎤=⎢⎥⎣⎦，其中E 、R 、O 、S分别为单位矩阵、随机矩阵、零矩阵和对角阵，试通过数值计算验证22E R RS A O S +⎡⎤=⎢⎥⎣⎦。

解: M 文件如下；5. 下面是一个线性方程组：123d4e56g91231112340.951110.673450.52111456x x x ⎡⎤⎢⎥⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎢⎥⎢⎥⎣⎦(1) 求方程的解。

(2) 将方程右边向量元素b 3改为0.53再求解，并比较b 3的变化和解的相对变化。

(3) 计算系数矩阵A 的条件数并分析结论。

解： M 文件如下：实验三选择结构程序设计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 值。

MATLAB 课后实验答案

实验一 MATLAB 运算基础1. 先求下列表达式的值，然后显示MATLAB 工作空间的使用情况并保存全部变量。

(1) 0122sin851z e =+ (2) 221ln(1)2z x x =++，其中2120.455i x +⎡⎤=⎢⎥-⎣⎦(3) 0.30.330.3sin(0.3)ln , 3.0, 2.9,,2.9,3.022a a e e a z a a --+=++=--(4) 2242011122123t t z t t t t t ⎧≤<⎪=-≤<⎨⎪-+≤<⎩，其中t =0:0.5:2.5 解： 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;z3=(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) 运算结果：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;z3=(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)z1 =0.2375z2 =0.7114 - 0.0253i 0.8968 + 0.3658i0.2139 + 0.9343i 1.1541 - 0.0044iz3 =Columns 1 through 40.7388 + 3.1416i 0.7696 + 3.1416i 0.7871 + 3.1416i 0.7913 + 3.1416iColumns 5 through 80.7822 + 3.1416i 0.7602 + 3.1416i 0.7254 + 3.1416i 0.6784 + 3.1416iColumns 9 through 120.6196 + 3.1416i 0.5496 + 3.1416i 0.4688 + 3.1416i 0.3780 + 3.1416iColumns 13 through 160.2775 + 3.1416i 0.1680 + 3.1416i 0.0497 + 3.1416i -0.0771 + 3.1416iColumns 17 through 20-0.2124 + 3.1416i -0.3566 + 3.1416i -0.5104 + 3.1416i -0.6752 + 3.1416iColumns 21 through 24-0.8536 + 3.1416i -1.0497 + 3.1416i -1.2701 + 3.1416i -1.5271 + 3.1416iColumns 25 through 28-1.8436 + 3.1416i -2.2727 + 3.1416i -2.9837 + 3.1416i -37.0245 Columns 29 through 32-3.0017 -2.3085 -1.8971 -1.5978Columns 33 through 36-1.3575 -1.1531 -0.9723 -0.8083Columns 37 through 40-0.6567 -0.5151 -0.3819 -0.2561Columns 41 through 44-0.1374 -0.0255 0.0792 0.1766Columns 45 through 480.2663 0.3478 0.4206 0.4841Columns 49 through 520.5379 0.5815 0.6145 0.6366Columns 53 through 560.6474 0.6470 0.6351 0.6119Columns 57 through 600.5777 0.5327 0.47740.4126Column 610.3388z4 =0 0.2500 0 1.2500 1.0000 2.2500 2. 已知：求下列表达式的值：(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]解：M 文件：A=[12 34 -4;34 7 87;3 65 7];B=[1 3 -1;2 0 3;3 -2 7]; A+6.*BA-B+eye(3)A*BA.*BA^3A.^3A/BB\A[A,B][A([1,3],:);B^2]运算结果：A=[12 34 -4;34 7 87;3 65 7];B=[1 3 -1;2 0 3;3 -2 7]; A+6.*BA-B+eye(3)A*BA.*BA^3A.^3A/BB\A[A,B][A([1,3],:);B^2]ans =18 52 -1046 7 10521 53 49ans =12 31 -332 8 840 67 1ans =68 44 62309 -72 596154 -5 241ans =12 102 468 0 2619 -130 49ans =37226 233824 48604247370 149188 60076678688 454142 118820 ans =1728 39304 -6439304 343 65850327 274625 343 ans =16.4000 -13.6000 7.600035.8000 -76.2000 50.200067.0000 -134.0000 68.0000109.4000 -131.2000 322.8000-53.0000 85.0000 -171.0000-61.6000 89.8000 -186.2000ans =12 34 -4 1 3 -134 7 87 2 0 33 65 7 3 -2 7ans =12 34 -43 65 74 5 111 0 1920 -5 403. 设有矩阵A和B(1) 求它们的乘积C。

Matlab实验指导书(含答案)

实验一：Matlab操作环境熟悉一、实验目的1．初步了解Matlab操作环境。

2．学习使用图形函数计算器命令funtool及其环境。

二、实验内容熟悉Matlab操作环境，认识命令窗口、内存工作区窗口、历史命令窗口；学会使用format命令调整命令窗口的数据显示格式；学会使用变量和矩阵的输入，并进行简单的计算；学会使用who和whos命令查看内存变量信息；学会使用图形函数计算器funtool，并进行下列计算：1．单函数运算操作。

求下列函数的符号导数(1) y=sin(x);(2) y=(1+x)^3*(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)^2+cos(x)^2;(3) x+sin(x)+2*x-3*cos(x)+4*x*sin(x);2．函数与参数的运算操作。

从y=x^2通过参数的选择去观察下列函数的图形变化(1) y1=(x+1)^2(2) y2=(x+2)^2(3) y3=2*x^2(4) y4=x^2+2(5) y5=x^4(6) 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);(2) x/(1+x^2);(3) 1/(x-1)/(x-2);求复合函数(1) y=exp(u) u=sin(x)(2) y=sqrt(u) u=1+exp(x^2)(3) y=sin(u) u=asin(x)(4) y=sinh(u) u=-x三、设计提示1．初次接触Matlab应该注意函数表达式的文本式描述。

Matlab实验第一次实验答案

实验一Matlab使用方法和程序设计一、实验目的1、掌握Matlab软件使用的基本方法；2、熟悉Matlab的数据表示、基本运算和程序控制语句3、熟悉Matlab绘图命令及基本绘图控制4、熟悉Matlab程序设计的基本方法二、实验内容：1、帮助命令使用help命令，查找 sqrt（开方）函数的使用方法；解：sqrtSquare rootSyntaxB = sqrt(X)DescriptionB = sqrt(X) returns the square root of each element of the array X. For the elements of X that are negative or complex, sqrt(X) produces complex results.RemarksSee sqrtm for the matrix square root.Examplessqrt((-2:2)')ans =0 + 1.4142i0 + 1.0000i1.00001.41422、矩阵运算（1）矩阵的乘法已知A=[1 2;3 4]; B=[5 5;7 8]; 求A^2*B解：A=[1 2;3 4 ];B=[5 5;7 8 ];A^2*B（2）矩阵除法已知 A=[1 2 3;4 5 6;7 8 9]; B=[1 0 0;0 2 0;0 0 3];A\B,A/B解：A=[1 2 3;4 5 6;7 8 9 ];B=[1 0 0;0 2 0;0 0 3 ];A\B,A/B（3）矩阵的转置及共轭转置已知A=[5+i,2-i,1;6*i,4,9-i];求A.', A'解：A=[5+1i,2-1i,1;6*1i,4,9-1i ];A1=A.',A2=A'（4）使用冒号表达式选出指定元素已知： A=[1 2 3;4 5 6;7 8 9];求A中第3列前2个元素；A中所有列第2，3行的元素；方括号[]解：A=[1 2 3;4 5 6;7 8 9 ];B=A([1,2],[3]),C=A(2:end, : )用magic 函数生成一个4阶魔术矩阵，删除该矩阵的第四列 3、多项式（1）求多项式 42)(3--=x x x p 的根解：A=[1 0 -2 -4];B=roots(A)（2）已知A=[1.2 3 5 0.9;5 1.7 5 6;3 9 0 1;1 2 3 4] ，求矩阵A 的特征多项式;解：A=[1.2 3 5 .9 ; 5 1.7 5 6 ;3 9 0 1 ;1 2 3 4]; A=poly(A); A=poly2sym(A)把矩阵A作为未知数代入到多项式中；4、基本绘图命令（1）绘制余弦曲线 y=cos(t)，t∈[0，2π]解：t=0:.1:2*pi;y=cos(t);plot(t,y),grid（2）在同一坐标系中绘制余弦曲线y=cos(t-0.25)和正弦曲线y=sin(t-0.5)，t∈[0，2π] t=0:.1:2*pi;y1=cos(t-0.25);y2=sin(t-0.25);plot(t,y1，t,y2)grid5、基本绘图控制绘制[0，4π]区间上的x1=10sint曲线，并要求：（1）线形为点划线、颜色为红色、数据点标记为加号；（2）坐标轴控制：显示范围、刻度线、比例、网络线（3）标注控制：坐标轴名称、标题、相应文本；t=0:.01:4*pi;y=10*sin(t);plot(t,y,'-.',t,y,'r')grid6、基本程序设计（1）编写命令文件：计算1+2+…+n<2000 时的最大n值；解法1：s=0;i=0;while(s<2000) i=i+1;s=s+i; ends=s-i,i=i-1解法2：s=0;for i=1:1000; s=s+i;if(s>2000) ,break;endends=s-i,i=i-1（2）编写函数文件：分别用for和while循环结构编写程序，求2的0到n次幂的和。

MATLAB)课后实验答案

A = 1.0000 0 0 0 0 0 1.0000 0 0 0 0 0 1.0000 0 0 0.5383 0.9961 0.0782 1.0000 0 0.4427 0.1067 0.9619 0 2.0000
a = 1.0000 0 0 0 0 0 1.0000 0 0 0 0 0 1.0000 0 0 1.0767 1.9923 0.1564 1.0000 0 1.3280 0.3200 2.8857 0 4.0000
P= 1 1 1 1 1 2 3 4 1 3 6 10 1 4 10 20 1 5 15 35
1
5
15
35Βιβλιοθήκη 70Hh = 3.7493e-012
Hp = 1
Th = 4.7661e+005
Tp = 8.5175e+003 因为它们的条件数 Th>>Tp,所以 pascal 矩阵性能更好。 3. 建立一个 5×5 矩阵，求它的行列式值、迹、秩和范数。解： M 文件如下：
实验一
MATLAB 运算基础
1. 先求下列表达式的值，然后显示 MATLAB 工作空间的使用情况并保存全部变量。 (1) z1
2sin 850 1 e2
1 2
1 2i 2 5 0.45
(2) z2 ln( x 1 x 2 ) ，其中 x (3) z3
用 if 语句实现，分别输出 x=-5.0,-3.0,1.0,2.0,2.5,3.0,5.0 时的 y 值。解：M 文件如下：
运算结果有： f(-5) y= 14 >> f(-3)
y= 11 >> f(1) y= 2 >> f(2) y= 1 >> f(2.5) y= -0.2500 >> f(3) y= 5 >> f(5) y= 19

东北大学MATLAB实验答案

东北大学《MATLAB语言与应用》实验东北大学MATLAB实验课习题答案第一部分MATLAB语言编程、科学绘图与基本数学问题求解2、>> A=[1,2,3,4;4,3,2,1;2,3,4,1;3,2,4,1]A =1 2 3 44 3 2 12 3 4 13 24 1B=[1+4j,2+3j,3+2j,4+1j;4+1j,3+2j,2+3j,1+4j;2+3j,3+2j,4+1j,1+4j;3+2j,2+3j,4+1j,1+4j]B =1.0000 + 4.0000i2.0000 +3.0000i 3.0000 + 2.0000i4.0000 + 1.0000i4.0000 + 1.0000i 3.0000 + 2.0000i 2.0000 + 3.0000i 1.0000 + 4.0000i2.0000 +3.0000i 3.0000 + 2.0000i4.0000 + 1.0000i 1.0000 + 4.0000i3.0000 + 2.0000i 2.0000 + 3.0000i4.0000 + 1.0000i 1.0000 + 4.0000i >> A(5,6)=5A =1 2 3 4 0 04 3 2 1 0 02 3 4 1 0 03 24 1 0 00 0 0 0 0 53、A=magic(8)A =64 2 3 61 60 6 7 579 55 54 12 13 51 50 1617 47 46 20 21 43 42 2440 26 27 37 36 30 31 3332 34 35 29 28 38 39 2541 23 22 44 45 19 18 4849 15 14 52 53 11 10 568 58 59 5 4 62 63 1 >> B=A(2:2:end,:)B =9 55 54 12 13 51 50 1640 26 27 37 36 30 31 3341 23 22 44 45 19 18 488 58 59 5 4 62 63 1 4．i=0:63;s=sum(2.^i)s =1.8447e+0195、(1) >> z=sin(1./t);Warning: Divide by zero.>> plot(t,z)-1-0.8-0.6-0.4-0.20.20.40.60.81-1-0.8-0.6-0.4-0.200.20.40.60.81(2)>> t=[-pi:0.05:-1.8,-1.799:.001:-1.2,-1.2:0.05:1.2,1.201:0.001:1.8,1.81:0.05:pi]; >> y=sin(tan(t))-tan(sin(t)); >> plot(t,y)-4-3-2-101234-3-2-11236.三视图>> surf(x,y,z),shading flat;zlim([0,15])>> xx=[-2:.1:-1.2,-1,1:0.02:-0.9,-0.8:0.1:0.8,0.9:0.02:1.1,1.2:0.1:2]; >> yy=[-1:0.1:-0.2,-0.1:0.02:0.1,0.2:.1:1];>> [x,y]=meshgrid(xx,yy);>> subplot(224),surf(x,y,z)>> subplot(221),surf(x,y,z),view(0,90);>> subplot(222),surf(x,y,z),view(90,0);>> subplot(223),surf(x,y,z),view(0,0);1-2-1012-1-0.500.51-1-0.500.510204060-2-11202040607.（1）>> syms x;f=(3.^x+9.^x)^(1./x);L=limit(f,x,inf) L = 9（2）>>syms x y;f=x*y/(sqrt(x*y+1)-1);L1=limit(limit(f,x,0),y,0) L1 = 2 （3）>> syms x y;>> f=(1-cos(x^2+y^2))/((x^2+y^2)*exp(x^2+y^2)); >> L=limit(limit(f,x,0),y,0) L = 08．先建立M文件：function result=paradiff(y,x,t,n)if mod(n,1)~=0,error('n should positive integer,please correct') elseif n==1,result=diff(y,t)/diff(x,t);else,result=diff(paradiff(y,x,t,n-1),t)/diff(x,t);end,end然后调用函数：>> syms t;x=log(cos(t));y=cos(t)-t*sin(t);>> f=paradiff(y,x,t,1);>> [n,d]=numden(f);>> F=simple(n)/simple(d)F =(2*sin(t)+t*cos(t))*cos(t)/sin(t)>> syms t;x=log(cos(t));y=cos(t)-t*sin(t);>> f=paradiff(y,x,t,1);>> syms t;x=log(cos(t));y=cos(t)-t*sin(t);>> f=paradiff(y,x,t,2);>> [n,d]=numden(f);>> F=simple(n)/simple(d)F =-cos(t)*(3*cos(t)^2*sin(t)+cos(t)^3*t-2*sin(t)-2*t*cos(t))/sin(t)^3 >> subs(F,t,pi/3)ans =1.53879.>> syms x y t;>> f=exp(-t^2);>> I=simple(int(f,t,0,x*y))I =1/2*pi^(1/2)*erf(x*y)>> F=x/y*diff(I,x,2)-2*diff(diff(I,x),y)+diff(I,y,2)F =2*x^2*y^2*exp(-x^2*y^2)-2*exp(-x^2*y^2)-2*x^3*y*exp(-x^2*y^2) >> K=simple(F)K =-2*exp(-x^2*y^2)*(-x^2*y^2+1+x^3*y)10.（1）>> syms n;>> S=symsum(1/((2*n)^2-1),n,1,inf)S =1/2（2）>> syms k n>> limit(n*symsum(1/(n^2+k*pi),k,1,n),n,inf)ans =111.（1）>> syms t;>> syms a positive;>> x=a*(cos(t)+t*sin(t));>> y=a*(sin(t)-t*cos(t));>> I=int((x^2+y^2)*sqrt(diff(x,t)^2+diff(y,t)^2),t,0,2*pi)I =2*a^3*pi^2+4*a^3*pi^4（2）>> syms t;>> syms a b c positive;>> x=c/a*cos(t);>> y=c/b*sin(t);>> F=[y*x^3+exp(y),x*y^3+x*exp(y)-2*y];>> ds=[diff(x,t);diff(y,t)];>> I=int(F*ds,t,pi,0)I =2/15*c*(-2*c^4+15*b^4)/a/b^412.首先编写M程序：function A=vander(v)n=length(v);v=v(:);A=sym(ones(n));for j=n-1:-1:1,A(:,j)=v.*A(:,j+1);end>> syms a b c d e;>> A=[a,b,c,d,e];>> V=vander(A)V =[ a^4, a^3, a^2, a, 1][ b^4, b^3, b^2, b, 1][ c^4, c^3, c^2, c, 1][ d^4, d^3, d^2, d, 1][ e^4, e^3, e^2, e, 1]>> det(V),simple(ans)ans =(c-d)*(b-d)*(b-c)*(a-d)*(a-c)*(a-b)*(-d+e)*(e-c)*(e-b)*(e-a) 13.>> A=[-2,0.5,-0.5,0.5;0,-1.5,0.5,-0.5;2,0.5,-4.5,0.5;2,1,-2,-2]; >> [V J]=jordan(sym(A))V =[ 0, 1/2, 1/2, -1/4][ 0, 0, 1/2, 1][ 1/4, 1/2, 1/2, -1/4][ 1/4, 1/2, 1, -1/4]J =[ -4, 0, 0, 0][ 0, -2, 1, 0][ 0, 0, -2, 1][ 0, 0, 0, -2]14.先编写M文件：function X=lyap(A,B,C)if nargin==2,C=B;B=A';end[nr,nc]=size(C);A0=kron(A,eye(nc))+kron(eye(nr),B');tryC1=C';x0=-inv(A0)*C1(:);X=reshape(x0,nc,nr)';catch,error('singular matrix found.'),end数值解为：>> A=[3,-6,-4,0,5;1,4,2,-2,4;-6,3,-6,7,3;-13,10,0,-11,0;0,4,0,3,4]; >> B=[3,-2,1;-2,-9,2;-2,-1,9];>> C=[-2,1,-1;4,1,2;5,-6,1;6,-4,-4;-6,6,-3];>> X=lyap(A,B,C)X =-4.0569 -14.5128 1.56530.0356 25.0743 -2.74089.4886 25.9323 -4.41772.6969 21.6450 -2.88517.7229 31.9100 -3.7634>> norm(A*X+X*B+C)ans =3.9870e-013解析解为：>> X=lyap(sym(A),B,C)X =[ -434641749950/107136516451, -4664546747350/321409549353, 503105815912/321409549353][ 3809507498/107136516451, 8059112319373/321409549353, -880921527508/321409549353][ 1016580400173/107136516451, 8334897743767/321409549353, -1419901706449/321409549353][ 288938859984/107136516451, 6956912657222/321409549353, -927293592476/321409549353][ 827401644798/107136516451, 10256166034813/321409549353, -1209595497577/321409549353]>> A*X+X*B+Cans =[ 0, 0, 0][ 0, 0, 0][ 0, 0, 0][ 0, 0, 0][ 0, 0, 0]15.（1）>> A=[-4.5,0,0.5,-1.5;-0.5,-4,0.5,-0.5;1.5,1,-2.5,1.5;0,-1,-1,-3];>> A=sym(A);syms t;>> expm(A*t)ans =[ 1/2*exp(-5*t)+1/2*exp(-3*t)-1/2*t*exp(-3*t)+1/2*t^2*exp(-3*t), t*exp(-3*t)-1/2*exp(-3*t)+1/2*exp(-5*t),1/2*t^2*exp(-3*t)+1/2*t*exp(-3*t),1/2*t^2*exp(-3*t)-1/2*t*exp(-3*t)-1/2*exp(-3*t)+1/2*exp(-5*t)][ -1/2*exp(-3*t)+1/2*exp(-5*t)+1/2*t*exp(-3*t), 1/2*exp(-5*t)+1/2*exp(-3*t),1/2*t*exp(-3*t),-1/2*exp(-3*t)+1/2*exp(-5*t)+1/2*t*exp(-3*t)][ 1/2*exp(-3*t)-1/2*exp(-5*t)+1/2*t*exp(-3*t), 1/2*exp(-3*t)-1/2*exp(-5*t),1/2*t*exp(-3*t)+exp(-3*t),1/2*exp(-3*t)-1/2*exp(-5*t)+1/2*t*exp(-3*t)][ -1/2*t^2*exp(-3*t), -t*exp(-3*t), -t*exp(-3*t)-1/2*t^2*exp(-3*t), exp(-3*t)-1/2*t^2*exp(-3*t)]（2）编写M程序function F=funm(A,fun,x)[V,J]=jordan(A);v1=[0,diag(J,1)'];v2=[find(v1==0),length(v1)+1];for i=1:length(v2)-1v_lambda(i)=J(v2(i),v2(i));v_n(i)=v2(i+1)-v2(i);endm=length(v_lambda);F=sym([]);for i=1:mJ1=J(v2(i):v2(i)+v_n(i)-1,v2(i):v2(i)+v_n(i)-1);fJ=funJ(J1,fun,x);F=diagm(F,fJ);endF=V*F*inv(V);function fJ=funJ(J,fun,x)lam=J(1,1);f1=fun;fJ=subs(fun,x,lam)*eye(size(J));H=diag(diag(J,1),1);H1=H;for i=2:length(J)f1=diff(f1,x);a1=subs(f1,x,lam);fJ=fJ+a1*H1;H1=H1*H/i;endfunction A=diagm(A1,A2)A=A1;[n,m]=size(A);[n1,m1]=size(A2);A(n+1:n+n1,m+1:m+m1)=A2;>> A=[-4.5,0,0.5,-1.5;-0.5,-4,0.5,-0.5;1.5,1,-2.5,1.5;0,-1,-1,-3];>> syms x t;>> A1=funm(sym(A),sin(x*t),x)A1 =[ -1/2*sin(5*t)+1/2*sin(3*t)*t^2-1/2*cos(3*t)*t-1/2*sin(3*t), -1/2*sin(5*t)+cos(3*t)*t+1/2*sin(3*t),1/2*cos(3*t)*t+1/2*sin(3*t)*t^2,-1/2*sin(5*t)+1/2*sin(3*t)+1/2*sin(3*t)*t^2-1/2*cos(3*t)*t][ -1/2*sin(5*t)+1/2*cos(3*t)*t+1/2*sin(3*t), -1/2*sin(5*t)-1/2*sin(3*t),1/2*cos(3*t)*t, -1/2*sin(5*t)+1/2*cos(3*t)*t+1/2*sin(3*t)] [ 1/2*sin(5*t)+1/2*cos(3*t)*t-1/2*sin(3*t), 1/2*sin(5*t)-1/2*sin(3*t),-sin(3*t)+1/2*cos(3*t)*t,1/2*sin(5*t)+1/2*cos(3*t)*t-1/2*sin(3*t)][ -1/2*sin(3*t)*t^2, -cos(3*t)*t, -cos(3*t)*t-1/2*sin(3*t)*t^2, -sin(3*t)-1/2*sin(3*t)*t^2]（3）>> A=[-4.5,0,0.5,-1.5;-0.5,-4,0.5,-0.5;1.5,1,-2.5,1.5;0,-1,-1,-3];>> syms x t;>> A1=funm(sym(A),exp(x*t)*sin(x^2*exp(x*t)*t),x)A1 =[ 1/2*exp(-5*t)*sin(25*exp(-5*t)*t)+1/2*t^2*exp(-3*t)*sin(9*exp(-3*t)*t)+t*exp(-3* t)*cos(9*exp(-3*t)*t)*(-6*exp(-3*t)*t+9*t^2*exp(-3*t))-1/2*exp(-3*t)*sin(9*exp(-3* t)*t)*(-6*exp(-3*t)*t+9*t^2*exp(-3*t))^2+1/2*exp(-3*t)*cos(9*exp(-3*t)*t)*(2*exp( -3*t)*t-12*t^2*exp(-3*t)+9*t^3*exp(-3*t))-1/2*t*exp(-3*t)*sin(9*exp(-3*t)*t)-1/2*e xp(-3*t)*cos(9*exp(-3*t)*t)*(-6*exp(-3*t)*t+9*t^2*exp(-3*t))+1/2*exp(-3*t)*sin(9* exp(-3*t)*t),1/2*exp(-5*t)*sin(25*exp(-5*t)*t)+t*exp(-3*t)*sin(9*exp(-3*t)*t)+exp(-3*t)*cos(9* exp(-3*t)*t)*(-6*exp(-3*t)*t+9*t^2*exp(-3*t))-1/2*exp(-3*t)*sin(9*exp(-3*t)*t),1/2*t*exp(-3*t)*sin(9*exp(-3*t)*t)+1/2*exp(-3*t)*cos(9*exp(-3*t)*t)*(-6*exp(-3*t)* t+9*t^2*exp(-3*t))+1/2*t^2*exp(-3*t)*sin(9*exp(-3*t)*t)+t*exp(-3*t)*cos(9*exp(-3 *t)*t)*(-6*exp(-3*t)*t+9*t^2*exp(-3*t))-1/2*exp(-3*t)*sin(9*exp(-3*t)*t)*(-6*exp(-3*t)*t+9*t^2*exp(-3*t))^2+1/2*exp(-3*t)*cos(9*exp(-3*t)*t)*(2*exp(-3*t)*t-12*t^2 *exp(-3*t)+9*t^3*exp(-3*t)),1/2*exp(-5*t)*sin(25*exp(-5*t)*t)-1/2*exp(-3*t)*sin(9*exp(-3*t)*t)+1/2*t^2*exp(-3 *t)*sin(9*exp(-3*t)*t)+t*exp(-3*t)*cos(9*exp(-3*t)*t)*(-6*exp(-3*t)*t+9*t^2*exp(-3 *t))-1/2*exp(-3*t)*sin(9*exp(-3*t)*t)*(-6*exp(-3*t)*t+9*t^2*exp(-3*t))^2+1/2*exp(-3*t)*cos(9*exp(-3*t)*t)*(2*exp(-3*t)*t-12*t^2*exp(-3*t)+9*t^3*exp(-3*t))-1/2*t*ex p(-3*t)*sin(9*exp(-3*t)*t)-1/2*exp(-3*t)*cos(9*exp(-3*t)*t)*(-6*exp(-3*t)*t+9*t^2* exp(-3*t))][1/2*exp(-5*t)*sin(25*exp(-5*t)*t)+1/2*t*exp(-3*t)*sin(9*exp(-3*t)*t)+1/2*exp(-3*t) *cos(9*exp(-3*t)*t)*(-6*exp(-3*t)*t+9*t^2*exp(-3*t))-1/2*exp(-3*t)*sin(9*exp(-3*t) *t),1/2*exp(-5*t)*sin(25*exp(-5*t)*t)+1/2*exp(-3*t)*sin(9*exp(-3*t)*t),1/2*t*exp(-3*t)*sin(9*exp(-3*t)*t)+1/2*exp(-3*t)*cos(9*exp(-3*t)*t)*(-6*exp(-3*t)* t+9*t^2*exp(-3*t)),1/2*exp(-5*t)*sin(25*exp(-5*t)*t)+1/2*t*exp(-3*t)*sin(9*exp(-3*t)*t)+1/2*exp(-3*t)*cos(9*exp(-3*t)*t)*(-6*exp(-3*t)*t+9*t^2*exp(-3*t))-1/2*exp(-3*t)*sin(9*exp(-3*t) *t)][-1/2*exp(-5*t)*sin(25*exp(-5*t)*t)+1/2*t*exp(-3*t)*sin(9*exp(-3*t)*t)+1/2*exp(-3*t )*cos(9*exp(-3*t)*t)*(-6*exp(-3*t)*t+9*t^2*exp(-3*t))+1/2*exp(-3*t)*sin(9*exp(-3* t)*t),-1/2*exp(-5*t)*sin(25*exp(-5*t)*t)+1/2*exp(-3*t)*sin(9*exp(-3*t)*t),exp(-3*t)*sin(9*exp(-3*t)*t)+1/2*t*exp(-3*t)*sin(9*exp(-3*t)*t)+1/2*exp(-3*t)*cos( 9*exp(-3*t)*t)*(-6*exp(-3*t)*t+9*t^2*exp(-3*t)),-1/2*exp(-5*t)*sin(25*exp(-5*t)*t)+1/2*t*exp(-3*t)*sin(9*exp(-3*t)*t)+1/2*exp(-3*t )*cos(9*exp(-3*t)*t)*(-6*exp(-3*t)*t+9*t^2*exp(-3*t))+1/2*exp(-3*t)*sin(9*exp(-3* t)*t)][-1/2*t^2*exp(-3*t)*sin(9*exp(-3*t)*t)-t*exp(-3*t)*cos(9*exp(-3*t)*t)*(-6*exp(-3*t) *t+9*t^2*exp(-3*t))+1/2*exp(-3*t)*sin(9*exp(-3*t)*t)*(-6*exp(-3*t)*t+9*t^2*exp(-3*t))^2-1/2*exp(-3*t)*cos(9*exp(-3*t)*t)*(2*exp(-3*t)*t-12*t^2*exp(-3*t)+9*t^3*e xp(-3*t)),-t*exp(-3*t)*sin(9*exp(-3*t)*t)-exp(-3*t)*cos(9*exp(-3*t)*t)*(-6*exp(-3*t)*t+9*t^2 *exp(-3*t)),-t*exp(-3*t)*sin(9*exp(-3*t)*t)-exp(-3*t)*cos(9*exp(-3*t)*t)*(-6*exp(-3*t)*t+9*t^2 *exp(-3*t))-1/2*t^2*exp(-3*t)*sin(9*exp(-3*t)*t)-t*exp(-3*t)*cos(9*exp(-3*t)*t)*(-6 *exp(-3*t)*t+9*t^2*exp(-3*t))+1/2*exp(-3*t)*sin(9*exp(-3*t)*t)*(-6*exp(-3*t)*t+9* t^2*exp(-3*t))^2-1/2*exp(-3*t)*cos(9*exp(-3*t)*t)*(2*exp(-3*t)*t-12*t^2*exp(-3*t) +9*t^3*exp(-3*t)),exp(-3*t)*sin(9*exp(-3*t)*t)-1/2*t^2*exp(-3*t)*sin(9*exp(-3*t)*t)-t*exp(-3*t)*cos(9 *exp(-3*t)*t)*(-6*exp(-3*t)*t+9*t^2*exp(-3*t))+1/2*exp(-3*t)*sin(9*exp(-3*t)*t)*(-6*exp(-3*t)*t+9*t^2*exp(-3*t))^2-1/2*exp(-3*t)*cos(9*exp(-3*t)*t)*(2*exp(-3*t)*t-12*t^2*exp(-3*t)+9*t^3*exp(-3*t))]第二部分数学问题求解与数据处理1.(1)>> syms a t;>> f=sin(a*t)/t;laplace(f)ans =atan(a/s)(2)>> syms t a;f=t^5*sin(a*t);laplace(f)ans =60*i*(-1/(s-i*a)^6+1/(s+i*a)^6)(3)>> syms t a;>> f=t^8*cos(a*t);laplace(f)ans =20160/(s-i*a)^9+20160/(s+i*a)^92.(1)>> syms s a b;F=1/(s^2*(s^2-a^2)*(s+b));ilaplace(F)ans =1/2/a^3/b^2/(a^2-b^2)*(2*t*a*b^3+2*(1-exp(-b*t)-b*t)*a^3+( -2*a+exp(a*t)*(a-b)+exp(-a*t)*(a+b))*b^2)（2）>> syms a b s;F=sqrt(s-a)-sqrt(s-b);ilaplace(F)ans =1/2/t^(3/2)/pi^(1/2)*(exp(b*t)-exp(a*t))（3）>> syms a b s; F=log((s-a)/(s-b));ilaplace(F)ans =(exp(b*t)-exp(a*t))/t3(1)>> syms x;f=x^2*(3*sym(pi)-2*abs(x));F=fourier(f) F =-6*(4+pi^2*dirac(2,w)*w^4)/w^4>> ifourier(F)ans =x^2*(-4*x*heaviside(x)+3*pi+2*x)(2)>> syms t;f=t^2*(t-2*sym(pi))^2;F=fourier(f)F =2*pi*(dirac(4,w)-4*pi^2*dirac(2,w)+4*i*pi*dirac(3,w)) >> ifourier(F)ans =x^2*(-2*pi+x)^24.(1)>> syms k a T;f=cos(k*a*T);F=ztrans(f)(z-cos(a*T))*z/(z^2-2*z*cos(a*T)+1)>> f1=iztrans(F)f1 =cos(a*T*n)(2)>> syms k T a;f=(k*T)^2*exp(-a*k*T);F=ztrans(f)F =T^2*z*exp(-a*T)*(z+exp(-a*T))/(z-exp(-a*T))^3>> f1=iztrans(F)f1 =T^2*(1/exp(a*T))^n*n^2(3)>> syms a k T;f=(a*k*T-1+exp(-a*k*T))/a;>> F=ztrans(f)F =1/a*(a*T*z/(z-1)^2-z/(z-1)+z/exp(-a*T)/(z/exp(-a*T)-1)) >> iztrans(F)ans =(-1+a*T*n+(1/exp(a*T))^n)/a5.用数值方法求解(1)>> syms x;>> x1=solve('exp(-(x+1)^2+pi/2)*sin(5*x+2)')x1 =-2/5验证过程>> subs('exp(-(x+1)^2+pi/2)*sin(5*x+2)',x,x1)ans =(2)>> syms x;>> y1=solve('(x^2+y^2+x*y)*exp(-x^2-y^2-x*y)=0','y')y1 =(-1/2+1/2*i*3^(1/2))*x(-1/2-1/2*i*3^(1/2))*x验证过程>>y2=simple(subs('(x^2+y^2+x*y)*exp(-x^2-y^2-x*y)=0','y',y1)) y2 =(x^2+(-1/2+1/2*i*3^(1/2))^2*x^2+x^2*(-1/2+1/2*i*3^(1/2)))* exp(-x^2-(-1/2+1/2*i*3^(1/2))^2*x^2-x^2*(-1/2+1/2*i*3^(1/2) )) = 0(x^2+(-1/2-1/2*i*3^(1/2))^2*x^2+x^2*(-1/2-1/2*i*3^(1/2)))*e xp(-x^2-(-1/2-1/2*i*3^(1/2))^2*x^2-x^2*(-1/2-1/2*i*3^(1/2))) = 06.首先求出积分：>> syms x c;y=int((exp(x)-c*x)^2,x,0,1)y =-1/2-2*c+1/2*exp(2)+1/3*c^2编写一个出M文件：function y=new(c)y=-1/2-2*c+1/2*exp(2)+1/3*c^2;>> x=fminsearch('new',0)x =3.00007.编写M文件：function [c,ce]=f2(x);ce=[];c=[x(1)+x(2);x(1)*x(2)-x(1)-x(2)+1.5;-10-x(1)*x(2)];>> f=@(x)exp(x(1))*(4*x(1)^2+2*x(2)^2+4*x(1)*x(2)+2*x(2)+1); >> A=[];B=[];Aeq=[];Beq=[];xm=[-10;-10];xM=[10;10];>> x0=(xm+xM)/2;>> ff=optimset;ff.TolX=1e-10;ff.TolFun=1e-20;>> x=fmincon(f,x0,A,B,Aeq,Beq,xm,xM,@f2,ff)Maximum number of function evaluations exceeded;increase OPTIONS.MaxFunEvals.x =0.41950.4195>> i=1;x=x0;>> while (1)[x,a,b]=fmincon('f',x,A,B,Aeq,Beq,xm,xM,'f2',ff);if b>0,break;endi=i+1;end>> x,ix =1.1825-1.7398i =58书上所描述的ipslv_mex下载地址已经失效了，其他网站上也没有这个函数的下载地址，所以这个题目没有找到该函数，运行失败。

MATLAB 课后实验答案

实验一 MATLAB 运算基础1. 先求下列表达式的值，然后显示MATLAB 工作空间的使用情况并保存全部变量。

(1) 0122sin851z e =+ (2) 221ln(1)2z x x =++，其中2120.455i x +⎡⎤=⎢⎥-⎣⎦(3) 0.30.330.3sin(0.3)ln , 3.0, 2.9,,2.9,3.022a a e e a z a a --+=++=--(4) 2242011122123t t z t t t t t ⎧≤<⎪=-≤<⎨⎪-+≤<⎩，其中t =0:0.5:2.5 解： 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;z3=(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) 运算结果：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;z3=(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)z1 =0.2375z2 =0.7114 - 0.0253i 0.8968 + 0.3658i0.2139 + 0.9343i 1.1541 - 0.0044iz3 =Columns 1 through 40.7388 + 3.1416i 0.7696 + 3.1416i 0.7871 + 3.1416i 0.7913 + 3.1416iColumns 5 through 80.7822 + 3.1416i 0.7602 + 3.1416i 0.7254 + 3.1416i 0.6784 + 3.1416iColumns 9 through 120.6196 + 3.1416i 0.5496 + 3.1416i 0.4688 + 3.1416i 0.3780 + 3.1416iColumns 13 through 160.2775 + 3.1416i 0.1680 + 3.1416i 0.0497 + 3.1416i -0.0771 + 3.1416iColumns 17 through 20-0.2124 + 3.1416i -0.3566 + 3.1416i -0.5104 + 3.1416i -0.6752 + 3.1416iColumns 21 through 24-0.8536 + 3.1416i -1.0497 + 3.1416i -1.2701 + 3.1416i -1.5271 + 3.1416iColumns 25 through 28-1.8436 + 3.1416i -2.2727 + 3.1416i -2.9837 + 3.1416i -37.0245 Columns 29 through 32-3.0017 -2.3085 -1.8971 -1.5978Columns 33 through 36-1.3575 -1.1531 -0.9723 -0.8083Columns 37 through 40-0.6567 -0.5151 -0.3819 -0.2561Columns 41 through 44-0.1374 -0.0255 0.0792 0.1766Columns 45 through 480.2663 0.3478 0.4206 0.4841Columns 49 through 520.5379 0.5815 0.6145 0.6366Columns 53 through 560.6474 0.6470 0.6351 0.6119Columns 57 through 600.5777 0.5327 0.47740.4126Column 610.3388z4 =0 0.2500 0 1.2500 1.0000 2.2500 2. 已知：求下列表达式的值：(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]解：M 文件：A=[12 34 -4;34 7 87;3 65 7];B=[1 3 -1;2 0 3;3 -2 7]; A+6.*BA-B+eye(3)A*BA.*BA^3A.^3A/BB\A[A,B][A([1,3],:);B^2]运算结果：A=[12 34 -4;34 7 87;3 65 7];B=[1 3 -1;2 0 3;3 -2 7]; A+6.*BA-B+eye(3)A*BA.*BA^3A.^3A/BB\A[A,B][A([1,3],:);B^2]ans =18 52 -1046 7 10521 53 49ans =12 31 -332 8 840 67 1ans =68 44 62309 -72 596154 -5 241ans =12 102 468 0 2619 -130 49ans =37226 233824 48604247370 149188 60076678688 454142 118820 ans =1728 39304 -6439304 343 65850327 274625 343 ans =16.4000 -13.6000 7.600035.8000 -76.2000 50.200067.0000 -134.0000 68.0000109.4000 -131.2000 322.8000-53.0000 85.0000 -171.0000-61.6000 89.8000 -186.2000ans =12 34 -4 1 3 -134 7 87 2 0 33 65 7 3 -2 7ans =12 34 -43 65 74 5 111 0 1920 -5 403. 设有矩阵A和B(1) 求它们的乘积C。

(完整版)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 a z 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 分别为单位矩阵、随机矩阵、零矩阵和对角阵，试通过数值计算验证22E 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 值。

参考答案Matlab实验报告

实验一 Matlab基础知识一、实验目的：1.熟悉启动和退出Matlab的方法。

2.熟悉Matlab命令窗口的组成。

3.掌握建立矩阵的方法。

4.掌握Matlab各种表达式的书写规则以及常用函数的使用。

二、实验内容：1.求[100，999]之间能被21整除的数的个数。

（rem）2.建立一个字符串向量，删除其中的大写字母。

（find）3.输入矩阵，并找出其中大于或等于5的元素。

（find）4.不采用循环的形式求出和式6312ii=∑的数值解。

(sum)三、实验步骤：●求[100，199]之间能被21整除的数的个数。

（rem）1.开始→程序→Matlab2.输入命令：»m=100:999;»p=rem(m,21);»q=sum(p==0)ans=43●建立一个字符串向量，删除其中的大写字母。

（find）1.输入命令：»k=input('’,’s’)；Eie48458DHUEI4778»f=find(k>=’A’&k<=’Z’)；f=9 10 11 12 13»k(f)=[ ]K=eie●输入矩阵，并找出其中大于或等于5的元素。

（find）1.输入命令：»h=[4 8 10;3 6 9; 5 7 3]；»[i,j]=find(h>=5)i=3 j=11 22 23 21 32 3●不采用循环的形式求出和式的数值解。

（sum）1.输入命令：»w=1:63；»q=s um(2.^w)q=1.8447e+019- 1 -实验二 Matlab 基本程序一、实验目的：1. 熟悉Matlab 的环境与工作空间。

2. 熟悉M 文件与M 函数的编写与应用。

3. 熟悉Matlab 的控制语句。

4. 掌握if,switch,for 等语句的使用。

二、实验内容：1. 根据y=1+1/3+1/5+……+1/（2n-1）,编程求：y<5时最大n 值以及对应的y 值。

MATLAB)课后实验答案[1]

试验【2 】一 MATLAB运算基本1. 先求下列表达式的值,然后显示MATLAB工作空间的应用情形并保存全体变量.(1)0 122sin851ze =+(2)21ln( 2z x=,个中2120.455i x+⎡⎤=⎢⎥-⎣⎦(3)0.30.330.3sin(0.3)ln, 3.0, 2.9,,2.9,3.0 22a ae e az a a--+=++=--(4)2242011122123t tz t tt t t⎧≤<⎪=-≤<⎨⎪-+≤<⎩,个中t=0:0.5:2.5解：4. 完成下列操作：(1) 求[100,999]之间能被21整除的数的个数.(2) 树立一个字符串向量,删除个中的大写字母. 解：(1) 成果：(2). 树立一个字符串向量例如：ch='ABC123d4e56Fg9';则请求成果是：试验二 MATLAB矩阵剖析与处理1. 设有分块矩阵33322322E RAO S⨯⨯⨯⨯⎡⎤=⎢⎥⎣⎦,个中E.R.O.S分离为单位矩阵.随机矩阵.零矩阵和对角阵,试经由过程数值盘算验证22E R RSAO S+⎡⎤=⎢⎥⎣⎦.解: M文件如下;5. 下面是一个线性方程组：1231112340.951110.673450.52111456x x x ⎡⎤⎢⎥⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎢⎥⎢⎥⎣⎦(1) 求方程的解.(2) 将方程右边向量元素b 3改为0.53再求解,并比较b 3的变化息争的相对变化. (3) 盘算系数矩阵A 的前提数并剖析结论. 解： M 文件如下：试验三选择构造程序设计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.个中90分~100分为A,80分~89分为B,79分~79分为C,60分~69分为D,60分以下为E.请求：(1) 分离用if语句和switch语句实现.(2) 输入百分制成绩后要断定该成绩的合理性,对不合理的成绩应输出出错信息.解：M文件如下3. 硅谷公司员工的工资盘算办法如下：(1) 工作时数超过120小时者,超过部分加发15%.(2) 工作时数低于60小时者,扣发700元.(3) 其余按每小时84元计发.试编程按输入的工号和该号员工的工时数,盘算应发工资. 解：M文件下试验四轮回构造程序设计1. 依据2222211116123nπ=++++,求π的近似值.当n分离取100.1000.10000时,成果是若干？请求：分离用轮回构造和向量运算（应用sum函数）来实现.解：M文件如下：运行成果如下：2. 依据11113521yn=++++-,求：(1) y<3时的最大n值.(2) 与(1)的n值对应的y值. 解：M—文件如下：3. 斟酌以下迭代公式：1n n a x b x +=+个中a.b 为正的学数.(1) 编写程序求迭代的成果,迭代的终止前提为|x n+1-x n |≤10-5,迭代初值x 0=1.0,迭代次数不超过500次.(2) 假如迭代进程收敛于r,那么r 的精确值是24b b a-±+,当(a,b)的值取(1,1).(8,3).(10,0.1)时,分离对迭代成果和精确值进行比较.解：M 文件如下：运算成果如下;5. 若两个持续天然数的乘积减1是素数,则称这两个边疆天然数是密切数对,该素数是密切素数.例如,2×3-1=5,因为5是素数,所以2和3是密切数,5是密切素数.求[2,50]区间内：(1) 密切数对的对数.(2) 与上述密切数对对应的所有密切素数之和.解：M文件：试验五函数文件4. 设2411()(2)0.1(3)0.01f xx x=+-+-+,编写一个MATLAB函数文件fx.m,使得挪用f(x)时,x可用矩阵代入,得出的f(x)为同阶矩阵.解：运算成果：5. 已知(40)(30)(20)fyf f=+(1) 当f(n)=n+10ln(n2+5)时,求y的值.(2) 当f(n)=1×2+2×3+3×4+...+n×(n+1)时,求y的值. 解：(1)(2).试验八数据处理与多项式盘算2. 将100个学生5门功课的成绩存入矩阵P中,进行如下处理：(1) 分离求每门课的最高分.最低分及响应学生序号.(2) 分离求每门课的平均分和标准方差.(3) 5门课总分的最高分.最低分及响应学生序号.(4) 将5门课总分按从大到小次序存入zcj中,响应学生序号存入xsxh.提醒：上机调试时,为避免输入学生成绩的麻烦,可用取值规模在[45,95]之间的随机矩阵来表示学生成绩.解：M文件：运行成果：3. 某气候不雅测得某日6:00~18:00之间每隔2h的室表里温度（0C）如试验表1所示.试验表1 室表里温度不雅测成果（0C）时光h 6 8 10 12 14 16 18室内温度t1 18.0 20.0 22.0 25.0 30.0 28.0 24.0室外温度t2 15.0 19.0 24.0 28.0 34.0 32.0 30.0试用三次样条插值分离求出该日室表里6:30~18:30之间每隔2h各点的近似温度（0C）.解：M 文件：运行成果：4. 已知lgx 在[1,101]区间10个整数采样点的函数值如试验表2所示.试验表2 lgx 在10个采样点的函数值x 1 11 21 31 41 51 61 71 81 91 101lgx 0 1.0414 1.3222 1.4914 1.6128 1.7076 1.7853 1.8513 1.9085 1.9510 2.0043 试求lgx 的5次拟合多项式p(x),并绘制出lgx 和p(x)在[1,101]区间的函数曲线. 解：M 文件：5. 有3个多项式P1(x)=x4+2x3+4x2+5,P2(x)=x+2,P3(x)=x2+2x+3,试进行下列操作：(1) 求P(x)=P 1(x)+P 2(x)P 3(x).(2) 求P(x)的根.(3) 当x 取矩阵A 的每一元素时,求P(x)的值.个中：1 1.2 1.40.752 3.505 2.5A --⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦(4) 当以矩阵A 为自变量时,求P(x)的值.个中A 的值与第(3)题雷同.解：M 文件：试验九数值微积分与方程数值求解1. 求函数在指定点的数值导数.试验六高层画图操作3. 已知01ln(02x y x x ≤=⎨⎪>⎪⎩在-5≤x ≤5区间绘制函数曲线.解：M 文件：2. 用数值办法求定积分.(1) 210I π=⎰的近似值. (2) 2220ln(1)1x I dt x π+=+⎰解：M 文件：运行成果：3. 分离用3种不同的数值办法解线性方程组.6525494133422139211x y z u x y z u x y z u x y u +-+=-⎧⎪-+-=⎪⎨++-=⎪⎪-+=⎩解：M 文件：运行成果：4. 求非齐次线性方程组的通解.1234123412342736352249472x x x x x x x x x x x x +++=⎧⎪+++=⎨⎪+++=⎩解：M文件：.5. 求代数方程的数值解.(1) 3x +sin x -e x =0在x 0=1.5邻近的根.(2) 在给定的初值x 0=1,y 0=1,z 0=1下,求方程组的数值解.23sin ln 70321050y x y z x z x y z ⎧++-=⎪+-+=⎨⎪++-=⎩解：M 文件：(2). M 文件：运行成果：6. 求函数在指定区间的极值.(1) 3cos log ()x x x x x f x e ++=在(0,1)内的最小值.(2) 33212112122(,)2410f x x x x x x x x =+-+在[0,0]邻近的最小值点和最小值. 解：M文件：8. 求微分方程组的数值解,并绘制解的曲线.123213312123'''0.51(0)0,(0)1,(0)1y y y y y y y y y y y y =⎧⎪=-⎪⎨=-⎪⎪===⎩解：令y1=x,y2=y,y3=z; 如许方程变为:'''0.51(0)0,(0)1,(0)1x yz y xz z xyx y z =⎧⎪=-⎪⎨=-⎪⎪===⎩,自变量是tM 文件：试验十符号盘算基本与符号微积分一.1. 已知x=6,y=5,应用符号表达式求z =提醒：界说符号常数x=sym(‘6’),y=sym(‘5’).解：M 文件：运行成果：2. 分化因式.(1) x 4-y 4(2) 5135解：M文件：运行成果：5. 用符号办法求下列极限或导数.22sin tan 3013222220,1(1)2(1)(1)lim (2)lim sin 1cos(2)(3),',''(4),,,cos ln (5)(,)(2),,x x x x x x y xy x y x e e xa t x dA d A d A y y y A xdx dt dxdt t x x y f f x y x x ex x y +→→----==+--⎡⎤-==⎢⎥⎣⎦∂∂=-∂∂∂求已知分别求已知求解：M 文件：运行成果：6. 用符号办法求下列积分.482ln 22400(1)(2)11(3)(4)(1)1x x dx x x x dx e e dx x +∞+++++⎰⎰⎰ 解：M 文件：运行成果：。

MATLAB)课后实验答案[1]

实验一 MATLAB运算基础之欧侯瑞魂创作1. 先求下列表达式的值，然后显示MATLAB工作空间的使用情况并保管全部变量。

(1)0 122sin851ze =+(2)21ln( 2z x=+，其中2120.455i x+⎡⎤=⎢⎥-⎣⎦(3)0.30.330.3sin(0.3)ln, 3.0, 2.9,,2.9,3.0 22a ae e az a a--+=++=--(4)2242011122123t tz t tt t t⎧≤<⎪=-≤<⎨⎪-+≤<⎩，其中t解：4. 完成下列操纵：(1) 求[100,999]之间能被21整除的数的个数。

(2) 建立一个字符串向量，删除其中的大写字母。

解：(1) 结果：(2). 建立一个字符串向量例如：ch='ABC123d4e56Fg9';则要求结果是：实验二 MATLAB 矩阵分析与处理1. 设有分块矩阵33322322E R A O S ⨯⨯⨯⨯⎡⎤=⎢⎥⎣⎦，其中E 、R 、O 、S 分别为单位矩阵、随机矩阵、零矩阵和对角阵，试通过数值计算验证22E R RS A O S +⎡⎤=⎢⎥⎣⎦。

解: M 文件如下；5. 下面是一个线性方程组：(1) 求方程的解。

(2) 将方程右边向量元素b 3改为0.53再求解，并比较b 3的变更和解的相对变更。

(3) 计算系数矩阵A的条件数并分析结论。

解： M文件如下：实验三选择结构程序设计1. 求分段函数的值。

用if语句实现，分别输出x=-5.0,-3.0,1.0,2.0,2.5,3.0,5.0时的y值。

解：M文件如下：2. 输入一个百分制成绩，要求输出成绩等级A、B、C、D、E。

其中90分~100分为A，80分~89分为B，79分~79分为C，60分~69分为D，60分以下为E。

要求：(1) 分别用if语句和switch语句实现。

(2) 输入百分制成绩后要判断该成绩的合理性，对分歧理的成绩应输出出错信息。