matlab15练习.docx
- 文档编号:3294213
- 上传时间:2022-11-21
- 格式:DOCX
- 页数:13
- 大小:238.55KB
matlab15练习.docx
《matlab15练习.docx》由会员分享,可在线阅读,更多相关《matlab15练习.docx(13页珍藏版)》请在冰豆网上搜索。
matlab15练习
1.请生成从0到10间隔为0.1的列向量。
A=0:
0.1:
10
2.矩阵A为5×5的魔方矩阵,请从A中提取8和16产生新的列向量B。
a=magic(5)
b=(a[12],[45])
已知A=
请通过矩阵变换,由A生成B=
A=[1234;5678;9123;4567];
B=A(:
[4231])
3.写出产生序列[
…
]的程序。
A=pi/100:
pi/100:
pi
4.下列变量名中,属于合法变量名的是____A_________。
A.flower2B.2flowerC._whatD.who's_it
5.矩阵M=[112733;295712;734537]
则M(2,3)=____A_________。
A.12B.29C.37D.45
6.表达式'
'的执行结果是_____B________。
A.a=38.18B.a=293.18C.a=292.18D.a=37.56
7.某同学设计了一个程序文件myprogram.m,并将其保存到了c:
\mydocument中,但在命令窗口中输入文件名>>myprogram后,MATLAB系统提示:
?
?
?
Undefinedfunctionorvariable‘myprogram’
试分析产生错误的原因并给出解决办法。
没有修改路径;file-setpath-addFolder
8.MATLAB语句后,加上___A______,则运行时不显示中间结果。
A.分号B.逗号C.句号D.空格
9.使用语句t=0:
110生成_____C________个元素的向量。
A.1100B.910C.111D.110
10.已知两个矩阵:
,
将矩阵A中所有等于-1的元素改为-3,将矩阵B中等于0的元素值改为新的矩
阵A中相应位置元素的值,最后将新的矩阵A的第三、四行的元素全改为1,
编写相应的程序。
A=[-3-76-5;-213-1;0013;26-10];
A([24],[34])=[-3]
B=[1468;2370;-2-157;11-10];
A([24],4)=[0]
11.struct函数创建空结构体person,结构体的域包括name、score、rank和class,给name域赋值字符串变量MATLAB,然后删除(参考rmfield命令)class域,编写程序实现。
person=struct('name','','score','','rank','','class','');
person.name='MATLAB'
rmfield(person,'class')
12.已知多项式
和
,编写程序计算:
1)多项式a(x)+b(x)
p1=[1201];
>>p2=[0101];
>>p1+p2
2)多项式a(x)的根。
roots(p1)
13.已知多项式
和
,编写程序计算:
1)多项式a(x)*b(x)
p1=[02035];
>>p2=[10301];
>>p=conv(p1,p2)
2)多项式a(x)/b(x)的部分分式展开。
[r,p,k]=residue(p1,p2)
14.求解下列方程:
a=[34-7-12;5-742;108-5;-65-210]
b=[4;-3;9;-8]
a\b
15.一天中各个时段温度观测结果如下:
表1室内外温度观测结果(℃)
时间/h
5
7
9
11
13
15
17
室内t1
15.0
18.0
20.0
23.0
27.0
28.0
26.0
室外t2
12.0
15.0
22.0
28.0
32.0
33.0
30.0
试用三次样条(spline)插值法分别求出室内外5:
30-16:
30之间每隔2h,各点的近似温度(℃)。
t=5:
2:
17;
>>y1=[15182023272826];
>>y2=[12152228323330];
>>tt=5.5:
2:
16.5;
>>yy1=spline(t,y1,tt)
yy2=spline(t,y2,tt)
16.求解在x=8时多项式(x-1)(x-2)(x-3)(x-4)的值。
p1=[1-1];
p2=[1,-2];
p3=[1-3];
p4=[1-4];
a=conv(p1,p2)
b=conv(a,p3)
c=conv(b,p4)
polyval(c,8)
17.关于稀疏矩阵,下列哪句话是正确的______C____________。
A.用sparse函数可将一个稀疏矩阵转化为一个满矩阵
B.若A为满矩阵,B为稀疏矩阵,那么A+B仍然为一个稀疏矩阵
C.满矩阵能够参与的运算稀疏矩阵同样能够参与
D.稀疏矩阵是指元素值为0的元素个数占矩阵元素总个数的90%以上
18.当x=-3.2,使用取整函数得出-4,则该取整函数名为______D_________。
A.fixB.roundC.ceilD.floor
19.使用命令LOGSPACE(X1,X2)创建X1=0,X2=4
的行向量,有20个元素,查看其元素分布情况。
(1)写出程序代码;logspace(0,4*pi,20)
logspace(0,4*pi,20)
(2)分析结果中“Columns1through12”为什么都是0.0000?
如何显示正常的结果?
没有选择合适的数值显示格式
File-preference-numberformat-long
20.已知函数
,求该函数在区间
的最小值。
x=-1:
0.01:
3;
p=[10-3];
y=polyval(p,x)
min(y)
21.符号表达式“f=sym('sin(a*b)+cos(c*d)')”中的自由符号变量是d___。
22.创建符号表达式
的命令是_______。
symsabtx
a=[2+3,a*b;sin(t),sqrt(x+1)]
23.用两种方法,求解下列方程:
(1)a=[34-7-12;5-742;108-5;-65-210]
b=[4;-3;9;-8]
a\b
(2)eq1=sym('3*x+4*y-7*z-12*w=4');
eq2=sym('5*x-7*y+4*z+2*w=-3');
eq3=sym('x+8*z-5*w=9');
eq4=sym('-6*x+5*y-2*z+10*w=-8');
s=solve(eq1,eq2,eq3,eq4)
24.请用62位数值解的精度表示‘pi’。
vpa(pi,62)
25.求方程组
关于
的解。
eq1=sym('u*y^2+vz+w=0');
eq2=sym('y+z+w=0');
[y,z]=solve(eq1,eq2)
26.利用符号积分法求
,并得到数值解。
27.symsx
f=1/sqrt(2*pi)*exp((-x^2)/2)
I=int(f,0,1)
28.计算定积分
的命令,给出计算结果。
f=(abs(x)+sin(x))*x^2
I=int(f,-1,1)
29.求
的3阶导数的命令。
y=atan((x+3)/(x-2))-log(1+exp(-2*x))
diff(y,3)
30.求符号极限
的命令。
symsxnm
f=tan(m*x)/n*x
limit(f)
31.用符号计算验证三角等式:
sin
(1)cos
(2)-cos
(1)sin
(2)=sin(1-2)
symsphi1phi2
y=sin(phi1)*cos(phi2)-cos(phi1)*sin(phi2)
y1=simple(y)
true=y1==sin(phi1-phi2)
32.在[0,2pi]范围内,用红线画sin(x),用绿圈画cos(x)。
并加标注x轴:
“自变量X”、y轴:
“函数Y”、图名:
“示意图”,并加格栅。
Grid
x=0:
0.1:
2*pi;
plot(x,sin(x),'r',x,cos(x),'g')
xlabel('自变量x')
ylabel('函数Y')
title('示意图')
grid
33.请用程序画出下图
x=0:
0.1:
2*pi;
plot(x,sin(x),'b',x,cos(x),'g')
text(sin(pi),'t=pi')
text(1.5,cos(1.5),'cos(x)')
legend('sin(x)','cos(x)',4)
34.将屏幕分割成四块,并分别画出y=sin(x),z=cos(x),a=sin(x)cos(x),b=sin(x)/cos(x),并加标注,如下图所示。
x=0:
0.1:
2*pi;
y1=sin(x);
subplot(2,2,1)
plot(y1);
title('sin(x)');
y2=cos(x);
subplot(2,2,2)
plot(y2);
title('cos(x)');
y3=sin(x).*cos(x)
subplot(2,2,3)
plot(y3);
title('sin(x)*cos(x)');
y4=sin(x)./cos(x)
subplot(2,2,4)
plot(y4);
title('sin(x)/cos(x)');
35.利用plotyy指令在同一张图上绘制
和
在
上的曲线。
x=0:
0.1:
4;
plotyy(x,sin(x),x,10.^x)
36.用曲面图表示函数
(如下图)。
x=-2:
0.1:
2;
y=x;
[xx,yy]=meshgrid(x,y);
z=xx.^2+yy.^2;
surf(xx,yy,z)
37.请画出如下3D网线图和3D曲面图,其中
,
x=-4:
0.1:
4;
y=x;
[xx,yy]=meshgrid(x,y);
z=9*(1-xx).^2*exp(-xx.^2/2-(yy+1).^2);
surf(xx,yy,z)
38.请用符号函数绘图法绘制函数x=sin(3t)cos(t),y=sin(3t)sin(t)的图形,t的变化范围为[0,2]。
ezplot('sin(3*t)*cos(t)','sin(3*t)*sin(t)',[0,2*pi])
39.在同一图形窗口绘制
其中,y1以黑色实线标识,y2以蓝色虚线标识,
,给横坐标加上标注“x轴”
t=0:
0.1:
3*pi;
y1=2*exp(-t).*sin(2*pi*t);
y2=exp(-2*t).*cos(2*pi*t+pi/3);
plot(t,y1,'k-')
hold
plot(t,y2,'b--')
40.请找出下列程序中出现错误的行____A_______。
A.x=0:
0.2:
10;y=sqrt(x^2+sqrt(x));
B.h=plot(x,y,'o',x,y);
C.set(h,'color','r','linestyle','--');
D.xlabel('时间');ylabel('温度');
41.命令plot(x,y,'m:
*'),其中字符串“m:
*”中,第一个字符“m”____B________。
A.表示曲线为m符号构成的线
B.表示曲线颜色为紫色
C.表示曲线x坐标的最大值为m变量的值
D.表示曲线y坐标的最大值为m变量的值
42.编写实现比较a,b两个值大小功能的函数:
result=ex01(a,b),a,b中的大值赋值给result;若a=b,则输出result“两值相等”并同时屏幕提示“请输入两个大小不同的数值”,要求写出完整的函数体。
functionresult=ex01(a,b)
if(a>b)
reult=a
elseifa
result=b
else
reult='两值相等'
end
43.Fibonacci数列各元素为:
1、1、2、3、5、8、…,满足以下关系:
F1=1,F2=1,Fn=Fn-1+Fn-2,用M函数文件实现:
找到Fibonacci数组中第一个大于10000的元素,即n=?
时,Fn>10000.
%ex01
F1=1;F2=1;
whileFn<10000
Fn=Fn-1+Fn-2
n=n+1
end
44.
45.A=randn(4,3)>0.5,编程实现查找矩阵A的每行中第一个0元素所在的列。
将结果存放着一个列向量中。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- matlab15 练习