matlab基础与应用教程课后答案

8. x=linspace(-3,3,100);y=linspace(-3,3,100); [x y]=meshgrid(x,y); fxy=-5./(1+x.^2+y.^2); i=find(abs(x)<=0.8 & abs(y)<=0.5);
fxy(i)=NaN; surf(x,y,fxy) 9. u=linspace(1,10,100);v=linspace(-pi,pi,100); [u,v]=meshgrid(u,v); x=3.*u.*sin(v);y=2.*u.*cos(v);z=4*u.^2; x=3*u.*sin(v);y=2*u.*cos(v);z=4*u.^2; surf(x,y,z);shading interp; light('position',[1,0,1]); 10.
(1) B=A(1:3,:)
C=A(:,1:2)
D=A(2:4,3:4)
E=B*C
(2)A(find(A>=50 & A<=100))
二、实验题
1.x=-74/180*pi;y=-27/180*pi;
sin(x*x+y*y)/sqrt(tan(abs(x+y)))+pi
2.a=-3.0:0.1:3.0;
stem(t,y) subplot(2,2,4); semilogy(t,y); 5. theta=linspace(-pi,pi,100); ro=5.*cos(theta)+4; polar(theta,ro); fi=linspace(0,2*pi,100);a=1 r=a.*(1+cos(fi)); polar(fi,r);
b=[1,0,0,0,0]'; x1=inv(A)*b; x2=A\b; [L,U]=lu(A); x3=U\(L\b); [Q,R]=qr(a); [Q,R]=qr(A);
x4=R\(Q\b)
R=chol(A);
x5=R\(R'\b)
5.
B=sparse(A);
x1=inv(B)*b;
x2=B\b;
gz=dj*120+1.15*dj*(h-120); elseif h<60
gz=dj*h-700; else
gz=dj*h; end format bank; display([gh,gz]) 3. 循环结构 n=input('input n:'); s=0;
for k=1:n s=s+1/k^2;
t=0:pi/100:2*pi; y=sin(t); comet(t,y) 第 4章 一、思考题
5. （1）A=eye(3); （2）C=100+(200-100)*rand(5,6); （3）D=1+sqrt(0.2)*randn(1,500); （4）E=ones(size(A)); （5）A=A+30*eye(size(A)); （6）B=diag(diag(A)) 二、实验题
nA=norm(A);nb=norm(B); cA=cond(A);cB=cond(B); 3. A=[31,1,0;-4,-1,0;4,-8,-2]； [V,D]=eig(A); 4. A=diag([-1,-1,-1,-1],-1)+diag([-1,-1,-1,-1],1)+diag([2,2,2,2,2])
plot(x,x-x.^3/3/2) （2）plot(x,exp(-x.*x/2)/2/pi) （3）x=-8:0.01:8; plot(x,sqrt((64-x.*x)/2)) （4）t=0:pi/100:6*pi; plot(t.*sin(t),t.*cos(t)) 2. （1） x1=linspace(0,1,100); y1=2*x1-0.5; t=linspace(0,pi,100); x=sin(3*t).*cos(t);y=sin(3*t).*sin(t); plot(x1,y1,'r-',x,y,'b:'); text(0.4,1,'y=2x-0.5'); text(-0.4,0.5,'x=sin(3t)cos(t)'); text(-0.4,0.3,'y=sin(3t)sin(t)'); （2）
x=3*exp(-t); % 求各采样点样本值x
dt=t(2)-t(1); f=1/dt;
% 采样周期 % 采样频率(Hz)
X=fft(x);
% 计算x的快速傅立叶变换X
F=X(1:N/2+1); f=f*(0:N/2)/N;
% F(k)=X(k)(k=1:N/2+1) % 使频率轴f从零开始
plot(f,abs(F),'-*')
6. t=0:pi/20:2*pi; x=exp(-t./20).*cos(t);y=exp(-t./20).*sin(t);z=t; plot3(x,y,z); t=0:0.01:1; x=t;y=t.^2;z=t.^3; plot3(x,y,z); 7. x=-30:0.1:0;y=0:0.1:30; [x,y]=meshgrid(x,y); z=10.*sin(sqrt(x.^2+y.^2))./sqrt(1+x.^2+y.^2); meshc(x,y,z);
第 1章
一、思考题
4. （1） B=A(2:5,1:2:5)
（2）A(7)=[]
（3） A=A+30
（4）size(A)
（5） t(find(t==0))=eps
（6）t=reshape(x,3,4)
（7）abs('matlab')
（8）char(93)
5. A=[97,67,34,10;-78,75,65,5;32,5,-23,-59;0,-12,54,7]
729000
3375
（3）A/B
ans =
-5.8845 5.3549 -0.2028
6.3554 -5.6596 4.3293
-6.6325 6.2078 9.0241
B/A
ans =
30.2855 19.2643 -8.5614
33.4394 21.1547 -9.3974 -0.7443 -0.3938 0.2830 （4）[A,B] ans =
% 绘制振幅-频率图
xlabel('Frequency');
3.
x=linspace(1,10,50);
y=log(x);
f=polyfit(x,y,5);
yy=polyval(f,x);
plot(x,y,'r-',x,yy,'g.')
4. N=64;
% 采样点数
T=5;
% 采样时间终点
t=linspace(0,T,N);
% 给出N个采样时间ti(I=1:N)
end display(sqrt(s*6)) 向量运算
n=input('input n:'); k=1:n; display(sqrt(sum(1./k.^2)*6)) 4. y=0;k=0; while y<3
k=k+1; y=y+1/(2*k-1); end display([k-1,y-1/(2*k-1)]) 5. x0=0;x=1;k=0; a=input('a='); b=input('b='); while abs(x-x0)>=1e-5 && k<500 x0=x; x=a/(b+x0); k=k+1; end display([k,x]); display([(-b+sqrt(b^2+4*a))/2,(-b-sqrt(b^2+4*a))/2]);
subplot(1,2,1); scatter(x1,y1,10)
title('y=2x-0.5'); subplot(1,2,2); scatter(x,y,10) 3.
subplot(1,21); x=1:1:100; y=sin(1./x); plot(x,y) subplot(1,2,2); fplot('sin(1/x)',[1,100]) 4. subplot(2,2,1); bar(t,y); subplot(2,2,2); stairs(t,y); subplot(2,2,3);
exp(-0.3*a).*sin(a+0.3)
3.x=[2,4;-0.45,5];
log(x+sqrt(1+x.*x))/2
4. A=[3,54,2;34,-45,7;87,90,15];B=[1,-2,67;2,8,74;9,3,0];
（1）A*B
ans =
129
432
4197
7
-407
-1052
3 54 2 1 -2 67 34 -45 7 2 8 74 87 90 15 9 3 0 [A([1,3],:);B^2] ans =
3 54 2 87 90 15 600 183 -81 684 282 726 15 6 825 5.a=1+2i;b=3+4i; c=exp(pi*i/6) c= 0.8660 + 0.5000i c+a*b/(a+b) ans = 1.6353 + 1.8462i 第 2章 一、思考题 3.s=0; for n=0:63 s=s+2^n; end disp(s)
1. P=pascal(5);H=hilb(5); Dp=det(P);Dh=det(H);
Kp=cond(P);Kh=cond(H); P 矩阵的性能更好，因为 Kp 较小 2. A=[1,-1,2,3;0,9,3,3;7,-5,0,2;23,6,8,3] B=[3,pi/2,45;32,-76,sqrt(37);5,72,4.5e-4;exp(2),0,97] A1=diag(A);B1=diag(B); A2=triu(A);B2=triu(B); A3=tril(A);B3=tril(B); rA=rank(A);rB=rank(B);
402
591
12489
A.*B
ans =
3 -108 134
68 -360 518
783 270 0
（2） A^3
ans =
-28917
240246
-4368
137883 -259101
27669
171333
252504
38673
A.^3
ans =
27
157464
8
39304
-91125
343
658503
[L,U]=lu(B);
x3=U\(L\b); 第 5章
一、思考题
3.
A=randn(10,5);
mean(A)
std(A)
max(max(A))
min(min(A))
sum(A,2)
sum(sum(A))
sort(A,1)
sort(A,2,'descend')
二、实验题
1. A=rand(1,Байду номын сангаас0000);
6.
y=fun(40)/(fun(30)+fun(20)) （1）函数文件fun.m function f=fun(n) f=n+log(n^2+5); （2）函数文件fun.m function f=fun(n) a=1:n; f=sum(a.*(a+1));
第 3章 一、思考题
4. t=0:0.01:1.5; y=sqrt(3)/2*exp(-4*t).*sin(4*sqrt(3)*t+pi/3); 5. x=-10:0.01:10; y=linspace(-6,6,size(x,2)) z=x.^3+3*x.*y.^2; plot3(x,y,z) 6. x=100:100:400; y=100:100:400; z=[636,697,624,478;698,712,630,478;680,674,598,412;662,626,552,334]; [X,Y]=meshgrid(x,y); mesh(X,Y,z) 二、实验题 1. （1）x=-10:0.01:10;
mean(A)
std(A)
max(A)
min(A)
size(find(A>0.5))/size(A)
2.
h=[466,715,950,1422,1635];
w=[7.04,4.28,3.40,2.52,2.13];
hh=[500,900,1500];
ww=interp1(h,w,hh,'spline')
n=0:63; s=sum(2.^n) 二、实验题 1. x=input('输入一个 四位整 数：'); y=[fix(x/1000),mod(fix(x/100),10),mod(fix(x/10),10),mod(x,10)] z=mod((y+7),10) x=z(3)*1000+z(4)*100+z(1)*10+z(2) 2. gh=input('输入工号'); h=input('输入工时'); dj=84; if h>120