2016年秋matlab期中大作业
MatLab考试题 期中3
MatLab 测试题（高难度）一．编程题（第一第二各6分，第六第七题各5分，其余每题5分，共56分）1. A 是一个維度m×n 的矩阵. 写一段程序, 算出A 中有多少个零元素A=input('请输入一个矩阵')[m,n]=size(A)sig=0for i=1:mfor j=1:nif A(i,j)==0sig=sig+1;endendend。
2. 对任意给定n 值求和：()()111,21231-,2n n n n n n -⎧+⎪⎪-+-+-=⎨⎪⎪⎩L clearn=input('输入n=')switch mod(n,2)case 0-n/2otherwise(n+1)/2end3.用MATLAB 语句表达“如果a 等于 b 且c 等于0就让d 等于3，否则如果a 大于b 且c=1让d 等于 0，其余情况让d 等于3”；if a==b&c==0d=3else if a>b&c==1d=0elsed=3endend4.产生7×9阶的正态分布随机数矩阵R1, 求其各列的平均值和均方差。
并求该矩阵全体数的平均值和均方差。
x=randn(7,9)m=mean(x)mm=mean(mean(x))偶数 奇数s=std(x)ss=std(x(:))5. f(x)=x^5-4x^4+3x^2-2x+6取x=[-3,13]之间函数的值，求它的零点。
fun1=@(x)x.^5-4*x.^4+3*x.^2-2*x+6fzero(fun1,[-3,13])fplot(fun1,[-3,13])grid on6. 设计一段程序：(1) 将0.10.10.10.11223,,sin ,cos x x x x y e y e y e x y e x ----==-==分别用实线、点线、点划线、虚线，且颜色分别用黑、红、绿、蓝，而线宽为2磅。
(2) 用黑色星号将点()()()()()()0,0,0,0.5,0,,0,0,0,1.5,0,2ππππ绘制在同图上。
matlab大作业
2014年春季学期MATLAB 课程考查题一.必答题(80分)1.如何设置当前目录和搜索路径，在当前目录上的文件和在搜索路径上的文件有何区别？(2)答：方法一：在MATLAB命令窗口中输入editpath或pathtool命令或通过【File】/【SetPath】菜单，进入“设置搜索路径”对话框，通过该对话框编辑搜索路径。
方法二：在命令窗口执行“path(path,…D:\Study ‟)”，然后通过“设置搜索路径”对话查看“D:\Study”是否在搜索路径中。
方法三：在命令窗口执行“addpathD:\Study- end”，将新的目录加到整个搜索路径的末尾。
如果将end改为begin，可以将新的目录加到整个搜索路径的开始。
区别：当前文件目录是正在运行的文件的目录，显示文件及文件夹的详细信息，且只有将文件设置为当前目录才能直接调用。
搜索路径中的文件可以来自多个不同目录，在调用时不用将其都设置为当前目录，为同时调用多个文件提供方便。
2.创建符号变量和符号表达式有哪几种方法？(4)答：定义符号变量：方法一：sym函数，可以定义单个符号变量，调用格式为：符号量名=sym（'符号字符串'）；方法二，syms函数，定义多个符号变量，调用格式：syms 符号变量名1 符号变量名2 符号变量名3 符号变量名4 。
空格隔开。
定义符号表达式：方法一：用单引号定义符号表达式；方法二：用sym函数定义符号表达式；方法三：用已经定义好的符号变量组成符号表达式；3.GUIDE提供哪些常用的控件工具，各有什么功能？（5分）答：一、控件风格和外观（1）BackgroundColor：设置控件背景颜色，使用[R G B]或颜色定义。
（2）CData：在控件上显示的真彩色图像，使用矩阵表示。
（3）ForegroundColor：文本颜色。
（4）String属性：控件上的文本，以及列表框和弹出菜单的选项。
（5）Visible：控件是否可见。
MatLab考试题期中3
MatLab考试题期中3MatLab考试题期中3考察的主要内容为MatLab软件的应用和具体操作技巧。
以下为详细解答。
1. 给定两个向量A和B，编写MatLab代码计算它们的内积。
解答：```MatLabA = [1, 2, 3]; % 输入向量AB = [4, 5, 6]; % 输入向量Bproduct = dot(A, B); % 计算内积disp(product); % 显示内积结果```2. 编写一个MatLab函数来计算n阶单位矩阵，其中n为输入参数。
解答：```MatLabfunction I = unitMatrix(n)I = eye(n); % 利用MatLab内置函数计算单位矩阵end```3. 给定一个矩阵A，编写MatLab代码来计算A的逆矩阵，如果不存在逆矩阵则显示错误提示。
解答：```MatLabA = [1, 2; 3, 4]; % 输入矩阵Aif det(A) == 0 % 判断A是否可逆，如果行列式为0则不可逆disp('Error: A is singular.'); % 显示错误提示elseinv_A = inv(A); % 计算逆矩阵disp(inv_A); % 显示逆矩阵end```4. 给定一个输入矩阵A，编写MatLab代码来计算矩阵A的特征值和特征向量。
解答：```MatLabA = [1, 2; 3, 4]; % 输入矩阵A[V, D] = eig(A); % 计算特征值和特征向量，V为特征向量矩阵，D 为特征值矩阵disp('Eigenvalues:'); % 显示特征值disp(D);disp('Eigenvectors:'); % 显示特征向量disp(V);```5. 编写MatLab代码来生成一个5行5列的随机矩阵，并计算矩阵的迹。
解答：```MatLabM = rand(5); % 生成一个5行5列的随机矩阵trace_M = trace(M); % 计算矩阵的迹disp(trace_M); % 显示迹的结果```通过以上解答，我们可以看到MatLab考试题期中3主要考察了MatLab软件的基本操作和常用函数的应用。
Matlab 期中试题
>> compose(f,g,x,y,z) ans = cos(z+2*t)^2/(1+u) 5. >> syms x y; >> z=y*x^2+3*y^2*x+2*y^3; diff(z,y) a =x^2+6*y*x+6*y^2
8．MATLAB 语言中用来设置曲面光源的图形光照控制函数是 sphere() 。
1
B ，须保持这些元素的相对位置不变。
f 2 =1
fn fn1 fn2 ,(n 2)
编写一个 matlab 程序，求 Fibonacci 数列的第 27 项。 for i=3:27; f(1)=1; f(2)=1; f(i)=f(i-1)+f(i-2); end >> f(27) ans = 196418
⒊ 计算级数 S ⒋ 设 f ( x)
1 及其前 100 项的部分和； 2 n 1 n
>> y1=sin(x1); >> y2=cos(x2); >> plot(x1,y1,'- g',x2,y2,'--r','linewidth',7); >> legend('正弦曲线','余弦曲线','Location','North') >>
x2 , g ( y ) cos( y 2t ), 求复合函数 f ( g ( z )) ； 1 u
⒌ 计算 z yx2 3 y 2 x 2 y3 的
得分 阅卷人
六、 （10 分）已知： x 5 2 1 6 3 7 ，求 （1） x 各元素的和； （2） x 各元素的积； （3）向量 x 的算术平均值和中值。
MATLAB作业6(2016)作业本
MATLAB 作业61、在图形绘制语句中，若函数值为不定式NaN ，则相应的部分不绘制出来，试利用该规律绘制sin()z xy =的表面图，并剪切下2220.5x y +≤的部分。
解：>> [x,y]=meshgrid(-1:.1:1); z=sin(x.*y);ii=find(x.^2+y.^2<=0.5^2); z(ii)=NaN; surf(x,y,z)2、试求解下面的无约束最优化问题。
222222211433222424min 100()(1)90()(1)10.1[(1)(1)]19.8(1)(1)xx x x x x x x x x x -+-+-+-+-+-+--解：>> f=inline(['100*(x(2)-x(1)^2)^2+(1-x(1))^2+',...'90*(x(4)-x(3)^2)+(1-x(3)^2)^2+',...'10.1*((x(2)-1)^2+(x(4)-1)^2)+',...'19.8*(x(2)-1)*(x(4)-1)'],'x');x=fminunc(f,ones(7,1))x =1.0e+002 *0.105464467987131.112320667672340.06782323911149-1.115047464577260.010000000000000.010000000000000.010000000000003、 试用图解法求解下面的非线性规划问题，并用数值求解算法验证结果。
321211221212min (44)20..100,0x x x x x s t x x x x +-+-+≥⎧⎪-+-≥⎨⎪≥≥⎩ 解：>> [x1,x2]=meshgrid(0:0.02:1,1:0.02:2);z=x1.^3+x2.^2+4*x1+4;ii=find(x1-x2+2<0); z(ii)=NaN;ii=find(-x1.^2+x2-1<0); z(ii)=NaN;ii=find(x1<0); z(ii)=NaN; ii=find(x2<0); z(ii)=NaN;surf(x1,x2,z)function [c,ce]=exc6f4(x)ce=[];c=[x(1)^2-x(2)+1];>> f_opt=inline('x(1)^3+x(2)^2+4*x(1)+4','x');A=[-1 1]; B=2; Aeq=[]; Beq=[]; xm=[0;0];x=fmincon(f_opt,[0;1],A,B,Aeq,Beq,xm,[],'exc6f4');4、 试求解此线性规划问题：6712341236723571,2,,7min 421..390x x x x x x x x x x x s t x x x x x ++++=⎧⎪-+--+=⎪⎨+++=⎪⎪≥⎩解：>> [x1,x2]=meshgrid(0:0.02:1,1:0.02:2);z=x1.^3+x2.^2+4*x1+4;ii=find(x1-x2+2<0); z(ii)=NaN;ii=find(-x1.^2+x2-1<0); z(ii)=NaN;ii=find(x1<0); z(ii)=NaN; ii=find(x2<0); z(ii)=NaN;surf(x1,x2,z)function [c,ce]=exc6f4(x)ce=[];c=[x(1)^2-x(2)+1];>> f_opt=inline('x(1)^3+x(2)^2+4*x(1)+4','x');A=[-1 1]; B=2; Aeq=[]; Beq=[]; xm=[0;0];x=fmincon(f_opt,[0;1],A,B,Aeq,Beq,xm,[],'exc6f4');5、 试求解下面的二次型规划问题，并用图示的形式解释结果。
Matlab期中考试
5.MATLAB在工作区如果某个较长的命令一行输入不下,可以在命令行后面用三点“...”表示续行。
6.MATLAB可以输入字母、汉字，但标点符号必须在英文状态下书写。
(2) A,B矩阵对应元素相乘，相除，左除；
（3）对比（1）和（2）不同之处；
（4）同时给A矩阵乘3次幂运算。
（5）结果中为什么会出现INF。
解(1): >> A=[1 2 3;4 5 6;7 8 9];
>> B=[1 0 1;5 1 0;7 2 100];
>> A+B
ans =
2 2 4
9 6 6
1.0000 4.0000 0.0900
>> A.\B
ans =
1.0000 0 0.3333
1.2500 0.2000 0
1.0000 0.2500 11.1111
(3)(1)是两个可乘除矩阵的乘除(2)是同阶矩阵内对应元素相乘相除并左除
(4) >> A=[1 2 3;4 5 6;7 8 9];
>> A^3
>> xor(A,B)
ans =
1 1 0 0 0 0 0 1 1 0 1 1
>> ~A
ans =
0 0 0 1 0 1 1 1 1 0 0 1
>> ~B
ans =
1 1 0 1 0 1 1 0 0 0 1 0
6、下面定义的用户变量是否正确，如果错误，请说明为什么，再修改.5分
Matlab期中复习题
一、单项选择题1、求矩阵特征值的函数是 【 A 】(A) eig (B) det (C) trace (D) eye2、MATLAB 中rand 函数的作用是 【 D 】(A) 矩阵分解 (B) 求矩阵的秩 (C) 矩阵转臵 (D) 产生随机矩阵3、下列标识符中可以用作Matlab 变量名的是 【 B 】(A) 123 (B) _abc (C) sin (D) end5、MATLAB 中获取函数句柄的符号是 【 C 】(A) & (B) ~ (C) @ (D) ^6、以下有关Matlab 说法中不正确的是 【 C 】(A) Matlab 区分字母的大小写；(B) Matlab 中数组下标是从1开始；(C) Matlab 中可以定义指针类型变量；(D) Matlab 中的函数可以返回可变数目的输出参数；7、根据一组给定的采样点x0=0:0.1:1; y0=[-.447,1.978,3.11,5.25,5.02, 4.66,4.01,4.58, 3.45,5.35,9.22];,求3阶拟合多项式的Matlab 命令是【 C 】(A) polyval(x0,y0,3) (B) interp1(x0,y0,3) (C) polyfit(x0,y0,3) (D)spline(x0,y0,3)8、以下有关Matlab 函数的说法中正确的是 【 B 】(A) 函数不在Matlab 的当前工作目录下时不可执行；(B) 用户定义函数被调用时，以所存储的m 文件的文件名作为函数名；(C) 一个函数文件中只允许定义一个函数；(D) 用户定义的函数不在Matlab 的搜索路径、当前目录中时，也可以执行。
二、填空题9、表示)cos(2x y π=的Matlab 表达式为_______sqrt(pi/2)*cos(x) ___。
10、Matlab 中与linspace(0,2*pi,101) 等价的命令是 0:pi/50:2*pi 。
matlab期中考试试题
1、在matlab环境下有6个变量a、b、c、x、y和z，现欲将这6个变量全部保存到c：盘下avar.mat文件中，应如何实现？若仅将指定的变量x、y和z保存到c:盘下变量文件mvar.mat中应如何实现？若仅将c:盘下mvar.mat文件中的变量z加入到工作空间中应如何实现？（6分）答：(1)save avar a b c x y z(2)save c:\avar x y z(3)load c:\avar z2、现有5行7列的矩阵A，如何取出矩阵A中第2、3、4行的全部元素构造矩阵B？如何将矩阵A中的第3行和第4列的元素删除？（6分）答：(1)B=A(2:4 , :);(2)A(3:4 , :)=[];3、如何将矩阵A的行数和列数分别保存到变量L和H中？如何生成向量1、3、5、……、97、99并保存到变量b中？（6分）答：(1)[L , H] =size(A);(2)b=1:2:99;4、生成全0矩阵、全1矩阵和单位矩阵的函数分别是什么？（4分）答：全0矩阵——zeros；全1矩阵——ones；单位矩阵——eye。
5、有矩阵A=[1，2，3；4，5，6]和矩阵B=[2，3，4；5，6，7]，求A.*B的结果是多少？如何将4行6列的矩阵A转换成3行8列（8分）答：(1)A=[1,2,3;4,5,6];B=[2,3,4;5,6,7]; A.*B（结果略去）(2)reshape(A , 3 , 8) ;6、如何计算矩阵A中值大于5的元素的个数？如何提取矩阵A中第2条对角线上的元素？矩阵A=[1，2，3；4，5，6]，则A的转置矩阵？（9分）答：k=find(A>5);length(k);答：triu(A , 2) 答：A’（结果略去）7、用求逆矩阵的方法求下面线性方程组：（5分）8、有两个矩阵A和B，如何实现并判断两个矩阵的性能好坏？（5分）x+2y+3z=5 答：可通过使用矩阵的条件数来判断矩阵的性能好坏。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
2016年秋matlab 期中大作业（40分，共5道题，每题8分）姓名：党朴成 学号： 1151200218 学院：理学院专业：数学系1. 解线性方程 （第二章线性方程）⎪⎩⎪⎨⎧=++=++=++105481272x 6963x 恰定方程z y x z y z y 和 ⎪⎪⎩⎪⎪⎨⎧=-+=++=++=++6321054812726963超定方程z y xz y xz y xz y x （1）分别求上两个系数矩阵的行列式det 、逆inv 、伪逆pinv （2）分别采用左除法、逆乘法和伪逆乘法求解；（3） 对比两个方程，解释逆乘法和伪逆乘法区别。
（8分）>> run('C:\Users\Administrator\Desktop\question_1.m') A1_det = 108 A1_inv =-0.1667 0.0833 0.0833 0.2407 -0.3981 0.2685 0.0741 0.1852 -0.1481A1_pinv =-0.1667 0.0833 0.0833 0.2407 -0.3981 0.26850.0741 0.1852 -0.1481A2_pinv =-0.1659 0.0724 0.0925 -0.00860.2207 -0.1232 0.0377 0.21680.0836 0.0547 -0.0386 -0.1028 x1_1 =0.33330.07411.4074x1_2 =0.33330.07411.4074x1_3 =0.33330.07411.4074x2_1 =0.24922.1857 0.4057 x2_3 =0.2492 2.1857 0.4057逆乘法和伪逆乘法的区别：如果矩阵A 不是一个方阵，或者A 是一个非满秩的方阵时，矩阵A 没有逆矩阵，伪逆函数在求解这种系数矩阵为严重病态问题时可避免“伪解”的产生，即找到一个与A 的转置矩阵A ‘同型的矩阵B ，使得：A ·B ·A=A ；B ·A ·B=B ；此时称矩阵B 为矩阵A 的伪逆，也称为广义逆矩阵。
从结果中的可看出逆矩阵和伪逆矩阵区别在于如果方阵行列式为零，得到的结果不同。
方阵行列式不为零，逆矩阵和伪逆矩阵的计算值相同，反之则不同2. 四个矩阵a 、b 、c 和d ，求matlab 运算 （第二章矩阵数组运算）2936a ⎛⎫= ⎪⎝⎭，8714b ⎛⎫= ⎪⎝⎭，275294c ⎛⎫ ⎪=- ⎪ ⎪-⎝⎭，5201031d -⎛⎫= ⎪-⎝⎭ （1） 哪些矩阵可以做加法运算（至少3个）Ab,c’d,cd’（2） 哪些可以点乘数组运算（至少3个）Ab,c’d,cd（3） 哪些可以矩阵乘法运算（至少3个）Ab,ad,bd（4） 说明数组运算和矩阵运算的特点。
（8分）>> a = [2,9;3,6];b = [8,7;1,4];c = [2,7;5,-2;-9,4];d = [5,-2,0;10,3,-1]d =5 -2 010 3 -1>> a_b = a+ba_b =10 164 10>> aabb = a.*baabb =16 633 24>> ab = a*bab =25 5030 45>> ad = a*dad =100 23 -975 12 -6>> bd = b*dbd =110 5 -745 10 -4只有维数相同的矩阵才能进行加、减、点乘、点除运算。
注意只有当两个矩阵中前一个矩阵的列数和后一个矩阵的行数相同时，才可以进行乘法运算。
a\b运算等效于求a*x=b的解；而a/b等效于求x*b=a的解。
只有方阵才可以求幂3. 已知计算分段函数220()x/sin(x)0111111x x xf x xx⎧-+≤⎪=<≤⎨⎪-<⎩（第三章编程，条件语句）(1)要求采用条件语句if编写脚本文件和函数文件计算分段函数文件；(2)x由input指令从键盘给出,数据的输出f(x)值采用fprintf('f=%e',f)(3)运行脚本文件计算x=-1的f值和运行函数文件计算x=5的f值(4)说明脚本文件和函数文件区别脚本文件x = input('x =');if x<=0f = x^2-x+2;elseif x>11f = -11;elsef = x/sin(x);endfprintf('f=%e',f);>> fxx =-1f=4.000000e+00>>>> fxx =5f=-5.214176e+00>>函数文件function f= fenduan(x)if x<=0f = x^2-x+2;elseif x>11f = -11;elsef = x/sin(x);endfprintf('f=%e',f);>> fenduan(-1); f=4.000000e+00>> >> fenduan(5);f=-5.214176e+00>>脚本式M 文件：无函数定义行；无输入和输出量，也不一定要返回结果。
在base workspace 中数据操作，运行后变量驻留其中；全局变量函数式M 文件：有函数定义行；可有输入和输出变量，并有返回结果；中间变量存在临时工作空间，它随函数结束而删除；局部变量，除特别声明4. 多项式运算 （第四章多项式） 已知：第一个多项式是 （1） 表示多项式p1=[] （2） 采用roots 求根（3） 采用polyval 求x=6时y 值 （4） 采用polyder 对多项式一阶求导 （5） 采用polyint 对多项式不定积分 已知第二个多项式为 （6） 求两个多项式的乘积conv （p1，p2） （7） 求两个多项式的除法deconv （p1，p2）>> p1 = [1,0,0,4,5]; >> r = roots(p1) r =1.0987 + 1.4551i 1.0987 - 1.4551i -1.0987 + 0.5449i -1.0987 - 0.5449i>> y = polyval(p1,6) y =y =x 4+ 4 x + 5 y=x 3+ 3 x 2 + 51325>> pd = polyder(p1)pd =4 0 0 4>> pn = polyint(p1)pn =0.2000 0 0 2.0000 5.0000 0>> p2 = [1,3,0,5];>> p1p2 = conv(p1,p2)p1p2 =1 3 0 9 17 15 20 25>> dp1p2 = deconv(p1,p2)dp1p2 =1 -35.第四章差值拟合某年美国旧车价格的调查资料如下表，其中xi表示轿车的使用年数，yi表示相应的平均价格。
(1)2次多项式拟合（polyfit）(2)采用2次拟合多项式预测每隔0.5年的价格(polyval)，再用插值法（interpl）预测一下；(3)绘制年数-价格曲线，要求用plot不同线型（LineSpec）绘制原始数据点、插值数据和2次多项式曲线，标注坐标轴为‘年份’和‘价格’（xlabel、ylabel）、图形标题（title）为‘插值’，以及标注图例legend，说明插值和拟合的特点x i 1 2 3 4 5 6 7 8 9 10 y i2615 1943 1494 1087 765 538 484 290 226 204>> x = 1:10;>> y = [2615,194,149,1087,76,53,48,29,22,20];>> p2 = polyfit(x,y,2)p2 =1.0e+03 *0.0460 -0.6775 2.3845>> xi = [1:0.5:10];>> yi1 = polyval(p2,xi)yi1 =1.0e+03 *1 至13 列1.7530 1.4718 1.2136 0.9783 0.7661 0.5768 0.4106 0.2673 0.1471 0.0498 -0.0244 -0.0757 -0.103914 至19 列-0.1092 -0.0915 -0.0507 0.0130 0.0997 0.2095>> yi2 = interp1(x,y,xi)yi2 =1.0e+03 *1 至13 列2.6150 1.4045 0.1940 0.1715 0.1490 0.6180 1.0870 0.5815 0.0760 0.0645 0.0530 0.0505 0.048014 至19 列0.0385 0.0290 0.0255 0.0220 0.0210 0.0200>> plot(x,y,'or',xi,yi1,'-*k', xi,yi2,'-dg');legend('原始数据','2阶多项式拟合曲线','插值曲线') ;xlabel('年份');ylabel('价格');title('插值')相同点：都需根据已知数据构造函数；可用得到函数计算未知点值。
不同点：插值构造函数必须通过已知点；拟合则不要求,只要均方差最小。