matlab实验五答案1
MATLAB)课后实验答案[1]
![MATLAB)课后实验答案[1]](https://img.taocdn.com/s3/m/7f76368c0b1c59eef9c7b44e.png)
实验一 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 RAO S⨯⨯⨯⨯⎡⎤=⎢⎥⎣⎦，其中E、R、O、S辨别为单位矩阵、随机矩阵、零矩阵和对角阵，试通过数值计较验证22E R RS AO S+⎡⎤=⎢⎥⎣⎦.解: M文件如下；5. 下面是一个线性方程组：(1) 求方程的解.(2) 将方程右边向量元素b3改成0.53再求解，并比较b3的变更息争的相对变更.(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) 输入百分制成绩后要判断该成绩的公道性，对不公道的成绩应输出出错信息.解：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. 考虑以下迭代公式：其中a、b为正的学数.(1) 编写程序求迭代的结果，迭代的终止条件为|x n+1-x n|≤10-5，迭代初值x0=1.0，迭代次数不超出500次.(2) 如果迭代进程收敛于r，那么r的准确值是(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 x x 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试用三次样条插值辨别求出该日室内外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 101试求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)=P1(x)+P2(x)P3(x).(2) 求P(x)的根.(3) 当x取矩阵A的每一元素时，求P(x)的值.其中：(4) 当以矩阵A为自变量时，求P(x)的值.其中A的值与第(3)题相同.解：M文件：实验九 数值微积分与方程数值求解1. 求函数在指定点的数值导数.实验六 高层画图操纵3. 已知在-5≤x ≤5区间绘制函数曲线.解：M 文件：2. 用数值办法求定积分.(1)210I π=⎰的近似值. (2) 2220ln(1)1x I dt x π+=+⎰ 解：M 文件：运行结果：3. 辨别用3种不合的数值办法解线性方程组.解：M文件：运行结果：4. 求非齐次线性方程组的通解.解：M文件：.5. 求代数方程的数值解.(1) 3x+sin x-e x=0在x0=1.5邻近的根.(2) 在给定的初值x0=1，y0=1，z0=1下，求方程组的数值解.解：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. 求微分方程组的数值解，并绘制解的曲线. 解： 令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，利用符号表达式求提示：定义符号常数x=sym(‘6’)，y=sym(‘5’).解：M文件：运行结果：2. 分化因式.(1) x4-y4(2) 5135解：M文件：运行结果：5. 用符号办法求下列极限或导数.解：M文件：运行结果：6. 用符号办法求下列积分.解：M文件：运行结果：。
MATLAB)课后实验答案[1]
![MATLAB)课后实验答案[1]](https://img.taocdn.com/s3/m/661460de43323968001c92d7.png)
实验一 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 RAO S⨯⨯⨯⨯⎡⎤=⎢⎥⎣⎦，其中E、R、O、S分别为单位矩阵、随机矩阵、零矩阵和对角阵，试通过数值计算验证22E R RS AO S+⎡⎤=⎢⎥⎣⎦.解: M文件如下；5. 上面是一个线性方程组：(1) 求方程的解.(2) 将方程右侧向量元素b3改为0.53再求解，并比较b3的变更和解的绝对变更.(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) 输入百分制成绩后要判断该成绩的合理性，对分歧理的成绩应输出出错信息.解：M文件如下3. 硅谷公司员工的工资计算方法如下：(1) 工作时数超出120小时者，超出部分加发15%.(2) 工作时数低于60小时者，扣发700元.(3) 其余按每小时84元计发.试编程按输入的工号和该号员工的工时数，计算应发工资.解：M文件下实验四 轮回结构程序设计1. 根据2222211116123n π=++++，求π的近似值.当n 分别取100、1000、10000时，结果是多少？ 请求：分别用轮回结构和向量运算（使用sum 函数）来实现.解：M 文件如下：运转结果如下：2. 根据11113521y n =++++-，求：(1) y<3时的最大n 值.(2) 与(1)的n 值对应的y 值.解：M —文件如下：3. 考虑以下迭代公式：其中a 、b 为正的学数.(1) 编写程序求迭代的结果，迭代的终止条件为|x n+1-x n |≤10-5，迭代初值x 0=1.0，迭代次数不超出500次.(2) 如果迭代过程收敛于r ，那么r 的精确值是2b -，当(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 x x 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试用三次样条插值分别求出该日室内外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 101试求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)=P1(x)+P2(x)P3(x).(2) 求P(x)的根.(3) 当x取矩阵A的每一元素时，求P(x)的值.其中：(4) 当以矩阵A为自变量时，求P(x)的值.其中A的值与第(3)题不异.解：M文件：实验九 数值微积分与方程数值求解1. 求函数在指定点的数值导数.实验六 高层绘图操纵3. 已知在-5≤x ≤5区间绘制函数曲线.解：M 文件：2. 用数值方法求定积分.(1)210I π=⎰的近似值. (2) 2220ln(1)1x I dt x π+=+⎰ 解：M 文件：运转结果：3. 分别用3种分歧的数值方法解线性方程组.解：M文件：运转结果：4. 求非齐次线性方程组的通解.解：M文件：.5. 求代数方程的数值解.(1) 3x+sin x-e x=0在x0=1.5附近的根.(2) 在给定的初值x0=1，y0=1，z0=1下，求方程组的数值解.解：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. 求微分方程组的数值解，并绘制解的曲线. 解： 令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，利用符号表达式求提示：定义符号常数x=sym(‘6’)，y=sym(‘5’).解：M文件：运转结果：2. 分解因式.(1) x4-y4(2) 5135解：M文件：运转结果：5. 用符号方法求以下极限或导数.解：M文件：运转结果：6. 用符号方法求以下积分.解：M文件：运转结果：。
(完整版)MATLAB)课后实验答案[1]
![(完整版)MATLAB)课后实验答案[1]](https://img.taocdn.com/s3/m/f6c08c5902d276a201292ece.png)
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 运算基础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 +⎡⎤=⎢⎥⎣⎦。
matlab实验指导答案详解（非常详细正确）
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教程课后实验报告题目及解答[第一至第五章]](https://img.taocdn.com/s3/m/512b4436f111f18583d05a96.png)
MATLAB教程实验报告实验项目名称实验一 Matlab基本操作学生姓名汪德旺专业班级 09数教（1）班学号 0301090131实验成绩日期一. 实验目的和要求1、了解MATLAB 的开发环境。
2、熟悉Matlab的基本操作。
3、掌握建立矩阵的方法。
4、掌握MATLAB各种表达式的书写规则以及常用函数的使用。
5、填写实验报告，实验报告文件取名为report1.doc。
6、于邮件附件形式将实验报告文件report1.doc 发到邮箱*******************，邮件主题为班级学号姓名，如：09数教1班15号张三。
二、实验内容1、先求下列表达式的值，然后显示MATLAB工作空间的使用情况并保存全部变量。
（1）0 122sin851ze =+（2）2212 1ln(0.4552i z x x+⎡⎤==⎢⎥-⎣⎦其中（3）0.30.330.3sin(0.3)ln,22a ae e az a--+=++a=-3.0,-2.9,-2.8,…, 2.8, 2.9,3.0（4）2242,011,12,0:0.5:2.521,23t tz t tt t t⎧≤<⎪=-≤<=⎨⎪-+≤<⎩其中t2.已知：1234413134787,2033657327 A 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 和B123453016678910A=,B=17-691112131415023-41617181920970212223242541311⎡⎤⎢⎥⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎢⎥⎢⎥⎣⎦（1）求它们的乘积C 。
（2）将矩阵C 的右下角3*2子矩阵赋给D 。
（3）查看MATLAB 工作空间的使用情况。
MATLAB)课后实验标准答案[1]
![MATLAB)课后实验标准答案[1]](https://img.taocdn.com/s3/m/29a8406be518964bcf847c75.png)
实验一 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 az 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)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 =实验二 MATLAB 矩阵分析与处理1. 设有分块矩阵33322322E R A O S ⨯⨯⨯⨯⎡⎤=⎢⎥⎣⎦，其中E 、R 、O 、S 分别为单位矩阵、随机矩阵、零矩阵和对角阵，试通过数值计算验证22E R RS A O S +⎡⎤=⎢⎥⎣⎦。
MATLAB全部实验及答案
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 值。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
实验五
1、编写程序，该程序在同一窗口中绘制函数在
[]0,2π之间的正弦曲线和余弦曲线，步长
为200/π，线宽为2个象素，正弦曲线设置为蓝色实线，余弦曲线颜色设置为红色虚线，两条曲线交点处，用红色星号标记并标注sin(x)=cos(x)；通过函数方式在生成的图形中添加注释，至少应包括：标题，文本注释，图例和坐标轴标注。
%x=linspace(0,2*pi,1000);x=[0:pi/200:2*pi]sinx =sin(x);cosx =cos(x);
k=find(abs(sinx-cosx)<1e-2);x1=x(k);
plot(x,sinx,'LineWidth',2)
hold on ,plot(x,cosx,'r:','LineWidth',2)hold on ,plot(x1,sin(x1),'r*')xlabel('x:(0-2\pi)');
ylabel('y:sin(x)/cos(x)');
title('正弦-余弦曲线');
text(x1+0.1,sin(x1),'sin(x)=cos(x)');legend('sin(x)','cos(x)');
1
2
34
5
6
7
x:(0-2π)
y :s i n (x )/c o s (x )
正弦-余弦曲线
2、绘制图像：双曲抛物面：22
x y z =-
，1616x -<<，44y -<<，并对绘制的双曲抛物面尝试进行视点控制。
[X,Y]=meshgrid(-16:0.4:16,-4:0.1:4);Z =X.^2/16-Y.^2/4;
subplot(1,3,1),plot3(X,Y,Z),view(0,180),title('azimuth =0,elevation =180');
subplot(1,3,2),plot3(X,Y,Z),view(-37.5,-30),title('azimuth =
-37.5,elevation =-30');
subplot(1,3,3),plot3(X,Y,Z),view([3,3,2]),title('viewpoint=[3,3,1]');
3、表中列出了4个观测点的6次测量数据，将数据绘制成为分组形式和堆叠形式的条形图。
>>y=[3696;6774;7323;4252;2487;8744];
>>bar(y)
>>bar(y,’stack’)
第1次
第2次第3次第4次第5次第6次观测点1367428观测点2673247观测点3972584观测点4
6
4
3
2
7
4
4、某次考试优秀、良好、中等、及格和不及格的人数分别为7、17、23、19和5，，绘制饼图进行成绩统计分析并将优秀成绩分离出来。
x=[71723195];L=[10000];colormap(hot);
pie(x,L),title('学生成绩统计饼图');legend('优秀','良好','中等','及格','不及
格');
24%
32%
27%
7%
学生成绩统计饼图 优秀良好中等及格不及格
5、2
2
y x
xe
z --=，当x 和y 的取值范围均为-2到2时，用建立子窗口的方法在同一个图形窗
口中绘制出三维线图、网格图、表面图和带渲染效果的表面图。
>>[x,y]=meshgrid([-2:.2:2]);>>z=x.*exp(-x.^2-y.^2);>>mesh(x,y,z)>>subplot(2,2,1),plot3(x,y,z)>>title('plot3(x,y,z)')>>subplot(2,2,2),mesh(x,y,z)>>title('mesh (x,y,z)')>>subplot(2,2,3),surf(x,y,z)>>title('surf (x,y,z)')>>subplot(2,2,4),surf(x,y,z),shading interp >>title('surf (x,y,z),shading interp')
6、用sphere函数产生球表面坐标，绘制不通明网格图、透明网格图、表面图和带剪孔的
表面图。
>>[x,y,z]=sphere(30);
>>mesh(x,y,z)
>>mesh(x,y,z),hidden off
>>surf(x,y,z)
>>z(18:30,1:5)=NaN*ones(13,5); >>surf(x,y,z)
axis square
axis off。