哈工大matlab选修课试题
matlab考试题库及答案大学
matlab考试题库及答案大学# MATLAB考试题库及答案大学一、选择题1. MATLAB中用于生成0到1之间均匀分布随机数的函数是： - A. rand- B. randi- C. randperm- D. randn答案：A2. 下列哪个命令可以用于绘制函数图像？- A. plot- B. text- C. title- D. xlabel答案：A3. MATLAB中，以下哪个不是矩阵的属性？- A. size- B. length- C. rank- D. transpose答案：D4. 以下哪个函数可以用于求解线性方程组？- A. solve- B. linsolve- C. linprog- D. fsolve答案：A5. MATLAB中，用于计算矩阵特征值的函数是：- A. eig- B. eign- C. eigen- D. eigenvalue答案：A二、填空题1. MATLAB的基本数据单位是________。
答案：矩阵2. 使用MATLAB进行数值计算时，可以利用________来存储数据。
答案：变量3. MATLAB中的向量可以是________或________。
- 答案：行向量；列向量4. 矩阵的转置可以通过________操作实现。
答案：.'5. MATLAB中，使用________函数可以计算矩阵的行列式。
答案：det三、简答题1. 简述MATLAB中使用循环结构的注意事项。
答案：在MATLAB中使用循环结构时，应注意以下几点：确保循环体内部逻辑正确，避免无限循环；使用for循环时，循环变量的初始化和步长设置要合理；使用while循环时，循环条件要明确，确保循环能够正常退出。
2. 描述MATLAB中如何实现数组的多维索引。
答案：在MATLAB中，多维数组的索引可以通过使用圆括号来实现，每个维度的索引用逗号分隔。
例如，对于一个三维数组A，可以使用A(i,j,k)来访问第i行、第j列、第k层的元素。
大学matlab考试题及答案
大学matlab考试题及答案大学MATLAB考试题及答案一、选择题1. MATLAB的全称是什么？A. Matrix LaboratoryB. Microprocessor Application Tool SetC. Microsoft Advanced Tool SetD. Microprocessor Application Test System答案：A2. 在MATLAB中，以下哪个命令用于绘制三维图形？A. plotB. scatterC. surfD. bar答案：C3. MATLAB中用于求解线性方程组的函数是？A. solveB. linsolveC. linprogD. fsolve答案：A二、简答题1. 简述MATLAB的基本数据类型有哪些，并给出至少两个每种类型的示例。
答案：MATLAB的基本数据类型包括数值数组、字符数组和单元数组。
数值数组可以是向量、矩阵或多维数组。
例如，向量 `v = [1 2 3]`，矩阵 `M = [1 2; 3 4]`。
字符数组是由单引号或双引号括起来的字符序列，如 `C = 'Hello'`。
单元数组可以包含不同类型的数据，如`{1, 'text', [1; 2; 3]}`。
2. 描述如何在MATLAB中实现矩阵的转置和翻转。
答案：矩阵的转置可以通过 `'T'` 来实现，例如 `A'` 表示矩阵A 的转置。
矩阵的翻转可以通过 `flip` 函数实现，例如 `flip(A)` 可以翻转矩阵A的所有行和列，`flipud(A)` 仅翻转矩阵A的行，而`fliplr(A)` 仅翻转矩阵A的列。
三、编程题1. 编写一个MATLAB函数，该函数接受一个向量作为输入，并返回向量中所有元素的和以及平均值。
```matlabfunction [sumVal, avgVal] = calculateSumAndAverage(V)sumVal = sum(V);avgVal = mean(V);end```2. 设计一个MATLAB脚本来解决以下问题：给定一个3x3的矩阵，找出其中最大的元素，并显示其位置。
哈工大MATLAB选修课第二次matlab作业
1. 表1 用三次样条方法插值计算0-90 度内整数点的sin 值和0-75 度内整数点的正切值，然后用5 次多项式拟合方法计算相同的函数。
a（度）0 15 30 45 60 75 90Sin（a）0 0.2588 0.5000 0.7071 0.8660 0.9659 1.0000tan（a）0 0.2679 0.5774 1.0000 1.7320 3.732解：分别对应的程序如下：正弦函数：x = pi*(0:90)/180;y = sin(x);xx = pi*(0:.25:90)/180;yy = spline(x,y,xx);plot(x,y,'o',xx,yy)正切函数：x = pi*(0:75)/180;y = tan(x);xx = pi*(0:.25:75)/180;yy = spline(x,y,xx);plot(x,y,'o',xx,yy)正弦拟合：figurex=pi*(0:15:90)/180;y=[0,0.2588,0.5,0.7071,0.866,0.9659,1.0]; xx=pi*(1:0.05:90)/180;p2=polyfit(x,y,5);yy=polyval(p2,xx);plot(x,y,'-ro',xx,yy);正切拟合：figurex=pi*(0:15:75)/180;y=[0,0.2679,0.5774,1,1.732,3.732];xx=pi*(1:0.05:75)/180;p2=polyfit(x,y,5);yy=polyval(p2,xx);plot(x,y,'-ro',xx,yy);legend('描点显示','五次拟合')2. 采用最近点法、线性法和3 次样条法插值计算1-100 整数间平方根n 1 4 9 16 25 36 49 64 81 100Sqtr（n）1 2 3 4 5 6 7 8 9 10解：程序如下：x=[1,4,9,16,25,36,49,64,81,100];y=[1,2,3,4,5,6,7,8,9,10];xx=1:100;yy=interp1(x,y,xx)subplot(2,2,1)plot(x,y,'-ro',xx,yy,'dr');title('线性法');subplot(2,2,2);y2=interp1(x,y,xx,'nearest');plot(x,y,'-ro',xx,y2,'dr');title('最近点法')subplot(2,2,3);y3=interp1(x,y,xx,'spline');plot(x,y,'-ro',xx,y3,'dr');title('3次样条法')仿真的结果：3. 已知p(x)=2x^4-3x^3+5x+13，求p（x）的全部根，由方程p（x）=0 的根构造一个多项式f（x），并和p（x）比较。
大学matlab考试题及答案
大学matlab考试题及答案一、选择题（每题2分，共10分）1. MATLAB中用于创建向量的命令是？A. vectorB. arrayC. linspaceD. [ ]答案：D2. 以下哪个函数用于计算矩阵的行列式？A. detB. rankC. invD. eig答案：A3. 在MATLAB中，以下哪个命令用于绘制三维曲面？A. plot3B. surfC. meshD. contour答案：B4. MATLAB中，以下哪个函数用于求解非线性方程？A. fsolveB. fzeroC. linsolveD. solve答案：A5. 如何在MATLAB中创建一个全为1的3x3矩阵？A. ones(3,3)B. zeros(3,3)C. eye(3,3)D. rand(3,3)答案：A二、填空题（每题3分，共15分）1. 在MATLAB中，使用________函数可以计算矩阵的逆。
答案：inv2. 要将向量[1, 2, 3]转置成行向量，可以使用________运算符。
答案：.'3. MATLAB中，________函数用于计算两个矩阵的点乘。
答案：dot4. 若要在MATLAB中生成一个从0到1的100个等间隔的数，可以使用________函数。
答案：linspace5. 使用________函数可以在MATLAB中创建一个随机数矩阵。
答案：rand三、简答题（每题5分，共20分）1. 描述MATLAB中矩阵的基本操作有哪些？答案：矩阵的基本操作包括矩阵的创建、转置、加法、减法、乘法、点乘、求逆、行列式计算、特征值和特征向量的计算等。
2. MATLAB中如何实现循环结构？答案：MATLAB中可以使用for循环、while循环和switch-case结构来实现循环结构。
3. 简述MATLAB中函数定义的基本语法。
答案：MATLAB中函数定义的基本语法是：function [输出参数列表] = 函数名(输入参数列表)函数体end4. MATLAB中如何实现条件判断？答案：MATLAB中可以使用if-else语句、switch-case语句来实现条件判断。
matlab试题及答案
matlab试题及答案# MATLAB试题及答案一、选择题1. MATLAB的基本数据单位是：A. 矩阵B. 向量C. 标量D. 数组答案：A2. 下列哪个命令可以用来绘制函数图形？A. `plot`B. `graph`C. `draw`D. `chart`答案：A3. MATLAB中，以下哪个是正确的矩阵转置操作？A. `transpose(A)`B. `A'`C. `A^T`D. `flip(A)`答案：B二、简答题1. 简述MATLAB中矩阵的基本操作。
答案：在MATLAB中，矩阵是最基本的数据结构，可以进行加、减、乘、除等基本运算。
矩阵的创建可以使用方括号`[]`，例如`A = [1 2;3 4]`。
矩阵的转置使用单引号`'`，例如`A'`。
矩阵的求逆使用`inv`函数，例如`inv(A)`。
2. MATLAB中如何实现循环结构？答案：MATLAB中实现循环结构主要有两种方式：`for`循环和`while`循环。
`for`循环用于已知迭代次数的情况，例如：```matlabfor i = 1:5disp(i);end````while`循环用于迭代次数未知的情况，例如：```matlabi = 1;while i <= 5disp(i);i = i + 1;end```三、计算题1. 给定矩阵A和B，请计算它们的乘积C，并求C的行列式。
A = [1 2; 3 4]B = [5 6; 7 8]答案：首先计算矩阵乘积C：```matlabC = A * B;```然后计算C的行列式：```matlabdetC = det(C);```结果为：```matlabC = [19 22; 43 50]detC = -16```2. 编写一个MATLAB函数，计算并返回一个向量的范数。
答案：```matlabfunction norm_value = vector_norm(v)norm_value = norm(v);end```四、编程题1. 编写一个MATLAB脚本，实现以下功能：- 随机生成一个3x3的矩阵。
哈工大matlab春季学期考查题答案
2012年春季学期 MATLAB 课程考查题一、必答题： 1，简述MATLAB组成部分？答：（1）开发环境（development environment）（2）matlab 数据函数库（Math Function Library）（3）matlab 语言(matlab language) （4）句柄图形（handle graphics）（5）应用程序接口(Applied Function Interface) 2，说明使用M文件编辑/调试器的方法和优点？答：（1）方法：a. 打开文件编辑器(Opening the file editor): 最简单的方法是在操作桌面的工具栏上选择新建文件键（New M- File）或打开已有文件键(Open File),也可以在命令窗口输入命令edit建立新文件或输入edit filename, 打开名为filename 的M文件；b. 编写程序内容(Compose the program content): 编写新的文件或修改已有文件；c. 保存文件(save file): 文件运行前必须完成保存操作，与一般的文件编辑保存操作相同；d. 运行文件(running file): 在命令窗口输入文件名即可运行。
如要在编辑器中直接完成运行，可在编辑器的Debug菜单下save and run选项，或按Run快捷键，最快捷的方法是直接按F5键执行运行。
（2）优点：a.功能强大的数值运算：在MATLAB环境中，有超过500种数学、统计、科学及工程方面的函数可使用，函数的标示自然，使得问题和解答像数学式子一般简单明了，让使用者可全力发挥在解题方面，而非浪费在电脑操作上。
b.先进的资料视觉化功能： MATLAB的物件导向图形架构让使用者可执行视觉数据分析，并制作高品质的图形，完成科学性或工程性图文并茂的文章。
c.高阶但简单的程序环境：做为一种直译式的程序语言，MATLAB让使用者在短时间内写完程序，所花的时间约为用FORTRAN 或C 的几分之一，而且不需要编译(compile)及联结(link) 即能执行，同时包含了更多及更容易使用的内建功能。
MATLAB选修课作业
目录第二章 MATLAB数据及其运算1.下列命令执行后 L2、L3、L4的值分别是多少？2.建立单元矩阵A并回答问题第三章 MATLAB矩阵分析与处理1.求(x+y)5的展开式：pascal(6)2.求矩阵A的特征值和相应的特征向量第四章 MATLAB程序设计1.已知：s=1+2+2^2+2^3+…2^63，分别用循环结构和调运MATLAB的SUM函数求S的值；2.分别用if和switch语句实现以下计算，其中a,b,c的值从键盘输入;第五章 MATLAB绘图1.绘制下列曲线2．绘制下列极坐标图第六章 MATLAB数值计算1.求数值积分2．绘制下列极坐标图第七章 MATLAB符号计算1.分解因式2.求非线性方程的符号1第二章 MATLAB数据及其运算1.下列命令执行后 L2、L3、L4的值分别是多少？A=1:9;B=10-A;L2=A<=5;L3=A>3&A<7:；L4=find(A>3&A<7);L2 =1 1 1 1 1 0 0 0 0L3 =0 0 1 1 1 0 0 0L4 =4 5 62.建立单元矩阵A并回答问题A{1,1}=1;A{1,2}='brenden';A{2,1}=reshape(1:9,3,3);A{2,2}={12,34,2;54,21,3;4,23,67};(2）A(2)和A(4)的值分别是多少？(3)A(3)=[ ]和A{3}=[ ]的值分别是多少？A(2)ans = [3x3 double]A(4)ans = {3x3 cell}A(3)=[ ]A= [1] [3x3 double] {3x3 cell}> A{3}=[ ]A = [1] [3x3 double] [ ]第三章 MATLAB矩阵分析与处理1.求(x+y)5的展开式：pascal(6)2.求矩阵A的特征值和相应的特征向量A=[1 1 0.5;1 1 0.25;0.5 0.25 2]; [V,D]=eig(A)V =0.7212 0.4443 0.5315-0.6863 0.5621 0.4615-0.0937 -0.6976 0.7103D =-0.0166 0 00 1.4801 00 0 2.53653第四章 MATLAB程序设计1.已知：s=1+2+2^2+2^3+…2^63，分别用循环结构和调运MATLAB的SUM函数求S的值；s=2^0;n=63;for i=1:63;s=s+2^i;endss =1.8447e+019n=63;i=0:n;f=2.^i;s=sum(f)s =1.8447e+0192.分别用if和switch语句实现以下计算，其中a,b,c的值从键盘输入;a=8;b=5;c=10;x=0.5:0.1:5.5;if x>=0.5&x<1.5y=a*x^2+b*x+c;endif x>=1.5&x<3.5y=a*sin(c*b)+x;endif x>=3.5&x<5.5y=log(b+c/x);endyx=input('请输入X的值')switch fix（x）case{1}y=a*x^2+b*x+c;case{2}y=a*sin(c*b)+x;case{3}y=log(b+c/x);end第五章 MATLAB绘图1.绘制下列曲线(2)x=-2*pi:pi/180:2*pi;y=(1/2*pi).*exp(-x.^2/2);plot(x,y)5(3)x=0:pi/180:2*pi;y=cos(x);z=sin(x);plot(y,z)2．绘制下列极坐标图(3)theta=-2*pi:0.01:2*pi;rho=5./cos(theta);polar(theta,rho ,'-')(4)theta=-2*pi:pi/180:2*pi;rho=(pi/3)*(theta.^2);polar(theta,rho,'k')7第六章 MATLAB 数值计算1.求数值积分1.50sin sin 5dx x x π⎰z=quad('(sin(x).^5).*sin(5*x)',0,pi) z =0.09822.20sin dx 1cos x x xπ+⎰ y=quad('(x.*sin(x))./(1+cos(x).^2)',0,pi) y =2.4674第七章 MATLAB 符号计算1.分解因式(1)x 9-1syms xa=x^9-1;factor(a)ans=(x-1)*(x^2+x+1)*(x^6+x^3+1)(4)x 2+y 2+z2+2(xy+yz+zx)syms x y z;d=x^2+y^2+z^2+2*(x*y+y*z+z*x);C=factor(d)C =(x+y+z)^22.求非线性方程的符号解(1)ax2+bx+c=0x=solve('a*x^2+b*x+c=0','x')x =-(b+(b^24*a*c)^(1/2))/(2*a) -(b-(b^2-4*a*c)^(1/2))/(2*a)9。
2019年Matlab选修课考试题目
2019年Matlab选修课考试题目一、填空题（每题4分，共40分）1. MATLAB的含义为_________________ 。
2. fix(-3.5)=__________,round(4.7)=________________.3. A=[4,5,6;3,1,2;7 8 9]; A(:,[1,3])=[ ]; A=__________________.4. A=rand(4,5);b=size(A);c=length(A)；则b和c的值分别为_________________.5. 为了使两个plot的图形在同一个坐标显示，可以使用__ _____命令进行图形保持；可以用___________命令为图形添加网格。
6. 已知A=[1 2 3；4 5 0；7 8 9]；B=[2 0 3；1 5 7；0 1 2]；写出下列各指令运行的结果。
(1) A==B 运行结果_______；(2) A.*B 运行结果_______.7．运行如下程序后, 命令窗口（command windows）显示的结果为_______.k=0;for n=400:600if rem(n,21)~=0 %R=rem(X,Y),求余数函数,X,Y应该为正数k=k+1;continueendbreak;endk8. 运行如下程序后, 命令窗口（command windows）显示的结果为_______.A=[13,-56,78; 25,63,-735; 56,25,563; 1,0,-1];y=max(max(A))9. 分析下面程序运行后s1,n1,s2,n2,m的值。
s1=0;s2=0;n1=0;n2=0；x=[1,-4,-8,3,10,-9,7,-3,10, 8,-5,-2,2,0 ]；m=length(x);for i=1:if x(i)<0s1=s1+x(i);n1=n1+1;elses2=s2+x(i);n2=n2+1;endends1,n1,s2,n2,m请回答s1,n1,s2,n2,m的结果: s1=_______, n1=_______, s2=_______, n2=_______,m=_______.10. 有以下程序：执行结果sum=_______.sum=0；i=1；While(i<=150)sum=sum+i；i=i+1；endsum二、编程题（共60分）1. （10分）（将程序保存为test1.m文件）编写程序：计算1~400中即能被3整除，又能被7整除的所有数之和。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
2015年春季学期《MATLAB语言及应用》课程试卷姓名：李清蔚学号：11408010304学院：机电学院专业：机械设计制造及其自动化(1)请概括Matlab中冒号(:)，逗号(,)，分号(;)的功能？(6分)：用于生成一位数组，表示一位数组的全部元素或多维数组的某一维的全部元素，用于要显示计算结果的命令之间的分隔符，用于输入变量之间的分隔符，用于数组行元素之间的分隔符。
;用于不显示计算结果命令行的结尾，用于不显示计算结果命令之间的分隔符，用于数组元素行之间的分隔符。
(2)a=[1 -6 0 0.45 0 100]，b=[0 1.5 4 0 0.5 -34]，请问在进行逻辑运算时，a、b分别相当于什么样的逻辑量？a与b的逻辑与、逻辑或都分别是什么？(4分)元素方式的逻辑变量a&b=[0 1 0 0 0 1], a|b=[1,1,1.1,1](3)有几种建立矩阵的方法？各有什么优点？(4分)直接输入矩阵，eye()生成单位矩阵，通过M文件建立矩阵，通过外部数据文件导入矩阵(4)生成一个7阶的魔方矩阵A，再采用相应的MATLAB命令，将其全部奇数行提取出来，赋给矩阵B，给出执行语句以及计算结果。
(4分)A=magic(7)table=[1 3 5 7];[n,len]=size(table);B=A(table(1),:);for i=2:lenB=[B;A(table(i),:)];EndB= 30 39 48 1 10 19 2846 6 8 17 26 35 3713 15 24 33 42 44 422 31 40 49 2 11 20(5)生成5阶魔方矩阵A，再分别计算A中元素的平方、A矩阵的平方，给出执行语句以及计算结果。
(5分)A=magic(5)B=A.^2C=A^2A =17 24 1 8 1523 5 7 14 164 6 13 20 2210 12 19 21 311 18 25 2 9B =289 576 1 64 225529 25 49 196 25616 36 169 400 484100 144 361 441 9121 324 625 4 81C =1090 900 725 690 820850 1075 815 720 765700 840 1145 840 700765 720 815 1075 850820 690 725 900 1090(6)生成矩阵84126115689a⎡⎤-⎢⎥=--⎢⎥⎢⎥-⎣⎦，计算a的行列式、逆矩阵，给出执行语句以及计算结果。
(5分)b =784c =0.0753 0.0765 0.14290.0306 0.1837 0.1429-0.0230 0.1122 0.1429(7)生成矩阵23459161514151196i iA i i ii⎡⎤+--⎢⎥=-++⎢⎥⎢⎥-⎣⎦，计算A的转置矩阵B、共轭转置矩阵C，给出执行语句以及计算结果。
(5分)A=[2+3i 4-5i -9;1-6i 1+i 5+14i;15 11 9-6i]B=A.'C=A'B=2.0000 +3.0000i 1.0000 - 6.0000i 15.0000 + 0.0000i4.0000 -5.0000i 1.0000 + 1.0000i 11.0000 + 0.0000i-9.0000 + 0.0000i 5.0000 +14.0000i 9.0000 - 6.0000iC =2.0000 -3.0000i 1.0000 + 6.0000i 15.0000 + 0.0000i4.0000 +5.0000i 1.0000 - 1.0000i 11.0000 + 0.0000i-9.0000 + 0.0000i 5.0000 -14.0000i 9.0000 + 6.0000i(8)对于BAX=，如果4927184357A⎡⎤-⎢⎥=-⎢⎥⎢⎥-⎣⎦，⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=183657B，求解X，给出执行语句以及计算结果。
(4分)A=[4 -9 2;7 18 -4;-3 5 7];B=[57;36;18];X=A\BX =10.0000-0.31517.0822(9)说明break语句、continue语句和return语句的用法及区别。
(3分)break语句用于终止for循环和while循环的执行。
当遇到break时，则退出循环体，继续执行循环体外的下一个语句。
在嵌套循环中，break往往存在于内层的循环中continue命令用于结束本次for或while循环，与break命令不同的是continue指结束本次循环而继续进行下次循环return语句用来终止被调用函数的运行，后面的程序代码将不再执行，直接返回到上一级调用函数。
(10)根据x= [23 71 12 8 66 16 55]绘制三维饼图，并将第五个切块分离出来，给出执行语句以及计算结果。
(4分)x= [23 71 12 8 66 16 55]explode=[0 0 0 0 1 0 0]pie(x,explode)(11)试求出如下极限，给出执行语句以及计算结果(6分)（1）xxxx163)(lim+∞→；（2）114lim0-+→→xyxyyx（1）syms xf=(x^3+x^6)^(1/x);y=limit(f,x,inf) y = 1 （2）syms x,y ;f=(4*x*y)/(sqrt(x*y+1)-1)z=limit(limit(f,x,0),y,0) z=8(12)假定某天的气温变化记录如下表，试作图描述这一天的气温变化规律，再用二次多项式拟合的方法得出其变化规律，给出执行语句以及二次多项式的系数。
(6分)时刻t(h) 0 1 2 3 4 5 6 7 8 9 10 11 12 温度o C(t) 15o 14o 14o 14o 14o 15o 16o 18o 20o 22o 23o 25o 28o 时刻t(h) 13 14 15 16 17 18 19 20 21 22 23 24温度o C(t) 31o 32o 31o 29o 27o 25o 24o 22o 20o 18o 17o 16ox=[0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24]y=[15 14 14 14 14 15 16 18 20 22 23 25 28 31 32 31 29 27 25 24 22 20 18 17 16] a=plot(x,y)p=polyfit(x,y,5)z=vpa(poly2sym(p),2)0.016(13)算以下两个多项式A、B的乘积A*B，以及A/B的商及余数，给出执行语句以及计算结果。
(4分)A=6x4-4x3+5x-9B=6x3+2x2-5x+7. A=[6 -4 5 -9]B=[6 2 -5 7]P1=conv(A,B)y=poly2sym(P1)[q,r]=deconv(A,B);y1=poly2sym(q)y2=poly2sym(r)y = 36*x^6 - 12*x^5 - 8*x^4 + 18*x^3 - 71*x^2 + 80*x - 63y1 =1y2 = - 6*x^2 + 10*x – 16(14)计算z=3x4-6xy+3y2x+2y3错误!未找到引用源。
+12的zy∂∂错误!未找到引用源。
和2zy x∂∂∂。
(4分)错误!未找到引用源。
syms x y;z=3*x^4-6*x*y+3*y^2*x+2*y^3+12; a=diff(z,y)b=diff(a,x)a = 6*y^2 + 6*x*y - 6*xb = 6*y – 6(15)用符号计算验证三角等式：11sin()cos()sin()sin()22=++-a b a b a b，给出执行语句以及计算结果。
(4分)a=sym('a','real')b=sym('b','real')x=sin(a)*cos(b)y=sin(a+b)/2+sin(a-b)/2simple(x-y)ans =(16)求矩阵sin()cos()cos()sin()a a A a a ⎡⎤=⎢⎥-⎣⎦的逆、行列式的值，给出执行语句以及计算结果。
(4分)A=sym('[sin(a) cos(a); -cos(a) sin(a)]')B=inv(A);c=det(A);disp(B)disp(c)[ sin(a)/(cos(a)^2 + sin(a)^2), -cos(a)/(cos(a)^2 + sin(a)^2)][ cos(a)/(cos(a)^2 + sin(a)^2), sin(a)/(cos(a)^2 + sin(a)^2)]cos(a)^2 + sin(a)^2(17)某仪器需安装一个电子元件，需要电子元件的使用寿命不低于1000小时即可。
现有甲乙两厂的电子元件可供选择，甲厂生产的电子元件的寿命服从正态分布年N(1100, 502)，乙厂生产的电子元件的寿命服从正态分布N(1150, 802)。
问应选哪个工厂的产品呢？给出执行语句以及计算结果。
(4分)r1 = 1 - normcdf(1000,1100,50);r2 = 1 - normcdf(1000,1150,80);delta = (r1 - r2) > 0;open delta;(18)用下列公式计算自然对数e 的近似值，分别用for 循环、while 循环来进行编程，直到最后一项的绝对值小于10-7，给出执行语句以及计算结果。
(8分)1111111!2!3!4!5!e ≈++++++while:n=1i=1sum=1while(n>=10^(-7))n=1/factorial(i)i=i+1sum=sum+nenddisp(sum)e=2.7183for: n=1i=1sum=1for i=(1:1:12)n=1/factorial(i)i=i+1sum=sum+nenddisp(sum)(19)房间里面有序号分别为1，2，…，100的电灯，初始时都关闭；又有100个序号分别为1，2，…，100的学生，当他们依次从这100盏灯前走过时分别将序号是本人序号的整数倍的所有灯开关切换一次（原来为开的变为关，原来为关的变为开）。
请问当所有的人都走过之后，那些灯是亮的？编程并给出亮着灯的序号。
(8分)a = -1 * ones(1,100);i = -1;n = 1;while (i >= -100)while( (n*i) >= -100)a(1,-n*i) = i * a(1,-n*i);n = n + 1;end;n = 1;i = i - 1;endb = (a > 0);find(b);(20)考虑将一枚均匀硬币掷N次，当N很大时，正面出现的机率接近0.5，采用matlab生成随机数列，再设计一个随机模拟程序验证这一现象。