matlab习题二答案
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习题,并提供相应的答案。
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运算基础%第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基础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课后习题解答第二章](https://img.taocdn.com/s3/m/e02c0759f6ec4afe04a1b0717fd5360cba1a8daf.png)
第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第二章习题答案

第一大题:(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参考答案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。
第一章:基本操作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提供了丰富的图形绘制函数,可以绘制线图、散点图、柱状图等。
matlab第二版课后习题答案

matlab第二版课后习题答案
《MATLAB第二版课后习题答案》
MATLAB是一种强大的数学软件,被广泛应用于工程、科学和金融等领域。
《MATLAB第二版》是一本经典的教材,为了帮助学生更好地掌握MATLAB的使用,书中提供了大量的课后习题。
下面我们将为大家总结一些MATLAB第二版课后习题的答案,希望能对大家的学习有所帮助。
1. 第一章课后习题答案
第一章主要介绍了MATLAB的基本操作,包括变量的定义、矩阵的运算、函数的使用等。
在课后习题中,有一道题目是要求计算一个矩阵的逆矩阵。
答案是使用MATLAB中的inv函数,将原矩阵作为参数传入即可得到逆矩阵。
2. 第二章课后习题答案
第二章介绍了MATLAB中的绘图功能,包括二维和三维图形的绘制。
有一道课后习题是要求绘制一个正弦曲线和余弦曲线,并在同一张图上显示。
答案是使用MATLAB中的plot函数,分别绘制正弦曲线和余弦曲线,并使用legend函数添加图例。
3. 第三章课后习题答案
第三章介绍了MATLAB中的控制流程,包括if语句、for循环和while循环等。
有一道课后习题是要求编写一个程序,计算1到100之间所有偶数的和。
答案是使用for循环遍历1到100之间的所有数,判断是否为偶数并累加。
通过以上几个例子,我们可以看到MATLAB第二版课后习题的答案涵盖了各种基本和高级的操作,对于学习MATLAB是非常有帮助的。
希望大家在学习MATLAB的过程中能够多加练习,掌握更多的技巧和方法。
matlab课后习题答案

习题二之阳早格格创做1.怎么样明白“矩阵是MATLAB最基础的数据对于象”?问:果为背量不妨瞅成是仅有一止或者一列的矩阵,单个数据(标量)不妨瞅成是仅含一个元素的矩阵,故背量战单个数据皆不妨动做矩阵的惯例去处理.果此,矩阵是MATLAB最基础、最要害的数据对于象.2.设A战B是二个共维共大小的矩阵,问:(1)A*B战A.*B的值是可相等?问:没有相等.(2)A./B战B.\A的值是可相等?问:相等.(3)A/B战B\A的值是可相等?问:没有相等.(4)A/B战B\A所代表的数教含意是什么?问:A/B等效于B的顺左乘A矩阵,即A*inv(B),而B\A等效于B矩阵的顺左乘A矩阵,即inv(B)*A.3.写出完毕下列支配的下令.(1)将矩阵A第2~5止中第1, 3, 5列元素赋给矩阵B.问:B=A(2:5,1:2:5); 或者B=A(2:5,[1 3 5])(2)简略矩阵A的第7号元素.问:A(7)=[](3)将矩阵A的每个元素值加30.问:A=A+30;(4)供矩阵A的大小战维数.问:size(A);ndims(A);(5)将背量t的0元素用呆板整去代替.问:t(find(t==0))=eps;(6)将含有12个元素的背量x变换成34 矩阵.问:reshape(x,3,4);(7)供一个字符串的ASCII码.问:abs(‘123’); 或者double(‘123’);(8)供一个ASCII码所对于应的字符.问:char(49);4.下列下令真止后,L1、L2、L3、L4的值分别是几?A=1:9;B=10-A;...L1=A==B;L2=A<=5;L3=A>3&A<7;L4=find(A>3&A<7);问:L1的值为[0, 0, 0, 0, 1, 0, 0, 0, 0]L2的值为[1, 1, 1, 1, 1, 0, 0, 0, 0]L3的值为[0, 0, 0, 1, 1, 1, 0, 0, 0]L4的值为[4, 5, 6]5.已知完毕下列支配:(1)与出A的前3止形成矩阵B,前二列形成矩阵C,左下角32⨯子矩阵形成矩阵D,B与C的乘积形成矩阵E.问:B=A(1:3,:);C=A(:,1:2);D=A(2:4,3:4);E=B*C;(2)分别供E<D、E&D、E|D、~E|~D战find(A>=10&A<25).问:E<D=010001⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦,E&D=110111⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦,E|D=111111⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦,~E|~D=001000⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦find(A>=10&A<25)=[1; 5].6.当A=[34, NaN, Inf, -Inf, -pi, eps, 0]时,分解下列函数的真止截止:all(A),any(A),isnan(A),isinf(A),isfinite(A).问:all(A)的值为0any(A)的值为1isnan(A)的值为[0, 1, 0, 0, 0, 0, 0]isinf(A)的值为[0, 0, 1, 1, 0, 0, 0]isfinite(A)的值为[1, 0, 0,0, 1, 1, 1]7.用结构体矩阵去保存5名教死的基础情况数据,每名教死的数据包罗教号、姓名、博业战6门课程的结果.问:student(1).id='0001';student(1).name='Tom';student(1).major='computer';student(1).grade=[89,78,67,90,86,85];8.修坐单元矩阵B并回问有闭问题.B{1,1}=1;B{1,2}='Brenden';B{2,1}=reshape(1:9,3,3);B{2,2}={12,34,2;54,21,3;4,23,67};(1)size(B)战ndims(B)的值分别是几?问:size(B)的值为2,2.ndims(B)的值为2.(2)B(2)战B(4)的值分别是几?问:B(2)=147258369⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦,B(4)=[12][34][2][54][21][3][4][23][67]⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦(3)B(3)=[]战B{3}=[]真止后,B的值分别是几?问:当真止B(3)=[]后,B={1, [1, 4, 7; 2, 5, 8; 3, 6, 9], {12, 34, 2; 54, 21, 3; 4, 23, 67}}当真止B{3}=[]后,B={1,[]; [1, 4, 7; 2, 5, 8; 3, 6, 9],{12, 34, 2; 54, 21, 3; 4, 23, 67}}习题三1.写出完毕下列支配的下令.(1)修坐3阶单位矩阵A.问:A=eye(3);(2)修坐5×6随机矩阵A,其元素为[100,200]范畴内的随机整数.问:round(100+(200-100)*rand(5,6));(3)爆收均值为1,圆好为0.2的500个正态分散的随机数.问:1+sqrt(0.2)*randn(5,100);(4)爆收战A共样大小的幺矩阵.问:ones(size(A));(5)将矩阵A对于角线的元素加30.问:A+eye(size(A))*30;(6)从矩阵A提与主对于角线元素,并以那些元素形成对于角阵B.问:B=diag(diag(A));2.使用函数,真止圆阵左旋90o或者左旋90o的功能.比圆,本矩阵为A,A左旋后得到B,左旋后得到C.问:B=rot90(A);C=rot90(A,-1);3.修坐一个圆阵A,供A的顺矩阵战A的止列式的值,并考证A与A-1是互顺的.问:A=rand(3)*10;B=inv(A);C=det(A);先估计B*A,再估计A*B,由估计可知B*A=A*B,即A·A-1=A-1·A是互顺.4.供底下线性圆程组的解.问:A=[4,2,-1;3,-1,2;12,3,0];b=[2;10;8];x=inv(A)*b圆程组的解为x=6.000026.666727.3333 -⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦5.供下列矩阵的主对于角线元素、上三角阵、下三角阵、秩、范数、条件数战迹.(1)112351423052111509A-⎡⎤⎢⎥-⎢⎥=⎢⎥⎢⎥⎣⎦(2)0.434328.9421B⎡⎤=⎢⎥-⎣⎦问:(1)与主对于角线元素:diag(A);上三角阵:triu(A);下三角阵:tril(A);秩:rank(A);范数:norm(A,1); 或者 norm(A);或者 norm(A,inf);条件数:cond(A,1); 或者 cond(A,2); 或者 cond(A,inf)迹:trace(A);(2)【请参照(1)】.6.供矩阵A的特性值战相映的特性背量.问:[V,D]=eig(A);习题四1.从键盘输进一个4位整数,按如下准则加稀后输出.加稀准则:每位数字皆加上7,而后用战除以10的余数与代该数字;再把第一位与第三位接换,第二位与第四位接换.问:a=input('请输进4位整数:');A=[a/1000,a/100,a/10,a];A=fix(rem(A,10));A=rem(A+7,10);b=A(3)*1000+A(4)*100+A(1)*10+A(2);disp(['加稀后的值为:',num2str(b)]);2.分别用if语句战switch语句真止以下估计,其中a、b、c 的值从键盘输进.问:(1) 用if语句真止估计:a=input('请输进a的值:');b=input('请输进b的值:');c=input('请输进c的值:');x=input('请输进x的值:');y=a*x^2+b*x+c;endy=a*((sin(b))^c)+x;endy=log(abs(b+c/x));enddisp(['y=',num2str(y)]);(2) 用switch语句真止估计:a=input('请输进a的值:');b=input('请输进b的值:');c=input('请输进c的值:');x=input('请输进x的值:');switch fix(x/0.5)case {1,2}y=a*x^2+b*x+c;case num2cell(3:6)y=a*((sin(b))^c)+x;case num2cell(7:10)y=log(abs(b+c/x));enddisp(['y=',num2str(y)]);3.爆收20个二位随机整数,输出其中小于仄衡值的奇数.问:A=fix(10+89*rand(1,20));sum=0;for i=1:20sum=sum+A(i);endB=A(find(A<(sum/20)));C=B(find(rem(B,2)==0));disp(C);4. 输进20个数,供其中最大数战最小数.央供分别用循环结媾战调用MATLAB 的max 函数、min 函数去真止. 问:(1) 用循环结构真止:v_max=0;v_min=0;for i=1:20x=input(['请输进第', num2str(i), '数:']);if x>v_maxv_max=x;end;if x<v_minv_min=x;end;enddisp(['最大数为:', num2str(v_max)]);disp(['最小数为:', num2str(v_min)]);(2) 用max 函数、min 函数真止:for i=1:5A(i)=input(['请输进第', num2str(i), '数:']);enddisp(['最大数为:', num2str(max(A))]);disp(['最小数为:', num2str(min(A))]);5. 已知:236312222s =+++++,分别用循环结媾战调用MATLAB 的sum 函数供s 的值.问:(1) 用循环结构真止:s=0;for i=0:63s=s+2^i;ends(2) 调用sum函数真止:s=0:63;s=2.^s;sum(s)6.当n分别与100、1000、10000时,供下列各式的值.(1)111111(1)(ln2) 234nn+-+-++-+=(2)1111() 3574π-+-+=(3) 11111() 4166443n+++++=(4)224466(2)(2)133557(21)(21)2n nn nπ⎛⎫⨯⨯⨯⨯⎛⎫⎛⎫⎛⎫⎛⎫=⎪⎪⎪⎪ ⎪⨯⨯⨯-⨯+⎝⎭⎝⎭⎝⎭⎝⎭⎝⎭央供分别用循环结媾战背量运算(使用sum或者prod函数)去真止.问:(1) 用循环结构真止:sum=0;for k=1:100sum=sum+(-1)^(k+1)/k;endsum使用sum函数:x=[];for k=1:10000x=[x, (-1)^(k+1)/k];endsum(x)(2)用循环结构真止:sum=0;for k=1:100sum=sum+(-1)^(k+1)/(2*k-1);endsum使用sum函数:x=[];for k=1:100x=[x, (-1)^(k+1)/(2*k-1)];endsum(x)(3)用循环结构真止:sum=0;for k=1:100sum=sum+1/(4^k);endsum使用sum函数真止:x=[];for k=1:100x=[x, 1/(4^k)];endsum(x)(4)用循环结构真止:t=1;for k=1:100t=t*(((2*k)*(2*k))/((2*k-1)*(2*k+1)));endt使用prod函数真止:x=[];for k=1:100x=[x, ((2*k)*(2*k))/((2*k-1)*(2*k+1))];endprod(x)7. 编写一个函数文献,供小于任性自然数n的斐波那契(Fibnacci)数列各项.斐波那契数列定义如下:问:function x=fibnacci(n)for i=1:nif i<=2x(i)=1;elsex(i)=x(i-1)+x(i-2);endend8.编写一个函数文献,用于供二个矩阵的乘积战面乘,而后正在下令文献中调用该函数.问:函数文献myfnc.m:function [x, y]=myfnc(A, B)tryx=A*B;catchx=[];endy=A.*B;下令文献myexe.m:A=input('请输进矩阵A:');B=input('请输进矩阵B:');[x, y]=myfnc(A, B);if length(x)==0display('二矩阵的维数没有匹配,无法举止乘积运算!');disp('矩阵A 战矩阵B 的乘积为:'); xenddisp('矩阵A 战矩阵B 的面乘为:'); y9. 先用函数的递归调用定义一个函数文献供1n m i i =∑,而后调用该函数文献供100501021111k k k k k k ===++∑∑∑. 问:函数文献myfnc.m :function sum=myfnc(n, m)if n<=1sum=1;elsesum= myfnc (n-1, m)+n^m;end正在下令窗心中调用myfnc.m 文献,估计100501021111k k k k k k ===++∑∑∑: sum=myfnc(100, 1)+myfnc(50, 2)+myfnc(10,-1) 10.写出下列步调的输出截止.① s=0;a=[12,13,14;15,16,17;18,19,20;21,22,23]; for k=afor j=1:4if rem(k(j),2)~=0s=s+k(j);endendends问:真止截止为②下令文献exe.m真止后的截止为:x =41220y=2 4 6。
matlab习题2答案

matlab习题2答案MATLAB习题2答案在MATLAB习题2中,我们学习了如何使用MATLAB编程语言解决各种数学问题和数据分析任务。
在本篇文章中,我们将探讨一些常见的习题,并给出相应的答案和解释。
第一道习题是关于矩阵运算的。
我们需要编写一个MATLAB程序,计算两个矩阵的乘积。
在MATLAB中,可以使用矩阵乘法运算符“*”来实现这个功能。
例如,如果我们有两个矩阵A和B,我们可以使用以下代码来计算它们的乘积:```matlabC = A * B;```接下来,我们需要解决一个关于插值的问题。
我们需要编写一个程序,使用插值方法来估算给定点的函数值。
在MATLAB中,可以使用interp1函数来实现这个功能。
例如,如果我们有一个已知的函数y和一组离散的点x,我们可以使用以下代码来进行插值计算:```matlabx_interp = 0:0.1:10;y_interp = interp1(x, y, x_interp, 'spline');```最后,我们需要解决一个关于数据可视化的问题。
我们需要编写一个程序,绘制给定数据的散点图和拟合曲线。
在MATLAB中,可以使用scatter和plot函数来实现这个功能。
例如,如果我们有一组数据点x和y,我们可以使用以下代码来进行数据可视化:```matlabscatter(x, y);hold on;p = polyfit(x, y, 1);y_fit = polyval(p, x);plot(x, y_fit, 'r');```通过以上习题的练习,我们可以更好地掌握MATLAB编程语言的基本语法和常用函数。
希望本篇文章对大家有所帮助,也希望大家能够在日常工作和学习中充分利用MATLAB的强大功能,提高工作效率和解决问题的能力。
(完整版)MATLAB作业2参考答案

(完整版)MATLAB作业2参考答案MATLAB 作业⼆参考答案1、试求出如下极限。
【求解】极限问题可以由下⾯语句直接求解。
>> syms x; f=(x+2)A(x+2)*(x+3)A(x+3)/(x+5F(2*x+5); limit(f,x,i nf) ans = exp(-5) >> syms x yfa=(x A 2*y+x*y A 3)/(x+y)A 3; limit(limit(fa,x,-1),y,2) ans = -6>> fc=(1-cos(xA2+yA2))*exp(xA2+yA2)/(xA2+yA2); limit(limit(fc,x,O),y,O) ans = 02、试求出下⾯函数的导数。
(1) y(x) . xsinx 、1e x , (2) atan yIn(x 2 y 2)x【求解】由求导函数diff() 可以直接得出如下结果,其中 (2)为隐函数,故需要⽤隐函数求导公式得出导数。
>> syms x;f=sqrt(x*si n( x)*sqrt(1-exp(x))); simple(diff(f)) ans =1/2/(x*s in (x)*(1-exp(x))A(1/2))A(1/2)*(si n(x)*(1-exp(x))A(1/2)+ x*cos(x)*(1-exp(x))A(1/2)-1/2*x*sin (x)/(1-exp(x))A(1/2)*exp(x)) >> syms x,y; f=ata n(y /x)_log(xA2+yA2); f1=simple(-diff(f,x)/diff(f,y)) f1 =(y+2*x)/(x-2*y)2 23、假设U cos 1 J *,试验证⼀———。
\y x y y x【求解】证明⼆者相等亦可以由⼆者之差为零来证明,故由下⾯的语句直接证明。
>> syms x y; u=acos(x/y); diff(diff(u,x),y)-diff(diff(u,y),x) ans =(1) limx (x 2)x 2(x 3)x3(x 5)2x 52..x y lim y 21(x xy 3y)3H xy2Xos2Xe2o o2xyt 2x2f2f2f4、假设 f (x, y) e dt ,试求 22 2y x x yy【求解】由下⾯的命令可以得出所需结果。
matlab习题二答案

matlab习题二答案Matlab习题二答案Matlab是一种强大的数值计算和科学计算软件,广泛应用于工程、科学、金融等领域。
在学习和使用Matlab的过程中,习题是一种很好的练习和巩固知识的方式。
本文将为大家提供一些Matlab习题二的答案,希望能对大家的学习有所帮助。
1. 编写一个函数,输入一个矩阵A,输出矩阵A的转置矩阵。
```matlabfunction B = transposeMatrix(A)B = A';end```2. 编写一个函数,输入一个矩阵A,输出矩阵A的每一行的平均值。
```matlabfunction avg = rowAverage(A)avg = mean(A, 2);end```3. 编写一个函数,输入一个矩阵A和一个标量k,输出矩阵A中大于k的元素个数。
```matlabfunction count = countGreaterThanK(A, k)count = sum(A(:) > k);end```4. 编写一个函数,输入一个矩阵A和一个标量k,输出矩阵A中大于k的元素的索引。
```matlabfunction indices = findGreaterThanK(A, k)indices = find(A > k);end```5. 编写一个函数,输入一个矩阵A,输出矩阵A的每一列的方差。
```matlabfunction variances = columnVariance(A)variances = var(A);end```6. 编写一个函数,输入一个矩阵A和一个标量k,将矩阵A中小于k的元素替换为0。
```matlabfunction B = replaceLessThanK(A, k)B = A;B(A < k) = 0;end```7. 编写一个函数,输入一个矩阵A和一个标量k,将矩阵A中小于k的元素替换为k。
```matlabfunction B = replaceLessThanK(A, k)B = A;B(A < k) = k;end```8. 编写一个函数,输入一个矩阵A和一个标量k,将矩阵A中大于k的元素替换为k。
matlab练习一、二答案

练习一答案1、一个三位整数各位数字的立方和等于该数本身,则称该数为水仙花数。
试编写程序找到,并输出全部水仙花数。
for n=100:999a=floor(n/100);b=floor(rem(n,100)/10);c=rem(n,10);if n==a^3+b^3+c^3s=nendend2、若一个数等于它的各个真因子之和,则称该数为完数,如6=1+2+3,所以6 是完数。
求[1 500]之间的全部完数。
for i=2:500s=1;for j=2:(i-1)if (rem(i,j)==0)s=s+j;endendif s==isendend3、写一个函数pifun 计算下列级数:F(n)=4*(1-1/3+1/5-1/7+...)并使用tic 和toc 指令来测量pifun(1000000)的计算时间。
n=input('请输入n值:');tic;f=0;for i=1:nf=f+4*((-1)^(i+1))/(2*i-1) ;endftoc4、输入一个百分制成绩,要求输出成绩等级为A~E,即90~100 为A,80~89 为B,70~79为C,60~69 为D,60 以下为E。
m=input('请输入一个百分制成绩:');m=floor(m/10);switch mcase {10,9}disp('A');case 8disp('B');case 7disp('C');case 6disp('D');disp('E');otherwisedisp('E');end5、输入x,y 的值,并将它们的值互换后输出。
x = input('x=');y = input('y=');disp ([x,y]);t=x;x=y;y=t;disp ([x,y]);6、某商场对顾客所购买的商品实行打折销售,标准如下(商品价格用price 来表示):price<200 没有折扣;200≤price<500 3%折扣;500≤price<1000 5%折扣;1000≤price<2500 8%折扣;2500≤price<5000 10%折扣;5000≤price 14%折扣输入所售商品的价格,求其实际销售价格。
Matlab程序设计教程(第二版)课后参考答案

Matlab程序设计教程(第二版)课后参考答案第一章实验1.实验一第1题2.自己验证总结Matlab的主要优点3.实验一第2题4.实验一第3题5.网站思考练习1.启动见书P5 退出见书P62.Matlab主要功能见书P2 4种功能3.分行输入行末尾加续行符,即三个点“…”4.见书P115.直接在命令窗口输入fac第二章实验1.(1) w=sqrt(2)*(1+0.34245*10^(-6))w = 1.4142(2) 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 =0.9829(3) a=3.32;b=-7.9;y=2*pi*a^2*((1-pi/4)*b-(0.8333-pi/4)*a)y = -128.4271(4) t=[2,1-3i;5,-0.65];z=0.5*exp(2*t)*log(t+sqrt(1+t.*t))z =1.0e+004 *0.0048 + 0.0002i 0.0048 - 0.0034i1.58992.0090 - 1.3580i2.实验二第1题3.实验二第2题4. H=hilb(5);P=pascal(5);Hh=det(H)Hh = 3.7493e-012Hp=det(P)Hp = 1Th=cond(H)Th = 4.7661e+005Tp=cond(P)Tp = 8.5175e+003条件数越趋近于1,矩阵的性能越好,所以帕斯卡矩阵性能更好。
5. A=[-29,6,18;20,5,12;-8,8,5]A =-29 6 1820 5 12-8 8 5[V,D]=eig(A)V =0.7130 0.2803 0.2733-0.6084 -0.7867 0.87250.3487 0.5501 0.4050D =-25.3169 0 00 -10.5182 00 0 16.8351V为A的特征向量,D为A的特征值。
Matlab上机练习二答案

Matlab上机练习⼆答案Matlab 上机练习⼆班级学号姓名按要求完成题⽬,并写下指令和运⾏结果。
(不需要画图)1、求??+-+-+-+-++=i 44i 93i 49i 67i 23i 57i 41i 72i 53i 84x 的共轭转置。
>> x=[4+8i 3+5i 2-7i 1+4i 7-5i;3+2i 7-6i 9+4i 3-9i 4+4i]; >> x’ans =4.0000 - 8.0000i 3.0000 - 2.0000i3.0000 - 5.0000i 7.0000 + 6.0000i2.0000 + 7.0000i 9.0000 - 4.0000i1.0000 - 4.0000i 3.0000 + 9.0000i7.0000 + 5.0000i 4.0000 - 4.0000i2、计算=572396a 与??=864142b 的数组乘积。
>> a=[6 9 3;2 7 5];>> b=[2 4 1;4 6 8];>> a.*bans =12 36 38 42 403、对于B AX =,如果=753467294A ,=282637B ,求解X 。
>> A=[4 9 2;7 6 4;3 5 7];>> B=[37 26 28]’;>> X=A\BX =-0.51184.04271.33184、 -=463521a ,-=263478b ,观察a 与b 之间的六种关系运算的结果。
>> a=[1 2 3;4 5 6];>> b=[8 –7 4;3 6 2];>> a>bans =0 1 01 0 1>> a>=bans =0 1 01 0 1>> aans =1 0 10 1 0>> a<=bans =1 0 10 1 0>> a==bans =0 0 00 0 0>> a~=bans =1 1 11 1 15、[]7.0802.05--=a ,在进⾏逻辑运算时,a 相当于什么样的逻辑量。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 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页。