程序设计语言第二版课后题答案

《MATLAB7.x》程序设计语言（第二版）部分课后习题答案（楼顺天，姚若玉，沈俊霞编著）

说明：所有答案均是本人在备考过程中亲自整理的，收录了一部分题目的答案。答案可能存在不足甚至谬误，很多算法也可能不是最优的，仅供参考。本人尽力整理的一点点心血，希望对大家有所帮助。

第二章

8、 a=rand(5,5);

[i,j]=find(a>0.5); for u=1:length(i); b(u)=a(i(u),j(u)); end b'

10、aa=any(a');

i=find(aa==0) a(i,:)=[]; 第四章

1、function flag=isprime(m) %m为素数，flag=1 %m不为素数，flag=0 k=sqrt(m); flag=1; for i=2:k;

if rem(m,i)==0; flag=0; break end end

m=input('input an integer\\n'); flag=isprime(m); if flag==1

disp([num2str(m),'是素数']) else

disp([num2str(m),'不是素数']) end

4、function [x1,x2]=jfc(a,b,c) d=b^2-4*a*c; if d>0;

x1=(-b-sqrt(d))/(2*a); x2=(-b+sqrt(d))/(2*a); elseif d==0;

x1=-b/(2*a);

/ 61 / 6

1 x2=x1; end

a=input('a='); b=input('b='); c=input('c='); [x1,x2]=jfc(a,b,c); if x1~=x2;

disp('原方程存在两个不同的根') disp(['x1=',num2str(x1)]) disp(['x2=',num2str(x2)]); elseif x1==x2;

disp('原方程存在两个相同的根') disp(['x1=x2=',num2str(x1)]) else

disp('原方程的根不存在') end

%MATLAB上机作业，184页习题4，求二次方程的实根 function [s1,s2]=solve(a,b,c) d=b^2-4*a*c; if d>0

s1=(-b-sqrt(d))/(2*a) s1=(-b-sqrt(d))/(2*a)

%disp(['原方程有两个不同的根'num2str(s1)'和'nums2tr(s2)]) elseif d==0

s1=(-b-sqrt(d))/(2*a); s1=(-b-sqrt(d))/(2*a);

%disp(['原方程有两个相同的根'num2str(s1)) else

disp('原方程无实根') end end

5、x=-3:0.01:3;

if x>=-3 & x<-1;

y=(-x.^2-4*x-3)/2; elseif x>=-1 & x<1; y=-x.^2+1; else

y=(-x.^2+4*x-3)/2; end

plot(x,y)

6、%MATLAB上机作业，184页习题6，点不同按钮产生不同分布的数

s=menu('请选择随机数类型','U[-10,10]','U[-5,5]','U[-1,1]','N(0,1)'); switch s

case 1,n=(rand(1)-0.5)*20;

2 / 62 / 6

case 2,n=(rand(1)-0.5)*10; case 3,n=(rand(1)-0.5)*2; case 4,n=randn(1);

otherwise disp('请选择！') end n

7、load a.txt [m,n]=size(a); w=zeros(m,1); for i=1:m;

w(i)=(2*a(i,1)+3*a(i,2)+2*a(i,3)+4*a(i,4)+2.5*a(i,5)+a(i,6))/14.5; end w

8、function k=jj(i)

if i<=10

k=0.1*i; elseif i<=20

k=1+(i-10)*0.05; elseif i<=40

k=1.5+(i-20)*0.02; else

k=1.9+(i-40)*0.01; end

s=0;a=2;b=1; for i=1:15; s=s+a/b; c=b; b=a; a=a+c; end s

9、a=ones(15,1);b=ones(15,1);c=ones(15,1);a(1)=2;c(1)=2;

for i=2:15;

a(i)=b(i-1)+a(i-1); b(i)=a(i-1); c(i)=a(i)./b(i) end sum(c) 第五章

3、t=1:1:10;

t=t';

y=[4.842,4.362,3.754,3.368,3.169,3.083,3.304,3.016,3.012,3.005]';

3 / 63 / 6

x1=[ones(size(t)),exp(-t)]; x2=[ones(size(t)),t.*exp(-t)]; p1=x1\\y; p2=x2\\y;

tt=linspace(1,10,1000);

plot(t,y,'*',tt,p1(1)+p1(2)*exp(-tt),'r-',tt,p2(1)+p2(2)*tt.*exp(-tt),'c-') grid on,legend('给定数据','y1拟合','y2拟合')

5、（1）function dy=cwf(t,y)

dy=[5*y(1)-5*y(2)-6*y(3);3*y(1)-2*y(2)+5*y(3);2*y(1)-y(2)-4*y(3)];

x0=[1,-4,5]; tspan=[13,16];

[t,y]=ode45('cwf',tspan,x0); plot(t,y)

（2）function dy=cwf(t,y)

dy=[y(1)+2*y(2)-3*y(3)+y(4);3*y(1)+y(3)-2*y(4);y(1)-2*y(2)+5*y(4);2*y(1)+3*y(2)+y(4)];

x0=[1,-1,2,1]; tspan=[15,16];

[t,y]=ode45('cwf',tspan,x0); plot(t,y) 7、t=1:10;tt=1:0.01:10;

y=[15,39.5,66,85.5,89,67.5,12,-86.4,-236.9,-448.4]; p1=polyfit(t,y,2); y1=polyval(p1,tt); p2=polyfit(t,y,3); y2=polyval(p2,tt);

plot(t,y,'b*',tt,y1,'r-',tt,y2,'c-');

legend('样本点','二次拟合','三次拟合'),grid on 10&11、pa=[2 3 -4];

pb=[4 -2 5]; pc=[3 -2 5 6]; d1=conv(pa,pb)

[q2,r2]=deconv(pc,pa) [q3,r3]=deconv(pc,pb)

dy1=polyder(pa,pb) [q2,r2]=polyder(pc,pa) [q3,r3]=polyder(pc,pb)

12、x=-5:1:5;xx=-5:0.01:5;

y=10.*exp(-abs(x));

y1=interp1(x,y,xx,'nearest'); y2=interp1(x,y,xx,'linear'); y3=interp1(x,y,xx,'spline');

4 / 64 / 6

y4=interp1(x,y,xx,'cubic');

plot(x,y,'*',xx,y1,'r-',xx,y2,'c-',xx,y3,'k-',xx,y4,'g-')

grid on,legend('样本点','最临近内插','线性内插','三次样条内插','三次曲线内插')

13、a=rand(1,50);

amax=max(a) amin=min(a) ua=mean(a) vara=(std(a))^2

b=randn(1,50); bmax=max(a) bmin=min(a) ub=mean(a) varb=(std(a))^2

14、t=[0,0.2,0.4,0.6,0.8,1,2,5]';

tt=0:0.01:5;

y=[1,1.51,1.88,2.13,2.29,2.4,2.6,-4]'; p1=polyfit(t,y,3); y1=polyval(p1,tt);

x=[ones(size(t)),exp(t)]; p2=x\\y;

y2=p2(1)+p2(2)*exp(tt); plot(t,y,'k*',tt,y1,'r-',tt,y2,'b-')

legend('样本点','三次多项式拟合','指数拟合'),grid on 17、function y=jifen1(x)

y=exp(-2*x);

function y=jifen2(x) y=exp(2*x);

function y=jifen3(x) y=x.^2-3*x+0.5;

z1=quadl('jifen1',0,2) z2=quadl('jifen2',0,2) z3=quadl('jifen3',-1,1) 19、function z=jifen(x,y)

z=exp(-x*y)-2*x*y;

q=dblquad('jifen',0,1,-1,1)

20、function dy=cwf(t,y)

dy=[0.5-y(1);y(1)-4*y(2)];

5 / 65 / 6

x0=[1,-0.5]; tspan=[0,25];

[t,y]=ode45('cwf',tspan,x0); plot(y(1,:),y(2,:));grid on

6 / 66 / 6