《MATLAB基础与应用教程习题答案 作者 蔡旭晖 刘卫国 蔡立燕 第1-8章答案.docx》由会员分享，可在线阅读，更多相关《MATLAB基础与应用教程习题答案 作者 蔡旭晖 刘卫国 蔡立燕 第1-8章答案.docx（9页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、第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(1) B=A(1:3,:)C=A(:,1:2)D=A(2:4,3:4)E=B*C(2)A(find(A=50 & A120gz=dj*120+1.15*dj*(h-120); elseif h60gz=dj*h-700; elsegz=dj*h;e
2、ndformat bank; display(gh,gz) 3.循环结构n=input(input n:); s=0;for k=1:ns=s+1/k2;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=1e-5 & k500 x0=x;x=a/(b+x0); k=k+1;end display(k,x);display(-b+sqrt(b2+4*a)/2,(-b-sqrt(b2+4*a)/2);6.y=fun(40)/(fun(30)+fun(20)（
3、1） 函数文件fun.m function f=fun(n) f=n+log(n2+5);（2） 函数文件fun.mfunction 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
4、;680,674,598,412;662,626,552,334;X,Y=meshgrid(x,y); mesh(X,Y,z)二、实验题1.（1）x=-10:0.01:10;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=s
5、in(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）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
6、,2,1);bar(t,y);subplot(2,2,2);stairs(t,y);subplot(2,2,3); 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);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
7、);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);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)/size(A) 2. h=466,715,950,1422,1635;w=7.04,4
8、.28,3.40,2.52,2.13; hh=500,900,1500;ww=interp1(h,w,hh,spline)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);x=3*exp(-t);% 求各采样点样本值x dt=t(2)-t(1);f=1/dt; X=fft(x);采样%点数采样时% 间终点给%出N个采样时间ti(I=1:N)采%样周期采样%频率(Hz)计算%x的快速傅立叶变换XF=X(1:N/2
9、+1);% F(k)=X(k)(k=1:N/2+1)f=f*(0:N/2)/N; plot(f,abs(F),-*) xlabel(Frequency);ylabel(|F(k)|)5.使%频率轴f从零开始 绘% 制振幅-频率图（1）p1=1 2 0 0 7;p2=1 -2;p3=1 0 0 5 1;p12=conv(p1,p2);p=p12+zeros(1,size(p12,2)-size(p3,2),p3; roots(p)（2）A=-1,4,3;2,1,5;0,5,6;Px=polyval(p,A) Pxm=polyvalm(p,A) 6.（1） z=fzero(3*x-sin(x)+1
10、,0)（2） 建立函数文件myfun.m function F=myfun(X) x=X(1);y=X(2);F(1)=x*x+y*y-9;F(2)=x+y-1;在 命 令 窗 口 中 输 入 以 下 命 令 ： x=fsolve(myfun,3,0,optimset(Display,off) 第6章一、思考题2.fx=inline(1./(1+x.2);I,n=quad(fx,-100000,100000,1e-3);I,n=quadl(fx,-100000,100000,1e-3); x=-100000:0.01:100000;y=1./(1+x.*x);trapz(x,y);3.(1)
11、fx=inline(-2*y+2*x*x+2*x); t,y=ode23(fx,0,0.5,1)(2) fx=inline(y-exp(x)*cos(x); t,y=ode23(fx,0,3,1)二、实验题1.for x=1:3fx=x,x2,x3;1,2*x,3*x;0,2,6*x; diff(fx)end2.（1）x=0:0.01:1;y=x.10+10.x+1./log10(x); dy=diff(y)/0.01;（2）x=0:0.01:1;y=log(1+x); dy=diff(y,2)/0.01; plot(x(1:99),dy) 3.（1）fx=inline(x.2.*sqrt(2
12、*x.*x+3); quad(fx,1,5)（2） fx=inline(x./sin(x).2); quad(fx,pi/4,pi/3)（3） fx=inline(abs(cos(x+y); dblquad(fx,0,pi,0,pi)（4） syms x y; fx=x*y;int(int(fx,y2,y+2),-1,2) x的积分区间为【0，2】时fx=inline(x.*y);dblquad(fx,0,2,-1,2)4.x=0.3:0.2:1.5; y=0.3895,0.6598,0.9147,1.1611,1.3971,1.6212,1.8325;trapz(x,y)5.（1）yp=in
13、line(-(1.2+sin(10*x)*y);t,y=ode23(yp,0,5,1);（2）令 x1= y, x2= y, x3 = y ，则可写出原方程的状态方程形式：x = x12x = x23x = cos t + 5 cos 2t x - x -1x3(t + 1)23010 x 0001 1 0 x 1 23 + sin t 1x-x = 2 15 cos 2t 2 + x - 1 x cos t 3 3 + sin t建立函数文件ztfun.m function y=ztfun(t,x) b=0;0;cos(t);(t + 1)2 3y=0,1,0;0,0,1;-1/(3+sin
14、(t),-1,5*cos(2*t)/(t+1)2*x+b;解微分方程组t,y=ode23(ztfun,0,5,1;0;2); 6.建立函数文件ztfun.mfunction yy=ztfun(t,y)yy=y(2)*y(3);-y(1)*y(3);-0.51*y(1)*y(2);解微分方程组t,y=ode23(ztfun,0,5,0;1;1) 第7章一、思考题3.（1） 数值积分fx=inline(exp(x).*(1+exp(x).2); quad(fx,0,log(2)符号积分f=sym(exp(x)*(1+exp(x)2); v=int(f,0,log(2);eval(v)（2） 略 二
15、、实验题1.A=sym(1,2,3;x,y,z;3,2,1)rank(A) inv(A)det(A)2.(1) y=sym(sqrt(x+sqrt(x+sqrt(x); y1=diff(y)y2=diff(y,x,2)(2) syms x y; fxy=sin(x2*y)*exp(-x2-y); diff(diff(fxy,x),y)3.(1) syms x int(1/(1+x4)(2) syms x tint(-2*x*x+1)/(2*x*x-3*x+1)2,x,cos(t),exp(2*t)4.syms n x symsum(1/(2*n+1)/(2*x+1)(2*n+1),n,0,in
16、f) symsum(1/(2*n+1)/(2*x+1)(2*n+1),n,0,5) 5.(1) syms x taylor(exp(x)+exp(-x)/2,5,0)(2) syms a xtaylor(exp(-5*x)*sin(3*x+pi/3),5,a)6.(1)x=solve(sym(x3+a*x+1=0)(2)x y=solve(sqrt(x2+y2)-100=0,3*x+5*y-8=0)7. 方程转化为： y = y1y1+ y = 1 - t 2p符号解y1,y11=dsolve(Dy=y1,Dy1+y=1-t2/pi,y(-2)=5,y1(-2)=5,t) 数值解编 写 函 数 文 件 ztfun.m function yy=ztfun(t,y) yy=y(2);1-t2/pi-y(1); 在命令窗口输入以下命令t,y=ode45(ztfun,-2,7,-5;5); t=linspace(-2,7,49)y2=y8.x,y=dsolve(Dx=3*x+4*y,Dy=-4*x+3*y,x(0)=0,y(0)=1)