MATLAB基础教程薛山第二版课后习题答案
MATLAB程序设计教程(第二版)课后答案(可编辑修改word版)
MATLAB 第二版课后答案unit3-8 unit3实验指导1、 n=input('请输入一个三位数：');a=fix(n/100);b=fix((n-a*100)/10);c=n-a*100-b*10;d=c*100+b*10+a2（1）n=input('请输入成绩');switch ncase num2cell(90:100)p='A';case num2cell(80:89)p='B';case num2cell(70:79)p='C';case num2cell(60:69)p='D';otherwisep='E';endprice=p（2）n=input('请输入成绩');if n>=90&n<=100p='A';elseif n>=80&n<=89p='B';elseif n>=70&n<=79p='C';elseif n>=60&n<=69p='D';elsep='E';endprice=p（3）tryn;catchprice='erroe'end3n=[1,5,56,4,3,476,45,6,3,76,45,6,4,3,6,4,23,76,908,6];b=n(1);for m=2:20if n(m)>aa=n(m);elseif n(m)<bb=n(m);endendmax=amin=b法2n=[1,5,56,4,3,476,45,6,3,76,45,6,4,3,6,4,23,76,908,6];min=min(n)max=max(n)4b=[-3.0:0.1:3.0];for n=1:61a=b(n);y(n)=(exp(0.3*a)-exp(-0.3*a))/2*sin(a+0.3)+log((0.3+a)/2);endy5y1=0;y2=1;n=input('请输入n 的值：');for i=1:ny1=y1+1/i^2;y2=y2*((4*i*i)/((2*i-1)*(2*i+1)));endy1y26A=[1,1,1,1,1,1;2,2,2,2,2,2;3,3,3,3,3,3;4,4,4,4,4,4;5,5,5,5,5,5;6,6,6,6,6,6]; n=input('请输入n 的值:');if n<=5&n>=0disp(A([n],:));elseif n<0disp(lasterr);else disp(A([6],:));disp(lasterr);end7（1）f=[];f(n)=n+10*log(n^2+5);endy=f(40)/(f(30)+f(20))（2）f=[];a=0;for n=1:40f(n)=a+n*(n+1);a=f(n);endy=f(40)/(f(30)+f(20))8y=0;m=input('输入m 的值:');n=input('输入n 值：');for i=1:ny=y+i^m;endy************************************************************ function s=shi8_1(n,m)s=0;for i=1:ns=s+i^m;end************************************************************ shi8_1(100,1)+shi8_1(50,2)+shi8_1(10,1/2)思考练习2N=[1,2,3,4,5];2.*NN./21./N1./N.^23s=fix(100*rand(1,20)*9/10+10)y=sum(s)/20j=0;for i=1:20if s(i)<y&rem(s(i),2)==0j=j+1;A(j)=s(i);else continue;endendA4y1=0;y2=0;n=input('请输入n 的值：'); for i=1:ny1=y1+-(-1)^i/(2*i-1);y2=y2+1/4^i;endy1y2unit4实验指导1（1）x=-10:0.05:10;y=x-x.^3./6;plot(x,y)（2）x=-10:0.5:10;ezplot('x^2+2*y^2-64',[-8,8]); grid on;2t=-pi:pi/10:pi;y=1./(1+exp(-t));subplot(2,2,1);bar(t,y);title('条形图(t,y)');axis([-pi,pi,0,1]);subplot(2,2,2);stairs(t,y,'b');title('阶梯图(t,y)');axis([-pi,pi,0,1]);subplot(2,2,3);stem(t,y,'k');title('杆图(t,y)');axis([-pi,pi,0,1]);subplot(2,2,4);loglog(t,y,'y');title('对数坐标图(t,y)');3（1）t=0:pi/50:2*pi;r=5.*cos(t)+4;polar(t,r);title('\rho=5*cos\theta+4'); （2）t=-pi/3:pi/50:pi/3;r=5.*((sin(t)).^2)./cos(t); polar(t,r);4（1）t=0:pi/50:2*pi;x=exp(-t./20).*cos(t);y=exp(-t./20).*sin(t);z=t;plot3(x,y,z);grid on;（2）[x,y]=meshgrid(-5:5);z=zeros(11)+5;mesh(x,y,z);shading interp;5[x,y,z]=sphere(20);surf(x,y,z);axis off;shading interp;m=moviein(20);for i=1:20axis([-i,i,-i,i,-i,i])m(:,i)=getframe;endmovie(m,4);思考练习2（1）x=-5:0.1:5;y=(1./(2*pi)).*exp((-(x.^2))/2); plot(x,y);（2）t=-2*pi:0.1:2*pi;x=t.*sin(t);y=t.*cos(t);plot(x,y);grid on;3t=0:pi/1000:pi;x=sin(3.*t).*cos(t);y1=sin(3.*t).*sin(t);y2=2.*x-0.5;plot(x,y1,'k',x,y2);hold on;k=find(abs(y1-y2)<1e-2);x1=x(k);y3=2.*x1-0.5;plot(x1,y3,'rp');4x=-2:0.01:2;y=sin(1./x);subplot(2,1,1);plot(x,y);subplot(2,1,2);fplot('sin(1./x)',[-2,2],1e-4);5（1）i=-4*pi:0.1:10;j=12./sqrt(i);polar(i,j);title('{\rho}=12/sqrt(\theta)')（2）a=-pi/6:0.01:pi/6;b=3.*asin(a).*cos(a)./((sin(a)).^3+(cos(a)).^3); polar(a,b);6（1）[u,v]=meshgrid(-4:0.1:4);x=3.*u.*sin(v);y=2.*u.*cos(v);z=4.*u.^2;subplot(2,1,1);mesh(x,y,z);subplot(2,1,2);surf(x,y,z);（2）[x,y]=meshgrid(-3:0.2:3);z=-5./(1+x.^2+y.^2);subplot(1,2,1);mesh(x,y,z);subplot(1,2,2);surf(x,y,z);unit5实验指导1A=randn(10,5)x=mean(A)y=std(A)Max=max(max(A))Min=min(min(A))Sumhang=sum(A,2)SumA=sum(Sumhang)B=sort(A);C=sort(B,2,'descend');C2（1）a=0:15:90;b=a./180.*pi;s=sin(b)c=0:15:75;d=c./180.*pi;t=tan(d)e=input('请输入想计算的值：'); S=sin(e/180*pi)T=tan(e/180*pi)S1=interp1(a,s,e,'spline')T1=interp1(c,t,e,'spline')P1=polyfit(a,s,5);P2=polyfit(c,t,5);S2=polyval(P1,e)T2=polyval(P2,e)（2）n=[1,9,16,25,36,49,64,81,100]; N=sqrt(n);x=input('ji suan zhi : ');interp1(n,N,x,'cubic')3N=64;T=5;t=linspace(0,T,N);h=exp(-t);dt=t(2)-t(1);f=1/dt;X=fft(t);F=X(1:N/2+1);f=f*(0:N/2)/N;plot(f,abs(F),'-*')4P=[2,-3,0,5,13];Q=[1,5,8];p=polyder(P)q=polyder(P,Q)[a,b]=polyder(P,Q)5P1=[1,2,4,0,5];P2=[0,1,2];P3=[1,2,3];P=P1+conv(P2,P3)X=roots(P)A=[-1,1.2,-1.4;0.75,2,3.5;0,5,2.5];p=polyval(P,A)思考练习4A=rand(1,30000);a=mean(A)b=std(A)Max=max(A)Min=min(A)n=0;for i=1:30000if(A(i)>0.5)n=n+1;endendny=n/300005p=[45,74,54,55,14;78,98,45,74,12;87,98,85,52,65][M,S]=max(p)[N,H]=min(p)junzhi=mean(p,1)fangcha=std(p,1,1)zong=sum(p,2);[Max,wei]=max(zong)[Min,wei]=min(zong)[zcj,xsxh]=sort(zong,'descend')6x=[1:10:101];y=[0,1.0414,1.3222,1.4914,1.6128,1.7076,1.7853,1.8513,1.9085,1.9590,2.0043]; [p,s]=polyfit(x,y,5)a=1:5:101;y1=polyval(p,a);plot(x,y,':o',a,y1,'-*')unit6实验指导1A=[1/2,1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/6];p=[0.95,0.67,0.52]';x=A\pA=[1/2,1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/6]; p=[0.95,0.67,0.53]';x=A\pcond(A)2（1）x1=fzero(@funx1,-1)function fx=funx1(x)fx=x^41+x^3+1;（2）x2=fzero(@funx2,0.5)function fx=funx2(x)fx=x-sin(x)/x;（3）options=optimset('Display','off');x=fsolve(@fun3,[1,1,1]',options)q=fun3(x)function q=fun3(p)x=p(1);y=p(2);z=p(3);q(1)=sin(x)+y^2+log(z)-7;q(2)=3*x+2^y-z^3+1;q(3)=x+y+z-5;3（1）t0=0;tf=5;y0=1;[t,y]=ode23(@fun4,[t0,tf],y0);t'y'function yp=fun4(t,y)yp=-(1.2+sin(10*t))*y;（2）t0=0;tf=5;y0=1;[t,y]=ode23(@fun5,[t0,tf],y0);t'y'function yp=fun5(t,y)yp=cos(t)-y/(1+t^2);4x=fminbnd(@mymin,0,2);-mymin(x)function fx=mymin(x)fx=-(1+x.^2)/(1+x.^4);5options=optimset('Display','off');[x,fval]=fmincon(@fun6,[0,0,0],[],[],a,b,lb,ub)-fvalfunction f=fun6(x)f=-(sqrt(x(1))+(400-x(1))*1.1+(sqrt(x(2))+(400-x(1))*1.1-x(2))*1.1+sqrt(3)+(((400- x(1))*1.1-x(2))*1.1-x(3))*1.1+sqrt(x(x4)));思考练习1（1）A=[2,3,5;3,7,4;1,-7,1];B=[10,3,5]';C1=inv(A)*BC2=A\B[L,U]=lu(A);x=U\(L\B)（2）A=[6,5,-2,5;9,-1,4,-1;3,4,2,-2;3,-9,0,2];B=[-4,13,1,11]';C1=inv(A)*BC2=A\B[L,U]=lu(A);x=U\(L\B)2（1）x1=fzero(@funx1,1.5)function fx=funx1(x)fx=3*x+sin(x)-exp(x);（2）x1=fzero(@funx2,1)function fx=funx2(x)fx=x-1/x+5;（3）options=optimset('Display','off');x=fsolve(@fun3,[3,0]',options)q=fun3(x)function q=fun3(p)x=p(1);y=p(2);q(1)=x^2+y^2-9;q(2)=x+y-1;3（1）t0=0;tf=5;y0=[0,1];[t,y]=ode45(@vdpol,[t0,tf],y0);[t,y]function ydot=vdpol(t,y);ydot(1)=(2-3*y(2)-2*t*y(1))./(1+t^2);ydot(2)=y(1);ydot=ydot';（2）t0=0;tf=5;y0=[1;0;2];[t,y]=ode45(@vdpoll,[t0,tf],y0);[t,y]function ydot=vdpoll(t,y);ydot(1)=cos(t)-y(3)./(3+sin(t))+5*y(1).*cos(2*t)/((t+1).^2)-y(2); ydot(2)=y(1);ydot(3)=y(2);ydot=ydot';4x=fminbnd(@mymin,0,pi);-mymin(x)function fx=mymin(x)fx=-sin(x)-cos(x.^2);5[x,y1]=fminbnd(@mymax,0,1.5);-y1function fx=mymax(x);fx=-(9*x+4*x.^3-12*x.^2);unit7实验指导1（1）format longfx=inline('sin(x)./x');[I,n]=quadl(fx,0,2,1e-10)（2）format longfx=inline('1./((x-0.3).^2+0.01)-1./((x-0.9).^2+0.04)-6');[I,n]=quad(fx,0,1,1e-10)2（1）global ki;ki=0;I=dblquad(@fxy,0,1,0,1)ki（2）f=inline('abs(cos(x+y))','x','y');I=dblquad(f,0,pi,0,pi)3X=0.3:0.2:1.5;F=[0.3895,0.6598,0.9147,1.1611,1.3971,1.6212,1.8325];trapz(X,F)4p=0:pi/5:2*pi;for n=1:3nDX=diff(sin(p),n)end5f=inline('sin(x)./(x+cos(2.*x))');g=inline('(cos(x).*(x+cos(2*x))-sin(x).*(1-2.*sin(2*x)))/(x+cos(2.*x)).^2');x=-pi:0.01:pi;p=polyfit(x,f(x),5);dp=polyder(p);dpx=polyval(dp,x); %求dp 在假设点的函数值dx=diff(f([x,3.01]))/0.01; %直接对f(x)求数值导数gx=g(x); %求函数f 的导函数g 在假设点的导数plot(x,dpx,x,dx,'.',x,gx,'-'); %作图思考练习2format longfx=inline('1./(1+x.^2)');[I,n]=quad(fx,-Inf,Inf,1e-10)[I,n]=quadl(fx,-Inf,Inf,1e-10)x=-100000:100000;y=1./(1+x.^2);trapz(x,y)format short3(1)format longfx=inline('log(1+x)./(1+x.^2)');[I,n]=quad(fx,0,1,1e-10)(2)format longfx=inline('sqrt(cos(t.^2)+4*sin((2*t).^2)+1)'); [I,n]=quad(fx,0,2*pi,1e-10)4f=inline('4.*x.*z.*exp(-z.^2.*y-x.^2)');I=triplequad(f,0,pi,0,pi,0,1)5f=inline('sin(x)');g=inline('cos(x)');x=0:0.01:2*pi;p=polyfit(x,f(x),5);dp=polyder(p);dpx=polyval(dp,x);dx=diff(f([x,2*pi+0.01]))/0.01;gx=g(x);plot(x,dpx,x,dx,'.',x,gx,'-')unit8实验指导1syms x y;s=x^4-y^4;factor(s)factor(5135)2syms x;f=(x-2)/(x^2-4);limit(f,x,2)sym x;f=(sqrt(pi)-sqrt(acos(x)))/sqrt(x+1);limit(f,x,-1,'right')3sym x;f=sin(1/x);diff(f,'x')diff(f,'x',2)sym x;f=(1-cos(2*x))/x;diff(f,'x')diff(f,'x',2)4sym x;f=sqrt(exp(x)+1);int(f,'x')syms x y;f=x/(x+y);int(f,'y')sym x;f=exp(x)*(1+exp(x))^2;int(f,'x',0,log(2))sym x;f=x*log(x);int(f,'x',1,exp(1))5sym x;s=symsum((-1)^(x+1)/x,1,Inf)sym y;z=symsum(y^(2*y-1)/(2*y-1),1,Inf)6sym x;f1=(exp(x)+exp(-x))/2;f2=sqrt(x^3-2*x+1);taylor(f1,x,5,0)taylor(f2,x,6,0)7syms x y a;x=solve('x^3+a*x+1=0','x')x=solve('sin(x)+2*cos(x)-sqrt(x)=0','x')[x y]=solve('log(x/y)=9','exp(x+y)=3','x','y')8syms n;[x,y]=dsolve('x*(D2y)+(1-n)*(Dy)+y=0','y(0)=0','Dy(0)=0','x') 思考练习2syms x B1 B2 a bs1=2*((cos(x))^2)-(sin(x))^2;s2=sin(B1)*cos(B2)-cos(B1)*sin(B2);s3=sqrt((a+sqrt(a^2-b))/2)+sqrt((a-sqrt(a^2-b))/2); s4=(4*x^2+8*x+3)/(2*x+1);h1=simplify(s1)h2=simplify(s2)h3=simplify(s3)h4=simplify(s4)3syms x a;f=abs(x)/x;limit(f,x,0,'left')f=(x+a/x)^x;limit(f,x,inf)4syms x y mf=sqrt(x+sqrt(x+sqrt(x)));m=diff(f,'x')diff(m,'x')syms x y;f=x+y-sqrt(x^2+y^2);z1=diff(f,x)z2=diff(f,y)5syms x y;f=x+y-sqrt(x^2+y^2);z1=diff(f,x)z2=diff(f,y)sym x;f=1/(asin(x)^2*(1-x^2)^(1/2));int(f)6syms xf=1/(1+x);int(f,0,4)sym x;f=x^3*sin(x)^2/(x^6+2*x^4+1);int(f,-1,1)sym x;f=x^3*sin(x)^2/(x^6+2*x^4+1);int(f,-1,1)7syms ns=symsum(1/4^n,1,inf)sym n;s=symsum(((n+1)/n)^(1/2),1,inf)eval(y)8syms xf=tan(x);taylor(f,x,3,0)syms xf=sin(x)^2;taylor(f,x,5,0)9syms xx=solve('log(1+x)-5/(1+sin(x))=2','x')syms x y z[x y z]=solve('4*x^2/(4*x^2+1)=y','4*y^2/(4*y^2+1)=z','4*z^2/(4*z^2+1)=x','x','y','z') 10[x ,y]=dsolve('Dx=3*x+4*y','Dy=5*x-7*y','x(0)=0','y(0)=1','t')。
MATLAB程序设计教程第二版课后详解
MATLAB第二版课后答案unit3-8 unit3实验指导1、n=input('请输入一个三位数：');a=fix(n/100);b=fix((n-a*100)/10);c=n-a*100-b*10;d=c*100+b*10+a2（1）n=input('请输入成绩');switch ncase num2cell(90:100)p='A';case num2cell(80:89)p='B';case num2cell(70:79)p='C';case num2cell(60:69)p='D';otherwisep='E';endprice=p（2）n=input('请输入成绩');if n>=90&n<=100p='A';elseif n>=80&n<=89p='B';elseif n>=70&n<=79p='C';elseif n>=60&n<=69p='D';elsep='E';endprice=p（3）tryn;catchprice='erroe'end3n=[1,5,56,4,3,476,45,6,3,76,45,6,4,3,6,4,23,76,908,6];a=n(1)；b=n(1);for m=2:20if n(m)>aa=n(m);elseif n(m)<bb=n(m);endendmax=amin=b法2n=[1,5,56,4,3,476,45,6,3,76,45,6,4,3,6,4,23,76,908,6];min=min(n)max=max(n)4b=[-3.0:0.1:3.0];for n=1:61a=b(n);y(n)=(exp(0.3*a)-exp(-0.3*a))/2*sin(a+0.3)+log((0.3+a)/2); endy5y1=0;y2=1;n=input('请输入n的值：');for i=1:ny1=y1+1/i^2;y2=y2*((4*i*i)/((2*i-1)*(2*i+1)));endy1y26A=[1,1,1,1,1,1;2,2,2,2,2,2;3,3,3,3,3,3;4,4,4,4,4,4;5,5,5,5,5,5;6,6,6,6,6,6]; n=input('请输入n的值:');if n<=5&n>=0disp(A([n],:));elseif n<0disp(lasterr);else disp(A([6],:));disp(lasterr);end7（1）f=[];for n=1:40f(n)=n+10*log(n^2+5);endy=f(40)/(f(30)+f(20))（2）f=[];a=0;for n=1:40f(n)=a+n*(n+1);a=f(n);endy=f(40)/(f(30)+f(20))8y=0;m=input('输入m的值:');n=input('输入n值：');for i=1:ny=y+i^m;endy************************************************************ function s=shi8_1(n,m)s=0;for i=1:ns=s+i^m;end************************************************************ shi8_1(100,1)+shi8_1(50,2)+shi8_1(10,1/2)思考练习2N=[1,2,3,4,5];2.*NN./21./N1./N.^23s=fix(100*rand(1,20)*9/10+10)y=sum(s)/20j=0;for i=1:20if s(i)<y&rem(s(i),2)==0j=j+1;A(j)=s(i);else continue;endendA4y1=0;y2=0;n=input('请输入n的值：');for i=1:ny1=y1+-(-1)^i/(2*i-1);y2=y2+1/4^i;endy1y2unit4实验指导1（1）x=-10:0.05:10;y=x-x.^3./6;plot(x,y)（2）x=-10:0.5:10;ezplot('x^2+2*y^2-64',[-8,8]); grid on;2t=-pi:pi/10:pi;y=1./(1+exp(-t));subplot(2,2,1);bar(t,y);title('条形图(t,y)');axis([-pi,pi,0,1]);subplot(2,2,2);stairs(t,y,'b');title('阶梯图(t,y)');axis([-pi,pi,0,1]);subplot(2,2,3);stem(t,y,'k');title('杆图(t,y)');axis([-pi,pi,0,1]);subplot(2,2,4);loglog(t,y,'y');title('对数坐标图(t,y)');3（1）t=0:pi/50:2*pi;r=5.*cos(t)+4;polar(t,r);title('\rho=5*cos\theta+4'); （2）t=-pi/3:pi/50:pi/3;r=5.*((sin(t)).^2)./cos(t); polar(t,r);4（1）t=0:pi/50:2*pi;x=exp(-t./20).*cos(t);y=exp(-t./20).*sin(t);z=t;plot3(x,y,z);grid on;（2）[x,y]=meshgrid(-5:5); z=zeros(11)+5;mesh(x,y,z);shading interp;5[x,y,z]=sphere(20);surf(x,y,z);axis off;shading interp;m=moviein(20);for i=1:20axis([-i,i,-i,i,-i,i])m(:,i)=getframe;endmovie(m,4);思考练习2（1）x=-5:0.1:5;y=(1./(2*pi)).*exp((-(x.^2))/2); plot(x,y);（2）t=-2*pi:0.1:2*pi;x=t.*sin(t);y=t.*cos(t);plot(x,y);grid on;3t=0:pi/1000:pi;x=sin(3.*t).*cos(t);y1=sin(3.*t).*sin(t);y2=2.*x-0.5;plot(x,y1,'k',x,y2);hold on;k=find(abs(y1-y2)<1e-2);x1=x(k);y3=2.*x1-0.5;plot(x1,y3,'rp');4x=-2:0.01:2;y=sin(1./x);subplot(2,1,1);plot(x,y);subplot(2,1,2);fplot('sin(1./x)',[-2,2],1e-4);5（1）i=-4*pi:0.1:10;j=12./sqrt(i);polar(i,j);title('{\rho}=12/sqrt(\theta)')（2）a=-pi/6:0.01:pi/6;b=3.*asin(a).*cos(a)./((sin(a)).^3+(cos(a)).^3); polar(a,b);6（1）[u,v]=meshgrid(-4:0.1:4);x=3.*u.*sin(v);y=2.*u.*cos(v);z=4.*u.^2;subplot(2,1,1);mesh(x,y,z);subplot(2,1,2);surf(x,y,z);（2）[x,y]=meshgrid(-3:0.2:3); z=-5./(1+x.^2+y.^2); subplot(1,2,1);mesh(x,y,z);subplot(1,2,2);surf(x,y,z);unit5实验指导1A=randn(10,5)x=mean(A)y=std(A)Max=max(max(A))Min=min(min(A))Sumhang=sum(A,2)SumA=sum(Sumhang)B=sort(A);C=sort(B,2,'descend');C2（1）a=0:15:90;b=a./180.*pi;s=sin(b)c=0:15:75;d=c./180.*pi;t=tan(d)e=input('请输入想计算的值：'); S=sin(e/180*pi)T=tan(e/180*pi)S1=interp1(a,s,e,'spline')T1=interp1(c,t,e,'spline')P1=polyfit(a,s,5);P2=polyfit(c,t,5);S2=polyval(P1,e)T2=polyval(P2,e)（2）n=[1,9,16,25,36,49,64,81,100]; N=sqrt(n);x=input('ji suan zhi : '); interp1(n,N,x,'cubic')3N=64;T=5;t=linspace(0,T,N);h=exp(-t);dt=t(2)-t(1);f=1/dt;X=fft(t);F=X(1:N/2+1);f=f*(0:N/2)/N;plot(f,abs(F),'-*')4P=[2,-3,0,5,13];Q=[1,5,8];p=polyder(P)q=polyder(P,Q)[a,b]=polyder(P,Q)5P1=[1,2,4,0,5];P2=[0,1,2];P3=[1,2,3];P=P1+conv(P2,P3)X=roots(P)A=[-1,1.2,-1.4;0.75,2,3.5;0,5,2.5]; p=polyval(P,A)思考练习4A=rand(1,30000);a=mean(A)b=std(A)Max=max(A)Min=min(A)n=0;for i=1:30000if(A(i)>0.5)n=n+1;endendny=n/300005p=[45,74,54,55,14;78,98,45,74,12;87,98,85,52,65][M,S]=max(p)[N,H]=min(p)junzhi=mean(p,1)fangcha=std(p,1,1)zong=sum(p,2);[Max,wei]=max(zong)[Min,wei]=min(zong)[zcj,xsxh]=sort(zong,'descend')6x=[1:10:101];y=[0,1.0414,1.3222,1.4914,1.6128,1.7076,1.7853,1.8513,1.9085,1.9590,2.0043]; [p,s]=polyfit(x,y,5)a=1:5:101;y1=polyval(p,a);plot(x,y,':o',a,y1,'-*')unit6实验指导1A=[1/2,1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/6];p=[0.95,0.67,0.52]';x=A\pA=[1/2,1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/6]; p=[0.95,0.67,0.53]';x=A\pcond(A)2（1）x1=fzero(@funx1,-1)function fx=funx1(x)fx=x^41+x^3+1;（2）x2=fzero(@funx2,0.5)function fx=funx2(x)fx=x-sin(x)/x;（3）options=optimset('Display','off');x=fsolve(@fun3,[1,1,1]',options)q=fun3(x)function q=fun3(p)x=p(1);y=p(2);z=p(3);q(1)=sin(x)+y^2+log(z)-7;q(2)=3*x+2^y-z^3+1;q(3)=x+y+z-5;3（1）t0=0;tf=5;y0=1;[t,y]=ode23(@fun4,[t0,tf],y0); t'y'function yp=fun4(t,y)yp=-(1.2+sin(10*t))*y; （2）t0=0;tf=5;y0=1;[t,y]=ode23(@fun5,[t0,tf],y0); t'y'function yp=fun5(t,y)yp=cos(t)-y/(1+t^2);4x=fminbnd(@mymin,0,2);-mymin(x)function fx=mymin(x)fx=-(1+x.^2)/(1+x.^4);5options=optimset('Display','off');[x,fval]=fmincon(@fun6,[0,0,0],[],[],a,b,lb,ub)-fvalfunction f=fun6(x)f=-(sqrt(x(1))+(400-x(1))*1.1+(sqrt(x(2))+(400-x(1))*1.1-x(2))*1.1+sqrt(3)+(((40 0-x(1))*1.1-x(2))*1.1-x(3))*1.1+sqrt(x(x4)));思考练习1（1）A=[2,3,5;3,7,4;1,-7,1];B=[10,3,5]';C1=inv(A)*BC2=A\B[L,U]=lu(A);x=U\(L\B)（2）A=[6,5,-2,5;9,-1,4,-1;3,4,2,-2;3,-9,0,2];B=[-4,13,1,11]';C1=inv(A)*BC2=A\B[L,U]=lu(A);x=U\(L\B)2（1）x1=fzero(@funx1,1.5)function fx=funx1(x)fx=3*x+sin(x)-exp(x);（2）x1=fzero(@funx2,1)function fx=funx2(x)fx=x-1/x+5;（3）options=optimset('Display','off'); x=fsolve(@fun3,[3,0]',options)q=fun3(x)function q=fun3(p)x=p(1);y=p(2);q(1)=x^2+y^2-9;q(2)=x+y-1;3（1）t0=0;tf=5;y0=[0,1];[t,y]=ode45(@vdpol,[t0,tf],y0);[t,y]function ydot=vdpol(t,y);ydot(1)=(2-3*y(2)-2*t*y(1))./(1+t^2);ydot(2)=y(1);ydot=ydot';（2）t0=0;tf=5;y0=[1;0;2];[t,y]=ode45(@vdpoll,[t0,tf],y0);[t,y]function ydot=vdpoll(t,y);ydot(1)=cos(t)-y(3)./(3+sin(t))+5*y(1).*cos(2*t)/((t+1).^2)-y(2); ydot(2)=y(1);ydot(3)=y(2);ydot=ydot';4x=fminbnd(@mymin,0,pi);-mymin(x)function fx=mymin(x)fx=-sin(x)-cos(x.^2);5[x,y1]=fminbnd(@mymax,0,1.5);-y1function fx=mymax(x);fx=-(9*x+4*x.^3-12*x.^2);unit7实验指导1（1）format longfx=inline('sin(x)./x');[I,n]=quadl(fx,0,2,1e-10)（2）format longfx=inline('1./((x-0.3).^2+0.01)-1./((x-0.9).^2+0.04)-6'); [I,n]=quad(fx,0,1,1e-10)2（1）global ki;ki=0;I=dblquad(@fxy,0,1,0,1)ki（2）f=inline('abs(cos(x+y))','x','y');I=dblquad(f,0,pi,0,pi)3X=0.3:0.2:1.5;F=[0.3895,0.6598,0.9147,1.1611,1.3971,1.6212,1.8325];trapz(X,F)4p=0:pi/5:2*pi;for n=1:3nDX=diff(sin(p),n)end5f=inline('sin(x)./(x+cos(2.*x))');g=inline('(cos(x).*(x+cos(2*x))-sin(x).*(1-2.*sin(2*x)))/(x+cos(2.*x)).^2'); x=-pi:0.01:pi;p=polyfit(x,f(x),5);dp=polyder(p);dpx=polyval(dp,x); %求dp在假设点的函数值dx=diff(f([x,3.01]))/0.01; %直接对f(x)求数值导数gx=g(x); %求函数f的导函数g在假设点的导数plot(x,dpx,x,dx,'.',x,gx,'-'); %作图思考练习2format longfx=inline('1./(1+x.^2)');[I,n]=quad(fx,-Inf,Inf,1e-10)[I,n]=quadl(fx,-Inf,Inf,1e-10)x=-100000:100000;y=1./(1+x.^2);trapz(x,y)format short3(1)format longfx=inline('log(1+x)./(1+x.^2)');[I,n]=quad(fx,0,1,1e-10)(2)format longfx=inline('sqrt(cos(t.^2)+4*sin((2*t).^2)+1)');[I,n]=quad(fx,0,2*pi,1e-10)4f=inline('4.*x.*z.*exp(-z.^2.*y-x.^2)'); I=triplequad(f,0,pi,0,pi,0,1)5f=inline('sin(x)');g=inline('cos(x)');x=0:0.01:2*pi;p=polyfit(x,f(x),5);dp=polyder(p);dpx=polyval(dp,x);dx=diff(f([x,2*pi+0.01]))/0.01;gx=g(x);plot(x,dpx,x,dx,'.',x,gx,'-')unit8实验指导1syms x y;s=x^4-y^4;factor(s)factor(5135)2syms x;f=(x-2)/(x^2-4);limit(f,x,2)sym x;f=(sqrt(pi)-sqrt(acos(x)))/sqrt(x+1); limit(f,x,-1,'right')3sym x;f=sin(1/x);diff(f,'x')diff(f,'x',2)sym x;f=(1-cos(2*x))/x;diff(f,'x')diff(f,'x',2)sym x;f=sqrt(exp(x)+1);int(f,'x')syms x y;f=x/(x+y);int(f,'y')sym x;f=exp(x)*(1+exp(x))^2;int(f,'x',0,log(2))sym x;f=x*log(x);int(f,'x',1,exp(1))5sym x;s=symsum((-1)^(x+1)/x,1,Inf)sym y;z=symsum(y^(2*y-1)/(2*y-1),1,Inf)sym x;f1=(exp(x)+exp(-x))/2;f2=sqrt(x^3-2*x+1);taylor(f1,x,5,0)taylor(f2,x,6,0)7syms x y a;x=solve('x^3+a*x+1=0','x')x=solve('sin(x)+2*cos(x)-sqrt(x)=0','x')[x y]=solve('log(x/y)=9','exp(x+y)=3','x','y')8syms n;[x,y]=dsolve('x*(D2y)+(1-n)*(Dy)+y=0','y(0)=0','Dy(0)=0','x') 思考练习2syms x B1 B2 a bs1=2*((cos(x))^2)-(sin(x))^2;s2=sin(B1)*cos(B2)-cos(B1)*sin(B2);s3=sqrt((a+sqrt(a^2-b))/2)+sqrt((a-sqrt(a^2-b))/2);s4=(4*x^2+8*x+3)/(2*x+1); h1=simplify(s1)h2=simplify(s2)h3=simplify(s3)h4=simplify(s4)3syms x a;f=abs(x)/x;limit(f,x,0,'left')f=(x+a/x)^x;limit(f,x,inf)4syms x y mf=sqrt(x+sqrt(x+sqrt(x))); m=diff(f,'x')diff(m,'x')syms x y;f=x+y-sqrt(x^2+y^2);z1=diff(f,x)z2=diff(f,y)5syms x y;f=x+y-sqrt(x^2+y^2);z1=diff(f,x)z2=diff(f,y)sym x;f=1/(asin(x)^2*(1-x^2)^(1/2)); int(f)6syms xf=1/(1+x);int(f,0,4)sym x;f=x^3*sin(x)^2/(x^6+2*x^4+1); int(f,-1,1)sym x;f=x^3*sin(x)^2/(x^6+2*x^4+1); int(f,-1,1)7syms ns=symsum(1/4^n,1,inf)sym n;s=symsum(((n+1)/n)^(1/2),1,inf) eval(y)8syms xf=tan(x);taylor(f,x,3,0)syms xf=sin(x)^2;taylor(f,x,5,0)9syms xx=solve('log(1+x)-5/(1+sin(x))=2','x')syms x y z[x y z]=solve('4*x^2/(4*x^2+1)=y','4*y^2/(4*y^2+1)=z','4*z^2/(4*z^2+1)=x','x','y','z')10[x ,y]=dsolve('Dx=3*x+4*y','Dy=5*x-7*y','x(0)=0','y(0)=1','t')。
最新MATLAB基础教程-薛山第二版-课后习题答案
上海市劳动和社会保障局所辖的“促进就业基金”，还专门为大学生创业提供担保，贷款最高上限达到５万元。
调研课题：4ຫໍສະໝຸດ WWW。google。com。cn。大学生政策2004年3月23日
为此，装潢美观，亮丽，富有个性化的店面环境，能引起消费者的注意，从而刺激顾客的消费欲望。这些问题在今后经营中我们将慎重考虑的。
7．编写脚本，计算上面第2题中的表达式。
clear,clc
disp('sin(60)=');
disp(sind(60))
disp('exp(3)=');
disp(exp(3))
disp('cos(3*pi/4)=');
disp(cos(3*pi/4))
8．编写脚本，输出上面第6题中的表达式的值。
clear,clc
a(4) = [];
vec1 = a==b;
vec2 = mod(a,2)==0;
c = sum(vec1);
vec3 = vec1+vec2;
d = vec3.*a;
vec4 = find(a > 5);
e = a(vec4) + 5;
vec5 = find(a < 5);
f = vec5.^2;
a=39;b=58;c=3;d=7;
disp('a>b');disp(a>b)
disp('a<c');disp(a<c)
disp('a>b&&b>c');,disp(a>b&&b>c)
MATLAB基础教程-薛山第二版-课后习题答案
clear,clc
a = [3, 7, 2, 7, 9, 3, 4, 1, 6];
b = [7];
a(4) = [];
vec1 = a==b;
vec2 = mod(a,2)==0;
c = sum(vec1);
vec3 = vec1+vec2;
d = vec3.*a;
vec4 = find(a > 5);
disp('vecN=[]');disp(vecN)
d=vec(2:2:end);
vecG=d(find(d~=2&d~=4));
e = a(vec4) + 5;
vec5 = find(a < 5);
f = vec5.^2;
disp('c=');disp(c)
disp('d=');disp(d)
disp('e=');disp(e)
disp('f=');disp(f)
3. 向量操作时MATLAB的主要部分，使用给出的向量来做下面的练习。
b=((exp(u)+v)^2)/(v^2-u)
c=(sqrt(u-3*v))/(u*v)
4．计算如下表达式：
(1)
(2)
clear,clc
(3-5*i)*(4+2*i)
sin(2-8*i)
5．判断下面语句的运算结果。
(1) 4 < 20
(2) 4 <= 20
(3) 4 == 20
(4) 4 ~= 20
(5) 'b'<'B'
clear,clc
matlab第二版习题答案
matlab第二版习题答案Matlab是一种强大的数学软件工具，被广泛应用于科学计算、数据分析和工程设计等领域。
对于学习和掌握Matlab的人来说，习题是不可或缺的一部分。
本文将为大家提供Matlab第二版习题的答案，帮助读者更好地理解和应用Matlab。
第一章：基本操作1.1 Matlab的启动和退出启动Matlab的方法有多种，可以通过桌面图标、命令行或者启动器来打开Matlab。
退出Matlab可以直接关闭窗口或者使用命令"exit"。
1.2 Matlab的基本语法Matlab的基本语法与其他编程语言相似，包括变量的定义、运算符的使用、条件语句和循环语句等。
例如，定义一个变量x并赋值为5可以使用语句"x = 5;"。
1.3 Matlab的数据类型Matlab支持多种数据类型，包括数值型、字符型和逻辑型等。
数值型可以是整数或者浮点数，字符型用单引号或双引号表示，逻辑型只有两个值true和false。
第二章：向量和矩阵操作2.1 向量的定义和运算向量是一维数组，可以通过一对方括号来定义。
Matlab提供了丰富的向量运算函数，如加法、减法、乘法和除法等。
2.2 矩阵的定义和运算矩阵是二维数组，可以通过方括号和分号来定义。
Matlab提供了矩阵的加法、减法、乘法、转置和求逆等运算。
2.3 矩阵的索引和切片可以使用索引和切片来访问矩阵中的元素。
索引从1开始，可以使用冒号表示全部元素。
切片可以用来选择矩阵的一部分。
第三章：函数和脚本文件3.1 函数的定义和调用函数是一段独立的代码块，可以接受输入参数并返回输出结果。
在Matlab中，函数的定义以关键字"function"开头，调用函数使用函数名和参数。
3.2 脚本文件的编写和运行脚本文件是一系列Matlab语句的集合，可以保存为.m文件。
通过运行脚本文件，可以一次性执行多个语句，提高效率。
第四章：图形绘制和数据可视化4.1 图形绘制函数Matlab提供了丰富的图形绘制函数，可以绘制线图、散点图、柱状图等。
MATLAB基础教程 薛山第二版 课后复习题答案
《MATLAB及应用》实验指导书《MATLAB及应用》实验指导书班级： T1243-7姓名：柏元强学号：总评成绩：汽车工程学院电测与汽车数字应用中心目录实验04051001 MATLAB语言基础 (1)实验04051002 MATLAB科学计算及绘图 (18)实验04051003 MATLAB综合实例编程 (31)实验04051001 MATLAB语言基础1实验目的1)熟悉MATLAB的运行环境2)掌握MATLAB的矩阵和数组的运算3)掌握MATLAB符号表达式的创建4)熟悉符号方程的求解2实验容第二章1.创建double的变量，并进行计算。
(1)a=87，b=190，计算 a+b、a-b、a*b。
clear,clca=double(87);b=double(190);a+b,a-b,a*b(2)创建 uint8 类型的变量，数值与(1)中相同，进行相同的计算。
clear,clca=uint8(87);b=uint8(190);a+b,a-b,a*b2．计算：(1)() sin60o(2) e3(3) 3cos 4⎛⎫π ⎪⎝⎭clear,clc a=sind(60) b=exp(3) c=cos(3*pi/4)3．设2u =，3v =，计算：(1) 4log uvv(2) ()22euv v u +-(3)clear,clc u=2;v=3; a=(4*u*v)/log(v) b=((exp(u)+v)^2)/(v^2-u) c=(sqrt(u-3*v))/(u*v)4．计算如下表达式：(1) ()()3542i i -+(2)()sin 28i -clear,clc (3-5*i)*(4+2*i) sin(2-8*i)5．判断下面语句的运算结果。
(1) 4 < 20 (2) 4 <= 20(3) 4 == 20 (4) 4 ~= 20 (5) 'b'<'B'clear,clc4 < 20 , 4 <= 20,4 == 20,4 ~= 20,'b'<'B'6．设39a =，58b =，3c =，7d =，判断下面表达式的值。
matlab第二版课后习题答案
matlab第二版课后习题答案
《MATLAB第二版课后习题答案》
MATLAB是一种强大的数学软件，被广泛应用于工程、科学和金融等领域。
《MATLAB第二版》是一本经典的教材，为了帮助学生更好地掌握MATLAB的使用，书中提供了大量的课后习题。
下面我们将为大家总结一些MATLAB第二版课后习题的答案，希望能对大家的学习有所帮助。
1. 第一章课后习题答案
第一章主要介绍了MATLAB的基本操作，包括变量的定义、矩阵的运算、函数的使用等。
在课后习题中，有一道题目是要求计算一个矩阵的逆矩阵。
答案是使用MATLAB中的inv函数，将原矩阵作为参数传入即可得到逆矩阵。
2. 第二章课后习题答案
第二章介绍了MATLAB中的绘图功能，包括二维和三维图形的绘制。
有一道课后习题是要求绘制一个正弦曲线和余弦曲线，并在同一张图上显示。
答案是使用MATLAB中的plot函数，分别绘制正弦曲线和余弦曲线，并使用legend函数添加图例。
3. 第三章课后习题答案
第三章介绍了MATLAB中的控制流程，包括if语句、for循环和while循环等。
有一道课后习题是要求编写一个程序，计算1到100之间所有偶数的和。
答案是使用for循环遍历1到100之间的所有数，判断是否为偶数并累加。
通过以上几个例子，我们可以看到MATLAB第二版课后习题的答案涵盖了各种基本和高级的操作，对于学习MATLAB是非常有帮助的。
希望大家在学习MATLAB的过程中能够多加练习，掌握更多的技巧和方法。
MATLAB程序设计教程(第二版)课后答案
MATLAB第二版课后答案unit3-8 unit3实验指导1、n=i nput('请输入一个三位数：');a=fix (n/100);b=fix(( n-a*100)/10);c=n-a*100-b*10;d=c*100+b*10+a2 (1)n=input('请输入成绩');switch ncase num2cell(90:100)p='A';case num2cell(80:89)p='B';case num2cell(70:79)p='C';case num2cell(60:69)p='D';otherwisep='E';end price二p(2)n=input('请输入成绩'); if n>=90&n <=100 p='A：elseif n>=80&n v=89p='B';elseif n>=70&n v=79p='C';elseif n>=60&n v=69p='D';elsep='E';endprice=p(3)tryn;catchprice二'erroe'end3n二[1,5,56,4,3,476,45,6,3,76,45,6,4,3,6,4,23,76,908,6]; a=n(1);b=n (1);for m=2:20if n(m )>aa=n( m);elseif n(m)<bb=n( m);endendmax=amin=b法2n二[1,5,56,4,3,476,45,6,3,76,45,6,4,3,6,4,23,76,908,6]; mi n=mi n(n)max二max (n)4b二[-3.0:0.130];for n=1:61a=b( n);y(n )=(exp(0.3*a)-exp(-0.3*a))/2*si n(a+0.3)+log((0.3+a)/2);5y1=0;y2=1;n=input('请输入n的值：');for i=1: ny1=y1+1/i A2;y2二y2*((4*i*i)/((2*i-1)*(2*i+1)));endyiy26A=[1,1,1,1,1,1;2,2,2,2,2,2;3,3,3,3,3,3;4,4,4,4,4,4;5,5,5,5,5,5;6,6,6,6,6,6]; n=input('请输入n的值:');if n<=5&n>=0disp(A([n],:));elseif n<0disp(lasterr);else disp(A([6],:));disp(lasterr);endf=[]；for n=1:40f(n )=n+10*log( n^2+5);endy=f(40)/(f(30)+f(20))(2)f=[];a=0;for n=1:40f(n )=a+n*( n+1);a=f(n);endy=f(40)/(f(30)+f(20))8y=0;m=input('输入m 的值:');n=input('输入n 值：');for i=1: ny=y+i A m;end************************************************************fun cti on s=shi8_1( n,m)s=0;for i=1: ns=s+i A m;end************************************************************ shi8_1(100,1)+shi8_1(50,2)+shi8_1(10,1/2)思考练习2N=[1,2,3,4,5];2.*NN./21./N1./N.A23s=fix(100*ra nd(1,20)*9/10+10)y=sum(s)/20j=0；for i=1:20if s(i)vy&rem(s(i),2)==0j=j+1；A(j)=s(i)；else con ti nue;endendA4y1=0;y2=0;n=input('请输入n的值：'); for i=1: n y1=y1+-(-1)A i/(2*i-1); y2=y2+1/4A i; endyiy2unit4实验指导1 (1)x=-10:0.05:10;y=x-x.A3./6;Plot(x,y)(2)x=-10:0.5:10;ezplot('x A2+2*y A2-64',[-8,8]);grid on;2t=-pi:pi/10:pi;y=1./(1+exp(-t)); subplot(2,2,1);bar(t,y); title('条形图(t,y)');axis([-pi,pi,0,1]); subplot(2,2,2);stairs(t,y,'b');title('阶梯图(t,y)');axis([-pi,pi,0,1]); subplot(2,2,3);stem(t,y,'k');title('杆图(t,y)');axis([-pi,pi,0,1]); subplot(2,2,4);loglog(t,y,'y');title('对数坐标图(t,y)');3 (1)t=0:pi/50:2*pi;r=5.*cos(t)+4;polar(t,r);title('\rho=5*cos\theta+4');(2)t=-pi/3:pi/50:pi/3;r=5.*((si n( t)).A2)./cos(t);polar(t,r);4 (1)t=0:pi/50:2*pi;x=exp(-t./20).*cos(t);y=exp(-t./20).*si n( t);z=t;plot3(x,y,z);grid on;(2)[x,y]=meshgrid(-5:5); z=zeros(11)+5; mesh(x,y,z); shad ing in terp;5[x,y,z]=sphere(20);surf(x,y,z);axis off;shad ing in terp;m=movie in( 20);axis([-i,i,-i,i,-i,i])m(:,i)=getframe;endmovie(m,4);思考练习2 (1)x=-5:0.1:5;y=(1./(2*pi)).*exp((-(x“2))/2); plot(x,y);(2)t=-2*pi:0.1:2*pi;x=t.*si n( t);y=t.*cos(t);grid on;3t=O:pi/1OOO:pi;x=si n( 3.*t).*cos(t);y1=si n(3.*t).*si n( t);y2=2.*x-0.5;plot(x,y1,'k',x,y2);hold on;k二fin d(abs(y1-y2)<1e-2); x1=x(k);y3=2.*x1-0.5;plot(x1,y3,'rp');4x=-2:0.01:2;y=s in (1./x);subplot(2,1,1);plot(x,y);subplot(2,1,2);fplot('si n(1./x)',[-2,2],1e-4);5 (1)i=-4*pi:0.1:10;j=12./sqrt(i);polar(ij);title('{\rho}=12/sqrt(\theta)')(2)a=-pi/6:0.01:pi/6;b=3.*asi n( a).*cos (a) ./((s in (a)).八3+(cos(a)).八3); polar(a,b);6 (1)[u,v]=meshgrid(-4:0.1:4);x=3.*u.*s in (v);y=2.*u.*cos(v);z=4.*u.A2;subplot(2,1,1);mesh(x,y,z);subplot(2,1,2);surf(x,y,z);(2)[x,y]=meshgrid(-3:0.2:3);z=-5./(1+x.A2+y.A2);subplot(1,2,1);mesh(x,y,z);subplot(1,2,2); surf(x,y,z);unit5实验指导1A=ra ndn (10,5)x=mea n(A)y=std(A)Max二max(max(A))Min二min( mi n( A))Sumha ng二sum(A,2) SumA=sum(Sumha ng) B=sort(A); C=sort(B,2,'desce nd'); C2 (1)a=0:15:90;b=a./180.*pi;s=s in(b)c=0:15:75;d二c./180.*pi;t=ta n(d)e=i nput('请输入想计算的值：');S=si n( e/180*pi)T=ta n(e/180*pi)S1=i nterp1(a,s,e,'spli ne')T1二in terp1(c,t,e,'spli ne')P1=polyfit(a,s,5);P2=polyfit(c,t,5);S2=polyval(P1,e)T2=polyval(P2,e)(2)n二[1,9,16,25,36,49,64,81,100]; N=sqrt (n);x=in put('ji sua n zhi :');in terp1( n,N,x,'cubic')3N=64;T=5;t=li nspace(0,T,N);h=exp(-t);dt=t( 2)-t(1);f=1/dt;X=fft(t);F=X(1:N/2+1);f=f*(0:N/2)/N;plot(f,abs(F),'-*')4P=[2,-3,0,5,13];Q=[1,5,8];p=polyder(P)q=polyder(P,Q)[a,b]=polyder(P,Q)5P仁[1,2,4,0,5];P2=[0,1,2];P3=[1,2,3];P=P1+co nv(P2,P3)X=roots(P)A二[-1,1.2,-1.4；0.75,2,3.5;0,5,2.5]; p=polyval(P,A)思考练习A=ra nd(1,30000);a=mea n(A)b二std(A)Max二max(A)Min 二mi n(A)n=0;for i=1:30000if(A(i)>0.5)n=n+1;endendny=n/300005p二[45,74,54,55,14;78,98,45,74,12;87,98,85,52,65] [M,S]=max(p)[N,H]=mi n(p)jun zhi=mea n( p,1)fan gcha=std(p,1,1)zon g=sum(p,2);[Max,wei]=max(z ong)[Mi n, wei]=mi n(zong)[zcj,xsxh]=sort(z on g,'desce nd')6x=[1:10:101];y二[0,1.0414,1.3222,1.4914,1.6128,1.7076,1.7853,1.8513,1.9085,1.9590,2.0043]; [p,s]=polyfit(x,y,5)a=1:5:101;y1=polyval(p,a);plot(x,y,':o',a,y1,'-*')unit6实验指导1A二[1/2,1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/6];p二[0.95,0.67,0.52]';x=A\pA=[1/2,1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/6];p二[0.95,0.67,0.53]';x=A\pcon d(A)2 (1)x1二fzero(@fu nx1,-1) fun cti on fx二fun x1(x) fx=x A41+x A3+1;(2)x2=fzero(@fu nx2,0.5)fun cti on fx二fun x2(x)fx=x-si n( x)/x;(3)optio ns=optimset('Display','off); x=fsolve(@fu n3,[1,1,1]',optio ns) q=fu n3(x) fun cti on q二fun 3(p)x=p(1);y=p (2);z=p (3)；q(1)=si n(x)+yA2+log (z)-7;q(2)=3*x+2Ay-zA3+1; q(3)=x+y+z-5;3 (1)t0=0;tf=5;yO=1;[t,y]=ode23(@fu n4 ,[tO,tf],yO); t'y'fun cti on yp二fu n4 (t,y)yp=-(1.2+si n(1O*t))*y;(2)t0=0;tf=5;y0=1;[t,y]=ode23(@fu n5,[t0,tf],y0); t'y'fun cti on yp二fu n5(t,y)yp=cos(t)-y/(1+t A2);4x=fm inbn d(@m ymin, 0,2);-mymi n(x)fun cti on fx=mymi n(x)fx=-(1+x.A2)/(1+x.A4);5optio ns=optimset('Display','off);[x,fval]=fmi neon (@fu n6,[0,0,0],[],[],a,b,lb,ub)-fvalfun eti on f二fun 6(x)f=-(sqrt(x(1))+(400-x(1))*1.1+(sqrt(x (2))+(400-x(1))*1.1-x (2))*1.1+sqrt(3)+(((400- x(1))*1.1-x(2))*1.1-x(3))*1.1+sqrt(x(x4)));思考练习1 (1)A=[2,3,5;3,7,4;1,-7,1];B=[10,3,5]';C1=i nv(A)*BC2=A\B[L,U]=lu(A);x=U\(L\B)(2)A=[6,5,-2,5;9,-1,4,-1;3,4,2,-2;3,-9,0,2]; B二[-4,13,1,11]';C1=i nv(A)*BC2=A\B[L,U]=lu(A);x=U\(L\B)2 (1)x1二fzero(@fu nx1,1.5)fun cti on fx二fun x1(x)fx=3*x+si n(x)-exp(x);(2)x1=fzero(@fu nx2,1)fun cti on fx二fun x2(x) fx=x-1/x+5;(3)optio ns=optimset('Display','off); x=fsolve(@fu n3,[3,0]',optio ns) q=fu n3(x) fun cti on q二fun 3(p)x=p(1);y=p (2)；q(1)=x A2+y A2-9;q(2)=x+y-1;3 (1)t0=0;tf=5;yo=[o,i];[t,y]=ode45(@vdpol,[t0,tf],y0);[t,y]fun cti on ydot二vdpol(t,y);ydot(1)=(2-3* y(2)-2*t*y(1))./(1+t A2);ydot(2)=y(1);ydot=ydot';(2)t0=0;tf=5;y0=[i;0;2];[t,y]=ode45(@vdpoll,[t0,tf],y0);[t,y]fun cti on ydot二vdpoll(t,y);ydot(1)=cos(t)-y(3)./(3+si n( t))+5*y(1).*cos(2*t)/((t+1)A2)-y(2); ydot(2)=y(1);ydot(3)=y(2);ydot=ydot';4x=fmi nbn d(@mymi n,0,pi);-mymi n(x)fun cti on fx=mymi n(x)fx=-si n(x)-cos(x.A2);5[x,y1]=fmi nbn d(@mymax,0,1.5);-yifun cti on fx=mymax(x);fx=-(9*x+4*x.A3-12*x.A2);unit7实验指导1 (1)format longfx二in li ne('s in( x)./x');[I, n ]=quadl(fx,0,2,1e-10)(2)format longfx=i nlin e('1./((x-0.3)A2+0.01)-1./((x-0.9)A2+0.04)-6'); [l,n ]=quad(fx,0,1,1e-10)2 (1)global ki;ki=0;l=dblquad(@fxy,0,1,0,1)ki(2)f=in li ne('abs(cos(x+y))','x','y');匸dblquad(f,O,pi,O,pi)3X=0.3:0.2:1.5;F二[0.3895,0.6598,0.9147,1.1611,1.3971,1.6212,1.8325]; trapz(X,F)4p=0:pi/5:2*pi;for n=1:3DX=diff(s in (p), n)end5f=inlin e('si n(x)./(x+cos(2.*x))');g=i nlin e('(cos(x).*(x+cos(2*x))-si n(x).*(1-2.*si n(2*x)))/(x+cos(2.*x)).八2'); x=-pi:0.01:pi;p=polyfit(x,f(x),5);dp=polyder(p);dpx=polyval(dp,x); %求dp在假设点的函数值dx=diff(f([x,3.01]))/0.01; %直接对f(x)求数值导数gx=g(x); %求函数f的导函数g在假设点的导数plot(x,dpx,x,dx,'.',x,gx,'-'); %作图思考练习2format longfx=i nlin e('1./(1+x.A2)');[I, n]二quad(fx,-l nf,ln f,1e-10)[I, n]=quadl(fx,-I nf,ln f,1e-10) x=-100000:100000;y=1./(1+x.A2);trapz(x,y)format short3(1)format longfx=i nli ne('log(1+x).心+x.A2)');[l,n ]=quad(fx,0,1,1e-10)⑵format longfx=i nli ne('sqrt(cos(t.A2)+4*s in ((2吐).八2)+1)'); [I, n]=quad(fx,0,2*pi,1e-10)4f=in li ne('4.*x.*z.*exp(-z.A2.*y-x.A2)');匸triplequad(f,0,pi,0,pi,0,1)5f=in li ne('s in( x)');g=i nli ne('cos(x)');x=0:0.01:2*pi;p=polyfit(x,f(x),5);dp二polyder(p);dpx=polyval(dp,x);dx二diff(f([x,2*pi+0.01]))/0.01;plot(x,dpx,x,dx,'.',x,gx,'-') unit8实验指导1syms x y;s=x A4-y A4;factor(s)factor(5135)2syms x;f=(x-2)/(xA2-4);Iimit(f,x,2)sym x;f=(sqrt(pi)-sqrt(acos(x)))/sqrt(x+1); limit(f,x,-1,'right')3sym x;f=si n(1/x);diff(f,'x')diff(f,'x',2)f=(1-cos(2*x))/x;diff(f,'x')diff(f,'x',2)4sym x;f=sqrt(exp(x)+1);in t(f,'x')syms x y;f=x/(x+y);in t(f,'y')sym x;f=exp(x)*(1+exp(x))八2;in t(f,'x',0,log(2))sym x;f=x*log(x);in t(f,'x',1,exp(1))5sym x;s=symsum((-1)八(x+1)/x,1,l nf) sym y;z=symsum(y A(2*y-1)/(2*y-1),1,I nf)6sym x;f1=(exp(x)+exp(-x))/2;f2=sqrt(xA3-2*x+1);(f1,x,5,0)(f2,x,6,0)7syms x y a;x=solve('xA3+a*x+1=0','x')x=solve('si n( x)+2*cos(x)-sqrt(x)=0','x')[x y]=solve('log(x/y)=9','exp(x+y)=3','x','y')syms n;[x,y]=dsolve('x*(D2y)+(1-n)*(Dy)+y=O','y(O)=O','Dy(O)=O','x') 思考练习2syms x B1 B2 a bs1=2*((cos(x))A2)-(si n(x))八2;s2=si n(B1)*cos(B2)-cos(B1)*si n(B2);s3二sqrt((a+sqrt(aA2-b))/2)+sqrt((a-sqrt(aA2-b))/2);s4=(4*xA2+8*x+3)/(2*x+1);h1=simplify(s1)h2=simplify(s2)h3=simplify(s3)h4=simplify(s4)3syms x a;f=abs(x)/x;limit(f,x,0,'left')f=(x+a/x)Ax;limit(f,x, inf)4syms x y m f=sqrt(x+sqrt(x+sqrt(x))); m=diff(f,'x') diff(m,'x') syms x y; f=x+y-sqrt(x A2+y A2);z1=diff(f,x)z2=diff(f,y)5syms x y; f=x+y-sqrt(xA2+yA2);z1=diff(f,x)z2=diff(f,y)sym x;f=1/(asi n(xF2*(1-xA2)A(1/2)); int(f)6syms xf=1/(1+x);in t(f,0,4)sym x;f=x A3*si n(xF2/(xA6+2*x%+1);in t(f,-1,1)sym x;f=xA3*si n(xF2/(xA6+2*xA4+1); in t(f,-1,1)7syms ns=symsum(1/4A n,1,i nf)sym n;s=symsum((( n+1)/n)A(1/2),1,i nf) eval(y)8syms xf=ta n( x);(f,x,3,0)syms xf=si n( x)A2;(f,x,5,0)9syms xx=solve('log(1+x)-5/(1+si n(x))=2','x')syms x y z[x yz]=solve('4*x A2/(4*x A2+1)=y','4*y A2/(4*y A2+1)=z','4*z A2/(4*z A2+1)=x','x','y','z')10[x ,y]=dsolve('Dx=3*x+4*y','Dy=5*x-7*y','x(0)=0','y(0)=1','t')。
MATLAB基础教程 薛山 课后答案[2-12章].khda
![MATLAB基础教程 薛山 课后答案[2-12章].khda](https://img.taocdn.com/s3/m/d96c4bde33d4b14e852468f5.png)
2. (1) y e x (2) y 1
x5
பைடு நூலகம்
x 8
x sin x
>> y=@(x)exp(x)-x^5; >> x = fzero(y,8) x= 12.7132 >> fplot(y,[x-1,x+1]); >> hold on >> plot(x,y(x),'r*');
2 >> y=@(x)x*sin(x); x= 0
5 1
% find the times and places of the specified letter in the string letter = 'a'; string = 'China'; % The specified letter to be searched for % The specified tring to be searched
a 39
b 58
c 3
d
7
(1) a b (2) a c (3) a b & &b c (4) a d (5) a | b c (6) ~~ d 1 >> a=39; >> b=58;
>> c=3; >> d=7; >> a>b ans = 0 2 >> a<c ans = 0 3 >> a>b&&b>c ans = 0 4 >> a==d ans = 0 5 >> a|b>c ans = 1 6 >> ~~d ans = 1 7
MATLAB程序设计教程(第二版)课后答案
h=exp(-t);
dt=t(2)-t(1);
f=1/dt;
X=fft(t);
F=X(1:N/2+1);
f=f*(0:N/2)/N;
plot(f,abs(F),'-*')
4
P=[2,-3,0,5,13];
Q=[1,5,8];
p=polyder(P)
q=polyder(P,Q)
end
y1
y2
6
A=[1,1,1,1,1,1;2,2,2,2,2,2;3,3,3,3,3,3;4,4,4,4,4,4;5,5,5,5,5,5;6,6,6,6,6,6];
n=input('请输入n的值:');
if n<=5&n>=0
disp(A([n],:));
elseif n<0
disp(lasterr);
3
t=0:pi/1000:pi;
x=sin(3.*t).*cos(t);
y1=sin(3.*t).*sin(t);
y2=2.*x-0.5;
plot(x,y1,'k',x,y2);
hold on;
k=find(abs(y1-y2)<1e-2);
x1=x(k);
y3=2.*x1-0.5;
plot(x1,y3,'rp');
4
b=[-3.0:0.1:3.0];
for n=1:61
a=b(n);
y(n)=(exp(0.3*a)-exp(-0.3*a))/2*sin(a+0.3)+log((0.3+a)/2);
end
MATLAB程序设计教程(第二版)课后答案
MATLAB第二版课后答案unit3-8 unit3实验指导1、 n=input('请输入一个三位数：');a=fix(n/100);b=fix((n-a*100)/10);c=n-a*100-b*10;d=c*100+b*10+a2（1）n=input('请输入成绩');switch ncase num2cell(90:100)p='A';case num2cell(80:89)p='B';case num2cell(70:79)p='C';case num2cell(60:69)p='D';otherwisep='E';endprice=p（2）n=input('请输入成绩');if n>=90&n<=100p='A';elseif n>=80&n<=89p='B';elseif n>=70&n<=79p='C';elseif n>=60&n<=69p='D';elsep='E';endprice=p（3）tryn;catchprice='erroe'end3n=[1,5,56,4,3,476,45,6,3,76,45,6,4,3,6,4,23,76,908,6];b=n(1);for m=2:20if n(m)>aa=n(m);elseif n(m)<bb=n(m);endendmax=amin=b法2n=[1,5,56,4,3,476,45,6,3,76,45,6,4,3,6,4,23,76,908,6];min=min(n)max=max(n)4b=[-3.0:0.1:3.0];for n=1:61a=b(n);y(n)=(exp(0.3*a)-exp(-0.3*a))/2*sin(a+0.3)+log((0.3+a)/2);endy5y1=0;y2=1;n=input('请输入n的值：');for i=1:ny1=y1+1/i^2;y2=y2*((4*i*i)/((2*i-1)*(2*i+1)));endy1y26A=[1,1,1,1,1,1;2,2,2,2,2,2;3,3,3,3,3,3;4,4,4,4,4,4;5,5,5,5,5,5;6,6,6,6,6,6]; n=input('请输入n的值:');if n<=5&n>=0disp(A([n],:));elseif n<0disp(lasterr);else disp(A([6],:));disp(lasterr);end7（1）f=[];f(n)=n+10*log(n^2+5);endy=f(40)/(f(30)+f(20))（2）f=[];a=0;for n=1:40f(n)=a+n*(n+1);a=f(n);endy=f(40)/(f(30)+f(20))8y=0;m=input('输入m的值:');n=input('输入n值：');for i=1:ny=y+i^m;endy************************************************************ function s=shi8_1(n,m)s=0;for i=1:ns=s+i^m;end************************************************************ shi8_1(100,1)+shi8_1(50,2)+shi8_1(10,1/2)思考练习2N=[1,2,3,4,5];2.*NN./21./N1./N.^23s=fix(100*rand(1,20)*9/10+10)y=sum(s)/20j=0;for i=1:20if s(i)<y&rem(s(i),2)==0j=j+1;A(j)=s(i);else continue;endendA4y1=0;y2=0;n=input('请输入n的值：'); for i=1:ny1=y1+-(-1)^i/(2*i-1);y2=y2+1/4^i;endy1y2unit4实验指导1（1）x=-10:0.05:10;y=x-x.^3./6;plot(x,y)（2）x=-10:0.5:10;ezplot('x^2+2*y^2-64',[-8,8]); grid on;2t=-pi:pi/10:pi;y=1./(1+exp(-t));subplot(2,2,1);bar(t,y);title('条形图(t,y)');axis([-pi,pi,0,1]);subplot(2,2,2);stairs(t,y,'b');title('阶梯图(t,y)');axis([-pi,pi,0,1]);subplot(2,2,3);stem(t,y,'k');title('杆图(t,y)');axis([-pi,pi,0,1]);subplot(2,2,4);loglog(t,y,'y');title('对数坐标图(t,y)');3（1）t=0:pi/50:2*pi;r=5.*cos(t)+4;polar(t,r);title('\rho=5*cos\theta+4');（2）t=-pi/3:pi/50:pi/3;r=5.*((sin(t)).^2)./cos(t); polar(t,r);4（1）t=0:pi/50:2*pi;x=exp(-t./20).*cos(t);y=exp(-t./20).*sin(t);z=t;plot3(x,y,z);grid on;（2）[x,y]=meshgrid(-5:5);z=zeros(11)+5;mesh(x,y,z);shading interp;5[x,y,z]=sphere(20);surf(x,y,z);axis off;shading interp;m=moviein(20);for i=1:20axis([-i,i,-i,i,-i,i])m(:,i)=getframe;endmovie(m,4);思考练习2（1）x=-5:0.1:5;y=(1./(2*pi)).*exp((-(x.^2))/2); plot(x,y);（2）t=-2*pi:0.1:2*pi;x=t.*sin(t);y=t.*cos(t);plot(x,y);grid on;3t=0:pi/1000:pi;x=sin(3.*t).*cos(t);y1=sin(3.*t).*sin(t);y2=2.*x-0.5;plot(x,y1,'k',x,y2);hold on;k=find(abs(y1-y2)<1e-2);x1=x(k);y3=2.*x1-0.5;plot(x1,y3,'rp');4x=-2:0.01:2;y=sin(1./x);subplot(2,1,1);plot(x,y);subplot(2,1,2);fplot('sin(1./x)',[-2,2],1e-4);5（1）i=-4*pi:0.1:10;j=12./sqrt(i);polar(i,j);title('{\rho}=12/sqrt(\theta)')（2）a=-pi/6:0.01:pi/6;b=3.*asin(a).*cos(a)./((sin(a)).^3+(cos(a)).^3); polar(a,b);6（1）[u,v]=meshgrid(-4:0.1:4);x=3.*u.*sin(v);y=2.*u.*cos(v);z=4.*u.^2;subplot(2,1,1);mesh(x,y,z);subplot(2,1,2);surf(x,y,z);（2）[x,y]=meshgrid(-3:0.2:3);z=-5./(1+x.^2+y.^2);subplot(1,2,1);mesh(x,y,z);subplot(1,2,2);surf(x,y,z);unit5实验指导1A=randn(10,5)x=mean(A)y=std(A)Max=max(max(A))Min=min(min(A))Sumhang=sum(A,2)SumA=sum(Sumhang)B=sort(A);C=sort(B,2,'descend');C2（1）a=0:15:90;b=a./180.*pi;s=sin(b)c=0:15:75;d=c./180.*pi;t=tan(d)e=input('请输入想计算的值：'); S=sin(e/180*pi)T=tan(e/180*pi)S1=interp1(a,s,e,'spline')T1=interp1(c,t,e,'spline')P1=polyfit(a,s,5);P2=polyfit(c,t,5);S2=polyval(P1,e)T2=polyval(P2,e)（2）n=[1,9,16,25,36,49,64,81,100]; N=sqrt(n);x=input('ji suan zhi : ');interp1(n,N,x,'cubic')3N=64;T=5;t=linspace(0,T,N);h=exp(-t);dt=t(2)-t(1);f=1/dt;X=fft(t);F=X(1:N/2+1);f=f*(0:N/2)/N;plot(f,abs(F),'-*')4P=[2,-3,0,5,13];Q=[1,5,8];p=polyder(P)q=polyder(P,Q)[a,b]=polyder(P,Q)5P1=[1,2,4,0,5];P2=[0,1,2];P3=[1,2,3];P=P1+conv(P2,P3)X=roots(P)A=[-1,1.2,-1.4;0.75,2,3.5;0,5,2.5];p=polyval(P,A)思考练习4A=rand(1,30000);a=mean(A)b=std(A)Max=max(A)Min=min(A)n=0;for i=1:30000if(A(i)>0.5)n=n+1;endendny=n/300005p=[45,74,54,55,14;78,98,45,74,12;87,98,85,52,65][M,S]=max(p)[N,H]=min(p)junzhi=mean(p,1)fangcha=std(p,1,1)zong=sum(p,2);[Max,wei]=max(zong)[Min,wei]=min(zong)[zcj,xsxh]=sort(zong,'descend')6x=[1:10:101];y=[0,1.0414,1.3222,1.4914,1.6128,1.7076,1.7853,1.8513,1.9085,1.9590,2.0043]; [p,s]=polyfit(x,y,5)a=1:5:101;y1=polyval(p,a);plot(x,y,':o',a,y1,'-*')unit6实验指导1A=[1/2,1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/6];p=[0.95,0.67,0.52]';x=A\pA=[1/2,1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/6]; p=[0.95,0.67,0.53]';x=A\pcond(A)2（1）x1=fzero(@funx1,-1)function fx=funx1(x)fx=x^41+x^3+1;（2）x2=fzero(@funx2,0.5)function fx=funx2(x)fx=x-sin(x)/x;（3）options=optimset('Display','off');x=fsolve(@fun3,[1,1,1]',options)q=fun3(x)function q=fun3(p)x=p(1);y=p(2);z=p(3);q(1)=sin(x)+y^2+log(z)-7;q(2)=3*x+2^y-z^3+1;q(3)=x+y+z-5;3（1）t0=0;tf=5;y0=1;[t,y]=ode23(@fun4,[t0,tf],y0);t'y'function yp=fun4(t,y)yp=-(1.2+sin(10*t))*y;（2）t0=0;tf=5;y0=1;[t,y]=ode23(@fun5,[t0,tf],y0);t'y'function yp=fun5(t,y)yp=cos(t)-y/(1+t^2);4x=fminbnd(@mymin,0,2);-mymin(x)function fx=mymin(x)fx=-(1+x.^2)/(1+x.^4);5options=optimset('Display','off');[x,fval]=fmincon(@fun6,[0,0,0],[],[],a,b,lb,ub)-fvalfunction f=fun6(x)f=-(sqrt(x(1))+(400-x(1))*1.1+(sqrt(x(2))+(400-x(1))*1.1-x(2))*1.1+sqrt(3)+(((400-x(1))*1.1-x(2))*1.1-x(3))*1.1+sqrt(x(x4)));思考练习1（1）A=[2,3,5;3,7,4;1,-7,1];B=[10,3,5]';C1=inv(A)*BC2=A\B[L,U]=lu(A);x=U\(L\B)（2）A=[6,5,-2,5;9,-1,4,-1;3,4,2,-2;3,-9,0,2];B=[-4,13,1,11]';C1=inv(A)*BC2=A\B[L,U]=lu(A);x=U\(L\B)2（1）x1=fzero(@funx1,1.5)function fx=funx1(x)fx=3*x+sin(x)-exp(x);（2）x1=fzero(@funx2,1)function fx=funx2(x)fx=x-1/x+5;（3）options=optimset('Display','off');x=fsolve(@fun3,[3,0]',options)q=fun3(x)function q=fun3(p)x=p(1);y=p(2);q(1)=x^2+y^2-9;q(2)=x+y-1;3（1）t0=0;tf=5;y0=[0,1];[t,y]=ode45(@vdpol,[t0,tf],y0);[t,y]function ydot=vdpol(t,y);ydot(1)=(2-3*y(2)-2*t*y(1))./(1+t^2);ydot(2)=y(1);ydot=ydot';（2）t0=0;tf=5;y0=[1;0;2];[t,y]=ode45(@vdpoll,[t0,tf],y0);[t,y]function ydot=vdpoll(t,y);ydot(1)=cos(t)-y(3)./(3+sin(t))+5*y(1).*cos(2*t)/((t+1).^2)-y(2); ydot(2)=y(1);ydot(3)=y(2);ydot=ydot';4x=fminbnd(@mymin,0,pi);-mymin(x)function fx=mymin(x)fx=-sin(x)-cos(x.^2);5[x,y1]=fminbnd(@mymax,0,1.5);-y1function fx=mymax(x);fx=-(9*x+4*x.^3-12*x.^2);unit7实验指导1（1）format longfx=inline('sin(x)./x');[I,n]=quadl(fx,0,2,1e-10)（2）format longfx=inline('1./((x-0.3).^2+0.01)-1./((x-0.9).^2+0.04)-6');[I,n]=quad(fx,0,1,1e-10)2（1）global ki;ki=0;I=dblquad(@fxy,0,1,0,1)ki（2）f=inline('abs(cos(x+y))','x','y');I=dblquad(f,0,pi,0,pi)3X=0.3:0.2:1.5;F=[0.3895,0.6598,0.9147,1.1611,1.3971,1.6212,1.8325];trapz(X,F)4p=0:pi/5:2*pi;for n=1:3nDX=diff(sin(p),n)end5f=inline('sin(x)./(x+cos(2.*x))');g=inline('(cos(x).*(x+cos(2*x))-sin(x).*(1-2.*sin(2*x)))/(x+cos(2.*x)).^2');x=-pi:0.01:pi;p=polyfit(x,f(x),5);dp=polyder(p);dpx=polyval(dp,x); %求dp在假设点的函数值dx=diff(f([x,3.01]))/0.01; %直接对f(x)求数值导数gx=g(x); %求函数f的导函数g在假设点的导数plot(x,dpx,x,dx,'.',x,gx,'-'); %作图思考练习2format longfx=inline('1./(1+x.^2)');[I,n]=quad(fx,-Inf,Inf,1e-10)[I,n]=quadl(fx,-Inf,Inf,1e-10)x=-100000:100000;y=1./(1+x.^2);trapz(x,y)format short3(1)format longfx=inline('log(1+x)./(1+x.^2)');[I,n]=quad(fx,0,1,1e-10)(2)format longfx=inline('sqrt(cos(t.^2)+4*sin((2*t).^2)+1)'); [I,n]=quad(fx,0,2*pi,1e-10)4f=inline('4.*x.*z.*exp(-z.^2.*y-x.^2)');I=triplequad(f,0,pi,0,pi,0,1)5f=inline('sin(x)');g=inline('cos(x)');x=0:0.01:2*pi;p=polyfit(x,f(x),5);dp=polyder(p);dpx=polyval(dp,x);dx=diff(f([x,2*pi+0.01]))/0.01;gx=g(x);plot(x,dpx,x,dx,'.',x,gx,'-')unit8实验指导1syms x y;s=x^4-y^4;factor(s)factor(5135)2syms x;f=(x-2)/(x^2-4);limit(f,x,2)sym x;f=(sqrt(pi)-sqrt(acos(x)))/sqrt(x+1);limit(f,x,-1,'right')3sym x;f=sin(1/x);diff(f,'x')diff(f,'x',2)sym x;f=(1-cos(2*x))/x;diff(f,'x')diff(f,'x',2)4sym x;f=sqrt(exp(x)+1);int(f,'x')syms x y;f=x/(x+y);int(f,'y')sym x;f=exp(x)*(1+exp(x))^2;int(f,'x',0,log(2))sym x;f=x*log(x);int(f,'x',1,exp(1))5sym x;s=symsum((-1)^(x+1)/x,1,Inf)sym y;z=symsum(y^(2*y-1)/(2*y-1),1,Inf)6sym x;f1=(exp(x)+exp(-x))/2;f2=sqrt(x^3-2*x+1);taylor(f1,x,5,0)taylor(f2,x,6,0)7syms x y a;x=solve('x^3+a*x+1=0','x')x=solve('sin(x)+2*cos(x)-sqrt(x)=0','x')[x y]=solve('log(x/y)=9','exp(x+y)=3','x','y')8syms n;[x,y]=dsolve('x*(D2y)+(1-n)*(Dy)+y=0','y(0)=0','Dy(0)=0','x')思考练习2syms x B1 B2 a bs1=2*((cos(x))^2)-(sin(x))^2;s2=sin(B1)*cos(B2)-cos(B1)*sin(B2);s3=sqrt((a+sqrt(a^2-b))/2)+sqrt((a-sqrt(a^2-b))/2); s4=(4*x^2+8*x+3)/(2*x+1);h1=simplify(s1)h2=simplify(s2)h3=simplify(s3)h4=simplify(s4)3syms x a;f=abs(x)/x;limit(f,x,0,'left')f=(x+a/x)^x;limit(f,x,inf)4syms x y mf=sqrt(x+sqrt(x+sqrt(x)));m=diff(f,'x')diff(m,'x')syms x y;f=x+y-sqrt(x^2+y^2);z1=diff(f,x)z2=diff(f,y)5syms x y;f=x+y-sqrt(x^2+y^2);z1=diff(f,x)z2=diff(f,y)sym x;f=1/(asin(x)^2*(1-x^2)^(1/2));int(f)6syms xf=1/(1+x);int(f,0,4)sym x;f=x^3*sin(x)^2/(x^6+2*x^4+1);int(f,-1,1)sym x;f=x^3*sin(x)^2/(x^6+2*x^4+1);int(f,-1,1)7syms ns=symsum(1/4^n,1,inf)sym n;s=symsum(((n+1)/n)^(1/2),1,inf)eval(y)8syms xf=tan(x);taylor(f,x,3,0)syms xf=sin(x)^2;taylor(f,x,5,0)9syms xx=solve('log(1+x)-5/(1+sin(x))=2','x')syms x y z[x y z]=solve('4*x^2/(4*x^2+1)=y','4*y^2/(4*y^2+1)=z','4*z^2/(4*z^2+1)=x','x','y','z') 10[x ,y]=dsolve('Dx=3*x+4*y','Dy=5*x-7*y','x(0)=0','y(0)=1','t')。
MATLAB基础教程-薛山第二版-课后习题答案
nums1=[7 1 3 5 32 12 1 99 10 24];
nums2=[54 1 456 9 20 45 48 72 61 32 10 94 11];
nums3=[44 11 25 41 84 77 998 85 2 3 15];
编写脚本文件创建相应的3个向量：newNums1、newNums2和newNums3，分别包含以上3个向量中从第一元素开始且间隔取值的元素。
矩阵是一个二维数组，所以矩阵的加、减、数乘等运算与数组运算是一致的。但有两点要注意：
(1)对于乘法、乘方和除法等三种运算，矩阵运算与数组运算的运算符及含义都不同：矩阵运算按线性变换定义，使用通常符号；数组运算按对应元素运算定义，使用点运算符；
(2)数与矩阵加减、矩阵除法在数学是没有意义的，在MATLAB中为简便起见，定义了这两类运算
6．设 ， ， ， ，判断下面表达式的值。
(1)
(2)
(3)
(4)
(5)
(6)
clear,clc
a=39;b=58;c=3;d=7;
a>b,a<c,a>b&&b>c,a==d,a|b>c,~~d
7．编写脚本，计算上面第2题中的表达式。
clear,clc
disp('sin(60)=');
disp(sind(60))
实验04051002MATLAB科学计算及绘图
操作成绩
报告成绩
1
1)熟悉MATLAB所提供的常用数值计算的函数(方程（组）的求解、插值、拟合)；
2)掌握MATLAB二维图形绘制命令及其图形控制（plot、loglog、contour、polar等）；
MATLAB基础教程薛山第二版课后习题答案讲解
《MATLAB及应用》实验指导书《MATLAB及应用》实验指导书班级： T1243-7姓名：柏元强学号： 20120430724总评成绩：汽车工程学院电测与汽车数字应用中心目录实验04051001 MATLAB语言基础 (1)实验04051002 MATLAB科学计算及绘图 (18)实验04051003 MATLAB综合实例编程 (31)实验04051001 MATLAB语言基础1实验目的1)熟悉MATLAB的运行环境2)掌握MATLAB的矩阵和数组的运算3)掌握MATLAB符号表达式的创建4)熟悉符号方程的求解2实验内容第二章1.创建double的变量，并进行计算。
(1)a=87，b=190，计算 a+b、a-b、a*b。
clear,clca=double(87);b=double(190);a+b,a-b,a*b(2)创建 uint8 类型的变量，数值与(1)中相同，进行相同的计算。
clear,clca=uint8(87);b=uint8(190);a+b,a-b,a*b2．计算：(1)() sin 60(2) e3(3) 3cos 4⎛⎫π ⎪⎝⎭clear,clc a=sind(60) b=exp(3) c=cos(3*pi/4)3．设2u =，3v =，计算：(1) 4log uvv(2) ()22euv v u +-(3)clear,clc u=2;v=3; a=(4*u*v)/log(v) b=((exp(u)+v)^2)/(v^2-u) c=(sqrt(u-3*v))/(u*v)4．计算如下表达式：(1) ()()3542i i -+(2)()sin 28i -clear,clc (3-5*i)*(4+2*i) sin(2-8*i)5．判断下面语句的运算结果。
(1) 4 < 20 (2) 4 <= 20(3) 4 == 20 (4) 4 ~= 20 (5) 'b'<'B'clear,clc4 < 20 , 4 <= 20,4 == 20,4 ~= 20,'b'<'B'6．设39a =，58b =，3c =，7d =，判断下面表达式的值。
matlab基础教程课后答案
matlab基础教程课后答案MATLAB是一种用于数值计算和数据可视化的高级编程语言和环境。
它广泛应用于科学、工程和金融等领域。
下面是MATLAB基础教程的课后答案，希望能帮助你更好地理解和掌握MATLAB的基本概念和操作。
1.MATLAB基础知识 (1) MATLAB是什么？它的特点和优势是什么？答：MATLAB是一种高级的编程语言和环境，用于数值计算和数据可视化。
它的特点和优势包括：•简单易学：MATLAB语法与数学表达式类似，易于理解和使用。
•功能强大：MATLAB提供了丰富的数值计算和数据处理函数，可以方便地进行各种科学计算和数据分析。
•可视化能力强：MATLAB内置了丰富的绘图函数，可以直观地展示数据和结果。
•扩展性好：MATLAB支持用户自定义函数和工具包，可以根据需要扩展功能。
•跨平台：MATLAB可以在Windows、Mac和Linux等操作系统上运行。
(2) 如何定义变量并赋值？答：在MATLAB中，可以使用等号（=）将一个值赋给一个变量。
例如，要将数值3赋给变量x，可以输入x=3。
(3) 如何进行算术运算？答：MATLAB支持常见的算术运算，包括加法、减法、乘法和除法。
可以使用+、-、*和/符号进行运算。
例如，要计算2加上3的结果，可以输入2+3。
(4) 如何进行矩阵运算？答：在MATLAB中，可以使用方括号（[]）来定义矩阵。
例如，要定义一个2×2的矩阵A，可以输入A=[1 2; 3 4]。
可以使用符号进行矩阵乘法运算。
例如，要计算矩阵A乘以矩阵B的结果，可以输入A B。
2.MATLAB的基本函数和操作 (1) 如何使用help命令查看函数的帮助文档？答：在MATLAB命令窗口中，可以使用help命令查看函数的帮助文档。
例如，要查看sin函数的帮助文档，可以输入help sin。
(2) 如何生成随机数？答：MATLAB提供了rand函数用于生成服从均匀分布的随机数。
MATLAB第二版课后答案Unit3-8
end
movie(m,4);
思考练习
2（1）
x=-5:0.1:5;
y=(1./(2*pi)).*exp((-(x.^2))/2);
（2）
t=-pi/3:pi/50:pi/3;
r=5.*((sin(t)).^2)./cos(t);
polar(t,r);
4（1）
t=0:pi/50:2*pi;
x=exp(-t./20).*cos(;
y=exp(-t./20).*sin(t);
z=t;
plot3(x,y,z);
grid on;
（2）
p='A';
case num2cell(80:89)
p='B';
case num2cell(70:79)
p='C';
case num2cell(60:69)
p='D';
otherwise
p='E';
end
price=p
（2）n=input('请输入成绩');
if n>=90&n<=100
p='A';
a=n(1)；
b=n(1);
for m=2:20
if n(m)>a
a=n(m);
elseif n(m)<b
b=n(m);
end
end
max=a
min=b
法2
n=[1,5,56,4,3,476,45,6,3,76,45,6,4,3,6,4,23,76,908,6];
min=min(n)
max=max(n)
Matlab程序设计教程(第二版)课后参考答案
Matlab程序设计教程（第二版）课后参考答案第一章实验1.实验一第1题2.自己验证总结Matlab的主要优点3.实验一第2题4.实验一第3题5.网站思考练习1.启动见书P5 退出见书P62.Matlab主要功能见书P2 4种功能3.分行输入行末尾加续行符，即三个点“…”4.见书P115.直接在命令窗口输入fac第二章实验1.(1) w=sqrt(2)*(1+0.34245*10^(-6))w = 1.4142(2) a=3.5;b=5;c=-9.8;x=(2*pi*a+(b+c)/(pi+a*b*c)-exp(2))/(tan(b+c)+a)x =0.9829(3) a=3.32;b=-7.9;y=2*pi*a^2*((1-pi/4)*b-(0.8333-pi/4)*a)y = -128.4271(4) t=[2,1-3i;5,-0.65];z=0.5*exp(2*t)*log(t+sqrt(1+t.*t))z =1.0e+004 *0.0048 + 0.0002i 0.0048 - 0.0034i1.58992.0090 - 1.3580i2．实验二第1题3．实验二第2题4． H=hilb(5);P=pascal(5);Hh=det(H)Hh = 3.7493e-012Hp=det(P)Hp = 1Th=cond(H)Th = 4.7661e+005Tp=cond(P)Tp = 8.5175e+003条件数越趋近于1，矩阵的性能越好，所以帕斯卡矩阵性能更好。
5． A=[-29,6,18;20,5,12;-8,8,5]A =-29 6 1820 5 12-8 8 5[V,D]=eig(A)V =0.7130 0.2803 0.2733-0.6084 -0.7867 0.87250.3487 0.5501 0.4050D =-25.3169 0 00 -10.5182 00 0 16.8351V为A的特征向量，D为A的特征值。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
5．判断下面语句的运算结果。
(1) 4 < 20
(2) 4 <= 20
(3) 4 == 20
(4) 4 ~= 20
(5) 'b'<'B'
clear,clc
4 < 20 , 4 <= 20,4 == 20,4 ~= 20,'b'<'B'
6．设 ， ， ， ，判断下面表达式的值。
(1)
vec3 = vec1+vec2;
d = vec3.*a;
vec4 = find(a > 5);
e = a(vec4) + 5;
vec5 = find(a < 5);
f = vec5.^2;
disp('c=');disp(c)
disp('d=');disp(d)
disp('e=');disp(e)
disp('f=');disp(f)
求y值为多少？
2. 在MATLAB中运行以下命令：
a = [3, 7, 2, 7, 9, 3, 4, 1, 6];
b = [7];
a(4) = [];
vec1 = a==b;
vec2 = mod(a,2)==0;
c = sum(vec1);
vec3 = vec1+vec2;
d = vec3.*a;
vec4 = find(a > 5);
a=double(87);
b=double(190);
a+b,a-b,a*b
(2)创建 uint8 类型的变量，数值与(1)中相同，进行相同的计算。
clear,clc
a=uint8(87);
b=uint8(190);
a+b,a-b,a*b
2．计算：
(1)
(2) e3
(3)
clear,clc
a=sind(60)
disp('a==d');disp(a==d)
disp('a|b>c');disp(a|b>c)
disp('~~d');disp(~~d)
第三章
1. 在命令提示符下输入以下两条命令：
>> x = [ 9 3 0 6 3]
>> y = mod((sqrt(length(((x+5).*[1 2 3 4 5]))*5)),3)
vecR=vec';
disp('vecR=[]');disp(vecR)
a=length(vec);
vecB=[vec(a/2+1:a) vec(1:a/2)];
disp('vecB=[]');disp(vecB)
C=find(vec<45);
vecS=vec(C);
disp('vecS=[]');disp(vecS)
vec3R=vec(end:-4:1);
disp('vec3R=[]');disp(vec3R)
vecN=find(vec==2|vec==4);
e = a(vec4) + 5;
vec5 = find(a < 5);
f = vec5.^2;
求c、d、e、f的值。
clear,clc
a = [3, 7, 2, 7, 9, 3, 4, 1, 6];
b = [7];
a(4) = [];
vec1 = a==b;
vec2 = mod(a,2)==0;
c = sum(vec1);
(2)创建一个新的向量vecB，使其为vec中的前半部分与后半部分对换的结果，这样vecB包含的元素为vec的后半部分紧接着vec的前半部分。
(3)创建一个新的向量vecS，使其包含vec中所有小于45的元素，且元素按照vec中的顺序排列。
(4)创建一个新的向量vec3R，使其从vec中从最后一个元素开始，并且间隔三个元素取一个元素，直到第一个元素为止。
《MATLAB及应用》实验指导书
《MATLAB及应用》实验指导书
班级：T1243-7
姓名：柏元强
学号：20120430724
总评成绩：
汽车工程学院
电测与汽车数字应用中心
实验04051001 MATLAB语言基础1
实验04051002 MATLAB科学计算及绘图13
实验04051003 MATLAB综合实例编程31
实验04051001MATLAB语言基础
操作成绩
报告成绩
1
1)熟悉MATLAB的运行环境
2)掌握MATLAB的矩阵和数组的运算
3)掌握MATLAB符号表达式的创建
4)熟悉符号方程的求解
2
第二章
1.创建double的变量，并进行计算。
(1)a=87，b=190，计算 a+b、a-b、a*b。
clear,clc
(2)
(3)
(4)
(5)
(6)
clear,clc
a=39;b=58;c=3;d=7;
a>b,a<c,a>b&&b>c,a==d,a|b>c,~~d
7．编写脚本，计算上面第2题中的表达式。
clear,clc
disp('sin(60)=');
disp(sind(60))
disp('exp(3)=');
disp(exp(3))
disp('cos(3*pi/4)=');
disp(cos(3*pi/4))
8．编写脚本，输出上面第6题中的表达式的值。
clear,clc
a=39;b=58;c=3;d=7;
disp('a>b');disp(a>b)
disp('a<c');disp(a<c)
disp('a>b&&b>c');,disp(a>b&&b>c)
3. 向量操作时MATLAB的主要部分，使用给出的向量来做下面的练习。
注意：
不要直接给出下列问题中任何一个的最终结果，不要在问题的任何部分使用迭代。
vec=[4 5 2 8 4 7 2 64 2 57 2 45 7 43 2 5 7 3 3 6253 3 4 3 0 -65 -343]
(1)创建一个新的向量vecR,使其为vec的转置。
b=exp(3)
c=cos(3*pi/4)
3．设 ， ，计算：
(1)
(2)
(3)
clear,clc
u=2;v=3;
a=(4*u*v)/log(v)
b=((exp(u)+v)^2)/(v^2-u)
c=(sqrt(u-3*v))/(u*v)
4．计算如下表达式：
(3-5*i)*(4+2*i)
(5)创建一个新的向量vecN，使其包含vec中所有等于2或4的元素的索引值。
(6)创建一个新的向量vecG，使其包含vec中去掉索引值为奇数且取值为2或4的元素后的所有元素。
clear,clc
vec=[4 5 2 8 4 7 2 64 2 57 2 45 7 43 2 5 7 3 3 6253 3 4 3 0 -65 -343];