matlab习题2

合集下载

matlab 部分答案

matlab 部分答案


D=A.*B


MATLAB: A=[1,4,8,13; -3,6,-5,-9; 2,-7,-12,-8];


8 13 ⎤ ⎡ 5 4 3 −2 ⎤ ⎡1 4 ⎥ ⎢ ⎥ 10. 设 A = ⎢ −3 6 −5 −9 ⎥ , B = ⎢ ⎢ 6 −2 3 − 8 ⎥ ⎢ ⎢ ⎣ −1 3 −9 7 ⎥ ⎦ ⎣ 2 −7 −12 −8⎥ ⎦
2
+ y2 )
,求定义域 x = [−2, 2] , y = [−2, 2] 内的 z 值(网格取 0.1 见方) ,并绘制
15.设 x = cos(t ) , y = sin( Nt + α ) ,若 N = 2 , α = 0, π / 3, π / 2, π ,在 4 个子图中分别画出其
MATLAB: t=0:0.01:3; N=2;
4. 用三点公式计算 y=f(x)在 x=1.0,1.2 处的导数值, f(x)值由下表给出; x f(x)
5
1.0 0.25
4
1.1
ww w
1.2
2
.k
1.3 0.1890 1.4 0.1736

7. 求解下列线性常微分方程的解析解。

6.设方程的根为 x = [−3, −5, −8, −9] ,求它们对应的 x 多项式的系数。
π


=∫
0
17.设 ( X , Y ) 的概率密度为

⎧ d 2 x dy + − x = et ⎪ ⎪ dt 2 dt ⎨ 2 ⎪ d y + dx + y = 0 ⎪ ⎩ dt 2 dt
dy 。 dx

第二次练习题 MATLAB应用基础 练习题

第二次练习题 MATLAB应用基础 练习题

第二章 MATLAB 应用基础 练习题1、下列变量名中____A _________是合法变量。

A 、pi,exe_01B 、x*y,x1C 、1a,ifD 、abs, b.m2、已知a=0:4, b=1:5, 下面的运算表达式出错的为_______D ______。

A 、a+bB 、a ./bC 、a ’*bD 、a*b3、将矩阵A=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡987654321用________D_____命令可以变为A=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡741852963。

A 、A ’ B 、det(A) C 、inv(A) D 、rot90(A)4、已知x 为一个向量,计算ln(x)的运算为______B______。

A 、ln(x))B 、log (x )C 、Ln(x)D 、log10(x)5、____A____产生均匀分布的随机矩阵,元素取值范围0.0~1.0。

A 、 rand(m,n)B 、eye(m,n)C 、magic(N)D 、randn(m,n)6、____D____产生正态分布的随机矩阵。

A 、 rand(m,n)B 、eye(m,n)C 、magic(N)D 、randn(m,n)7、已知x 为一个向量,计算其余弦函数的运算为______C ______。

A 、COS (X )B 、COS (x )C 、cos(x)D 、cosx8、用“from :step :to ”方式得到从0到4π步长为0.4π的变量;使用linspace 函数生成向量从0到2*pi 等分成100个点的变量;使用logspace 函数生成向量从1到100对数等分成3个点的变量。

9、求矩阵⎥⎦⎤⎢⎣⎡4321的转置矩阵、逆矩阵、矩阵的秩、矩阵的行列式值、矩阵的三次幂、矩阵的特征值和特征向量。

10.我国人口按照2000年第五次全国人口普查的结果为12.9533亿,如果年增长率为1.07%,求公元2010年末的人口数。

计算人口的公式为:P1=P0(1+r)^n ,其中:P1为几年后的人口,P0为人口初值,r 为年增长率,n 为年数。

[整理版]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课后习题

一行元素,并给出出错信息。
clear
a=rand(5,6);
n=input('请输入一个数:');
if
n>5
n=5
a(n,:) disp('输入超出范围'); else a(n,:) end 请输入一个数:4 ans =
请输入一个数:7 n=
5
ans =
输入超出范围 7.已知 (1)当 f(n)=n+10ln(n^2+5)时,y 的值是多少。 function a=n2a(n) a=n+10*log(n*n+5) y=n2a(40)/(n2a(30)+n2a(20)) a=
数来实现。
a=input('请输入二十个数');
%用循环结构实现
max=a(1);
min=a(1);
for i=1:20; if max<a(i); max=a(i); end if min>a(i); min=a(i); end
end max min 请输入二十个数[ 3 4 6 5 3 44 5 5 5 6 6 5 5 4 3 32 2 2 2 6]
case{9,10}
m='A';
case{8}
m='B';
case{7}
m='C';
case{6}
m='D';
casenum2cell(0:5)
m='E';
otherwise
m='输入出错';
end
disp(m);
请输入一个百分制分数:98
A
3.输入 20 个数,求其中最大数和最小数。要求分别用循环结构和调用 MATLAB 的 max 函数、min 函

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。

