数学实验matlab练习题
- 格式:doc
- 大小:190.50 KB
- 文档页数:12
数学实验matlab练习题模板1. 题目一:矩阵的基本运算请使用MATLAB编写一个程序,实现以下矩阵运算:- 创建一个3x3的单位矩阵A。
- 创建一个3x3的对角矩阵B,对角线上的元素依次为1, 2, 3。
- 计算矩阵A和B的和C。
- 计算矩阵A和B的乘积D。
- 计算矩阵B的转置矩阵E。
- 计算矩阵A的逆矩阵F(如果存在)。
2. 题目二:向量的线性组合给定向量a = [1, 2, 3]和b = [4, 5, 6],请使用MATLAB完成以下任务:- 计算向量a和b的点积。
- 计算向量a和b的叉积。
- 验证向量a和b是否线性相关。
- 找出一组实数x和y,使得x*a + y*b = [7, 8, 9]。
3. 题目三:多项式的运算定义两个多项式p(x) = 2x^3 - 3x^2 + 4x - 5和q(x) = x^3 + 2x^2 + 3x + 4,使用MATLAB完成以下操作:- 计算多项式p(x)和q(x)的和。
- 计算多项式p(x)和q(x)的差。
- 计算多项式p(x)和q(x)的乘积。
- 计算多项式p(x)的导数。
- 在x = 1处求多项式p(x)的值。
4. 题目四:数值积分使用MATLAB的数值积分函数计算以下定积分的近似值:- ∫(0 to 1) (x^2 + 3x - 2) dx- ∫(0 to π/2) (sin(x) + cos(x)) dx5. 题目五:微分方程的数值解给定一阶微分方程dy/dx = x^2 - y^2,初始条件y(0) = 1,使用MATLAB的数值方法求解该微分方程在区间[0, 1]上的解,并绘制解的图形。
6. 题目六:数据拟合给定一组数据点:(1, 2), (2, 3), (3, 5), (4, 7), (5, 11),请使用MATLAB完成以下任务:- 使用线性拟合方法拟合这组数据,并绘制拟合曲线。
- 使用多项式拟合方法拟合这组数据,并绘制拟合曲线。
- 比较两种拟合方法的效果,并讨论其适用性。
一、求下列方程的根1、5510x x ++=7、42254x dx x+⎰8.0sin 2x e xdx +∞-⎰二、综合题1.求由参数方程arctan x y t⎧⎪=⎨=⎪⎩dy dx 与二阶导数22d y dx 。
2.设函数y =f (x )由方程xy +e y = e 所确定,求y ′(x )。
3. 08x =展开(最高次幂为)4.求变上限函数x x ⎰对变量x 的导数。
5.求点(1,1,4)到直线L :31102x y z --==- 的距离 6、求矩阵211020413A -⎛⎫ ⎪= ⎪ ⎪-⎝⎭的逆矩阵1A - 及特征值和特征向量。
7、求二重极限1y x y →→《Matlab 数学实验》一、求下列方程的根1、5510x x ++=a=solve('x^5+5*x+1',0);a=vpa(a,6)a =1.10447+1.05983*i-1.00450+1.06095*i-.199936-1.00450-1.06095*i1.10447-1.05983*i7、42254x dx x+⎰ >> sym x;>> int(x^4/(25+x^2),x)ans =125*atan(x/5) - 25*x + x^3/38.0sin 2x e xdx +∞-⎰>> syms x;>> y=exp(-x)*sin(2*x);>> int(y,0,inf)ans =2/5二、综合题1求由参数方程arctan x y t ⎧⎪=⎨=⎪⎩dy dx 与二阶导数22d y dx 。
>> syms t>> x=log(sqrt(1+t^2));y=atan(t);>> diff(y,t)/diff(x,t)ans =1/t2.设函数y =f (x )由方程xy +e y = e 所确定,求y ′(x )。
数学实验练习题(MATLAB)注意:在下⾯的题⽬中m 为你的学号的后3位(1-9班)或4位(10班以上).第⼀次练习题1.求解下列各题: 1)30sin lim x mx mx x->- 2)(4)cos ,1000.0=x mx y e y 求 3)21/20mx e dx ?(求近似值,可以先⽤inline 定义被积函数,然后⽤quad 命令)4)4224x dx m x+? 50x =展开(最⾼次幂为8). 2.对矩阵21102041A m -?? ?= ? ?-??,分别求逆矩阵,特征值,特征向量,⾏列式,并求矩阵,P D (D 是对⾓矩阵),使得1A PDP -=。
3.已知21(),()2f x e x µσ=--分别在下列条件下画出)(x f 的图形: (1)/600m σ=,µ分别为0,1,1-(在同⼀坐标系上作图);(2)0µ=,σ分别为1,2,4,/100m (在同⼀坐标系上作图).4.画(1)sin 020cos 02100x u t t y u t u t z m ??=≤≤?=?≤≤??=?(2) sin()03,03z mxy x y =≤≤≤≤(3)sin()(/100cos )02cos()(/100cos )02sin x t m u t y t m u u z u ππ=+?≤≤?=+?≤≤?=?的图(第4题只要写出程序).5.对于⽅程50.10200m x x --=,先画出左边的函数在合适的区间上的图形,借助于软件中的⽅程求根的命令求出所有的实根,找出函数的单调区间,结合⾼等数学的知识说明函数为什么在这些区间上是单调的,以及该⽅程确实只有你求出的这些实根。
最后写出你做此题的体会.第⼆次练习题判断迭代收敛速度的程序x0=1;stopc=1;eps=10^(-8);a=1;c=1;b=2*c;d=a;k=0;f=inline('(a*x+b)/(c*x+d)');kmax=100;while stopc>eps&kk=k+1x0=f(a,b,c,d,x0)stopc=abs(x0^2-2);end1.设 ,131211p p p n nx ++++= }{n x 是否收敛?若收敛,其值为多少?精确到8位有效数字。
2015-2016数学实验练习题一、选择题1.清除Matlab工作空间(wordspace)变量的命令是(B )A. clcB. clearC. clfD.delete2. 清除当前屏幕上显示的所有内容,但不清除工作空间中的数据的命令是( A )A. clcB. clearC. clfD.delete3. 用来清除图形的命令( C )A. clcB. clearC. clfD.delete4. 在MATLAB程序中,使命令行不显示运算结果的符号是( A )A. ;B. %C. #D. &5. 在MATLAB程序中,可以将某行表示为注释行的符号是( B )A. ;B. %C. #D. &6.在循环结构中跳出循环,执行循环后面代码的命令为 ( B )A. returnB. breakC. continueD. Keyboard7.在循环结构中跳出循环,但继续下次循环的命令为( C )A. returnB. breakC. continueD. Keyboard8. MATLAB中用于声明全局变量的关键字是( C )A. infB. symsC. globalD. function9. 用户可以通过下面哪项获得指令的使用说明( A )A. helpB. loadC. demoD. lookfor10.在MATLAB命令窗口中键入命令S=zoros(3);可生成一个三行三列的零矩阵,如果省略了变量名S,MATLAB表现计算结果将用下面的哪一变量名做缺省变量名( A )A. ans;B. pi;C. NaN;D. Eps.11. 9/0的结果是( B )A. NAN;B. Inf;C. eps;D. 012.在MATLAB中程序或语句的执行结果都可以用不同格式显示,将数据结果显示为分数形式,用下面哪一条命令语句( D )A. format long;B. format long e;C. format bank;D. fromat rat13. 下列MATLAB命令中是构造1行3列的(-1,1)均匀分布随机矩阵的命令的是(D)A. randn(1,3);B. rand(1,3);C. ones(3);D. 以上都不对14. 产生四维元素都为1矩阵的语句为( A )A. ones(4)B. eye(4)C. zeros(4)D. rand(4)15. 用round 函数对数组[2.48 6.39 3.93 8.52]取整,结果为 ( C )A. [2 6 3 8]B. [2 6 4 8]C. [2 6 4 9]D. [3 7 4 9]16. y=dsolve(‘Dy=1/(1+x^2)-2*y^2’,’y(0)=0’,’x ’); ezplot(y)的功能是( A )A. 求微分方程特解并绘图;B. 解代数方程;C. 求定积分;D.求微分方程通解.17. MATLAB 命令roots([1,0,0,-1])的功能是 ( D )A. 产生向量[1,0,0,1];B. 求方程310x +=的根;C. 求多项式31x -的值;D. 求方程310x -=的根。
数学实验(概率论)题目一.用MATLAB 计算随机变量的分布1.用MA TLAB 计算二项分布在一级品率为0.2的大批产品中,随机地抽取20个产品,求其中有2个一级品的概率。
1. 用MA TLAB 计算泊松分布用MATLAB 计算:保险公司售出某种寿险保单2500份.已知此项寿险每单需交保费120元,当被保人一年内死亡时,其家属可以从保险公司获得2万元的赔偿(即保额为2万元).若此类被保人一年内死亡的概率0.002,试求:(1)保险公司的此项寿险亏损的概率;(2)保险公司从此项寿险获利不少于10万元的概率; (3)获利不少于20万元的概率. 3.用MA TLAB 计算均匀分布乘客到车站候车时间ξ()0,6U ,计算()13P ξ<≤。
4.用MA TLAB 计算指数分布用MA TLAB 计算:某元件寿命ξ服从参数为λ(λ=11000-)的指数分布.3个这样的元件使用1000小时后,都没有损坏的概率是多少? 5。
用MATLAB 计算正态分布 某厂生产一种设备,其平均寿命为10年,标准差为2年.如该设备的寿命服从正态分布,求寿命不低于9年的设备占整批设备的比例? 二.用MATLAB 计算随机变量的期望和方差 1.用MA TLAB 计算数学期望(1)用MATLAB 计算离散型随机变量的期望 1)。
一批产品中有一、二、三等品、等外品及废品5种,相应的概率分别为0.7、0.1、0.1、0.06及0.04,若其产值分别为6元、5.4元、5元、4元及0元.求产值的平均值 2)。
已知随机变量X 的分布列如下:{}kk X p 21== ,,2,1n k =计算.EX (2)用MATLAB 计算连续型随机变量的数学期望假定国际市场上对我国某种商品的年需求量是一个随机变量ξ(单位:吨),服从区间[],a b 上的均匀分布,其概率密度为: 1()0a x bx b aϕ⎧≤≤⎪=-⎨⎪⎩其它计算我国该种商品在国际市场上年销售量的期望.ξE .(3)用MATLAB 计算随机变量函数的数学期望假定国际市场每年对我国某种商品的需求量是随机变量X (单位:吨),服从[20,40]上的均匀分布,已知该商品每售出1吨,可获利3万美元,若销售不出去,则每吨要损失1万美元,如何组织货源,才可使收益最大? 2. 用MA TLAB 计算方差(1)利用MATLAB 计算:设有甲、乙两种股票,今年的价格都是10元,一年后它们的试比较购买这两种股票时的投资风险.。
MATLAB数学实验答案(全)第⼀次练习教学要求:熟练掌握Matlab 软件的基本命令和操作,会作⼆维、三维⼏何图形,能够⽤Matlab 软件解决微积分、线性代数与解析⼏何中的计算问题。
补充命令vpa(x,n) 显⽰x 的n 位有效数字,教材102页fplot(‘f(x)’,[a,b]) 函数作图命令,画出f(x)在区间[a,b]上的图形在下⾯的题⽬中m 为你的学号的后3位(1-9班)或4位(10班以上) 1.1 计算30sin limx mx mx x →-与3sin lim x mx mxx →∞-syms xlimit((902*x-sin(902*x))/x^3) ans =366935404/3limit((902*x-sin(902*x))/x^3,inf)//inf 的意思 ans = 0 1.2 cos1000xmxy e =,求''y syms xdiff(exp(x)*cos(902*x/1000),2)//diff 及其后的2的意思 ans =(46599*cos((451*x)/500)*exp(x))/250000 - (451*sin((451*x)/500)*exp(x))/250 1.3 计算221100x y edxdy +??dblquad(@(x,y) exp(x.^2+y.^2),0,1,0,1)//双重积分 ans = 2.13941.4 计算4224x dx m x +? syms xint(x^4/(902^2+4*x^2))//不定积分 ans =(91733851*atan(x/451))/4 - (203401*x)/4 + x^3/12 1.5 (10)cos ,x y e mx y =求//⾼阶导数syms xdiff(exp(x)*cos(902*x),10) ans =-356485076957717053044344387763*cos(902*x)*exp(x)-3952323024277642494822005884*sin(902*x)*exp(x)1.6 0x =的泰勒展式(最⾼次幂为4).syms xtaylor(sqrt(902/1000+x),5,x)//泰勒展式 ans =-(9765625*451^(1/2)*500^(1/2)*x^4)/82743933602 +(15625*451^(1/2)*500^(1/2)*x^3)/91733851-(125*451^(1/2)*500^(1/2)*x^2)/406802 + (451^(1/2)*500^(1/2)*x)/902 +(451^(1/2)*500^(1/2))/500 1.7 Fibonacci 数列{}n x 的定义是121,1x x ==12,(3,4,)n n n x x x n --=+=⽤循环语句编程给出该数列的前20项(要求将结果⽤向量的形式给出)。
数学软件实验测试题Matlab作业电子版姓名:**学号:**一、选择语句和循环语句编程1. 用if 语句实现以下的计算,其中,a b 的值从键盘输入。
2lg cos cos ,52sin ln ,219,14x b ax b x x a y a b b x x e x ⎧+-≤<-⎪⎪=++-≤<⎨⎪⎪+≤<⎩请计算当1,2a b 时,==5,2,3x y 时的的值=--2. 当n 取50时,求1111416644n +++++L L代码如下:1.function y=f(x,a,b)if(x>=-2&x<1)y=a*sin(b)+log(abs(b+a/x)); else if ( -5<=x&x<-2)y=a*x^2+b*cos(x); else if(1<=x&x<=4)y=9^(log10(x))+exp(cos(b)); end end end 2. function y=f(n)y=0; for i=1:n y=y+1/i.^4; end二、函数和调用函数22(,)(1,1)(2,2)(3,3)(100,100)f a b a b y f f f f =+=+++1. 已知,求的值。
L (,)f a b y 要求:定义为函数,利用循环语句求。
代码如下:1. function m=f(a,b)m1=0;m2=0; for a=1:100 for b=1:100 m2=m2+b*b endm1=m1+m2+a*a; end三、画图和拟合1. 已知1cos(4),2sin(2)1,11y x y x y x ==⋅-≤≤,完成下列操作:(1)在同一坐标系下用不同的颜色和线型绘制两条曲线。
(2)以子图形式绘制两条曲线。
(以一行两列形式绘图) 2. 已知函数f (x )在[1,101]区间上10个整数采样点的函数值如实验表1所示:实验表1 10个采样点的函数值 x1112131415161718191101f(x) 0 1.041 1.322 1.491 1.612 1.707 1.785 1.851 1.908 1.959 2.004先利用10个采样点利用plot函数绘制f(x),在同一坐标系下利用4次拟合多项式绘制出f(x)的近似曲线p(x),并利用4次拟合多项式求出平p(75)和p(111)。
MATLAB数学实验第二版课后练习题含答案课后练习题MATLAB数学实验第二版的课后练习题如下:第一章课后练习题1.编写MATLAB程序,计算并输出下列公式的结果:y = \\frac{1}{\\sqrt{2\\pi\\sigma^2}} e^{-\\frac{(x-\\mu)^2}{2\\sigma^2}}其中,x, $\\mu$, $\\sigma$ 分别由用户输入。
要求输出结果精确至小数点后两位。
答案如下:x=input('请输入 x 的值:');mu=input('请输入 mu 的值:');sigma=input('请输入 sigma 的值:');y=1/sqrt(2*pi*sigma^2) *exp(-(x-mu)^2/ (2*sigma^2));fprintf('y = %.2f\', y);2.编写MATLAB程序,求解下列方程的解:4x + y = 11\\\\x + 2y = 7答案如下:A= [4,1;1,2];B= [11;7];X=inv(A) *B;fprintf('x = %.2f, y = %.2f\', X(1), X(2));第二章课后练习题1.编写MATLAB程序,计算下列多项式的值:P(x) = x^4 - 2x^3 + 3x^2 - x + 1其中,x 由用户输入。
要求输出结果精确至小数点后两位。
答案如下:x=input('请输入 x 的值:');y=x^4-2*x^3+3*x^2-x+1;fprintf('P(%.2f) = %.2f\', x, y);2.编写MATLAB程序,绘制下列函数的图像:f(x) = \\begin{cases} x + 1, & x < 0 \\\\ x^2, & 0 \\leq x < 1 \\\\ 2x - 1, & x \\geq 1 \\end{cases}答案如下:x=-2:0.01:2;y1=x+1;y2=x.^2.* ((x>=0) & (x<1));y3=2*x-1;plot(x,y1,x,y2,x,y3);legend('y1 = x + 1','y2 = x^2','y3 = 2x - 1');总结本文提供了《MATLAB数学实验第二版》的部分课后练习题及其答案。
一元函数微分学实验1 一元函数的图形(基础实验)实验目的 通过图形加深对函数及其性质的认识与理解, 掌握运用函数的图形来观察和分析 函数的有关特性与变化趋势的方法,建立数形结合的思想; 掌握用Matlab 作平面曲线图性的方法与技巧.初等函数的图形2 作出函数x y tan =和x y cot =的图形观察其周期性和变化趋势. 解:程序代码:>> x=linspace(0,2*pi,600); t=sin(x)./(cos(x)+eps);plot(x,t);title('tan(x)');axis ([0,2*pi,-50,50]); 图象:程序代码:>> x=linspace(0,2*pi,100); ct=cos(x)./(sin(x)+eps);plot(x,ct);title('cot(x)');axis ([0,2*pi,-50,50]); 图象:4在区间]1,1[-画出函数xy 1sin =的图形. 解:程序代码:>> x=linspace(-1,1,10000);y=sin(1./x); plot(x,y);axis([-1,1,-2,2]) 图象:二维参数方程作图6画出参数方程⎩⎨⎧==t t t y tt t x 3cos sin )(5cos cos )(的图形:解:程序代码:>> t=linspace(0,2*pi,100);plot(cos(t).*cos(5*t),sin(t).*cos(3*t)); 图象:极坐标方程作图8 作出极坐标方程为10/t e r =的对数螺线的图形. 解:程序代码:>> t=0:0.01:2*pi; r=exp(t/10);polar(log(t+eps),log(r+eps)); 图象:90270分段函数作图10 作出符号函数x y sgn =的图形. 解:程序代码:>> x=linspace(-100,100,10000); y=sign(x); plot(x,y);axis([-100 100 -2 2]);函数性质的研究12研究函数)3(log 3)(35x e x x f x -++=在区间]2,2[-上图形的特征. 解:程序代码:>> x=linspace(-2,2,10000);y=x.^5+3*exp(x)+log(3-x)/log(3); plot(x,y); 图象:实验2 极限与连续(基础实验)实验目的 通过计算与作图, 从直观上揭示极限的本质,加深对极限概念的理解. 掌握用 Matlab 画散点图, 以及计算极限的方法. 深入理解函数连续的概念,熟悉几种间断点的图形 特征,理解闭区间上连续函数的几个重要性质.作散点图14分别画出坐标为)10,,2,1(),4,(),,(3222 =+i i i i i i 的散点图, 并画出折线图. 解:散点图程序代码: >> i=1:10; plot(i,i.^2,'.')或:>> x=1:10;y=x.^2;for i=1:10;plot(x(i),y(i),'r')hold onend折线图程序代码:>> i=1:10;plot(i,i.^2,'-x')程序代码:>> i=1:10;plot(i.^2,4*(i.^2)+i.^3,'.')>> i=1:10;plot(i.^2,4*(i.^2)+i.^3,'-x')数列极限的概念16通过动画观察当∞→n 时数列21n a n =的变化趋势.解:程序代码: >> n=1:100; an=(n.^2); n=1:100; an=1./(n.^2); n=1:100; an=1./(n.^2); for i=1:100plot(n(1:i),an(1:i)),axis([0,100,0,1])pause(0.1) %在这个for 循环内每画一个点就暂停0.1s 因此有动画效果 end 图象:函数的极限18在区间]4,4[-上作出函数xx xx x f --=339)(的图形, 并研究 )(lim x f x ∞→ 和 ).(lim 1x f x →解:作出函数x x x x x f --=339)(在区间]4,4[-上的图形 >> x=-4:0.01:4;y=(x.^3-9*x)./(x.^3-x+eps); plot(x,y)从图上看,()f x 在x →1与x →∞时极限为0两个重要极限 20计算极限⎪⎭⎫⎝⎛+→x x x x x sin 11sin lim )1(0 x x e x 2lim )2(+∞→30sin tan lim )3(xx x x -→ x x x 0lim )4(+→ x xx ln cot ln lim )5(0+→ x x x ln lim )6(20+→ xx x x x x sin cos sin lim)7(20-→ 125523lim)8(323+++-∞→x x x x x x x x e e x x x sin 2lim )9(0----→ xx x x cos 110sin lim )10(-→⎪⎭⎫ ⎝⎛ 解:(1)>> limit(x*sin(1/x)+1/x*sin(x))ans =1(2) >> limit(x^2/exp(x),inf) ans = 0(3) >> limit((tan(x)-sin(8))/x^3) ans =NaN(4) >> limit(x^x,x,0,'right') ans =1(5) >> limit(log(cot(x))/log(x),x,0,'right') ans =-1(6) >> limit(x^2*log(x),x,0,'right') ans =016(7) >> limit((sin(x)-x.*cos(x))./(x.^2.*sin(x)),x,0) ans =1/3(8) >> limit((3*x.^3-2*x.^2+5)/(5*x.^3+2*+1),x,inf) ans =3/5(9) >> limit((exp(x)-exp(-x)-2*x)./(x-sin(x))) ans =2(10) >> limit((sin(x)/x).^(1/(1-cos(x)))) ans =exp(-1/3)实验3 导数(基础实验)实验目的 深入理解导数与微分的概念, 导数的几何意义. 掌握用Matlab 求导数与高 阶导数的方法. 深入理解和掌握求隐函数的导数, 以及求由参数方程定义的函数的导数的方法. 导数概念与导数的几何意义22作函数71232)(23+-+=x x x x f 的图形和在1-=x 处的切线. 解:作函数71232)(23+-+=x x x x f 的图形程序代码: >> syms x;>> y=2*x^3+3*x^2-12*x+7; >> diff(y) ans =6*x^2+6*x-12 >> syms x;y=2*x^3+3*x^2-12*x+7; >> f=diff(y) f =6*x^2+6*x-12 >> x=-1;f1=6*x^2+6*x-12 f1 = -12>> f2=2*x^3+3*x^2-12*x+7 f2 = 20>> x=linspace(-10,10,1000);y1=2*x.^3+3*x.^2-12*x+7; y2=-12*(x+1)+20; plot(x,y1,'r',x,y2,'g')求函数的导数与微分24求函数bx ax x f cos sin )(=的一阶导数. 并求.1⎪⎭⎫⎝⎛+'b a f解:求函数bx ax x f cos sin )(=的一阶导数程序代码: >> syms a b x y;y= sin(a*x)*cos(b*x); D1=diff(y,x,1) 答案:D1 =cos(a*x)*a*cos(b*x)-sin(a*x)*sin(b*x)*b求.1⎪⎭⎫ ⎝⎛+'b a f程序代码:>> x=1/(a+b); %直接输入表达式直接带入 >> cos(a*x)*a*cos(b*x)-sin(a*x)*sin(b*x)*b 答案:ans =cos(a/(a+b))*a*cos(b/(a+b))-sin(a/(a+b))*sin(b/(a+b))*b 拉格朗日中值定理26对函数),2)(1()(--=x x x x f 观察罗尔定理的几何意义. (1) 画出)(x f y =与)(x f '的图形, 并求出1x 与.2x 解:程序代码:>> syms x;f=x*(x-1)*(x-2); f1=diff(f) f1 =(x-1)*(x-2)+x*(x-2)+x*(x-1) >> solve(f1) ans =1+1/3*3^(1/2) 1-1/3*3^(1/2)>> x=linspace(-10,10,1000); y1=x.*(x-1).*(x-2);y2 =(x-1).*(x-2)+x.*(x-2)+x.*(x-1); plot(x,y1,x,y2)(2)画出)(x f y =及其在点))(,(11x f x 与))(,(22x f x 处的切线.程序代码:>> syms x; >> f=x*(x-1)*(x-2); >> f1=diff(f) f1 =(x-1)*(x-2)+x*(x-2)+x*(x-1) %本题意在讨论罗尔定理首先回忆一下罗尔定理 >> solve(f1) %(a,b )连续;[a,b]可导;f(a)=f(b)=0;存在ƞ ans = %使得f ’(ƞ)=01+1/3*3^(1/2) %程序就要先求出f (x )的零点 1-1/3*3^(1/2) %再画出f ’(a)f ’(b)的切线 >> x=linspace(-3,3,1000); >> y1=x.*(x-1).*(x-2);>> y2 =(x-1).*(x-2)+x.*(x-2)+x.*(x-1); >> plot(x,y1,x,y2) >> hold on>> x=1+1/3*3^(1/2); >> yx1=x*(x-1)*(x-2) yx1 =-0.3849>> x=1-1/3*3^(1/2);>> yx2=x*(x-1)*(x-2) yx2 =0.3849x=linspace(-3,3,1000); yx1 =-0.3849*x.^0; yx2 =0.3849*x.^0; plot(x,yx1,x,yx2)28求下列函数的导数: (1) 31+=x e y ; 解:程序代码:>> syms x y; y=exp((x+1)^3); D1=diff(y,1) 答案:D1 =3*(x+1)^2*exp((x+1)^3)(2) )]42ln[tan(π+=x y ;解:程序代码:>> syms x;y=log(tan(x/2+pi/4)); D1=diff(y,1) 答案:D1 =(1/2+1/2*tan(1/2*x+1/4*pi)^2)/tan(1/2*x+1/4*pi)(3) x x y sin ln cot 212+=;解:程序代码:>> syms x;y=1/2*(cot(x))^2+log(sin(x)); D1=diff(y,1) 答案:D1 =cot(x)*(-1-cot(x)^2)+cos(x)/sin(x) (4) xy 2arctan21=. 解:程序代码:>> syms x;>> y=sqrt(2)*atan(sqrt(2)/x); >> D1=diff(y,1) 答案:D1 =-2/x^2/(1+2/x^2)一元函数积分学与空间图形的画法实验4 一元函数积分学(基础实验)实验目的 掌握用Matlab 计算不定积分与定积分的方法. 通过作图和观察, 深入理解定积分的概念和思想方法. 初步了解定积分的近似计算方法. 理解变上限积分的概念. 提高应用 定积分解决各种问题的能力.不定积分计算30求.)1(532⎰-dx x x解:程序代码:>> syms x y;>> y=x^2*(1-x^3)^5; >> R=int(y,x) 答案:R =-1/18*x^18+1/3*x^15-5/6*x^12+10/9*x^9-5/6*x^6+1/3*x^332求.arctan 2⎰xdx x 解:程序代码:>> syms x y;>> y=x^2*atan(x); >> R=int(y,x) 答案:R =1/3*x^3*atan(x)-1/6*x^2+1/6*log(x^2+1)定积分计算34 求.)(102⎰-dx x x解:程序代码:>> syms x y; >> y=x-x^2;>> R=int(y,x,0,1) 答案: R =1/6变上限积分 36 画出变上限函数⎰x dt t t 02sin 及其导函数的图形.解:程序代码:>> syms x y t; >> y=t*sin(t^2); >> R=int(y,x,0,x) 答案:R =t*sin(t^2)*x 再求导函数 程序代码:>> DR=diff(R,x,1) 答案:DR =t*sin(t^2)实验5 空间图形的画法(基础实验)实验目的 掌握用Matlab 绘制空间曲面和曲线的方法. 熟悉常用空间曲线和空间曲面 的图形特征,通过作图和观察, 提高空间想像能力. 深入理解二次曲面方程及其图形.一般二元函数作图38作出函数2214y x z ++=的图形.解:程序代码:>> x=linspace(-5,5,500);[x,y]=meshgrid( %mesh 和meshgrid 是生成网格的命令 z=4./(1+x.^2+y.^2); %先命令出一个x 向量组,再生二位网格 mesh(x,y,z); %用meshgrid 再用mesh 命令出各个点的高 xlabel('x-axis'),ylabel('y-axis'),zlabel('z-axis');title('function')40作出函数)94cos(22y x z +=的图形. 解:程序代码:>> x=-10:0.1:10;[x,y]=meshgrid(x);z=cos(4*x.^2+9*y.^2); mesh(x,y,z);xlabel('x-axis'),ylabel('y-axis'),zlabel('z-axis');title('function')讨论:坐标轴选取范围不同时,图形差异很大,对本题尤为明显,如右图为坐标轴[-1,1]二次曲面42作出单叶双曲面1941222=-+z y x 的图形.(曲面的参数方程为 ,tan 3,cos sec 2,sin sec u z v u y v u x === (.20,2/2/πππ≤≤<<-v u ))解:程序代码:>> v=0:pi/100:2*pi; >> u=-pi/2:pi/100:pi/2; >> [U,V]=meshgrid(u,v); >> x=sec(U).*sin(V); >> y=2*sec(U).*cos(V); >> z=3*tan(U); >> surf(x,y,z)44 可以证明: 函数xy z =的图形是双曲抛物面. 在区域22,22≤≤-≤≤-y x 上作出它的图形.解:程序代码:>> x=-2:0.01:2;[x,y]=meshgrid(x); >> z=x.*y;>> mesh(x,y,z);46 画出参数曲面]2,001.0[],4,0[)5/2/ln(tan cos sin sin sin cos ∈∈⎪⎩⎪⎨⎧++===v u u v v z vu y v u x π 的图形.解:程序代码:>> v=0.001:0.001:2;>> u=0:pi/100:4*pi;>> [U,V]=meshgrid(u,v); >> x=cos(U).*sin(V); >> y=sin(U).*sin(V);>> z=cos(V)+log(tan(V/2)+U/5); >> mesh(x,y,z);空间曲线48 作出空间曲线)60(2,sin ,cos π≤≤===t t z t t y t t x 的图形. 解:程序代码:>> syms t;ezplot3(t*cos(t),t*sin(t),2*t,[0,6*pi])-1010-100100xx = t cos(t), y = t sin(t), z = 2 tz50绘制参数曲线 ⎪⎪⎩⎪⎪⎨⎧=+==t z t y t x arctan 211cos 2的图形.解:程序代码:>> t=-2*pi:pi/100:2*pi;x=cos(t).*cos(t);y=1./(1+2*t);z=atan(t); plot3(x,y,z);grid;xlabel('x'),ylabel('y'),zlabel('z')xyz多元函数微积分实验6 多元函数微分学(基础实验)实验目的 掌握利用Matlab 计算多元函数偏导数和全微分的方法, 掌握计算二元函数极值和条件极值的方法. 理解和掌握曲面的切平面的作法. 通过作图和观察, 理解二元 函数的性质、方向导数、梯度和等高线的概念.求多元函数的偏导数与全微分52设),(cos )sin(2xy xy z +=求.,,,222yx z x z y z x z ∂∂∂∂∂∂∂∂∂ 解:程序代码:>> syms x y;S=sin(x*y)+(cos(x*y))^2; D1=diff(S,'x',1); D2=diff(S,'y',1); D3=diff(S,'x',2);D4=diff(S,'y',2); %这是错误的! D1,D2,D3,D4答案: D1 = cos(x*y)*y-2*cos(x*y)*sin(x*y)*yD2 = cos(x*y)*x-2*cos(x*y)*sin(x*y)*xD3 =-sin(x*y)*y^2+2*sin(x*y)^2*y^2-2*cos(x*y)^2*y^2 D4 = -sin(x*y)*x^2+2*sin(x*y)^2*x^2-2*cos(x*y)^2*x^2实验7 多元函数积分学(基础实验)实验目的掌握用Matlab 计算二重积分与三重积分的方法; 深入理解曲线积分、曲面积分的 概念和计算方法. 提高应用重积分和曲线、曲面积分解决各种问题的能力.计算重积分54计算,2dxdy xy D⎰⎰ 其中D 为由,,2y x y x ==+ 2=y 所围成的有界区域.解:程序代码:>> syms x y;int(int(x*y^2,x,2-y,sqrt(y)),y,1,2) 答案:ans =193/120 重积分的应用56求旋转抛物面224y x z --=在Oxy 平面上部的面积.S 解:程序代码:>> int(2*pi*r,r,0,2) 答案: ans =4*pi无穷级数与微分方程实验8 无穷级数(基础实验) 实验目的观察无穷级数部分和的变化趋势,进一步理解级数的审敛法以及幂级数部分和对函数的 逼近. 掌握用Matlab 求无穷级数的和, 求幂级数的收敛域, 展开函数为幂级数以及展 开周期函数为傅里叶级数的方法.数项级数58(1) 观察级数∑∞=121n n的部分和序列的变化趋势.解:程序代码:for i=1:100 s=0; for n=1:i s=s+1/n^2; endplot(i,s,'.');hold on; end(2) 观察级数∑∞=11n n 的部分和序列的变化趋势.>> for i=1:100 s=0; for n=1:i s=s+1/n; endplot(i,s,'.'); hold on; end60 求∑∞=++123841n n n 的值. 解:程序代码:>> syms n;score=symsum(1/(4*n^2+8*n+3),1,inf) %这个用法我比较不熟悉,记一下 答案: score =1/6函数的幂级数展开62求x arctan 的5阶泰勒展开式. >> syms x;>> T5=taylor(atan(x),6) 答案:T5 =x-1/3*x^3+1/5*x^5实验9 微分方程(基础实验)实验目的 理解常微分方程解的概念以及积分曲线和方向场的概念,掌握利用 Matlab 求微分方程及方程组解的常用命令和方法.求解微分方程64求微分方程 22x xe xy y -=+'的通解. 解:程序代码:>> y=dsolve('Dy+2*x*y=x*exp(-x^2)','x') %还记得solve 命令吗?格式solve (‘’) 答案:y = %将可以解一元多次方程而且不用写出等号(1/2*x^2+C1)*exp(-x^2) %本例用于解微分方程,需要声明什么是变量66求微分方程x e y y y x 2cos 52=+'-''的通解. 解:程序代码:>> y=dsolve('D2y-2*Dy+5*y=exp(x)*cos(2*x)','x') %二次倒数‘D2y ’ 答案: y =exp(x)*sin(2*x)*C2+exp(x)*cos(2*x)*C1+1/4*exp(x)*sin(2*x)*x68求微分方程组⎪⎪⎩⎪⎨⎧=--=++02y x dtdy ey x dt dx t 在初始条件0,100====t t y x 下的特解.解:程序代码: %你还别说我真不会接微分方程组>> [x,y]=dsolve('Dx+x+2*y-exp(t)','Dy-x-y','x(0)=1','y(0)=0','t') 答案: x = cos(t)y = 1/2*sin(t)-1/2*cos(t)+1/2*exp(t)70求解微分方程,)1(122/5+=+-x x y dx dy 并作出积分曲线. 解:程序代码:>> syms x yy=dsolve('Dy-2*y/(x+1)-(x+1)^(5/2)','x') 答案:y =(2/3*(x+1)^(3/2)+C1)*(x+1)^2 做积分曲线 由>> syms x yx=linspace(-5,5,100); C=input('请输入C 的值:'); y=(2/3*(x+1).^(3/2)+C).*(x+1).^2; plot(x,y)例如对应有: 请输入C 的值:2 请输入C 的值:20矩阵运算与方程组求解实验10 行列式与矩阵实验目的掌握矩阵的输入方法. 掌握利用Matlab 对矩阵进行转置、加、减、数乘、相乘、乘方等运算, 并能求矩阵的逆矩阵和计算方阵的行列式.矩阵A 的转置函数Transpose[A]72 求矩阵⎪⎪⎪⎪⎪⎭⎫⎝⎛411365243271的转置. 解:程序代码:>> A=[1,7,2;3,4,2;5,6,3;1,1,4]; >> Sove=A' 答案:Sove =1 3 5 1 7 4 6 12 234 矩阵线性运算 73设,291724,624543⎪⎪⎭⎫⎝⎛=⎪⎪⎭⎫⎝⎛=B A 求.24,A B B A -+ 解:程序代码:>> A=[3,4,5;4,2,6]; B=[4,2,7;1,9,2]; S1=A+B S2=4*B-2*A 答案:S1 =7 6 12 5 11 8 S2 =10 0 18 -4 32 -474设,148530291724,36242543⎪⎪⎪⎪⎪⎭⎫ ⎝⎛=⎪⎪⎭⎫ ⎝⎛=mb ma 求矩阵ma 与mb 的乘积. 解:程序代码:>> ma=[3,4,5,2;4,2,6,3];>> mb=[4,2,7;1,9,2;0,3,5;8,4,1];>> Sove=ma*mb答案:Sove =32 65 5642 56 65矩阵的乘法运算 75设,101,530291724⎪⎪⎪⎭⎫ ⎝⎛=⎪⎪⎪⎭⎫ ⎝⎛=B A 求AB 与,A B T 并求.3A解:程序代码:>> A=[4 2 7;1 9 2;0 3 5];B=[1;0;1];>> AB=A*BAB =1135>> BTA=B'*ABTA =4 5 12>> A3=A^3A3 =119 660 555141 932 44454 477 260求方阵的逆76 设,5123641033252312⎪⎪⎪⎪⎪⎭⎫⎝⎛=A 求.1-A解:程序代码:>> A=[2,1,3,2;5,2,3,3;0,1,4,6;3,2,1,5];Y=inv(A)答案:Y =-1.7500 1.3125 0.5000 -0.68755.5000 -3.6250 -2.0000 2.37500.5000 -0.1250 0.0000 -0.1250-1.2500 0.6875 0.5000 -0.312577 设,221331317230,5121435133124403⎪⎪⎪⎪⎪⎭⎫⎝⎛=⎪⎪⎪⎪⎪⎭⎫ ⎝⎛=B A 求.1B A -解:程序代码:>> A=[3 0 4 4 ;2 1 3 3 ;1 5 3 4;1 2 1 5]; B=[0 3 2 ;7 1 3;1 3 3 ;1 2 2];Solve=A'*B %这一步和题目要求不符,应是Solve=inv(A)*B答案:Solve =16 16 1714 20 2225 26 2830 37 3978 解方程组⎪⎩⎪⎨⎧-=-+=+-=++.2442,63,723z y x z y x z y x解:程序代码:>> A=[3 2 1;1 -1 3;2 4 -4]; b=[7 6 -2];>> A\b' %也可以用inv (A )*b ’表达答案:ans =1.00001.00002.0000求方阵的行列式79 求行列式 .3351110243152113------=D 解:程序代码:>> A=[3,1,-1,2;-5,1,3,-4;2,0,1,-1;1,-5,3,-3];D=det(A)答案:D =4080求.11111111111122222222d d d d c cc c b bb b a a a a D ++++= 解:程序代码:>> syms a b c d;D=[a^2+1/a^2 a 1/a 1;b^2+1/b^2 b 1/b 1;c^2+1/c^2 c 1/c 1;d^2+1/d^2 d 1/d 1];det(D)答案:ans =-(-c*d^2*b^3+c^2*d*b^3-c^3*d^2*a+c^3*d*a^2*b^4+c*d^2*a^3-c^3*d^2*a*b^4-c^2*d*a^3-c*d^2*b^3*a^4+c^2*d*b^3*a^4+c^3*d^2*b*a^4-c^3*d*b^2*a^4-c^2*d^3*b*a^4+c*d^3*b^2*a^4+c*d ^2*a^3*b^4-c^2*d*a^3*b^4+c^3*d^2*b-c^3*d*b^2-c^2*d^3*b+c*d^3*b^2+c^3*d*a^2+c^2*d^3*a-c *d^3*a^2-b*d^2*a^3+b^2*d*a^3+b^3*d^2*a-b^3*d*a^2-b^2*d^3*a+b*d^3*a^2+b*c^2*a^3-b^2*c*a ^3-b^3*c^2*a+b^3*c*a^2+b^2*c^3*a-b*c^3*a^2+c^2*d^3*a*b^4-c*d^3*a^2*b^4-b*d^2*a^3*c^4+b ^2*d*a^3*c^4+b^3*d^2*a*c^4-b^3*d*a^2*c^4-b^2*d^3*a*c^4+b*d^3*a^2*c^4+b*c^2*a^3*d^4-b^2*c*a^3*d^4-b^3*c^2*a*d^4+b^3*c*a^2*d^4+b^2*c^3*a*d^4-b*c^3*a^2*d^4)/a^2/c^2/d^2/b^281 计算范德蒙行列式.1111145444342413534333231252423222154321x x x x x x x x x x x x x x x x x x x x解:程序代码:>> syms x1 x2 x3 x4 x5;>> A=[1,1,1,1,1;x1,x2,x3,x4,x5;x1^2,x2^2,x3^2,x4^2,x5^2;x1^3,x2^3,x3^3,x4^3,x5^3;x1^4,x2^4,x3^4,x4^4,x5^4];>> DC=det(A);>> DS=simple(DC)答案:DS =(-x5+x4)*(x3-x5)*(x3-x4)*(-x5+x2)*(x2-x4)*(x2-x3)*(-x5+x1)*(x1-x4)*(x1-x3)*(x1-x2)82 设矩阵 ,60975738723965110249746273⎪⎪⎪⎪⎪⎪⎭⎫⎝⎛----=A 求.),(|,|3A A tr A解:程序代码:>> A=[3,7,2,6,-4;7,9,4,2,0;11,5,-6,9,3;2,7,-8,3,7;5,7,9,0,-6];>> D=det(A),T=trace(A),A3=A^3答案:D =11592T =3A3=726 2062 944 294 -3581848 3150 26 1516 2281713 2218 31 1006 4041743 984 -451 1222 384801 2666 477 745 -125向量的内积83 求向量}3,2,1{=u 与}0,1,1{-=v 的内积.解:程序代码:>> u=[1 2 3];v=[1 -1 0];solve=dot(u,v)答案:solve =-184设,001001⎪⎪⎪⎭⎫ ⎝⎛=λλλA 求.10A 一般地?=k A (k 是正整数).解:程序代码:>> syms r;>> A=[r,1,0;0,r,1;0,0,r];>> A^10答案:ans =[ r^10, 10*r^9, 45*r^8][ 0, r^10, 10*r^9][ 0, 0, r^10]85.求⎪⎪⎪⎪⎪⎪⎭⎫ ⎝⎛+++++a a a aa 1111111111111111111111111的逆. 解:程序代码:>> syms aA=[1+a,1,1,1,1;1,1+a,1,1,1;1,1,1+a,1,1;1,1,1,1+a,1;1,1,1,1,1+a];solve=inv(A)答案:solve =[ 1/a*(a+4)/(a+5), -1/a/(a+5), -1/a/(a+5), -1/a/(a+5), -1/a/(a+5)][ -1/a/(a+5), 1/a*(a+4)/(a+5), -1/a/(a+5), -1/a/(a+5), -1/a/(a+5)][ -1/a/(a+5), -1/a/(a+5), 1/a*(a+4)/(a+5), -1/a/(a+5), -1/a/(a+5)][ -1/a/(a+5), -1/a/(a+5), -1/a/(a+5), 1/a*(a+4)/(a+5), -1/a/(a+5)][ -1/a/(a+5), -1/a/(a+5), -1/a/(a+5), -1/a/(a+5), 1/a*(a+4)/(a+5)] 实验11 矩阵的秩与向量组的极大无关组实验目的 学习利用Matlab 求矩阵的秩,作矩阵的初等行变换; 求向量组的秩与极大无关组. 求矩阵的秩86 设,815073*********⎪⎪⎪⎭⎫⎝⎛-------=M 求矩阵M 的秩.解:程序代码:>> M=[3,2,-1,-3,-2;2,-1,3,1,-3;7,0,5,-1,-8];R=rank(M)答案:R=2向量组的秩87求向量组)0,3,0,2(),2,5,4,0(),1,1,2,1(231=--=-=ααα的秩.解:程序代码:>> A=[1,2,-1,1;0,-4,5,-2;2,0,3,0];R=rank(A)答案:R =288向量组)7,5,1,3(),5,4,3,1(),1,1,1,1(),3,2,1,1(4321==-==αααα是否线性相关?解:由>> A=[1 1 2 3;1 -1 1 1;1 3 4 5;3 1 5 7];rank(A)ans = 3即rank(A)=3 小于阶数489向量组)3,1,1(),2,1,3(),7,2,2(321=-==ααα是否线性相关?解:由>> A3=[2,2,7;3,-1,2;1,1,3];R=rank(A3)得 R = 3即rank(A3)=3 等于阶数3故向量组线性无关。
2015-2016数学实验练习题一、选择题1、清除Matlab工作空间(wordspace)变量的命令就是( B )A、clcB、clearC、clfD、delete2、清除当前屏幕上显示的所有内容,但不清除工作空间中的数据的命令就是( A )A、clcB、clearC、clfD、delete3、用来清除图形的命令( C )A、clcB、clearC、clfD、delete4、在MATLAB程序中,使命令行不显示运算结果的符号就是 ( A )A、;B、%C、#D、&5、在MATLAB程序中,可以将某行表示为注释行的符号就是 ( B )A、;B、%C、#D、&6、在循环结构中跳出循环,执行循环后面代码的命令为 ( B )A、returnB、breakC、continueD、Keyboard7、在循环结构中跳出循环,但继续下次循环的命令为( C )A、returnB、breakC、continueD、Keyboard8、MATLAB中用于声明全局变量的关键字就是( C )A、infB、symsC、globalD、function9、用户可以通过下面哪项获得指令的使用说明( A )A、helpB、loadC、demoD、lookfor10.在MATLAB命令窗口中键入命令S=zoros(3);可生成一个三行三列的零矩阵,如果省略了变量名S,MATLAB表现计算结果将用下面的哪一变量名做缺省变量名( A )A、ans;B、pi;C、NaN;D、Eps、11、9/0的结果就是( B )A、NAN;B、Inf;C、eps;D、012.在MATLAB中程序或语句的执行结果都可以用不同格式显示,将数据结果显示为分数形式,用下面哪一条命令语句( D )A、format long;B、format long e;C、format bank;D、fromat rat13、下列MATLAB命令中就是构造1行3列的(-1,1)均匀分布随机矩阵的命令的就是( D )A 、 randn(1,3);B 、 rand(1,3);C 、 ones(3);D 、 以上都不对14、 产生四维元素都为1矩阵的语句为( A )A 、 ones(4)B 、 eye(4)C 、 zeros(4)D 、 rand(4)15、 用round 函数对数组[2、48 6、39 3、93 8、52]取整,结果为 ( C )A 、 [2 6 3 8]B 、 [2 6 4 8]C 、 [2 6 4 9]D 、 [3 7 4 9]16、 y=dsolve(‘Dy=1/(1+x^2)-2*y^2’,’y(0)=0’,’x ’); ezplot(y)的功能就是( A )A 、 求微分方程特解并绘图;B 、 解代数方程;C 、 求定积分;D 、求微分方程通解、17、 MATLAB 命令roots([1,0,0,-1])的功能就是 ( D )A 、 产生向量[1,0,0,1];B 、 求方程310x +=的根;C 、 求多项式31x -的值;D 、 求方程310x -=的根。
18.在MATLAB 命令窗口中,键入命令A=[1,2,3;4,5,6;7,8,0];A(2,:)*A(:,2);将出现的结果就是( B )A 、 15;B 、 81;C 、 36;D 、 30、19、 在MATLAB 命令窗口中键入命令A=[1 4 2;3 1 2;6 1 5]; 则det(A(1:2,2:3)、*A(2:3,2:3))的结果就是 ( C )A 、 -143;B 、 60;C 、 36;D 、 -1920.i=2; a=2i; b=2*i; c=2*sqrt(-1); 程序执行后a, b, c 的值分别就是 ( C ) A 、 a=4, b=4, c=2、0000i B 、 a=4, b=2、0000i, c=2、0000iC 、 a=2、0000i, b=4, c=2、0000iD 、 a=2、0000i, b=2、0000i, c=2、0000i21、 在MATLAB 命令窗口中键入命令A=[1 4;3 5]; B=[2 3;3 8]; 则A 、*B 的结果就是( A )A 、 ⎪⎪⎭⎫ ⎝⎛409122;B 、⎪⎪⎭⎫⎝⎛49213514; C 、 ⎪⎪⎭⎫⎝⎛4921122; D 、 ⎪⎪⎭⎫⎝⎛4093514、 22、 在MATLAB 命令窗口中键入命令A=[1 4;3 5]; B=[2 3;3 8]; 则A*B 的结果就是( B )A 、 ⎪⎪⎭⎫⎝⎛409122; B 、⎪⎪⎭⎫⎝⎛49213514; C 、 ⎪⎪⎭⎫⎝⎛4921122; D 、 ⎪⎪⎭⎫⎝⎛4093514、 23、 下面有关MATLAB 变量名与函数名的说法错误的就是( D )A 、 变量名的第一个字符必须就是一个英文字母B 、 变量名可由英文字母、数字与下划线混合组成C、变量名不得包含空格与标点,但可以有下连字符D、变量名与函数名对于英文的大小使用没有区别24、正确表达命题A或B都大于C的逻辑表达式应该用下面哪一行( )A、 A > C;B、B>C;C、 A >C & B >C;D、 A >C | B >C。
25、输入方阵A的数据,在MATLAB中计算A的逆矩阵用下面哪一命令( D )A、det(A) B.eig(A) C.diag(A) D.inv(A)26、输入方阵A的数据,在MATLAB中计算A的特征值用下面哪一命令( D )A、inv(A)B、diag(A)C、det(A)D、eig(A)27、输入方阵A的数据,在MATLAB中计算A的行列式的值用下面哪一命令( C )A、inv(A)B、diag(A)C、det(A)D、eig(A)28、输入方阵A的数据,在MATLAB中计算A的秩用下面哪一命令( C )A、inv(A)B、diag(A)C、rank(A)D、eig(A)29、命令A=magic(3)创建3阶魔方矩阵,求A的特征值绝对值最小用( A )A、min(abs(eig(A)));B、min(eig(abs(A)));C、min(eig(A));D、min(abs(A));30、下列哪个变量的定义就是不合法的( A )A、abcd-3B、xyz_3C、abcdefD、x3yz31、下列变量名中合法的就是( A )A、char_1B、x*yC、x\yD、end32、在MATLAB中下列数值的表示不正确的就是( B)A.+99B.1、3e-5C.2-3*e^2D.3-2*pi33、图形窗口(Figure)显示网格的命令就是( B )A、axis onB、grid onC、box onD、hold on34、为了使两个plot的图形在同一个坐标显示,可以使用哪个命令进行图形保持( D )A、axis onB、grid onC、box onD、hold on35.假设在当前工作目录中有一名为exa1的M文件,其中有些错误。
现在为了修改它,用键盘命令将这一文件打开编辑,使用下面的命令不对的就是( )A、edit exa1;B、edit exa1、m;C、type exa1;D、open exa136、如果x=1:-3:-18,则x(3)与x(7)分别就是( D )A.-6,-14B. -5,-14C. -8,-17D. -5,-1737、已知a=2:2:8, b=2:5,下面的运算表达式中,出错的为( C )A 、 a' *bB 、 a 、*bC 、 a*bD 、 a-b38、 设A=[2 4 3; 5 3 1; 3 6 7],则sum(A),length(A)与size(A)的结果( D )A 、 [10 13 11] 9 [3 3]B 、 [9 9 16] 3 [3 3]C 、 [9 9 16] 9 [3 3]D 、 [10 13 11] 3 [3 3]39、 输入矩阵A=[1 3 2;3 7 -5;5 6 9]; 使用全下标方式用A(3,2)取出元素“6”,使用单下标方式用哪个命令可以取得 ( D )A 、 A(4)B 、 A(7)C 、 A(8)D 、 A(6)40、 MATLAB 命令A=rand(5,5);创建55)(A ⨯=ij a ,求∑=51||max i ij ja 用 ( B )A 、 max(sum(abs(A)))B 、 max(sum(abs(A ’)))C 、 max(sum(A)))D 、sum(max(A))41、 MATLAB 命令x=[1,2,4,5,9]; mean(x)的计算结果就是 ( B )A 、 4B 、 4、2C 、 4、5D 、 2142、 运行如下程序后, 命令窗口(command windows)显示的结果为( D )A=[13,-56,78; 25,63,-735; 78,25,563; 1,0,-1]; y=max(max(A))A 、 y=564B 、y=9C 、y=-735D 、y=56343、 运行如下程序后, 命令窗口(command windows)显示的结果为( A )A=[13,0,-5;8,-6,2;15,3,7]; y=max(A)A 、 y=15 3 7B 、y=13 8 15C 、y=15 8 7D 、y=13 3 744、 命令factor()用于分解因式,syms x; f=4*x^3+9*x^2-30*x; factor(diff(f))的结果就是( B )A 、 (x-1)*(2*x-5)B 、 6*(x-1)*(2*x+5)C 、 6*(x+1)*(2*x+5)D 、 (x+1)*(2*x-5)45、 MATLAB 表达式2*2^3^2的结果就是 ( A ) A 、 128 B 、 4096 C 、 262144 D 、 256 46、 已知函数文件如下,则factor(4)=( C )function f=factor(n) if n<=1 f=1; elsef=factor(n-1)*n;endA 、 4B 、 12C 、 24D 、4847、 下面的程序执行后array 的值为 ( A )for k=1:10 if k>6 break; else array(k) = k; end endA 、 array = [1, 2, 3, 4, 5, 6]B 、 array = [1, 2, 3, 4, 5, 6, 7, 8, 9, 10]C 、 array =6D 、 array =10、48、 运行如下程序后, 命令窗口(command windows)显示的结果为 ( B )s=0;a=[12,13,14;15,16,17;18,19,20]; for k=a s=s+k; end disp(s');A 、 144B 、 39 48 57C 、145D 、 45 48 5149、 十二属相为“鼠牛虎兔龙蛇马羊猴鸡狗猪”,命令k=mod(2008,12)+1的结果就是 ( D )A 、 k 指向第二动物牛;B 、 k 指向第三动物虎;C 、 k 指向第四动物兔;D 、 k 指向第五动物龙。