刘卫国matlab课后习题答案Word下载.docx
- 文档编号:16337054
- 上传时间:2022-11-23
- 格式:DOCX
- 页数:12
- 大小:19.67KB
刘卫国matlab课后习题答案Word下载.docx
《刘卫国matlab课后习题答案Word下载.docx》由会员分享,可在线阅读,更多相关《刘卫国matlab课后习题答案Word下载.docx(12页珍藏版)》请在冰豆网上搜索。
a=100:
999;
b=find(rem(a,21)==0);
c=length(b)t4-2:
a=input(请输入一个字符串:
s);
b=find(a=aa=z);
a(b)=[];
disp(a);
实验二:
t1:
e=eye(3),r=rand(3,2),o=zeros(2,3),s=diag([1,2]);
a=[e,r;
o,s]disp(a^2=);
disp(a^2);
disp([e,r+rs;
o,s^2]);
b=[e,r+r*s;
o,s^2]t2:
h=hilb(5)p=pascal(5)hh=det(h)hp=det(p)th=cond(h)tp=cond(p)a=abs(th-1);
b=abs(tp-1);
ifab
disp(帕萨卡矩阵p性能更好);
elseifab
disp(希尔伯特矩阵h性能更好);
else
disp(两个矩阵性能相同);
endt3:
a=1:
a=reshape(a,5,5)
disp(行列式的值:
disp(det(a));
disp(矩阵的秩:
disp(rank(a));
disp(矩阵的迹:
disp(trace(a));
disp(矩阵的范数:
disp(norm(a));
t4:
a=[-29,6,18;
20,5,12;
-8,8,5][v,d]=eig(a)t5:
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]x1=a\bb(3)=0.53x2=a\b
disp(系数矩阵a的条件数:
disp(cond(a));
t6:
a=reshape(a,5,5)disp(b1=sqrtm(a));
b1=sqrtm(a)
disp(b2=sqrt(a));
b2=sqrt(a)disp(b1*b1);
b1*b1
disp(b2.*b2);
b2.*b2
实验三:
t1:
x=-5.0:
2:
5.0fori=1:
length(x);
if(x(i)0)(x(i)~=-3)y(i)=x(i)^2+x(i)-6;
elseif
(x(i)=0)(x(i)5)(x(i)~=2)(x(i)~=3);
y(i)=x(i)^2-5*x(i)+6;
y(i)=x(i)^2-x(i)-1;
endendyt2-if:
s=input(pleaseenterthescore:
while(s0||s100)
disp(thescoreisnotreasonable);
end
if(s=90s=100);
elseif(s=80s90);
disp(b);
elseif(s=70s80);
disp(c);
elseif(s=60s70);
disp(d);
disp(e);
t2-switch:
switchfix(s/10)case{9,10}disp(a);
case{8}disp(b);
case{7}disp(c);
case{6}disp(d);
case{0,1,2,3,4,5}disp(e);
t=input(请输入工时time=);
ift120
w=120*84+(t-120)*1.15*84;
elseift60
w=t*84-700;
w=84*t;
disp(应发工资为:
disp(w);
a=10+floor(rand
(1)*89)b=10+floor(rand
(1)*89)
2
s=input(请输入+或-或*或/,s);
while(s~=+s~=-s~=*s~=/)disp(输入的符号错误,请重新输入);
endswitchs
case{+}c=a+b;
case{-}c=a-b;
case{*}c=a*b;
case{/}c=a/b;
endct5:
a=rand(5,6)
n=input(请输入n:
while(n1)
disp(输入的n有误,请重新输入);
endifn5
b=a(5,:
)else
b=a(n,:
)end
实验四:
m=0;
fori=1:
n;
m=m+1/i^2;
pii=sqrt(6*m)t1-2:
m=1./(1:
n).^2;
s=sum(m);
pii=sqrt(6*s)t2:
n=1000;
y=0;
y=y+1/(2*i-1);
ify=3
disp(最大n值为:
i-1
disp(对应的y值为:
y-1/(2*i-1)break;
endendt3:
formatlong
a=input(请输入a:
b=input(请输入b:
x
(1)=a/(1.0+b);
i=1;
x
(2)=a/(x(i)+b);
while(abs(x(i+1)-x(i))0.00001i500)i=i+1;
x(i+1)=a/(x(i)+b);
r1=(-b+sqrt(b^2+4*a))/2;
r2=(-b-sqrt(b^2+4*a))/2;
disp([x,num2str(i+1),=,num2str(x(i+1),8)]);
disp([x,num2str(i),=,num2str(x(i),8)]);
disp([r1=,num2str(r1,8),,r2=,num2str(r2,8)]);
formatshortt4:
f
(1)=1;
f
(2)=0;
f(3)=1;
i=4;
while(i=100)
f(i)=f(i-1)-2*f(i-2)+f(i-3);
i=i+1;
fmax=max(f)fmin=min(f)fsum=sum(f)
fp=length(find(f0))fn=length(find(f0))f0=length(find(f==0))t5:
j=0;
s=0;
fori=3:
50
m=i*(i-1)-1;
ifisprime(m)==1m
3
s=s+m;
j=j+1;
endend
disp([亲密数的对数:
num2str(j)]);
disp([所有亲密素数之和:
num2str(s)]);
实验五:
fn1:
functionf=fn1(n)f=n+10*log(n.^2+5);
fn2:
functionf=fn2(n)f=0;
n
f=f+i*(i+1);
endfx:
functionf=fx(x)
f=1./((x-2).^2+0.1)+1./((x-3).^4+0.01);
x=input(请输入一个复数:
fe=exp(x)flg=log(x)fsin=sin(x)fcos=cos(x)t2:
m1=input(请输入m1:
m2=input(请输入m2:
t=input(请输入theta:
a=[m1*cos(t),-m1,-sin(t),0;
m1*sin(t),0,...cos(t),0;
0,m2,-sin(t),0;
0,0,-cos(t),1]b=[0,m1*9.8,0,m2*9.8]disp([a1,a2,n1,n2]=);
disp([a\b]);
j=1;
n=1;
k=1;
fori=10:
99%挑选出10~99中所有的素数
ifisprime(i)==1m(j)=i;
j=j+1;
fort=1:
length(m);
4
s(n)=10*rem(m(t),10)+fix(m(t)/10);
%挑选出的素数进行位置交换
ifisprime(s(n))==1%再挑选出交换绝对素数
r(k)=m(t);
k=k+1;
endn=n+1;
disp([所有两位绝对素数:
num2str(r)])t5:
y1=fn1(40)/(fn1(30)+fn1(20))%函数调用y2=fn2(40)/(fn2(30)+fn2(20))实验六:
x=linspace(0,2*pi,101);
y=(0.5+3*sin(x)./(1+x.^2)).*cos(x);
plot(x,y);
t2:
x=linspace(0,2*pi,100);
y1=x.^2;
y2=cos(2*x);
y3=y1.*y2;
subplot(2,2,1);
plot(x,y1,r:
x,y2,y-,x,y3,b-.);
subplot(2,2,2);
subplot(2,2,3);
plot(x,y2,y-);
subplot(2,2,4);
plot(x,y3,b-.);
t2-3:
x=linspace(0,2*pi,10);
subplot(3,4,1);
bar(x,y1,r);
axis([0,7,0,40]);
subplot(3,4,2);
stairs(x,y1,y);
subplot(3,4,3);
stem(x,y1,b);
subplot(3,4,4);
fill(x,y1,g);
subplot(3,4,5);
bar(x,y2,r);
axis([-1,7,-1,1]);
subplot(3,4,6);
stairs(x,y2,y);
subplot(3,4,7);
stem(x,y2,b);
subplot(3,4,8);
fill(x,y2,g);
subplot(3,4,9);
bar(x,y3,r);
axis([0,7,-30,50]);
subplot(3,4,10);
stairs(x,y3,y);
subplot(3,4,11);
stem(x,y3,b);
subplot(3,4,12);
fill(x,y3,g);
x=linspace(-5,5,100);
y=[];
forx0=x;
ifx0=0;
y=[y,(x0+sqrt(pi))/exp
(2)];
y=[y,1/2*log(x0+sqrt(1+x0^2))];
subplot(1,2,1);
length(x)ifx(i)=0
y(i)=(x(i)+sqrt(pi))/exp
(2);
y(i)=1/2*log(x(i)+sqrt(1+x(i)^2));
subplot(1,2,2);
theta=0:
0.001:
2*pi;
y=a.*sin(b+n.*theta);
polar(theta,y);
t5:
x=linspace(-5,5,21);
y=linspace(0,10,31);
[x,y]=meshgrid(x,y);
z=cos(x).*cos(y).*exp(-sqrt(x.^2+y.^2)/4);
surf(x,y,z);
contour3(x,y,z);
s=0:
pi/2;
t=0:
3*pi/2;
[s,t]=meshgrid(s,t);
x=cos(s).*cos(t);
y=cos(s).*sin(t);
z=sin(s);
shadinginterp;
实验八:
a=rand(1,30000);
disp(均值为:
mean(a)
disp(标准方差为:
std(a)
disp(最大元素为:
max(a)
disp(最小元素为:
min(a)
k=find(a0.5);
b=length(k);
disp(大于0.5的随机数个数占总数的百分比:
sprintf(%f%%,b/300)t2:
p=45+50*rand(100,5)disp(最高分及序号:
[y,u]=max(p,[],1)
disp(最低分及序号:
[y1,u1]=min(p,[],1)
disp(每门课平均分:
5
【篇二:
matlab程序设计与应用(刘卫国编)课后实验答案】
xt>
1.先求下列表达式的值,然后显示matlab工作空间的使用情况并保存全部变量。
2sin850
(1)z1?
1?
e2
(2)z2?
ln(x?
,其中x?
?
1
21?
2i?
2?
5?
0.45
2.9,3.0e0.3a?
e?
0.3a0.3?
asin(a?
0.3)?
ln,a?
3.0,?
2.9,(3)z3?
22
?
t20?
t?
2,其中t=0:
2.5(4)z4?
t2?
2t?
12?
3?
解:
【篇三:
matlab课后习题】
txt>
1matlab的主要功能包括数值计算和符号计算功能,绘图功能,编程语言以及应用工具箱。
2matlab的优点在于可以在窗口直接输入运算命令,并的到计算结果;
有许多现成的函数可以直接调用;
可以直接绘制图形对结果更直接的显示。
3matlab中help命令在窗口中输入后将会显示当前帮助系统中的所包含的所有项目,即搜索路径中所有的目录项目,所以可以搜索到。
4主要有控制系统工具箱,信号处理工具箱,神经网络工具箱,最优化工具箱,金融工具箱,统计学工具箱,功能性工具箱。
5inv为求方阵的逆矩阵,plot为对所给函数作图,max为取最大值,round为四舍五入到最临近的整数。
第2章matlab矩阵及其运算
2-1合法的变量有globalmy_expabcdab_c_
2-2直接输入法;
利用m文件建立;
利用冒号表达式建立一个向量;
建立大矩阵。
2-3a*b要求a阵的列数必须与b阵的行数相等,a.*b表示对应元素相乘,矩阵运算左除与右除值一般不同,点除也有类似的规律,但对于标量则都相同。
2-4
(1)a=3.5;
b=5;
c=-9.8;
x=(2*pi+(b+c)/(pi+a*b*c)-exp
(2))/tan(b+c)+a
x=3.4054
(2)w=sqrt
(2)*(1+0.34245e-6)
w=1.4142
(3)t=[21-3i;
5-0.65];
z=(exp(2*t)*log(t+sqrt(1+t^2)))/2
z=
1.0e+004*
0.0057-0.0007i0.0049-0.0027i
1.9884-0.3696i1.7706-1.0539i
2-5a=[-1,5,-4;
0,7,8;
3,61,7];
b=[8,3,-1;
2,5,3;
-3,2,0];
a+6*b
ans=
4723-10
123726
-15737
i=eye;
a^2-b+i
-18-21618
23533110
22868526
a*b
141416
-105121
125328180
a.*b
-8154
03524
-91220
b*a
-110-15
722853
3-128
a/b
1.2234-0.92552.9787
-0.94682.3511-0.9574
4.61703.872313.8936
b\a
-0.5106-8.6170-1.1277
0.734017.57451.8085
-0.8830-21.21280.4043
[a,b]
-15-483-1
078253
3617-320
2-6a=[23,10,-0.778,0;
41,-45,65,5;
32,5,0,32;
6,-9.54,54,3.14];
b=a(1:
3,:
)
b=
23.000010.0000-0.77800
41.0000-45.000065.00005.0000
32.00005.0000032.0000
c=a(:
1:
2)
c=
23.000010.0000
41.0000-45.0000
32.00005.0000
6.0000-9.5400
d=a(2:
end,3:
end)
d=
65.00005.0000
032.0000
54.00003.1400
e=b*c
e=
1.0e+003*
0.9141-0.2239
1.20802.7123
1.1330-0.2103
ed
01
00
11
e|d,~e|
~d
10
find(a=10a=25)
2-7a=[34,nan,inf,-inf,-pi,eps,0]
a=
34.0000naninf-inf-3.14160.00000
ans=0
any(a)
ans=1
isnan(a)
ans=0100000isinf(a)
ans=0011000isfinite(a)
ans=100012-8a=[1-123;
51-42;
3052;
111509];
diag(a)
9
triu(a)
1-123
01-42
0052
0009
tril(a)
1000
5100
3050
111509
inv(a)
-0.17580.16410.2016-0.0
227-0.1055-0.1016-0.03910.0664
-0.0508-0.08590.15160.00230.3906-0.0313-0.18130.0281det(a)
ans=1280
rank(a)11
trace(a)
ans=16
norm(a,2)
ans=21.3005
norm(a,1)
ans=20
norm(a,inf)
ans=35
cond(a,1)
ans=14.4531
cond(a,2)
ans=11.1739
cond(a,inf)
ans=22.0938
b=[0.43432;
-8.9421];
diag(b)
ans=0.4300
4.0000
triu(b)
0.430043.00002.000004.000021.0000tril(b)
0.430000-8.90004.00000pinv(b)
0.0022-0.0175
0.0234-0.0017
-0.00350.0405
rank(b)
ans=2
trace(b)
ans=4.4300
norm(b,1)
ans=47
norm(b,2)
ans=43.4271
norm(b,inf)
ans=45.4300
2-9formatrat
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 卫国 matlab 课后 习题 答案
![提示](https://static.bdocx.com/images/bang_tan.gif)