下载文档原格式
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 的工作环境，学会使用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 运算基础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）课后实验答案[1]](https://img.taocdn.com/s1/m/306ff29b690203d8ce2f0066f5335a8102d266cb.png)
（完整版）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 az 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 分别为单位矩阵、随机矩阵、零矩阵和对⾓阵，试通过数值计算验证2 2E 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 值。
解：M ⽂件如下：2. 输⼊⼀个百分制成绩，要求输出成绩等级A、B、C、D、E。
实验一 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)=[]实验二 MATLAB 矩阵分析与处理1. 设有分块矩阵33322322E R A O S ⨯⨯⨯⨯⎡⎤=⎢⎥⎣⎦，其中E 、R 、O 、S 分别为单位矩阵、随机矩阵、零矩阵和对角阵，试通过数值计算验证22E R RS A O S +⎡⎤=⎢⎥⎣⎦。
![MATLAB)课后实验答案[1]](https://img.taocdn.com/s1/m/adaedec65727a5e9846a61d9.png)
实验一 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 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 OS +⎡⎤=⎢⎥⎣⎦。
解： M 文件如下；5。
下面是一个线性方程组：1231112340.951110.673450.52111456x x x ⎡⎤⎢⎥⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎢⎥⎢⎥⎣⎦ch =123d4e56g9（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。
实验一 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 运算基础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。
2,(1)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(3) 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 ];B=A([1,2],[3]),C=A(2:end, : )(5) M=magic(4),M(:,4)=[]3,(1) p1=[1 0 2 4];PS1=poly2str(p1,'x')r=roots(p1)(2) p=poly(A)polyval(p,20)poly2str(p,’A’)4,(1) t=0:.1:2*pi;y=cos(t);plot(t,y),grid(2) t=0:.1:2*pi;y1=cos(t-0.25);y2=sin(t-0.25);plot(t,y1，t,y2)gridt=0:.01:4*pi;y=10*sin(t);plot(t,y,'r:+')gridaxis([0,10,-15,15])xlabel('x','FontSize',16)ylabel('y','FontSize',16)title('正弦函数')6，(1)sum=0;i=1;while(sum<2000)sum=sum+i;i=i+1;endn=i-2{ n=0;s=0;while s<2000n=n+1;s=s+n;endn=n-1}(2)解法1：function [ s i ] = mitifang1(n)s=0;for i=0:ns=s+2^i; if(i>n), break;endends,n解法2：function [ s i ] = mitifang2(n)s=0;i=0;while(i<=n) ,s=s+2^i;i=i+1;ends,n(3)reply=input('ÇëÊäÈë:','s');while reply=='y'||reply=='Y'||reply=='n'||reply=='N'if (reply=='y'||reply=='Y')x=1;disp('x=')disp(x)elseif (reply=='n'||reply=='N')x=0;disp('x=')disp(x)elsedisp('shu ru you wu')endreply=input('ÇëÊäÈë:','s');enddisp('shu ru you wu ')试验二1,(1)num=5*[1 5 6];den=[1 6 10 8];G=tf(num,den)step(G)impulse(G)num=5*[1 5 6];den=[1 6 10 8];G=tf(num,den)[A B C D]=tf2ss(num,den)x=[1 1 0]'initial(A,B,C,D,x)(2)wn=[2:2:12];z=0.7;t=0:0.1:12;hold onfor i=1:length(wn)Gc=tf(wn(i)^2,[1,2*z*wn(i),wn(i)^2]); step(Gc,t)endgrid onhold on(3)z=[0.2:0.2:2.0];wn=6;t=0:0.1:12;hold onfor i=1:length(z)Gc=tf(wn^2,[1,2*z(i)*wn,wn^2]); step(Gc,t)endgrid onhold on2,(1)wn=[2:2:12];z=0.7;hold onfor i=1:length(wn)num=wn(i)^2;den=[1,2*z*wn(i),wn(i)^2]; bode(num,den);endgrid onhold on(2)z=[0.2:0.2:2.0];wn=6;hold onfor i=1:length(z)num=wn^2;den=[1,2*z(i)*wn,wn^2];bode(num,den);endgrid onhold on3,num=[1];den=conv([1 0],conv([1 1],[1 2])); sys=tf(num,den);rlocus(sys)[x,y]=ginput(3);p=x+i*yK=rlocfind(sys,p)实验三1，（1）k=100;p=[0,-1,-20];z=-2;sys=zpk(z,p,k)sys=tf(sys)G_c=feedback(sys,1);%转换为闭环传递函数roots(G_c.den{1})%求闭环特征方程的根[rtab,msg]=routh(G_c.den{1})(2)sys=tf([1,3],conv(conv([1,0],[1,5]),conv([1,6],[1,2,2]))) rlocus(sys)[K,poles]=rlocfind(sys)(3)sys1=tf(2.7,[1 5 4 0])margin(sys1);%画出波特图[Gm,Pm,Wcg,Wcp]=margin(sys1)G_c=feedback(sys,1)step(G_c)%验证阶跃响应是否稳定Grid2,[A,B,C,D]=tf2ss([1,-1],[1,10,27,18])%取0,1时同理Tc=ctrb(A,B)rank(Tc)To=obsv(A,C)rank(To)3, b=[2 3 4 0];a=[1 3 3 2];n=2^3;Fs=40;[h,f]=freqz(b,a,n,Fs);plot(f,abs(h))grid ontitle('离散系统幅频特性曲线')p=angle(h);ph=p*180/pi;plot(f,ph)grid ontitle('离散系统相频特性曲线')t=1;dbode(b,a,t)[mag,phase,w]=dbode(b,a,t)figure(7)[gm,pm,wg,wp]=margin(mag,phase,w) dnyquist(b,a,t)。
实验一 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 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 OS +⎡⎤=⎢⎥⎣⎦。
解： M 文件如下；5。
下面是一个线性方程组：1231112340.951110.673450.52111456x x x ⎡⎤⎢⎥⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎢⎥⎢⎥⎣⎦ch =123d4e56g9（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。
实验一1.（1）z1=2*sin(85*pi/180)/(1+exp(2))（2）>> x=[2,1+2i;-0.45,5];>> z2=(1/2)*log(x+sqrt(1+x*x))(3)>> 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)(4)t=0:0.5:2.5;>>z4=(t.*t).*(t>=0&t<1)+(t.*t-1).*(t>=1&t<2)+(t.*t-2.*t+1).*(t>=2& t<3)2.(1)>> A=[12,34,-4;34,7,87;3,65,7];>> B=[1,3,-1;2,0,3;3,-2,7];>> A+6*B(2)>> A*B(3)>> A^3(4)>> A/B(5)>> [A,B]3.(1)>>A=[1,2,3,4,5;6,7,8,9,10;11,12,13,14,15;16,17,18,19,20;21,22,23,2 4,25];>> B=[3,0,16;17,-6,9;0,23,-4;9,7,0;4,13,11];>> C=A*B(2)>> D=C(3:end,2:end)Name value classA <5×5double> doubleB <5×3double> doubleC <5×3double> doubleD [520 397;705 557;…] double4.(1)>> X=100:999;>> length(find(~rem(X,21)))(2)>> y='ahHGJjskjJjsdsaSa';>> M=y(find(y<='z'&y>='a'))实验三1 x=input('请输入x的值')if x<0&x~=-3y=x^2+x-6;elseif x>=0&x<5&x~=2&x~=3y=x^2-5*x+6;y=x^2-x-1;endy2（1）x=input('请输入成绩') if x>=90&x<=100disp('等级为A')elseif x>=80&x<=89 disp('等级为B')elseif x>=70&x<=79 disp('等级为C')elseif x>=60&x<=69 disp('等级为D')elseif x>=0&x<=59 disp('等级为E')elsedisp('成绩有误')end(2)x=input('请输入成绩') switch(x)case{90,100}disp('等级为A')case{80,89}disp('等级为B')case{70,79}disp('等级为C')case{60,69}disp('等级为D')case{0,59}disp('等级为E') otherwisedisp('成绩错误')end>>3Num=input('输入员工的工号:');Time=input('该工号员工的工作时数:'); salary=input('该员工底薪:');if Time>=120salary=salary*(1+0.15);elseif Time<=60salary=salary-700;elseif Time>60&Time<120salary=Time*84;endsalary=salary4:x=input('请输入x');y=input('请输入y');n=input('请输入运算符','s');if n=='+'m=x+y;elseif n=='-'m=x-y;elseif n=='*'m=x*y;elsem=x/y;endm5:a=rand(5,6)n=input('请输入要输出的第二行数'); if n>5a(5,:)elsea(n,:)end实验四1：程序：n=input('输入n:'); for a=1:n;b(a)=1/(a^2);endpi=sqrt(6*sum(b))2：.程序： sum=0;for n=1:200;b(n)=1/(2*n-1);sum=sum+b(n);if(sum>=3)break;endendn-14.程序：f(1)=1;max=0;f(2)=0;min=0;f(3)=1;c=0;zhengshu=0;fushu=0;zero=0;for n=4:100f(n)=f(n-1)-2*f(n-2)+f(n-3); endfor m=1:100if(f(m)>max)max=f(m);endif(f(m)<min)min=f(m);endc=sum(f);if(f>0)zhengshu=zhengshu+1; endif(f(m)<0)fushu=fushu+1;endif(f(m)==0)zero=zero+1;endif(f(m)>0)zhengshu=zhengshu+1;endif(f(m)<0)fushu=fushu+1;endif(f(m)==0)zero=zero+1;enddisp('最大值') maxdisp('最小值')mindisp('和')cdisp('正数的个数') zhengshudisp('负数的个数') fushudisp('零的个数')zero5p=0;i=0;h=0;for m=2:49a=m;b=m+1;c=a*b-1;for k=1:cif rem(c,k)==0i=i+1;endendif i==2disp(['亲密数对是','(',num2str(a),',',num2str(b),')']) p=p+1;h=[h,c];endi=0;endsum=sum(h);disp(['亲密数对的个数是：',num2str(p)])disp(['亲密数对的和是：',num2str(sum)])实验五1，函数文件：文件名：Yunsuan程序：function[z,d,s,c]=yunsuan(a,b)z=exp(a+b*i);d=log(a+b*i);s=sin(a+b*i);c=cos(a+b*i);disp(['复数',num2str(a),'+',num2str(b),'i的指数为:',num2str(z),'，对数为：',num2str(d),'，正弦为：',num2str(s),'，余弦为：',num2str(c)])函数调用：程序：clearx=input('请输入复数实部的值');y=input('请输入复数虚部的值');[x,y]=yunsuan(x,y);2函数文件：文件名：zhu程序：function[p,q,w]=zhu(p,q,w)A=[p*cos(w),-p,-sin(w),0;p*sin(w),0,cos(w),0;0,q,-sin(w),0;0 ,0,-cos(w),1];B=[0;9.8*p;0;9.8*q];X=A\B函数调用：程序：clearm1=input('m1=');m2=input('m2=');th=input('theta=');theta=pi*th/180;[m1,m2,theta]=zhu(m1,m2,theta);5(1)function[f]=fn(n);f=n+10*log(n*n+5);[f40]=fn(40);[f30]=fn(30);[f20]=fn(20);y=f40/(f30+f20)(2)function[f]=fn(n);f=0;for i=1:nf=f+i*(i+1);end[f40]=fn(40);[f30]=fn(30);[f20]=fn(20);y=f40/(f30+f20)实验六2_1clear;x=linspace(0,2*pi,101);y1=x.^2;y2=cos(2.*x);y3=y1.*y2;plot(x,y1,'b-',x,y2,'r:',x,y3,'g-.'); %y1蓝色实线，y2红色虚线，y3绿色点画线2_2subplot(2,2,1); %分四个子图(先画2行2列第1块)plot(x,y1);subplot(2,2,2);plot(x,y2),subplot(2,2,3);plot(x,y3);2_3 ()subplot(3,4,1); %y1的四种图形bar(x,y1);subplot(3,4,2);stairs(x,y1),subplot(3,4,3);stem(x,y1);subplot(3,4,4);fill(x,y1,'b');subplot(3,4,5); %y2bar(x,y2); %条形图subplot(3,4,6);stairs(x,y2), %阶梯图subplot(3,4,7);stem(x,y2); %杆图subplot(3,4,8);fill(x,y2,'b'); %填充图，注意必须加填充颜色subplot(3,4,9); %y3bar(x,y3);subplot(3,4,10);stairs(x,y3),subplot(3,4,11);stem(x,y3);subplot(3,4,12);fill(x,y3,'b');3clear;x=-5:0.1:5;if x<=0y=(x+sqrt(pi)/exp(2));elsey=0.5.*log(x+sqrt(1+x.^2));endplot(x,y);4M文件，假设文件名为Untitled6 a=input('a=');b=input('b='); %b单位为pi/4 b=b*pi/4;n=input('n=');q=linspace(-2*pi,2*pi,100); p=a*sin(b+n*q);plot(q,p);hold on; %保持图形命令窗口调用情况>> Untitled6a=1b=1n=1>> Untitled6a=2b=2n=25.程序：x=linspace(-5,5,21);y=linspace(0,10,31);[x,y]=meshgrid(x,y);z=cos(x).*cos(y).*exp(-0.25*sqrt(x.^2+y.^2));subplot(1,2,1);title('surf(x,y,z)');surf(x,y,z);subplot(1,2,2);title('surfc(x,y,z)');surfc(x,y,z);实验七701clear allhf=figure('Color',[1,0,0],'KeyPressFcn','disp(''Left Button Pressed'')');702clear allx=linspace(-2*pi,2*pi,500);y=x.^2.*exp(2*x);h=plot(x,y);set(h,'Color','r','LineStyle',':','LineWidth',5)title('下图是y=x^2e^{2x}曲线的图像');703clear alla=linspace(-2*pi,2*pi,40);b=linspace(-2*pi,2*pi,40);[x,t]=meshgrid(a,b);v=10*exp(-0.01*x).*sin(2000*pi*t-0.2*x+pi);axes('view',[-37.5,30]);h=surface(x,t,v,'FaceColor','w','EdgeColor','flat'); grid on;title('函数图像如下');set(h,'FaceColor','flat');705%实验七底层绘图操作clear all;[x,y,z]=cylinder(3,500);%cylinder是生成柱体的函数surf(x,y,z)title('实验七第五题圆柱体的光照和材料处理');xlabel('X-axis');ylabel('Y-axis');zlabel('Z-axis');axis([-5,5,-5,5,0,1])grid off;light('Color','r','Position',[-4,0,0],'Style','infinite');shading interp; material shiny; view(0,10); lighting phong; axis off;。
