matlab习题与答案文档格式.docx
- 文档编号:18633343
- 上传时间:2022-12-30
- 格式:DOCX
- 页数:24
- 大小:286.74KB
matlab习题与答案文档格式.docx
《matlab习题与答案文档格式.docx》由会员分享,可在线阅读,更多相关《matlab习题与答案文档格式.docx(24页珍藏版)》请在冰豆网上搜索。
Columns37through42
Columns43through48
Columns49through54
Columns55through60
Column61
0
A=[1234-4;
34787;
3657];
B=[13-1;
203;
3-27];
A+6*B
ans=
1852-10
467105
215349
A-B+eye(3)
1231-3
32884
0671
A*B
684462
309-72596
154-5241
A.*B
121024
680261
9-13049
A^3
3722623382448604
247370149188600766
78688454142118820
A.^3
172839304-64
39304343658503
27274625343
A/B
16.4000-13.60007.6000
35.8000-76.200050.2000
67.0000-134.000068.0000
B\A
109.4000-131.2000322.8000
-53.000085.0000-171.0000
-61.600089.8000-186.2000
[A,B]
1234-413-1
34787203
36573-27
[A([1,3],:
);
B^2]
1234-4
3657
451
11019
20-540
(1)求它们的乘积C;
(2)将矩阵C的右下角3x2子矩阵赋给D;
(3)查看MATLAB工作空间的使用情况。
A=[12345;
678910;
1112131415;
1617181920;
2122232425];
B=[3016;
17-69;
023-4;
970;
41311];
C=A*B
C=
9315077
258335237
423520397
588705557
753890717
D=A(3:
5,4:
end)
D=
1415
1920
2425
who
Yourvariablesare:
ABCD
whos
NameSizeBytesClass
A5x5200doublearray
B5x3120doublearray
C5x3120doublearray
D3x248doublearray
Grandtotalis61elementsusing488bytes
4.完成下列操作:
(1)求[100,999]之间能被21整除的数的个数。
提示;
先利用冒号表达式,再利用find和Length函数。
a=100:
999;
c=rem(a/21,1);
d=find(c==0);
n=length(d)
n=
43
习题二 程序设计及函数和数据的可视化
1、掌握建立和执行M文件的方法;
2、掌握利用if、switch、for、While语句实现的方法;
3、掌握定义和调用MATLAB函数的方法;
4、掌握绘制二维、三维图形的常用函数以及绘制图形的辅助操作。
1.求下列分段函数的值
要求:
(1)用if语句实现,分别输出x=-5.0,-3.0,1.0,2.0,2.5,3.0,5.0时的y值。
x=input('
pleaseinputx:
'
ifx<
0&
x~=-3
y=x^2+x-6
elseif0<
=x<
10,x~=2&
x~=3
y=x^2-5*x+6
else
y=x^2-x-1
end
-5.0
y=
14
-3.0
30
1.0
2
2.0
2.5
-0.2500
3.0
5.0
6
2.输入一个百分制成绩,要求输出成绩等级A,B,C,D,E.其中90~99分为A,80~89
分为B,70~79分为C,60~69分为D,60分以下为E。
分别用if语句和switch语句实现。
functionswitch_case(number);
g=input('
pleaseinputnumber:
'
s'
switch(g)
case('
90~99'
)
disp('
A'
80~89'
B'
70~79'
C'
60~69'
D'
<
60'
E'
otherwise
输入错误!
end
3.根据
,求π的近似值,当n分别取100、1000、10000时,结果是多少?
n=input('
pleaseinputn:
s=0;
form=1:
n
s=s+1.0/(m*m);
t=sqrt(6*s);
t
100
t=
3.1321
1000
3.1406
10000
3.1415
4.一个三位整数各位数字的立方和等于该数本身则称该数为水仙花数。
试输出全部水仙花数。
fora=1:
9;
forb=0:
forc=0:
i=100*a+10*b+c;
ifa^3+b^3+c^3==i;
disp(i);
end
153
370
371
407
5.定义三个函数文件,分别求给定双曲正弦、余弦和正切函数(如下),并在命令文件中调用该函数文件。
并用你定义的函数分别绘制上述三条曲线。
function[imeanstdev]=stat(x)
mean=1/2*(exp(x)-exp(-x));
stdev=sinh(x);
mean=stdev;
function[imeanstdev]=stat(x)
mean=1/2*(exp(x)+exp(-x));
mean=(exp(x)-exp(-x))/(exp(x)+exp(-x));
6.已知y1=x2,y2=cos(2x),y3=y1*y2,完成下列操作:
(1)在同——坐标系下用不同的颜色和线型绘制3条曲线;
x=0:
2*pi;
y1=x.^2;
plot(x,y1,'
r--'
holdon;
y2=cos(2*x);
plot(x,y2,'
g-.'
y3=y1.*y2;
plot(x,y3);
holdoff;
(2)以子图形式绘制3条曲线。
x=0:
plot(x,y1);
plot(x,y2)
(3)分别用条形图、阶梯图、杆图和填充图绘制3条曲线。
stairs(x,y1,'
stairs(x,y2,'
stairs(x,y3);
stem(x,y1,'
stem(x,y2,'
stem(x,y3);
hist(x,y1,'
hist(x,y2,'
hist(x,y3);
7.绘制极坐标曲线P=asin(b+nθ),并分析参数a,b,n对曲线形状的影响.
t=0:
.01:
2*pi;
polar(t,2*sin(4+5*t),'
--r'
)
polar(t,22*sin(44+55*t),'
8.分别用plot和fplot函数绘制函数y=sin(1/x)的曲线,分析两曲线的差别。
alpha=1./x;
y=sin(alpha);
plot(x,y)
1:
50;
y=sin(1./x);
figure
(2)
fplot(@(x)sin(1./x),[0,50])
习题三 Matlab在工程中的应用
一、习题目的
1、掌握生成特殊矩阵和矩阵分析的方法;
2、掌握线性方程组的求解方法;
3、了解统计和分析的方法以及多项式的常用运算;
4、掌握求数值导数、数值积分、常微分方程数值求解、非线性代数方程数组求解的;
5、掌握符号表达式的运算法则、符号矩阵运算、符号函数极限及导数、符号函数定积分和不定积分的方法。
二、习题内容
1.分别用3种不同的数值方法解线性方程组。
A=[65-25;
9-14-1;
342-2;
3-902];
B=[-413111]'
;
x=A\B
x=
0.6667
-1.0000
1.5000
-0.0000
C=inv(A);
x=C*B
[x,y,z,u]=solve('
6*x+5*y-2*z+5*u+4'
9*x-y+4*z-u-13'
3*x+4*y+2*z-2*u-1'
3*x-9*y+2*u-11'
2/3
z=
-1
u=
3/2
2.将20个学生5门功课的成绩存入矩阵P中,进行如下处理:
(1)分别求每门课的最高分、最低分及相应学生序号。
(2)分别求每门课的平均分和标推方差。
(3)5门课总分的最高分、最低分及相应学生序号。
(4)将5门课总分按从大到小顺序存入zcj中,相应学生序号存入xsxh.
p=[1,67,78,65,87,89;
2,98,78,67,89,78;
...
3,99,78,98,67,78;
4,78,98,45,67,89;
5,64,65,78,99,86;
6,70,71,73,87,88;
7,82,73,67,56,89;
8,61,47,87,56,0;
9,46,97,67,87,56;
10,34,68,56,87,88;
11,56,78,78,74,75;
12,56,67,86,76,98;
13,78,44,56,43,56;
14,98,32,56,65,16;
15,67,76,57,87,78;
16,45,56,56,67,87;
17,67,67,56,83,84;
18,87,89,89,90,69;
19,9,78,78,66,87;
20,56,78,78,78,89];
a=p(:
2);
b=p(:
3);
c=p(:
4);
d=p(:
5);
e=p(:
6);
a1=max(a)
a2=min(a)
b1=max(b)
b2=min(b)
c1=max(c)
c2=min(c)
d1=max(d)
d2=min(d)
e1=max(e)
e2=min(e)
mean(p)
std(p)
a1=
99
a2=
9
b1=
98
b2=
32
c1=
c2=
45
d1=
d2=
e1=
e2=
10.500065.900070.900069.650075.550074.0000
5.916122.459116.511214.102214.336125.1438
3.有3个多项式Pl(x)=x4十2x3十4x2十5,P2(x)=x十2,P3(x)=x2十2x十3,试进行下列操作:
(1)求P(x)=P1(x)*P2(x)*P3(x)。
symsxp1p2p3p
p1=x^4+2*x^3+4*x^2+5;
p2=x+2;
p3=x^2+2*x+3;
p=p1*p2*p3
p=
(x^4+2*x^3+4*x^2+5)*(x+2)*(x^2+2*x+3)
(2)求P(x)的根。
x=solve('
(x^4+2*x^3+4*x^2+5)*(x+2)*(x^2+2*x+3)=0'
-2
0.31557922940312053903691057468041-0.9977863722378007037546503988143*i
0.9977863722378007037546503988143*i+0.31557922940312053903691057468041
-1.6836749431657242695009186065958*i-1.3155792294031205390369105746804
1.6836749431657242695009186065958*i-1.3155792294031205390369105746804
1-2^(1/2)*i
1+2^(1/2)*i
4.用数值方法求定积分。
symsx
I2=int(log(1+x)/(1+x^2),0,1)
I2=
(pi*log
(2))/8
5.求微分方程的数值解。
y=dsolve('
x*D2y-5*Dy+y=0'
y(0)=0,Dy(0)=0'
x'
C2*x^3*besselj(6,2*x^(1/2))
6.求代数方程的数值解。
3*x+sin(x)-exp(x)'
0.36042170296032440136932951583028
7.用符号方法术下列积分。
int('
1/(1+x^4+x^8)'
1/6*3^(1/2)*atan(1/3*(2*x-1)*3^(1/2))+1/6*3^(1/2)*atan(1/3*(1+2*x)*3^(1/2))-1/12*3^(1/2)*log(x^2-3^(1/2)*x+1)+1/12*3^(1/2)*log(x^2+3^(1/2)*x+1)
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- matlab 习题 答案