matlab习题二答案

合集下载

MATLAB作业2参考答案(2018)

MATLAB作业2参考答案(2018)

MATLAB 作业二参考答案1、试求下面线性代数方程的解析解与数值解,并检验解的正确性。

2932114010110503848246303356684953X -----⎡⎤⎡⎤⎢⎥⎢⎥----⎢⎥⎢⎥=⎢⎥⎢⎥---⎢⎥⎢⎥------⎣⎦⎣⎦【求解】求出A , [A;B ] 两个矩阵的秩,可见二者相同,所以方程不是矛盾方程,应该有无穷多解。

>> A=[2,-9,3,-2,-1; 10,-1,10,5,0; 8,-2,-4,-6,3; -5,-6,-6,-8,-4]; B=[-1,-4,0; -3,-8,-4; 0,3,3; 9,-5,3]; [rank(A), rank([A B])] ans =4 4用下面的语句可以求出方程的解析解,并可以验证该解没有误差。

>> x0=null(sym(A));x_analytical=sym(A)\B; syms a; x=a*[x0 x0 x0]+x_analytical x =[ a+967/1535, a-943/1535, a-159/1535] [ -1535/1524*a, -1535/1524*a, -1535/1524*a][ -3659/1524*a-1807/1535,-3659/1524*a-257/1535,-3659/1524*a-141/1535] [ 1321/508*a+759/1535, 1321/508*a-56/1535, 1321/508*a-628/1535] [ -170/127*a-694/307, -170/127*a+719/307, -170/127*a+103/307] >> A*x-B ans =[ 0, 0, 0][ 0, 0, 0][ 0, 0, 0][ 0, 0, 0]用数值解方法也可以求出方程的解,但会存在误差,且与任意常数a 的值有关。

>> x0=null(A); x_numerical=A\B; syms a;x=a*[x0 x0 x0]+x_numerical; vpa(x,10)ans =[ .2474402553*a+.1396556436, .2474402553*a-.6840666849, .2474402553*a-.1418 420333][-.2492262414*a+.4938507789,-.2492262414*a+.7023776988e-1,-.2492262414*a+ .3853511888e-1][ -.5940839201*a, -.5940839201*a, -.5940839201*a][ .6434420813*a-.7805411315, .6434420813*a-.2178190763,.6434420813*a-.50860 89095][-.3312192394*a-1.604263460, -.3312192394*a+2.435364854,-.3312192394*a+.3867176824]>> A*x-Bans =[ 1/18014398509481984*a, 1/18014398509481984*a, 1/18014398509481984*a] [ -5/4503599627370496*a, -5/4503599627370496*a, -5/4503599627370496*a] [ -25/18014398509481984*a, -25/18014398509481984*a,-25/18014398509481984*a][ 13/18014398509481984*a, 13/18014398509481984*a, 13/18014398509481984*a]2、求解下面的联立方程,并检验得出的高精度数值解(准解析解)的精度。

matlab习题二带答案

matlab习题二带答案

matlab习题二带答案Matlab习题二带答案Matlab是一种强大的数值计算和科学编程软件,广泛应用于各个领域的科学研究和工程实践中。

掌握Matlab的使用对于科学家和工程师来说至关重要。

为了帮助大家更好地掌握Matlab,下面将介绍一些常见的Matlab习题,并提供相应的答案。

1. 习题一:计算平均值编写一个Matlab函数,输入一个向量,输出该向量的平均值。

答案:```matlabfunction avg = computeAverage(vector)avg = sum(vector) / length(vector);end```2. 习题二:矩阵操作编写一个Matlab函数,输入两个矩阵A和B,输出它们的乘积C。

答案:```matlabfunction C = matrixMultiplication(A, B)[m, n] = size(A);[~, p] = size(B);C = zeros(m, p);for i = 1:mfor j = 1:pfor k = 1:nC(i, j) = C(i, j) + A(i, k) * B(k, j);endendendend```3. 习题三:图像处理编写一个Matlab函数,输入一张彩色图像,输出该图像的灰度图像。

