Matlab实验题目与答案
MATLAB实验答案
MATLAB实验答案实验一：变量和矩阵的赋值、矩阵的初等运算实验目的：1. 熟悉MATLAB的工作环境2. 掌握变量的赋值、矩阵的创建。
3. 掌握矩阵和数组的初等运算。
4. 熟悉和掌握inv、mean、var、randn、rand、ones、zeros、magic、eye函数的使用。
实验内容：1. 菜单栏File→Preferences→Command Window→Font and Colors修改字体，选择Use custom font改为24。
同样方法File→Preferences→Editor/Debugger→Font and Colors 修改字体，选择Use custom font改为24。
点击Apply，OK即可。
2. 在内产生均匀分布的200个点，形成1×200向量赋给变量X。
（结果不必记录）X=linspace(0,2*pi,200);3. 在内每间隔3产生向量Y。
Y=[0:3:10];4. 给矩阵赋值，，打开workspace查看矩阵的赋值。
使用save data A B语句存储到data.mat数据文件中。
输入Clear all命令清空所有变量，再使用load data加载矩阵A 和B。
A=[1 4 8 13;-3 6 -5 -9;2 -7 -12 18];B=[5 4 3 -2;6 -2 3 -8;-1 3 -9 7];save data A B;clear all;load data;5. 产生8*6阶的正态分布随机数矩阵R，求其各列的平均值和方差，并求全体的平均值和方差。
R=randn(8,6);a=mean(R);b= var(R);c=mean(R(:));d=var(R(:));6. 模拟选号程序，现有10000人，按顺序编号为1,2,…10000号，一次随机选出一个编号，要求随机数均匀分布。
R=round(ceil(rand(10000,1)*10000));7. 产生4*6阶的均匀分布随机数矩阵R,要求其元素在1～16之间取整数值，并求此矩阵前四列组成的方阵的逆阵。
信号与系统MATLAB实验课后习题答案
a=[1,1,1]; b=[1,1]; sys=tf(b,a); t=[0:0.01:10]; figure;subplot(2,2,1); step(sys);subplot(2,2,2);x_step=zeros(size(t)); x_step(t>0)=1; x_step(t==0)=1/2; lsim(sys,x_step,t); subplot(2,2,3); impulse(sys,t);title('Impulse Response'); xlabel('Time(sec)'); ylabel('Amplitude'); subplot(2,2,4);x_delta=zeros(size(t)); x_delta(t==0)=100;[y1,t]=lsim(sys,x_delta,t); y2=y1;plot(t,y2);title('Impulse Response'); xlabel('Time(sec)'); ylabel('Amplitude');00.511.5Step ResponseTime (sec)A m p l i t u d e0.511.5Linear Simulation ResultsTime (sec)A m p l i t u d e510-0.500.51Impulse ResponseTime(sec) (sec)A m p l i t u d e510-0.50.51Impulse ResponseTime(sec)A m p l i t u d e函数int1如下：function [F,tF]=int1(f,tf,a)T=tf(2)-tf(1);F=zeros(size(tf)); tF=zeros(size(tf)); tF=tf; for n=1:length(tf)-1;F(n+1)=F(n)+T*f(n); end验证如下：t=[-1:0.01:4]; e=zeros(size(t)); e=(t>-1/2&t<1);[z,zz]=int1(e,t,-1); figure;plot(zz,z);-1-0.500.51 1.52 2.53 3.5400.511.5T1=1;N1=10000; t1=linspace(0,T1-T1/N1,N1)';f1=1-2*t1;OMG=32*pi;K1=100;omg=linspace(-OMG/2,OMG/2-OMG/K1,K1)';X1=T1/N1*exp(-j*kron(omg,t1.'))*f1;fs1=OMG/2/pi/K1*exp(j*kron(t1,omg.'))*X1;T2=5;N2=10000;t2=linspace(0,T2-T2/N2,N2)';fs2=0*t2;f2=sawtooth(t2*2*pi,0);X2=T2/N2*exp(-j*kron(omg,t2.'))*f2;fs2=fs2+OMG/2/pi/K1*exp(j*kron(t2,omg.'))*X2;figure;subplot(2,2,1);plot(omg,abs(X1),'r');xlabel('Frequency'),ylabel('Amplitude')title('单个锯齿周期幅频特性曲线');subplot(2,2,2);plot(t1,fs1,'r');xlabel('Time'),ylabel('Amplitude')title('Function after recovered');subplot(2,2,3);plot(omg,abs(X2),'r');xlabel('Frequency'),ylabel('Amplitude')title('五个锯齿周期幅频特性曲线');subplot(2,2,4);plot(t2,fs2,'r');xlabel('Time'),ylabel('Function after recovered')title('Function after recovered');-100-5005000.20.40.60.8FrequencyA m p l i t u d e单个锯齿周期幅频特性曲线00.51-1-0.500.51TimeA m p l i t u d eFunction after recovered-100-5005000.511.52FrequencyA m p l i t u d e五个锯齿周期幅频特性曲线246-2-1012TimeF u n c t i o n a f t e r r e c o v e r e dFunction after recovered4-2fsana 函数如下：function F=fsana(t,f,N)omg1=2*pi/(max(t)-min(t)); k=[-N:N]';F=1/length(t)*exp(-j*kron(k*omg1,t.'))*f; fssyn 函数如下：function f=fssyn(F,t)omg1=2*pi/(max(t)-min(t)); N=floor(length(F)/2); k=[-N:N];f=exp(j*kron(t,k*omg1))*F; 验证如下: clc clearclose allT1=1;N1=256; t=linspace(0,T1-T1/N1,N1)'; f=1-2*t;subplot(3,1,1); plot(t,f);title('验证原函数') N=25;F1=fsana(t,f,N); subplot(3,1,2); stem(abs(F1),'s');title('前N 项傅立叶级数系数幅度曲线') f2=fssyn(F1,t) ;subplot(3,1,3); plot(t,f2);xlabel('time[s]'),ylabel('Amplitude'); title('傅立叶逆变换后时域函数');00.10.20.30.40.50.60.70.80.91-101验证原函数00.20.4前N 项傅立叶级数系数幅度曲线00.10.20.30.40.50.60.70.80.91-202time[s]A m p l i t u d e傅立叶逆变换后时域函数。
matlab上机实验答案-整理版
第一次实验答案1． 设要求以0.01秒为间隔，求出y 的151个点，并求出其导数的值和曲线。
clcclearx=0:0.01:1.5;y=sqrt(3)/2*exp(-4*x).*sin(4*sqrt(3)*x+pi/3)y1=diff(y)subplot(2,1,1)plot(x,y)subplot(2,1,2)plot(x(1:150),y1)2绘制极坐标系下曲线（a,b,n 自定数据）clccleara=10;b=pi/2;n=5;theta=0:pi/100:2*pi;rho=a*cos(b+n*theta);polar(theta,rho)3. 列出求下列空间曲面交线的程序clcclearx=[-5:0.5:5];[X,Y]=meshgrid(x);z1=X.^2-2*Y.^2;z2=X.*2-Y.*3;xlabel('x')ylabel('y')zlabel('z')surf(X,Y,z1)hold onsurf(X,Y,z2)k=find(abs(z1-z2)<0.5);x1=X(k)y1=Y(k)z3=x1.^2-2*y1.^2hold onplot3(x1,y1,z3,'*')⎪⎭⎫ ⎝⎛+⋅=-334sin 234πt e y t ()θρn b a +=cos 2212y x z -=y x z 322-=4、设 ⎥⎦⎤⎢⎣⎡++=)1(sin 35.0cos 2x x x y 把x=0~2π间分为101点，画出以x 为横坐标，y 为纵坐标的曲线，要求有图形标注。
clcclearx=-2*pi:0.1: 2*pi;y=cos(x).*(0.5+sin(x)*3./(1+x.^2));plot(x,y,'b*-');title('绘图');xlabel('x 坐标');ylabel('y 坐标');legend('原函数')gtext('y=cos(x)(0.5+3*sin(x)/(1+x^2))')5、求下列联立方程的解81025695832475412743-=+-+-=-+-=++-=--+w z y x w z x w z y x w z y x clccleara=[3,4,-7,-12;5,-7,4,2;1,0,8,-5;-6,5,-2,10];b=[4,-3,9,-8];c=b/a;x=c(1,1)y=c(1,2)z=c(1,3)w=c(1,4)6. 假设一曲线数据点为x = 0:2:4*pi;y = sin(x).*exp(-x/5);试将x 的间距调成 0.1，采用不同插值方法进行插值，并通过子图的形式将不同插值结果和原始数据点绘制在同一图形窗口。
matlab考试题目及答案
matlab考试题目及答案1. 题目：编写一个MATLAB函数，实现计算并返回一个向量中所有元素的平方和。
答案：以下是一个简单的MATLAB函数，用于计算并返回一个向量中所有元素的平方和。
```matlabfunction sumOfSquares = calculateSumOfSquares(vector)sumOfSquares = sum(vector.^2);end```2. 题目：给定一个3x3的矩阵A，使用MATLAB编写代码，求出矩阵A 的转置。
答案：可以通过简单的转置操作来求得矩阵A的转置。
```matlabA = [1 2 3; 4 5 6; 7 8 9];A_transpose = A';```3. 题目：编写一个MATLAB脚本，实现对一个二维数组进行排序，并返回排序后的数组。
答案：以下是一个MATLAB脚本，用于对一个二维数组进行排序，并返回排序后的数组。
```matlabfunction sortedArray = sort2DArray(array)sortedArray = sort(array(:));end```4. 题目：给定一个向量x，使用MATLAB编写代码，计算并返回向量x的元素个数。
答案：可以通过内置函数`numel`来计算向量x的元素个数。
```matlabx = [1, 2, 3, 4, 5];numElements = numel(x);```5. 题目：编写一个MATLAB函数，实现计算并返回两个向量元素的点积。
答案：以下是一个简单的MATLAB函数，用于计算两个向量的点积。
```matlabfunction dotProduct = calculateDotProduct(vector1, vector2)dotProduct = dot(vector1, vector2);end```6. 题目：给定一个矩阵B，使用MATLAB编写代码，求出矩阵B的行列式。
MATLAB实验题答案
1求以下变量的值，并在MATLAB^验证。
( 1 ) a = 1 : 2 : 5 >> result8=a.cAttempt to reference field of non-structurea =1 3 5 (2 ) b = [ a' , a' ,a' ；a ] b =1 1 13 3 35 5 51 3 5( 3 ) c = a + b ( 2 , : )c =4 6 82 、下列运算是否合法，为什么？如合法，结果是多少？>> result2=a*bError using *Inner matrix dimensions must agree. >> result3=a+barray.>> result9=a.\b result9 =2.0000 2.0000 -0.33330.2500 0.6000 0.8333 >> result10=a.A2result10 =1 4 916 25 36>> result11=2.Aaresult11 =2 4 816 32 643、用MATLAB求解下面的的方程组。
result3 =3 6 25 >> result4=b*d result4 = 8 112 x19 15 3 2 x2 7(1)2 2 11 5 x3 11 32 13 x4 031 22 22 >> B=[4 7 -1 0]40 49 13 >> B=B'>> result5=[b;c']*d >> x=inv(A)*Bresult5 = xyz1x 2y z w 831 22 22 (2)2x y 3w 340 49 13 J-5 -8 73x 3y 5z 6w 5>> result6=a.*b>> A1=[1 1 1 0;1 2 1 -1;2 -1 0 - 3;3 3 5 -6]result6 = >> B2=[1;8;3;5]>> x2=inv(A1)*B22 8 -34 15 307 2 1 2>> result7=a./b result7 =9A 15 3 2I \2 2 11 50.5000 0.5000 -3.0000 4.0000 1.6667 1.20001 32 13>> A=[7 2 1 -2;9 15 3 -2;-2 -2 11 5;1 3 2 13] 4、已知(1)求矩阵A的秩(rank)⑵求矩阵A的行列式(determinant)⑶求矩阵A的逆(inverse)⑷求矩阵A的特征值及特征向量(eige nvalue and eige nvector) >> A3=[7 2 1 -2;9 15 3 -2;-2 -2 11 5;1 3 2 13]>> r=ran k(A3)>> b=i nv(A3)>> a=det(A3)>> [V,D]=eig(A3)7、分别用if和switch语句实现，将百分制成绩转换为成绩等级A、B C、D E。
Matlab实验习题集答案.docx
1：用以上两种形式计算56+sin^ + e 3算术运算结果。
» 5A 6+sin(pi)+exp(3)ans =1.5645e+004» x=5A 6+sin(pi)+exp(3)»A=[1 2;1 2];» B=[l 1;2 2];» C=(A<B)&(A==B)C =0 0 03：对数d = 5 + sin7用五位定点、十五位定点以及有理数形式表示出来。
»a=5+sin(7); format short,a 5.6570» a=5+sin(7);» format long,aa =» A=[l 2 3;4 15 60;7 8 9]A =12 3 415 60 7 8 9 n i n i i i5：输入矩阵U 1 1丿。
％利用MATLAB 命令直接输入矩阵OneMatrix=ones(33J)OneMatrix =8 10 9 101 1 11 1 11 1 1‘0 0 0 0 0、c养」0 0 0 0 0,» OneMatrix=ones(2,5；l);ZeroMatrix=zeros(size(OneMatrix)) ZeroMatrix =0 0 0 0 00 0 0 0 07：生成3阶魔方矩阵。
» magic(3) ans =8 1 63 5 749 28：操作符冒号”：”的应用a）步长为1的等差数列b）步长为2的等差数列c）步长为-2的等差、递减数列» 0:1:10 ans =» 0:2:10 ans =» 10:(-2):0C）取矩阵A 的A （2）, A （3）, A （4）d）取矩阵A的第一行e）取矩阵A的第三列D把矩阵A的第一行第三列元素赋值给变吊: g）把矩阵A的第二行第列元素修改为100» A=[l 23 56;sin(3) 7 9;log(2) 6 1]» A(2,3)fl 3 100) c 二,利用A 与B 生成矩阵V 6 9 0 A ans = ans =23 » A ⑵,A ⑶,A(4)ans =0.141120008059867ans =0.693147180559945 ans =23»A(1,:)ans =1 23 56» A(:,3)ans =5691» x=A(l,3)56 » A(2,1)=100L0000 100.0000 0.6931 23.0000 7.0000 6.000056.00009.00001.00001 3 A =6 9 io ：己知矩阵 L° R _(A 0) D = (A B)曲=〔0 B 丿»A=[13;6 9];B=[1 5;0 8]; » C=[AB(:,1)];» C(l,3)=100c =1 3 1006 9 0» A=[l 3;6 9];B=[1 5;0 8];D=[A B]D =公比为10°" » A=[l 23 56;sin(3) 7 9;log(2) 6 1];»A(1,:)=[]A =0.1411 7.0000 9.00000.6931 6.0000 1.000012： 1 inspace(0, 1, 6) / %给出区间[0, 1]的6个等分点数据logspace(0, 1,6) / %给出区间的6个等比点数据,»linspace(0J,6) ans =» a=l:5 b=3:2:ll c=a>*a a = 1 2b = 3 5c = 1 4» a.*bans = 3 103 4 57 9 119 16 2521 36 55■兀 . sm(k ——),k = ±2, ±1,014：计算 2 的值。
MATLAB全部实验及答案
实验项目序号1 MATLAB基本操作一、实验目的熟悉MATLAB软件环境，掌握命令窗口的使用。
二、实验基本知识1、熟悉MATLAB的桌面工具和环境：2、Start菜单&~三、主要仪器及耗材计算机四、实验内容和步骤1、MATLAB的启动及操作界面认识说明有哪几种启动的方法；操作界面有哪些窗口；对常用的窗口作些简单的介绍。
使用Windows的“开始”菜单。
运行MATLAB系统的启动程序。
{利用桌面上的快捷方式。
MATLAB主窗口命令窗口(Command Window)当前目录窗口(Current Directory)工作空间窗口(Workspace)命令历史窗口(Command History)2、下列符号中哪些是MATLAB的合法变量名用给变量赋值的方法在机器上验证你的答案（3vars、global、help、My_exp、sin、X+Y、_input、E-4、AbCd、AB_C_。
）help，My_exp，sin，AbCd，AB_C_；3、命令窗口的简单使用1）简单矩阵的输入￥1 2 3A = 4 5 67 8 92）求[12+2×(7-4)]÷32的算术运算结果4、有关向量、矩阵或数组的一些运算1）设A=15；B=20；求C=A+B与c=a+b2）设A=[1 2 3;4 5 6;7 8 9]，B=[9 8 7;6 5 4;3 2 1]；求A*B与A.*BA*B遵循矩阵运算，A.*B遵循每个对应元素之间进行相乘3）设A=10，B=20；求C=A/B与D=A\B :4）设a=[1 -2 3;4 5 -4;5 -6 7]请设计出程序，分别找出小于0的矩阵元素及其位置（单下标、全下标的形式），并将其单下标转换成全下标>> a=[1 -2 3;4 5 -4;5 -6 7];/%定义矩阵>> b = a(find(a<0));%找到a中小于零的元素b =-2-6-4>> c = find(a<0) %找到小于零元素的单下标~c =468>> [x,y] = find(a<0) %找到小于零元素的双下标>> [x,y] = ind2sub(size(a),find(a<0)) %单下标转换为双下标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=[1 2 3;3 4 2;5 2 3];>> a^2 % 矩阵运算>> a.^2 %数组与矩阵的乘法7）有一段程序设计如下，请思考并说明运行结果及其原因clearX=[1 2;8 9;3 6];X( : )解：转化X为列向量8）使用两种方法，创建一稀疏矩阵？解：1）>>a = eye(5); %创建5阶单位矩阵>>sparse(a) %将a构造为稀疏矩阵ans =(1,1) 1(2,2) 1(3,3) 1(4,4) 1(5,5) 12）>> data = [1 2 2 2];]>> ir = [1 3 2 4];>> jc = [1 2 4 1];>> s=sparse(ir,jc,data,4,4)s =(1,1) 1(4,1) 2 (3,2) 2 (2,4) 29）写出下列指令的运行结果%>> A = [ 1 2 3 ]; B = [ 4 5 6 ]; >> C = 3.^A C =3 9 27 >> D = A.^BD =1 32 7295、已知⎪⎭⎫⎝⎛+⋅=-334sin 234πt e y t 若需要计算t ∈[-1,1]，取间隔为，试计算出相对应的y 值。
MATLAB练习题和答案
MATLAB练习题和答案控制系统仿真实验Matlab部分实验结果complexC 5x6 480 double complexD 2x3 96 double complex2 选择合适的步距绘制出下面的图形(,)t∈-sin(/)t，其中111t=[-1:0.1:1];y=sin(1./t);plot(t,y)3 对下面给出的各个矩阵求取矩阵的行列式、秩、特征多项式、范数、特征根、特征向量和逆矩阵。
75350083341009103150037193......A ⎡⎤⎢⎥⎢⎥=⎢⎥-⎢⎥⎣⎦，5765710876810957910B ⎡⎤⎢⎥⎢⎥=⎢⎥⎢⎥⎣⎦12345678910111213141516C ⎡⎤⎢⎥⎢⎥=⎢⎥⎢⎥⎣⎦，33245518118575131D --⎡⎤⎢⎥-⎢⎥=⎢⎥-⎢⎥---⎣⎦A=[7.5,3.5,0,0;8,33,4.1,0;0,9,103,-1.5;0,0,3.7,19.3];B=[5,7,6,5;7,10,8,7;6,8,10,9;5,7,9,10];C=[1:4;5:8;9:12;13:1rtf6];D=[3,-3,-2,4;5,-5,1,8;11,8,5,-7;5,-1,-3,-1];det(A);det(B);det(C);det(D);rank(A);rank(B);rank(C);rank(D);a=poly(A);b=poly(B);c=poly(C);d=poly(D);norm(A);norm(B);norm(C);norm(D);[v,d]=eig(A,'nobalance');[v,d]=eig(B,'nobalance');[v,d]=eig(C,'nobalance');[v,d]=eig(D,'nobalance');m=inv(A);n=inv(B);p=inv(C);q=inv(D);4 求解下面的线性代数方程，并验证得出的解真正满足原方程。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
（温馨提示：实验课结束后，请将所有作业（题目、代码、结果）利用word 整理成一个完整的实验报告，加上封面，打印，纸质档于18周周一交）第一次上机作业目的：1. 掌握MATLAB 各种表达式的书写规则2. 运行课堂上讲过的例子，熟悉矩阵、表达式的基本操作和运算。
作业：1. 熟悉matlab 集成环境界面。
回答以下问题，并操作相关的指令：（1） 分别写出清除命令窗口和清除变量的指令。
答： clc 和clear（2）在命令行输入命令后，matlab 的搜索过程是怎样的？答： （1）检查该命令对象是不是一个变量。
（2）检查该命令对象是不是一个内部函数。
（3）检查该命令对象是否为当前目录下的程序文件。
（4）检查该命令对象是否为MATLAB 搜索路径中其他目录下的M 文件。
（3）什么是matlab 的当前工作目录？写出两种设置当前工作目录的方法？ 答： 就是matlab 当前文件读取和存储的默认路径（1）在当前目录窗口中更改（2）在MATLAB 桌面工具栏中更改（3）使用cd 命令:cd c:\mydir---将c ：\mydir 设置为当前目录（4）什么是matlab 的搜索路径？写出两种设置搜索路径的方法？答： 指Matlab 运行文件时进行搜索的目录。
（1）用path 命令设置：（2）用Set Path 对话框设置（5）help 命令和doc 命令有什么作用，它们有什么区别？答： help 命令：最基本的帮助命令，查询信息直接显示在命令窗口。
doc 命令：在帮助窗口中显示HTML 帮助文档，显示函数的详细用法及 例子，比help 命令更详细。
2. 在matlab 中输入下列表达式，并求各表达式的值，显示MATLAB 工作空间的使用情况并用两种方式保存全部变量，变量保存的文件名必须包含自己的学号后四位数：（1）)1034245.01(26-⨯+⨯=w w=sqrt(2)*(1+0.34245*10^-6)w =1.4142（2）,)tan(22ac b e abc c b a x ++-+++=ππ 其中a=3.5,b=5,c=9.8。
a=3.5;b=5;c=9.8;x=(2*pi*a+(b+c)/(pi+a*b*c)-exp(2))/(tan(b+c)+a)x =（2）])48333.0()41[(22απβππα---=y ，其中32.3=α，9.7-=β alpha=3.32;beta=-7.9;y=2*pi*alpha^2*[(1-pi/4)*beta-(0.8333-pi/4)*alpha]y =-128.4271（3）)1ln(2122t t e z t ++=，其中⎥⎦⎤⎢⎣⎡--=65.05312i tt=[2,1-3i;5,-0.65]t =2.0000 1.0000 -3.0000i5.0000 -0.6500>> z=0.5*exp(2*t)*log(t+sqrt(1+t^2))z =1.0e+004 *0.0057 - 0.0007i 0.0049 - 0.0027i1.9884 - 0.3696i 1.7706 - 1.0539i3. 已知13100.545735069.544A ⎡⎤-⎢⎥-⎢⎥=⎢⎥⎢⎥-⎣⎦1） 取出A 的第2、4、8、12、10、9个元素构成一个新的矩阵A1，其中A 的第2、4个矩阵构成新矩阵的第一行；第8、12个矩阵构成新矩阵的第二行；第10、9个矩阵构成新矩阵的第三行；A=[13,10,-0.5;4,-5,7;3,5,0;6,-9.54,4]A =13.0000 10.0000 -0.50004.0000 -5.0000 7.00003.0000 5.0000 06.0000 -9.5400 4.0000A1=A([2,4;8,12;10,9])A1 =4.0000 6.0000-9.5400 4.00007.0000 -0.5000将A1的所有元素构成一个行向量；A2=A1(:)A2 =4.0000-9.54007.00006.00004.0000-0.5000>> A2'ans =4.0000 -9.5400 7.0000 6.0000 4.0000 -0.50004. 已知⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡---=14.35454.9632053256545410778.01023A 1） 取出A 的前三行构成矩阵B ，其前两列构成矩阵C ，其右下角3×2子矩阵构成的矩阵D ，B 与C 的乘积构成矩阵E 。
A=[23,10,-0.778,0;41,-45,65,5;32,5,0,32;6,-9.54,54,3.14]A =23.0000 10.0000 -0.7780 041.0000 -45.0000 65.0000 5.000032.0000 5.0000 0 32.00006.0000 -9.5400 54.0000 3.1400B=A(1:3,:)B =23.0000 10.0000 -0.7780 041.0000 -45.0000 65.0000 5.000032.0000 5.0000 0 32.0000C=A(1:4,1:2)C =23.0000 10.000041.0000 -45.000032.0000 5.00006.0000 -9.5400D=A(2:4,3:4)D =65.0000 5.00000 32.000054.0000 3.1400E=B*CE =1.0e+003 *0.9141 -0.22391.20802.71231.1330 -0.21032）分别求E<D、E&D、E|D、~E|~D。
E<Dans =0 10 00 1E&Dans =1 10 11 1E|Dans =1 11 11 1ans =0 01 00 05. 使用函数，实现方阵A 左旋 90、右旋 270的功能，⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=131074129631185210741A 。
A=[1,4,7,10;2,5,8,11;3,6,9,12;4,7,10,13]A =1 4 7 102 5 8 113 6 9 124 7 10 13B=rot90(A)B =10 11 12 137 8 9 104 5 6 71 2 3 4C=rot90(A,-3)C =10 11 12 137 8 9 104 5 6 71 2 3 46. 请将矩阵A 的第4列所有元素删除，并使用函数将A 矩阵转变为一个行向量B ，其中135243179732A ⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦。
A=[1,3,5,2;4,3,1,7;9,3,7,2]A =1 3 5 24 3 1 79 3 7 2A(:,[4])=[]A =1 3 54 3 19 3 7第二次上机作业准备&目的：1.掌握MA TLAB 各种表达式的书写规则2.运行第二章课堂上讲过的例子，熟悉矩阵、表达式的基本操作和运算，掌握字符串、单元矩阵和结构矩阵的相关内容。
作业：1.继续完成第一次上机实验未完成的作业。
2.已知⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--=7613870451A ，⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--=023352138B ，求下列表达式的值，并注意第（2）（3）题表达式的结果有何特点：（1）B A 6+、I B A +-2（其中I 为单位阵）；A=[-1,5,-4;0,7,8;3,61,7]A =-1 5 -40 7 83 61 7>> B-[8,3,-1;2,5,3;-3,2,0]B=[8,3,-1;2,5,3;-3,2,0]B =8 3 -12 5 3-3 2 0(1)A+6*Bans =47 23 -1012 37 26-15 73 7I=eye(3,3)I =1 0 00 1 00 0 1A^2-B+Ians =-18 -217 1722 533 10921 867 526（2）A*B、A.*B、B*A、B.*A；A*Bans =14 14 16-10 51 21125 328 180A.*Bans =-8 15 40 35 24-9 122 0B*Aans =7 228 533 -1 28B.*Aans =-8 15 40 35 24-9 122 0（3）A/B、B\A、A./B、B.\A；A/Bans =1.2234 -0.92552.9787-0.9468 2.3511 -0.95744.6170 3.8723 13.8936B\Aans =-0.5106 -8.6170 -1.12770.7340 17.5745 1.8085-0.8830 -21.2128 0.4043A./BWarning: Divide by zero.ans =-0.1250 1.6667 4.00000 1.4000 2.6667-1.0000 30.5000 InfB.\AWarning: Divide by zero.ans =-0.1250 1.6667 4.00000 1.4000 2.6667-1.0000 30.5000 Inf （4）[A, B]、[A([1 3],:)；B^2]。
-1 5 -4 8 3 -10 7 8 2 5 33 61 7 -3 2 0[A([1 3],:);B^2]ans =-1 5 -43 61 773 37 117 37 13-20 1 93.产生均值为3，方差为1 的5阶正态分布的随机方阵。
A=3+sqrt(1)*randn(5)A =2.5674 4.1909 2.81333.1139 3.29441.3344 4.1892 3.7258 4.0668 1.66383.1253 2.9624 2.4117 3.0593 3.71433.2877 3.3273 5.1832 2.90444.62361.8535 3.17462.8636 2.1677 2.30824.建立一个字符串向量（要求字符串向量中必须包含自己的姓名，大小写均可），删除其中的大写字母。
并统计小写字母的个数。
A='Ma Feiyue'A =Ma Feiyue>> B=find(A>='A'&A<='Z')B =1 4>> A(:,[1,4])=[]C=find(A>='a'&A<='z')C =2 5 6 7 8 9>> length(C)ans =65.分别用两种方法建立一个结构矩阵和一个单元矩阵来存储自己班上5名同学的基本情况数据（其中1名同学必须是自己，另外一名同学必须是女生），每名学生的数据包括学号、姓名、专业和6门课程的成绩。