完整word版,Matlab课后习题

完整word版,Matlab课后习题

习题 11. 执行下列指令,观察其运算结果, 理解其意义: (1) [1 2;3 4]+10-2i(2) [1 2; 3 4].*[0.1 0.2; 0.3 0.4] (3) [1 2; 3 4].\[20 10;9 2] (4) [1 2; 3 4].^2 (5) exp([1 2; 3 4]) (6)log([1 10 100]) (7)prod([1 2;3 4])(8)[a,b]=min([10 20;30 40]) (9)abs([1 2;3 4]-pi)(10) [1 2;3 4]>=[4,3;2 1](11)find([10 20;30 40]>=[40,30;20 10])(12) [a,b]=find([10 20;30 40]>=[40,30;20 10]) (提示:a 为行号,b 为列号) (13) all([1 2;3 4]>1) (14) any([1 2;3 4]>1) (15) linspace(3,4,5) (16) A=[1 2;3 4];A(:,2)2. 执行下列指令,观察其运算结果、变量类型和字节数,理解其意义: (1) clear; a=1,b=num2str(a),c=a>0, a= =b, a= =c, b= =c (2) clear; fun='abs(x)',x=-2,eval(fun),double(fun)3. 本金K 以每年n 次,每次p %的增值率(n 与p 的乘积为每年增值额的百分比)增加,当增加到rK 时所花费的时间为)01.01ln(ln p n rT +=(单位:年)用MA TLAB 表达式写出该公式并用下列数据计算:r =2, p =0.5, n =12.4.已知函数f (x )=x 4-2x 在(-2, 2)内有两个根。

取步长h =0.05, 通过计算函数值求得函数的最小值点和两个根的近似解。

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习题——2

matlab习题——2

matlab习题——2第二章1。

%3X3矩阵a1=eye(3)a2=ones(3)a3=zeros(3)a4=1-2*rand(3)a5=2*randn(3)+1%15X8矩阵b1=eye(15,8)b2=ones(15,8)b3=zeros(15,8)b4=1-2*rand(15,8) b5=2*randn(15,8)+1 2。

%生成a矩阵x=diag([2,-7,8]); xx=diag([3,5],-1); a=rot90(x+xx)%生成b矩阵x=diag([2 5 8]);xx=rot90(diag([7 0 4])); b=x+xx%变成行向量ah=reshape(a,1,9)bh=reshape(b,1,9)3。

a=5-10*rand(50,2);b=round(10*a)/10;c=num2str(b,'%+10.1f') 4。

alpha=-pi:pi/180:pi;x=sin(alpha);y=cos(alpha);%画曲线plot(alpha,x,'b-')hold onlot(alpha,y,'r-')5。

a=rand(10);b=a>=0.6;total=sum(sum(b))6。

a=randn(10);b=and(a<.5,a>-.5);total=sum(sum(b))7。

%输入a和ba=input('please input a:'); b=input('please input b:'); %分支结构if and(a<1,b<=.5)disp('语句1')elseif and(a<1,b>.5)disp('语句2')elseif and(a>1,b<=.5)disp('语句3')elsedisp('语句4')end8。

数字图像处理及应用(MATLAB)第2章习题答案

数字图像处理及应用(MATLAB)第2章习题答案

7.平均值说明f (x ,y )的平均值等于其傅里叶变换F (u ,v )在频率原点的值F (0,0)。

2-3证明离散傅里叶变换的频率位移和空间位移性质。

证明:)(2101),(1),(NvyM ux j M x N y e y x f MN v u F +--=-=∑∑=π),(),(1),(100)(21010)(2)(21010000v v u u F dxdy ey x f MNe ey x f MN y Nv v x M u u j M x N y N yv M x u j Nvy M ux j M x N y --==-+---=-=++--=-=∑∑∑∑πππ因为()()v u F y x f ,,⇔ 所以 ),(),(00)(200v v u u F e y x f N y v M x u j --⇔+π2-4小波变换是如何定义的?小波分析的主要优点是什么?小波之所以小,是因为它有衰减性,即是局部非零的;而称为波,则是因为它有波动性,即其取值呈正负相间的振荡形式,将)(2R L 空间的任意函数f (t )在小波基下展开,称其为函数f (t )的连续小波变换。

小波变换是时间(空间)频率的局部化分析,它通过伸缩平移运算对信号(函数)逐步进行多尺度细化,最终达到高频处时间细分,低频处频率细分,能自动适应时频信号的要求从而可聚焦到信号的任意细节,解决了Fourier 变换的困难问题。

2-5 在图像缩放中,采用最近邻域法进行放大时,如果放大倍数太大,可能会出现马赛克效应,这个问题有没有办法解决,或者有所改善。

可以利用线性插值法,当求出的分数地址与像素点不一致时,求出周围四个像素点的距离比,根据该比率, 由四个邻域的像素灰度值进行线性插值。

2-6 复合变换的矩阵等于基本变换的矩阵按顺序依次相乘得到的组合矩阵。

即,T=T N T N-1…T 1。

问矩阵顺序的改变能否影响变换的结果。

矩阵顺序的改变不会影响变换的结果。

matlab习题及答案

matlab习题及答案

matlab习题及答案《Matlab习题及答案:提升编程技能,解决实际问题》Matlab是一种强大的数值计算和数据分析工具,它被广泛应用于工程、科学和其他领域的计算和模拟。

为了帮助大家提升编程技能,解决实际问题,我们为大家准备了一些Matlab习题及答案,希望能够帮助大家更好地掌握Matlab的使用。

1. 习题一:编写一个Matlab程序,计算斐波那契数列的前20个数字,并将结果打印出来。

答案:下面是一个简单的Matlab程序,用于计算斐波那契数列的前20个数字。

```matlaba = 0;b = 1;fib = zeros(1, 20);fib(1) = a;fib(2) = b;for i = 3:20fib(i) = fib(i-1) + fib(i-2);enddisp(fib);```2. 习题二:编写一个Matlab程序,求解一个二次方程ax^2 + bx + c = 0的根。

答案:下面是一个简单的Matlab程序,用于求解二次方程的根。

```matlaba = 1;b = -3;c = 2;delta = b^2 - 4*a*c;if delta > 0x1 = (-b + sqrt(delta))/(2*a);x2 = (-b - sqrt(delta))/(2*a);disp(['The roots are ', num2str(x1), ' and ', num2str(x2)]);elseif delta == 0x = -b/(2*a);disp(['The root is ', num2str(x)]);elsedisp('The equation has no real roots');end```通过以上两个习题及答案的示例,我们可以看到Matlab的强大功能和灵活性。

通过练习这些习题,我们可以更好地掌握Matlab的基本语法和常用函数,从而在实际问题中更快更准确地解决数值计算和数据分析的挑战。

北航Matlab教程(R2016a)习题2解答 2-推荐下载

北航Matlab教程(R2016a)习题2解答 2-推荐下载

习题21.说出以下四条指令产生的结果各属于哪种数据类型,是“双精度”对象,还是“符号”对象?3/7+0.1, sym(3/7+0.1), vpa(sym(3/7+0.1))a=class(3/7+0.1)%双精度b=class(sym(3/7+0.1))%符号c=class(vpa(sym(3/7+0.1),4))%符号d=class(vpa(sym(3/7+0.1)))%符号2.在不加专门指定的情况下,以下符号表达式中的哪一个变量被认为是独立自由变量。

sym('sin(w*t)') , sym('a*exp(-X)' ) , sym('z*exp(j*th)')a=sym('sin(w*t)');symvar(a)b=sym('a*exp(-X)');symvar(b)c=sym('z*exp(j*th)');symvar(c)3.求以下两个方程的解:(提示:关于符号变量的假设要注意)(1)试写出求三阶方程05.443=-x 正实根的程序。

注意:只要正实根,不要出现其他根。

x=sym('x','positive');f=x^3-44.5;x=solve(f,x)(2)试求二阶方程在时的根。

022=+-a ax x 0>a a=sym('a','positive');syms x;f=x^2-a*x+a^a;x=solve(f,x)4.观察一个数(在此用@记述)在以下四条不同指令作用下的异同:a = @ ,b = sym( @ ),c = sym( @ ,'d ' ), d = sym( '@ ' )在此,@ 分别代表具体数值 7/3 , pi/3 , pi*3^(1/3) ;而异同通过vpa(abs(a-d)) , vpa(abs(b-d)) , vpa(abs(c-d))等来观察。

matlab试题与答案

matlab试题与答案

例8 绘制三维曲面图z=sin(x+sin(y))-x/10。 程序如下: [x,y]=meshgrid(0:0.25:4*pi); z=sin(x+sin(y))-x/10; mesh(x,y,z); axis([0 4*pi 0 4*pi -2.5 1]);
1.编程实现当α∈[-π,π],间隔为1°时求解 正弦和余弦的值。 2.利用rand函数产生(0,1)间均匀分布的 10×10随机矩阵A,然后统计A中大于等于 0.6的元素个数。 3.利用randn函数产生均值为0,方差为1的 10×10正态分布随机阵A,然后统计A中大 于-0.5且小于0.5的元素个数。
调用tran.m的命令文件main1.m: x=input('Please input x=:'); y=input('Please input y=:'); [rho,the]=tran(x,y); rho the
例1 在0≤x≤2p区间内,绘制曲线 y=2e-0.5xcos(4πx) 程序如下: x=0:pi/100:2*pi; y=2*exp(-0.5*x).*cos(4*pi*x); plot(x,y),grid on;
4.编程实现下表功能
b a A<1 a≥1 b≤0.5 语句 1 语句 3 B>0.5 语句 2 语句 4
5.有一矩阵A,找出矩阵中其值大于1的元素,并 将它们重新排列成列向量B。 6.在一测量矩阵A(100×3)中,存在有奇异值(假设 大于100的值认为是奇异值),编程实现删去奇异 值所在的行。 7.在给定的100×100矩阵中,删去整行内容全为0 的行,删去整列内容全为0的列。
0 0 8 a 0 7 5 2 3 0

MATLAB数据及其运算 习题2

MATLAB数据及其运算 习题2

第2章 MATLAB数据及其运算习题2一、选择题1.下列可作为MATLAB合法变量名的是()。

DA.合计 B.123 C.@h D.xyz_2a 2.下列数值数据表示中错误的是()。

CA.+10 B.1.2e-5 C.2e D.2i 3.使用语句t=0:7生成的是()个元素的向量。

AA.8 B.7 C.6 D.54.执行语句A=[1,2,3;4,5,6]后,A(3)的值是()。

BA.1 B.2 C.3 D.4 5.已知a为3×3矩阵,则a(:,end)是指()。

DA.所有元素 B.第一行元素 C.第三行元素 D.第三列元素6.已知a为3×3矩阵,则运行a (1)=[]后()。

AA.a变成行向量 B.a变为2行2列C.a变为3行2列 D.a变为2行3列7.在命令行窗口输入下列命令后,x的值是()。

B>> clear>> x=i*jA.不确定B.-1 C. 1 D.i*j 8.fix(354/100)+mod(354,10)*10的值是()。

DA.34 B.354 C.453 D.43 9.9.下列语句中错误的是()。

BA.x==y==3 B.x=y=3 C.x=y==3 D.y=3,x=y10.find(1:2:20>15)的结果是()。

CA.19 20 B.17 19C.9 10 D.8 911.输入字符串时,要用()将字符括起来。

CA.[ ] B.{ } C.' ' D." " 12.已知s='显示"hello"',则s的元素个数是()。

A A.9 B.11 C.7 D.18 13.eval('sqrt(4)+2')的值是()。

BA.sqrt(4)+2 B.4 C.2 D.2, 2 14.有3×4的结构矩阵student,每个结构有name(姓名)、scores (分数)两个成员,其中scores是以1×5矩阵表示的5门课的成绩,那么要删除第4个学生的第2门课成绩,应采用的正确命令是()。

MATLAB作业2参考答案(2018年)

MATLAB作业2参考答案(2018年)

MA TLAB 作业二参考答案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_analyticalx =[ 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-.1418420333][-.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-.5086089095][-.3312192394*a-1.604263460, -.3312192394*a+2.435364854, -.3312192394*a+.3867176824] >> A*x-B[ 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课后习题集附标准答案

第2章MATLAB概论1、与其他计算机语言相比较,MA TLAB 语言突出的特点是什么?答:起点高、人机界面适合科技人员、强大而简易的作图功能、智能化程度高、功能丰富,可扩展性强.2、MA TLAB 系统由那些部分组成?答:开发环境、MATLAB数学函数库、MATLAB语言、图形功能、应用程序接口3、安装MATLAB 时,在选择组件窗口中哪些部分必须勾选,没有勾选的部分以后如何补安装?答:在安装MATLAB时,安装内容由选择组件窗口中各复选框是否被勾选来决定,可以根据自己的需要选择安装内容,但基本平台(即MATLAB选项)必须安装. 第一次安装没有选择的内容在补安装时只需按照安装的过程进行,只是在选择组件时只勾选要补装的组件或工具箱即可. 矚慫润厲钐瘗睞枥庑赖。

4、MATLAB 操作桌面有几个窗口?如何使某个窗口脱离桌面成为独立窗口?又如何将脱离出去的窗口重新放置到桌面上?聞創沟燴鐺險爱氇谴净。

答:在MATLAB 操作桌面上有五个窗口,在每个窗口的右下角有两个小按钮,一个是关闭窗口的Close 按钮,一个是可以使窗口称为独立的Undock 按钮,点击Undock 按钮就可以使该窗口脱离桌面称为独立窗口,在独立窗口的view 菜单中选择Dock,菜单项就可以将独立的窗口重新防止的桌面上.残骛楼諍锩瀨濟溆塹籟。

5、如何启动M 文件编辑/调试器?答:在操作桌面上选择“建立新文件”或“打开文件”操作时,M 文件编辑/调试器将被启动.在命令窗口中键入edit 命令时也可以启动M 文件编辑/调试器.酽锕极額閉镇桧猪訣锥。

6、存储在工作空间中的数组能编辑吗?如何操作?答:存储在工作空间的数组可以通过数组编辑器进行编辑:在工作空间浏览器中双击要编辑的数组名打开数组编辑器,再选中要修改的数据单元,输入修改内容即可.彈贸摄尔霁毙攬砖卤庑。

7、命令历史窗口除了可以观察前面键入的命令外,还有什么用途?答:命令历史窗口除了用于查询以前键入的命令外,还可以直接执行命令历史窗口中选定的内容、将选定的内容拷贝到剪贴板中、将选定内容直接拷贝到M文件中. 謀荞抟箧飆鐸怼类蒋薔。

Matlab编程与应用习题和一些参考答案

Matlab编程与应用习题和一些参考答案

Matlab 上机实验一、二1.安装Matlab 软件。

2.验证所学内容和教材上的例子。

3.求下列联立方程的解⎪⎪⎩⎪⎪⎨⎧=+-+-=-+=++-=--+41025695842475412743w z y x w z x w z y x w z y x>> a=[3 4 -7 -12;5 -7 4 2;1 0 8 -5;-6 5 -2 10];>> b=[4;4;9;4];>> c=a\bc =5.22264.45701.47181.59944.设⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡------=81272956313841A ,⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-----=793183262345B ,求C1=A*B’;C2=A’*B;C3=A.*B,并求上述所有方阵的逆阵。

>> A=[1 4 8 13;-3 6 -5 -9;2 -7 -12 -8];>> B=[5 4 3 -2;6 -2 3 -8;-1 3 -9 7];>> C1=A*B'C1 =19 -82 3012 27 3-38 54 29>> C2=A'*BC2 =-15 16 -24 3663 -17 93 -10522 6 117 -6019 46 84 -10>> C3=A.*BC3 =5 16 24 -26-18 -12 -15 72-2 -21 108 -56>> inv(C1)ans =0.0062 0.0400 -0.0106-0.0046 0.0169 0.00300.0168 0.0209 0.0150>> inv(C2)Warning: Matrix is close to singular or badly scaled.Results may be inaccurate. RCOND = 8.997019e-019.ans =1.0e+015 *-0.9553 -0.2391 -0.1997 0.27000.9667 0.2420 0.2021 -0.2732-0.4473 -0.1120 -0.0935 0.1264-1.1259 -0.2818 -0.2353 0.3182>> inv(C3)Error using ==> invMatrix must be square.5.设 ⎥⎦⎤⎢⎣⎡++=)1(sin 35.0cos 2x x x y ,把x=0~2π间分为101点,画出以x 为横坐标,y 为纵坐标的曲线。

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

习题3
1. 要求在闭区间]2,0[π上产生具有10个等距采样点的一维数组。

试用两种不同
的指令实现。

〖提示〗用“:”产生10个等距采样点一维数组;用linspace(a,b,n)产生10个
等距采样点一维数组
2. 由指令rand('state',0),A=rand(3,5)生成二维数组A ,试求该数组中所有大于0.5
的元素的位置,分别求出它们的“全下标”和“单下标”。

〖提示〗rand('state',0)---将随机发生器置为0状态
〖答案〗
大于0.5的元素的全下标
行号 1 3 2 3 3 2 3 1 2 列号 1 1 2 2 3 4 4 5 5
大于0.5的元素的单下标
1 3 5 6 9 11 1
2 1
3 14
3. 下列命令执行后,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
4. 已知矩阵⎥⎦


⎣⎡=43
21
A ,运行指令B1=A.^(2), B2=A^(2), 可以观察到不同运算
方法所得结果不同。

(1)请分别写出根据B1, B2恢复原矩阵A 的程序。

(2)
用指令检验所得的两个恢复矩阵是否相等。

5. 在时间区间 [0,10]中,绘制
t
e
y t
2cos 15.0--=曲线。

要求分别采取“标量循环
运算法”和“数组运算法”编写两段程序绘图。

〖答案〗
510
00.5
11.5
t
y 1
510
00.5
1
1.5
t
y 2
6. 先运行clear,format long,rand('state',1),A=rand(3,3),然后根据A 写出两个矩阵:
一个对角阵B ,其相应元素由A 的对角元素构成;另一个矩阵C ,其对角元素全为0,而其余元素与对应的A 阵元素相同。

7. 利用基本矩阵产生3×3和15×8的单位矩阵、全1矩阵、全0矩阵、正态
分布矩阵(均值为1,方差为4) 〖提示〗参见教材P40页例题3.2
8.输入矩阵⎥⎥⎥


⎢⎢⎢⎣⎡=98
7
654321A ,使用全下标方式取出元素“3”,使用单下标方式 取出元素“8”,取出后2行子矩阵块,使用逻辑矩阵方式取出⎥⎦


⎣⎡97
31。

9.输入A 为3×3的魔方阵,B 为3×3单位矩阵,由A 、B 小矩阵组成3×6的
矩阵C 和6×3的大矩阵D ,将D 矩阵的最后一行构成小矩阵E 。

10.将矩阵⎥⎥⎥


⎢⎢⎢
⎣⎡=98
7
654
321A 用flipud 、fliplr 、rot90、diag 、tril 和triu 函数操作。

11.求矩阵⎥⎦


⎣⎡4321的转置矩阵、逆矩阵、矩阵的秩、矩阵的行列式值、矩阵的三
次幂、矩阵的特征值和特征向量。

12.求解方程组⎪⎪⎩


⎨⎧=+-+=++-=++=++-522778638232432143214214321x x x x x x x x x x x x x x x
13.计算数组⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=98
7
654321A ,⎥⎥⎥


⎢⎢
⎢⎣⎡=33
3
222111
B 的左除、右除以及点乘和点除
14.输入a=[1.6 -2,4 5.2 -0.2],分别使用数学函数ceil 、fix 、round 查看各种取整
地运算结果。

15.计算函数)4sin(1110
)(22
t e
z
t f t
--=的值,其中t 的范围为0~2π,步长取
0.1π,z 为0.707;0)()(1>=t f t f 为的部分,计算)(1t f 的值。

16.x 在[0,20]范围内,计算多项式12345234++++=x x x x y 的值,并根据x
和y 进行二阶、三阶和四阶拟合
17.已知⎥⎥
⎥⎥⎦
⎤⎢⎢⎢
⎢⎣⎡---=14.354
54
.96
3205325
654541
0778.01023A 完成下列操作 取A 的前三行构成B 矩阵,前两列构成C 矩阵,右下角3×2子矩阵D ,B
与C 相乘构成矩阵E 18.已知多项式1342)(234-+-+=x x x x x f 和1)(2-=x x g ,求两个多项式的和、差、积、商及f (x )的导函数。

相关文档
最新文档