答案:```matlabfunction grayImage = convertToGray(image)grayImage = rgb2gray(image);end```4. 习题四:数据拟合给定一组数据点(x, y),编写一个Matlab函数,拟合这些数据点为一条直线,并返回拟合直线的斜率和截距。

答案:```matlabfunction [slope, intercept] = fitLine(x, y)n = length(x);sx = sum(x);sy = sum(y);sxy = sum(x .* y);sxx = sum(x .* x);slope = (n * sxy - sx * sy) / (n * sxx - sx^2);intercept = (sy - slope * sx) / n;end```5. 习题五:数值积分编写一个Matlab函数,输入一个函数f(x)和积分区间[a, b],输出该函数在该区间上的数值积分结果。

matlab程序设计与应用第二版习题答案

matlab程序设计与应用第二版习题答案

matlab程序设计与应用第二版习题答案matlab程序设计与应用第二版习题答案【篇一:matlab程序设计与应用(第二版)实验答案】%实验一 matlab运算基础%第1题%(1)z1=2*sin(85*pi/180)/(1+exp(2))%(2)x=[2,1+2i;-0.45,5];z2=0.5*log(x+sqrt(1+x.^2))%(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.^2.*(t=0t1)+(t.^2-1).*(t=1t2)+(t.^2-2*t+1).*(t=2t3)%第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(size(a))a*ba.*ba^3a.^3a/bb\a[a,b][a([1,3],:);b^2]%第3题a=[1 2 3 4 5;6 7 8 9 10;11 12 13 14 15;16 17 18 19 20;21 22 23 24 25] b=[3 0 16;17 -6 9;0 23 -4;9 7 0;4 13 11]c=a*bf=size(c)d=c(f(1)-2:f(1),f(2)-1:f(2))whos%第4题%(1):a=100:999;b=rem(a,21);c=length(find(b==0))%(2):a=lsdhksdlkklsdkl;k=find(a=aa=z);a(k)=[]%实验二 matlab矩阵分析与处理 %第1题e=eye(3);r=rand(3,2);o=zeros(2,3);s=diag([2,3]);a=[e,r;o,s];a^2b=[e,(r+r*s);o,s^2]%第2题h=hilb(5)p=pascal(5)hh=det(h)hp=det(p)th=cond(h)tp=cond(p)%第3题a=fix(10*rand(5))h=det(a)trace=trace(a)rank=rank(a)norm=norm(a)%第4题a=[-29,6,18;20,5,12;-8,8,5][v,d]=eig(a)%数学意义略%第5题方法一%(1):a=[1/2,1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/6]; b=[0.95,0.67,0.52]; x=inv(a)*b%(2):b=[0.95,0.67,0.53];x=inv(a)*b%(3):cond(a)%第5题方法二a=hilb(4)a(:,1)=[]a(4,:)=[]b=[0.95,0.67,0.52];x=inv(a)*bb1=[0.95,0.67,0.53];x1=inv(a)*b1n=cond(b)n1=cond(b1)na=cond(a) %矩阵a为病态矩阵%第6题a=[1,4,9;16,25,36;49,64,81]b=sqrtm(a)c=sqrt(a) %sqrtm函数是以矩阵为单位进行计算,sqrt函数是以矩阵中的元素进行计算%实验三选择程序结构设计%第1题程序一x=[-5.0,-3.0,1.0,2.0,2.5,3.0,5.0];y=[]; %建立存放所有y值的矩阵for x0=xif x00x0~=-3y=[y,x0*x0+x0-6];elseif x0=0x05x0~=2x0~=3y=[y,x0*x0-5*x0+6];elsey=[y,x0*x0-x0-1];endendx%输出所有xy%输出所有y%第1题程序二x=[-5,-3,1,2,2.5,3,5];y=[];for a=1:7if x(a)0x(a)~=-3y=[y,(x(a))^2+x(a)-6];elseif x(a)=0x(a)5x(a)~=2x(a)~=3y=[y,(x(a))^2-5*x(a)+6];elsey=[y,x(a)*x(a)-x(a)-1];endend%第2题程序一x=input(请输入一个百分制成绩:);if x100|x0disp(您输入的成绩不是百分制成绩,请重新输入。

matlab程序设计与应用第二版习题答案

matlab程序设计与应用第二版习题答案

matlab程序设计与应用第二版习题答案Matlab程序设计与应用第二版习题答案Matlab是一种强大的数学软件,广泛应用于科学计算、数据分析和工程设计等领域。

《Matlab程序设计与应用》是一本经典的教材,对于学习和掌握Matlab编程语言具有重要的意义。

本文将为大家提供《Matlab程序设计与应用第二版》中部分习题的答案,帮助读者更好地理解和应用Matlab。

第一章:Matlab基础1.1 基本操作1. a = 3; b = 4; c = sqrt(a^2 + b^2); disp(c);2. x = linspace(-pi, pi, 100); y = sin(x); plot(x, y);3. A = [1 2 3; 4 5 6; 7 8 9]; B = [9 8 7; 6 5 4; 3 2 1]; C = A + B; disp(C);1.2 控制结构1. for i = 1:10disp(i);end2. n = 0; sum = 0; while sum < 100n = n + 1;sum = sum + n;enddisp(n);3. x = 5; if x > 0disp('x is positive');elseif x < 0disp('x is negative');elsedisp('x is zero');end第二章:向量和矩阵运算2.1 向量运算1. A = [1 2 3]; B = [4 5 6]; C = A .* B; disp(C);2. A = [1 2 3]; B = [4 5 6]; C = A ./ B; disp(C);3. A = [1 2 3]; B = [4 5 6]; C = dot(A, B); disp(C);2.2 矩阵运算1. A = [1 2 3; 4 5 6]; B = [7 8; 9 10; 11 12]; C = A * B; disp(C);2. A = [1 2 3; 4 5 6]; B = [7 8; 9 10; 11 12]; C = B * A; disp(C);3. A = [1 2 3; 4 5 6]; B = [7 8; 9 10; 11 12]; C = A .* B; disp(C); 第三章:函数和脚本文件3.1 函数1. function y = myfunc(x)y = x^2 + 3*x + 2;end2. function [y1, y2] = myfunc(x1, x2)y1 = x1^2 + 3*x1 + 2;y2 = x2^2 + 3*x2 + 2;end3. function [y1, y2] = myfunc(x)y1 = x^2 + 3*x + 2;y2 = sin(x);end3.2 脚本文件1. x = linspace(0, 2*pi, 100); y = sin(x); plot(x, y);2. x = linspace(-10, 10, 100); y = x.^2 + 3*x + 2; plot(x, y);3. x = linspace(0, 2*pi, 100); y1 = sin(x); y2 = cos(x); plot(x, y1, x, y2);通过以上习题的答案,读者可以对Matlab程序设计的基本语法和常用函数有一个初步的了解。

[整理版]matlab课后习题解答第二章

[整理版]matlab课后习题解答第二章

第2章0000符号运算0000习题2及解答00001 说出以下四条指令产生的结果各属于哪种数据类型,是“双精度”对象,还是“符号”符号对象?00003/7+0.1; sym(3/7+0.1); sym('3/7+0.1'); vpa(sym(3/7+0.1))0000〖目的〗0000不能从显示形式判断数据类型,而必须依靠class指令。

0000〖解答〗0000c1=3/7+0.10000c2=sym(3/7+0.1)0000c3=sym('3/7+0.1')0000c4=vpa(sym(3/7+0.1))0000Cs1=class(c1)0000Cs2=class(c2)0000Cs3=class(c3)0000Cs4=class(c4) 0000c1 =00000.52860000c2 =000037/700000c3 =00000.528571428571428571428571428571430000c4 =00000.528571428571428571428571428571430000Cs1 =0000double0000Cs2 =0000sym0000Cs3 =0000sym0000Cs4 =0000sym00002 在不加专门指定的情况下,以下符号表达式中的哪一个变量被认为是自由符号变量.0000sym('sin(w*t)'),sym('a*exp(-X)'),sym('z*exp(j*th)')0000〖目的〗0000● 理解自由符号变量的确认规则。

0000〖解答〗0000symvar(sym('sin(w*t)'),1) 0 0 0 0ans = 0 0 0 0 w 0 0 0 0symvar(sym('a*exp(-X)'),1) 0 0 0 0ans = 0 0 0 0 a 0 0 0 0symvar(sym('z*exp(j*th)'),1) 0 0 0 0ans = 0 0 0 0 z 0 0 0 05求符号矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=333231232221131211a a a a a a a a a A 的行列式值和逆,所得结果应采用“子表达式置换”简洁化。

matlab第二章习题答案

matlab第二章习题答案

第一大题:(1)a = 7/3b = sym(7/3)c = sym(7/3,'d')d = sym('7/3')v1=vpa(abs(a-d))v2=vpa(abs(b-d))v3=vpa(abs(c-d))a =2.3333b =7/3c =2.3333333333333334813630699500209d =7/3v1 =0.0v2 =0.0v3 =0.00000000000000014802973661668756666666667788716(2)a = pi/3b = sym(pi/3)c = sym(pi/3,'d')d = sym('pi/3')v1=vpa(abs(a-d))v2=vpa(abs(b-d))v3=vpa(abs(c-d))a =1.0472b =pi/3c =1.047197551196597631317786181171d =pi/3v1 =0.0v2 =0.0v3 =0.00000000000000011483642827992216762806615818554(3)a = pi*3^(1/3)b = sym(pi*3^(1/3))c = sym(pi*3^(1/3),'d')d = sym('pi*3^(1/3)')v1=vpa(abs(a-d))v2=vpa(abs(b-d))v3=vpa(abs(c-d))a =4.5310b =1275352044764433/281474976710656c =4.5309606547207899041040946030989d =pi*3^(1/3)v1 =0.00000000000000026601114166290944374842393221638 v2 =0.00000000000000026601114166290944374842393221638 v3 =0.0000000000000002660111416629094726767991785515第二大题:(1)c1=3/7+0.1c1 =0.5286双精度(2)c2=sym(3/7+0.1)c2 =37/70符号(3)c3=vpa(sym(3/7+0.1))c3 =0.52857142857142857142857142857143完整显示精度第三大题:(1)findsym(sym('sin(w*t)'),1)ans =w(2)findsym(sym('a*exp(-X)' ) ,1)ans =a(3)findsym(sym('z*exp(j*theta)'),1)ans =z第四大题:A=sym('[a11 a12 a13;a21 a22 a23;a31 a32 a33]')A =[ a11, a12, a13][ a21, a22, a23][ a31, a32, a33]DA=det(A)DA =a11*a22*a33 - a11*a23*a32 - a12*a21*a33 + a12*a23*a31 + a13*a21*a32 - a13*a22*a31w=inv(A)w =[ (a22*a33 - a23*a32)/(a11*a22*a33 - a11*a23*a32 - a12*a21*a33 +a12*a23*a31 + a13*a21*a32 - a13*a22*a31), -(a12*a33 -a13*a32)/(a11*a22*a33 - a11*a23*a32 - a12*a21*a33 + a12*a23*a31 +a13*a21*a32 - a13*a22*a31), (a12*a23 - a13*a22)/(a11*a22*a33 -a11*a23*a32 - a12*a21*a33 + a12*a23*a31 + a13*a21*a32 - a13*a22*a31)] [ -(a21*a33 - a23*a31)/(a11*a22*a33 - a11*a23*a32 - a12*a21*a33 +a12*a23*a31 + a13*a21*a32 - a13*a22*a31), (a11*a33 -a13*a31)/(a11*a22*a33 - a11*a23*a32 - a12*a21*a33 + a12*a23*a31 +a13*a21*a32 - a13*a22*a31), -(a11*a23 - a13*a21)/(a11*a22*a33 -a11*a23*a32 - a12*a21*a33 + a12*a23*a31 + a13*a21*a32 - a13*a22*a31)] [ (a21*a32 - a22*a31)/(a11*a22*a33 - a11*a23*a32 - a12*a21*a33 +a12*a23*a31 + a13*a21*a32 - a13*a22*a31), -(a11*a32 -a12*a31)/(a11*a22*a33 - a11*a23*a32 - a12*a21*a33 + a12*a23*a31 +a13*a21*a32 - a13*a22*a31), (a11*a22 - a12*a21)/(a11*a22*a33 -a11*a23*a32 - a12*a21*a33 + a12*a23*a31 + a13*a21*a32 - a13*a22*a31)] IAs=subexpr(w,'d')d =1/(a11*a22*a33 - a11*a23*a32 - a12*a21*a33 + a12*a23*a31 + a13*a21*a32 - a13*a22*a31)IAs =[ d*(a22*a33 - a23*a32), -d*(a12*a33 - a13*a32), d*(a12*a23 -a13*a22)][ -d*(a21*a33 - a23*a31), d*(a11*a33 - a13*a31), -d*(a11*a23 -a13*a21)][ d*(a21*a32 - a22*a31), -d*(a11*a32 - a12*a31), d*(a11*a22 -a12*a21)]第六大题:syms ksyms x positives_s=2/(2*k+1)*((x-1)/(x+1))^(2*k+1)s_ss=simple(symsum(s_s,k,0,inf))s_s =(2*((x - 1)/(x + 1))^(2*k + 1))/(2*k + 1)警告: simple will be removed in a future release. Use simplify instead. [> In sym.simple at 41]s_ss =log(x)第八大题:syms x clearsyms xh=exp(-abs(x))*abs(sin(x))si=vpa(int(h,-5*pi,1.7*pi),64)h =abs(sin(x))*exp(-abs(x))si =1.087849417255503701102633764498941389696991336803454392428439159 第九大题:syms x y clearsyms x yr=int(int(x^2+y^2,y,1,x^2),x,1,2)r =1006/105第十大题:syms t x;f=sin(t)/t;y=int(f,t,0,x)y1=subs(y,x,sym('4.5'))ezplot(y,[0,2*pi])y =sinint(x)y1 =syms x clearsyms x ny=sin(x)^nyn=int(y,0,1/2*pi)y31=vpa(subs(yn,n,sym('1/3')))y32=vpa(subs(yn,n,1/3))y =sin(x)^nyn =piecewise([-1 < real(n), beta(1/2, n/2 + 1/2)/2], [real(n) <= -1, int(x^n/(1 - x^2)^(1/2), x, 0, 1)])y31 =1.2935547796148952674767575125656y32 =1.2935547796148952674767575125656第二十题:clearsyms y xy=dsolve('(Dy*y)/5+x/4=0','x')y =2^(1/2)*(C6 - (5*x^2)/8)^(1/2)-2^(1/2)*(C6 - (5*x^2)/8)^(1/2)y1=subs(y,'C6',1)y1 =2^(1/2)*(1 - (5*x^2)/8)^(1/2)-2^(1/2)*(1 - (5*x^2)/8)^(1/2)clfhy1=ezplot(y1(1),[-2,2,-2,2],1)set(hy1,'Color','r')grid onhold onhy2=ezplot(y1(2),[-2,2,-2,2],1)set(hy2,'Color','b')grid onxlabel('Y')ylabel('X')hold offbox onlegend('y(1)','y(2)','Location','Best')hy1 =174.0155hy2 =177.0145。

matelab作业2参考答案

matelab作业2参考答案

matelab作业2参考答案Matlab作业2参考答案Matlab作业2是一项综合性的任务,要求学生运用Matlab编程语言解决一系列数学问题。

本文将为大家提供一份参考答案,帮助学生更好地理解和完成这项作业。

首先,我们将讨论作业的第一个问题,即给定一个矩阵A,求解其特征值和特征向量。

在Matlab中,可以使用eig函数来实现这一功能。

例如,假设我们有一个3×3的矩阵A,可以按照以下方式计算其特征值和特征向量:```A = [1 2 3; 4 5 6; 7 8 9];[eigenvectors, eigenvalues] = eig(A);```在上述代码中,变量eigenvectors将存储A的特征向量,而变量eigenvalues 将存储A的特征值。

通过打印这两个变量的值,我们可以得到矩阵A的特征值和特征向量。

接下来,我们将探讨作业的第二个问题,即求解线性方程组。

假设我们有一个3×3的系数矩阵A和一个3×1的常数向量b,我们需要求解方程组Ax=b。

在Matlab中,可以使用backslash运算符来求解线性方程组。

例如,假设我们有以下方程组:```A = [1 2 3; 4 5 6; 7 8 9];b = [10; 20; 30];x = A \ b;```在上述代码中,变量x将存储方程组的解。

通过打印变量x的值,我们可以得到方程组的解。

此外,作业的第三个问题要求学生使用Matlab绘制函数图像。

在Matlab中,可以使用plot函数来实现这一功能。

例如,假设我们要绘制函数y=sin(x),其中x的取值范围为0到2π,可以按照以下方式绘制函数图像:```x = 0:0.1:2*pi;y = sin(x);plot(x, y);```在上述代码中,变量x将存储x的取值范围,变量y将存储对应的函数值。

通过调用plot函数,我们可以将函数y=sin(x)的图像绘制出来。

matlab第二版习题答案

matlab第二版习题答案

matlab第二版习题答案Matlab是一种强大的数学软件工具,被广泛应用于科学计算、数据分析和工程设计等领域。

对于学习和掌握Matlab的人来说,习题是不可或缺的一部分。

本文将为大家提供Matlab第二版习题的答案,帮助读者更好地理解和应用Matlab。

第一章:基本操作1.1 Matlab的启动和退出启动Matlab的方法有多种,可以通过桌面图标、命令行或者启动器来打开Matlab。

退出Matlab可以直接关闭窗口或者使用命令"exit"。

1.2 Matlab的基本语法Matlab的基本语法与其他编程语言相似,包括变量的定义、运算符的使用、条件语句和循环语句等。

例如,定义一个变量x并赋值为5可以使用语句"x = 5;"。

1.3 Matlab的数据类型Matlab支持多种数据类型,包括数值型、字符型和逻辑型等。

数值型可以是整数或者浮点数,字符型用单引号或双引号表示,逻辑型只有两个值true和false。

第二章:向量和矩阵操作2.1 向量的定义和运算向量是一维数组,可以通过一对方括号来定义。

Matlab提供了丰富的向量运算函数,如加法、减法、乘法和除法等。

2.2 矩阵的定义和运算矩阵是二维数组,可以通过方括号和分号来定义。

Matlab提供了矩阵的加法、减法、乘法、转置和求逆等运算。

2.3 矩阵的索引和切片可以使用索引和切片来访问矩阵中的元素。

索引从1开始,可以使用冒号表示全部元素。

切片可以用来选择矩阵的一部分。

第三章:函数和脚本文件3.1 函数的定义和调用函数是一段独立的代码块,可以接受输入参数并返回输出结果。

在Matlab中,函数的定义以关键字"function"开头,调用函数使用函数名和参数。

3.2 脚本文件的编写和运行脚本文件是一系列Matlab语句的集合,可以保存为.m文件。

通过运行脚本文件,可以一次性执行多个语句,提高效率。

第四章:图形绘制和数据可视化4.1 图形绘制函数Matlab提供了丰富的图形绘制函数,可以绘制线图、散点图、柱状图等。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

习题二1.选择题(1)下列变量名中A是合法的。

A. char_1,i,jB. x*y,a.1C. x\y,a1234D. end,lbcx说明:end是关键字,变量名中不能有运算符、标点符号,可以有下划线。

内置常量可以作为变量名,但一般不提倡这样做。

(2)下列C是合法常量。

A. 3*e10B. 1e500C. -1.85e-56D. 10-2说明:10-2是表达式,1e500超过了realmax,3*e10也是表达式。

(3)x=uint8(2.3e10),则x所占的字节是D个。

A. 1B. 2C. 4D. 8(4)已知x=0:10,则x有B个元素。

A. 10B. 11C. 9D. 12(5)产生对角线上为全1其余为0的2行3列矩阵的命令是C。

A. ones(2,3)B. ones(3,2)C. eye(2,3)D. eye(3,2)(6)已知数组123456789a⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦,则a(:,end)是指C 。

A. 所有元素B. 第一行元素C. 第三列元素D. 第三行元素(7)已知数组123456789a⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦,则运行a(:,1)=[]命令后C。

A. a变成行向量B. a数组为2行2列C. a数组为3行2列D. a数组中没有元素3(8)已知数组123456789a⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦,则运行mean(a)命令是B。

A. 计算a每行的平均值B. 计算a每列的平均值C. a数组增加一行平均值D. a数组增加一列平均值(9)已知x为一个向量,计算ln(x)的MATLA B命令是计算B。

A. ln(x)B. log(x)C. Ln(x)D. lg10(x)(10)当a=2.4,使用取整函数计算得出3,则该取整函数名为C。

A. fixB. roundC. ceilD. floor(11)已知a=0:4,b=1:5,下面的运算表达式出错的为D。

A. a+bB. a./bC. a'*bD. a*b(12)命令day(now)是指C。

A. 按照日期字符串格式提取当前时间B. 提取当前时间C. 提取当前时间的日期D. 按照日期字符串格式提取当前日期(13)以下运算符中哪个的优先级最高B。

A. *B. ^C. ~=D. |(14)运行命令bitand(20,15)的结果是C。

A. 15B. 20C. 4D. 5bitand(20,15)ans =4(15)使用检测函数isinteger(15)的结果是B。

A. 1B. 0C. falseD. trueisinteger(15)ans =(16)计算三个多项式s1、s2和s3的乘积,则算式为C。

A. conv(s1,s2,s3)B. s1*s2*s3C. conv(conv(s1,s2),s3)D. conv(s1*s2*s3)2.复数变量a=2+3i,b=3-4i,计算a+b,a-b,c=a*b,d=a/b,并计算变量c的实部、虚部、模和相角。

a=2+3i;b=3-4i;a+ba-bc=a*bd=a/breal(c)imag(c)abs(c)angle(c)ans =5.0000 - 1.0000ians =-1.0000 + 7.0000ic =18.0000 + 1.0000id =-0.2400 + 0.6800ians =18ans =1ans =18.0278ans =0.05553.用“from:step:to”方式和linspace函数分别得到从0~4π步长为0.4π的变量x1和从0~4π分成10点的变量x2。

x1=0:0.4*pi:4*pix2=linspace(0,4*pi,10)x1 =Columns 1 through 100 1.2566 2.5133 3.7699 5.0265 6.2832 7.5398 8.7965 10.0531 11.3097Column 1112.5664x2 =0 1.3963 2.7925 4.1888 5.5851 6.9813 8.3776 9.7738 11.1701 12.56644.输入矩阵123456789a⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦,使用全下标方式取出元素“3”,使用单下标方式取出元素“8”,取出后两行子矩阵块,使用逻辑矩阵方式取出1379⎡⎤⎢⎥⎣⎦。

a=[1 2 3;4 5 6;7 8 9]a(1,3)a(6)a(2:3,:)a([1,3],[1,3])a =1 2 34 5 67 8 9ans =3ans =8ans =4 5 67 8 9ans =1 37 95.输入a为3×3的魔方阵,b为3×3的单位阵,并将a、b小矩阵组成3×6的大矩阵c和6×3的大矩阵d,将d矩阵的最后一行取出构成小矩阵e。

a=magic(3)b=eye(3)c=[a,b]d=[a;b]d(end,:)a =8 1 63 5 74 9 2b =1 0 00 1 00 0 1c =8 1 6 1 0 03 5 7 0 1 04 9 2 0 0 1d =8 1 63 5 74 9 21 0 00 1 00 0 1ans =0 0 16.将矩阵123456789a⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦用flipud、fliplr、rot90、diag、triu和tril函数进行操作。

a=reshape(1:9,3,3); a=a'flipud(a)fliplr(a)rot90(a)diag(a)triu(a)tril(a)a =1 2 3 4 5 6 7 8 9 ans =7 8 9 4 5 6 1 2 3 ans =3 2 1 6 54 9 8 7 ans =3 6 9 2 5 8 14 7 ans =159ans =1 2 3 0 5 6 0 0 9 ans =1 0 0 4 5 0 7 8 97.求矩阵1358⎡⎤⎢⎥⎣⎦的转置、秩、逆矩阵、矩阵的行列式值和矩阵的三次幂。

a=[1 3;5 8];a'rank(a)inv(a)det(a)a^3ans =1 53 8ans =2ans =-1.1429 0.4286 0.7143 -0.1429 ans =-7ans =151 264440 7678. 输入a=[1.6 -2.4 5.2 -0.2],分别使用数学函数ceil、fix、floor、round查看各种取整运算的结果。

a=[1.6 -2.4 5.2 -0.2];ceil(a)fix(a)floor(a)round(a)ans =2 -2 6 0ans =1 -2 5 0ans =1 -3 5 -1ans =2 -2 5 09.求解方程组12341241234123423283687 7225x x x xx x xx x x xx x x x-++=⎧⎪++=⎪⎨-++=⎪⎪+-+=⎩。

A=[2 -3 1 2;1 3 0 1;1 -1 1 8;7 1 -2 2] b=[8;6;7;5]A =2 -3 1 21 3 0 11 -1 1 87 1 -2 2b =8675解法一:x=A\bx =2.53531.20336.8299-0.1452解法二:x=inv(A)*bx =2.53531.20336.8299-0.1452解法三:Cram法则D=det(A);D1=det([b,A(:,2:4)]);D2=det([A(:,1),b,A(:,3:4)]);D3=det([A(:,1:2),b,A(:,4)]);D4=det([A(:,1:3),b]);x1=D1/D,x2=D2/D,x3=D3/D,x4=D4/Dx1 =2.5353x2 =1.2033x3 =6.8299x4 =-0.145210.计算数组123456789A⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦,111222333B⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦的左除、右除以及点乘和点除。

A=reshape(1:9,3,3);A=A' B=repmat([1;2;3],1,3)A/BA\BA.*BA./BA =1 2 34 5 67 8 9 B =1 1 12 2 23 3 3Warning: Matrix is singular to working precision. ans =NaN NaN Inf NaN NaN Inf NaN NaN InfWarning: Matrix is close to singular or badly scaled. Results may be inaccurate. RCOND = 1.541976e-018. ans =-0.3333 -0.3333 -0.3333 0.6667 0.6667 0.6667 0 0 0 ans =1 2 3 8 10 12 21 24 27 ans =1.00002.00003.0000 2.0000 2.5000 3.0000 2.3333 2.6667 3.000011.计算函数2()10s i n (4)tf t e t =-的值,其中t 的范围从0~20步长取0.2;f1(t)为f(t)≥0的部分,计算f1(t)的值。

t=0:0.2:20;f=10*exp(2*t)-sin(4*t) f1=f(f>=0) f =1.0e+018 *Columns 1 through 100.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000Columns 11 through 200.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000Columns 21 through 300.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000Columns 31 through 400.0000 0.0000 0.0000 0.0000Columns 41 through 500.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000Columns 51 through 600.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000Columns 61 through 700.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000Columns 71 through 800.0000 0.0000 0.0000 0.0000 0.0001 0.0001 0.0002 0.0002 0.0004 0.0005Columns 81 through 900.0008 0.0012 0.0018 0.0026 0.0039 0.0058 0.0087 0.0130 0.0194 0.0289Columns 91 through 1000.0431 0.0643 0.0959 0.1431 0.2135 0.3186 0.4752 0.7090 1.0577 1.5778Column 1012.3539f1 =1.0e+018 *Columns 1 through 100.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000Columns 11 through 200.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000Columns 21 through 300.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000Columns 31 through 400.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000Columns 41 through 500.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000Columns 51 through 600.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000Columns 61 through 700.0000 0.0000 0.0000 0.0000Columns 71 through 800.0000 0.0000 0.0000 0.0000 0.0001 0.0001 0.0002 0.0002 0.0004 0.0005Columns 81 through 900.0008 0.0012 0.0018 0.0026 0.0039 0.0058 0.0087 0.0130 0.0194 0.0289Columns 91 through 1000.0431 0.0643 0.0959 0.1431 0.2135 0.3186 0.4752 0.7090 1.0577 1.5778Column 1012.353912.创建三维数组a,第一页为1234⎡⎤⎢⎥⎣⎦,第二页为1221⎡⎤⎢⎥⎣⎦,第三页为1222⎡⎤⎢⎥⎣⎦,重排生成数组b为3行、2列、2页。

相关文档
最新文档