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+a
2（1）
n=input('请输入成绩');
switch n
case num2cell(90:100)
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';
elseif n>=80&n<=89
p='B';
elseif n>=70&n<=79
p='C';
elseif n>=60&n<=69
p='D';
else
p='E';
end
price=p
（3）try
n;
catch
price='erroe'
end
3
n=[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:20
if n(m)>a
a=n(m);
elseif n(m) 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) 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 y 5 y1=0; y2=1; n=input('请输入n的值：'); for i=1:n y1=y1+1/i^2; y2=y2*((4*i*i)/((2*i-1)*(2*i+1))); 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); else disp(A([6],:)); disp(lasterr); end 7（1） f=[]; f(n)=n+10*log(n^2+5); end y=f(40)/(f(30)+f(20)) （2） f=[];a=0; for n=1:40 f(n)=a+n*(n+1); a=f(n); end y=f(40)/(f(30)+f(20)) 8 y=0; m=input('输入m的值:'); n=input('输入n值：'); for i=1:n y=y+i^m; end y ************************************************************ function s=shi8_1(n,m) s=0; for i=1:n s=s+i^m; end ************************************************************ shi8_1(100,1)+shi8_1(50,2)+shi8_1(10,1/2) 思考练习 2 N=[1,2,3,4,5]; 2.*N N./2 1./N 1./N.^2 3 s=fix(100*rand(1,20)*9/10+10) y=sum(s)/20 j=0; for i=1:20 if s(i) j=j+1; A(j)=s(i); else continue; end end A 4 y1=0;y2=0; n=input('请输入n的值：'); for i=1:n y1=y1+-(-1)^i/(2*i-1); y2=y2+1/4^i; end y1 y2 unit4 实验指导 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; 2 t=-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:20 axis([-i,i,-i,i,-i,i]) m(:,i)=getframe; end movie(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; 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 x=-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 实验指导 1 A=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'); C 2（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') 3 N=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),'-*') 4 P=[2,-3,0,5,13]; Q=[1,5,8]; p=polyder(P) q=polyder(P,Q) [a,b]=polyder(P,Q) 5 P1=[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) 思考练习 4 A=rand(1,30000); a=mean(A) b=std(A) Max=max(A) Min=min(A) n=0; for i=1:30000 if(A(i)>0.5) n=n+1; end end n y=n/30000 5 p=[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') 6 x=[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 实验指导 1 A=[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\p A=[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\p cond(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); 4 x=fminbnd(@mymin,0,2); -mymin(x) function fx=mymin(x) fx=-(1+x.^2)/(1+x.^4); 5 options=optimset('Display','off'); [x,fval]=fmincon(@fun6,[0,0,0],[],[],a,b,lb,ub) -fval function 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)*B C2=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)*B C2=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'; 4 x=fminbnd(@mymin,0,pi); -mymin(x) function fx=mymin(x) fx=-sin(x)-cos(x.^2); 5 [x,y1]=fminbnd(@mymax,0,1.5); -y1 function fx=mymax(x); fx=-(9*x+4*x.^3-12*x.^2); unit7 实验指导 1（1） format long fx=inline('sin(x)./x'); [I,n]=quadl(fx,0,2,1e-10) （2） format long fx=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) 3 X=0.3:0.2:1.5; F=[0.3895,0.6598,0.9147,1.1611,1.3971,1.6212,1.8325]; trapz(X,F) 4 p=0:pi/5:2*pi; for n=1:3 n DX=diff(sin(p),n) end 5 f=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,'-'); %作图 思考练习 2 format long fx=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 short 3(1) format long fx=inline('log(1+x)./(1+x.^2)'); [I,n]=quad(fx,0,1,1e-10) (2) format long fx=inline('sqrt(cos(t.^2)+4*sin((2*t).^2)+1)'); [I,n]=quad(fx,0,2*pi,1e-10) 4 f=inline('4.*x.*z.*exp(-z.^2.*y-x.^2)'); I=triplequad(f,0,pi,0,pi,0,1) 5 f=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 实验指导 1 syms x y; s=x^4-y^4; factor(s) factor(5135) 2 syms 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') 3 sym 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) 4 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)) 5 sym x; s=symsum((-1)^(x+1)/x,1,Inf) sym y; z=symsum(y^(2*y-1)/(2*y-1),1,Inf) 6 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) 7 syms 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') 8 syms n; [x,y]=dsolve('x*(D2y)+(1-n)*(Dy)+y=0','y(0)=0','Dy(0)=0','x') 思考练习 syms x B1 B2 a b s1=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) 3 syms x a; f=abs(x)/x; limit(f,x,0,'left') f=(x+a/x)^x; limit(f,x,inf) 4 syms 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^2+y^2); z1=diff(f,x) z2=diff(f,y) 5 syms 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) 6 syms x f=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) 7 syms n s=symsum(1/4^n,1,inf) sym n; s=symsum(((n+1)/n)^(1/2),1,inf) eval(y) 8 syms x f=tan(x); taylor(f,x,3,0) syms x f=sin(x)^2; taylor(f,x,5,0) 9 syms x x=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 第二版课后答案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+a 2（1） n=input('请输入成绩'); switch n case num2cell(90:100) 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'; elseif n>=80&n<=89 p='B'; elseif n>=70&n<=79 p='C'; elseif n>=60&n<=69 p='D'; else p='E'; end price=p （3）try n; catch price='erroe' end 3 n=[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:20 if n(m)>a a=n(m); elseif n(m)=0 disp(A([n],:)); elseif n<0 disp(lasterr); else disp(A([6],:)); disp(lasterr); end 7（1） f=[]; 实验指导 1、 n=input('请输入一个三位数：'); a=fix(n/100); b=fix((n-a*100)/10); c=n-a*100-b*10; d=c*100+b*10+a 2（1） n=input('请输入成绩'); switch n case num2cell(90:100) 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'; elseif n>=80&n<=89 p='B'; elseif n>=70&n<=79 p='C'; elseif n>=60&n<=69 p='D'; else p='E'; end price=p （3）try n; catch price='erroe' end 3 n=[1,5,56,4,3,476,45,6,3,76,45,6,4,3,6,4,23,76,908,6]; a=n(1)； for m=2:20 if n(m)>a a=n(m); elseif n(m)=0 disp(A([n],:)); elseif n<0 disp(lasterr); else disp(A([6],:)); disp(lasterr); end 7（1） f=[]; for n=1:40 第2章符号运算 习题2及解答 1 说出以下四条指令产生的结果各属于哪种数据类型，是“双精度” 对象，还是“符号”符号对象 3/7+; sym(3/7+; sym('3/7+'); vpa(sym(3/7+) 〖目的〗 不能从显示形式判断数据类型，而必须依靠class指令。 〖解答〗 c1=3/7+ c2=sym(3/7+ c3=sym('3/7+') c4=vpa(sym(3/7+) Cs1=class(c1) Cs2=class(c2) Cs3=class(c3) Cs4=class(c4) c1 = c2 = 37/70 c3 = c4 = Cs1 = double Cs2 = sym Cs3 = sym Cs4 = sym 2 在不加专门指定的情况下，以下符号表达式中的哪一个变量被认 为是自由符号变量. sym('sin(w*t)'),sym('a*exp(-X)'),sym('z*exp(j*th)') 〖目的〗 理解自由符号变量的确认规则。 〖解答〗 symvar(sym('sin(w*t)'),1) ans = w symvar(sym('a*exp(-X)'),1) ans = a symvar(sym('z*exp(j*th)'),1) ans = z 5求符号矩阵???? ??????=3332 31 232221 131211 a a a a a a a a a A 的行列式值和逆，所得结果应采用“子表达式置换”简洁化。 〖目的〗 理解subexpr 指令。 〖解答〗 A=sym('[a11 a12 a13;a21 a22 a23;a31 a32 a33]') DA=det(A) IA=inv(A); [IAs,d]=subexpr(IA,d) A = [ a11, a12, a13] [ a21, a22, a23] [ a31, a32, a33] DA = a11*a22*a33 - a11*a23*a32 - a12*a21*a33 + a12*a23*a31 + a13*a21*a32 - a13*a22*a31 IAs = [ d*(a22*a33 - a23*a32), -d*(a12*a33 - a13*a32), d*(a12*a23 - a13*a22)] [ -d*(a21*a33 - a23*a31), d*(a11*a33 - a13*a31), -d*(a11*a23 - a13*a21)] [ d*(a21*a32 - a22*a31), -d*(a11*a32 - a12*a31), d*(a11*a22 - a12*a21)] d = 1/(a11*a22*a33 - a11*a23*a32 - a12*a21*a33 + a12*a23*a31 + a13*a21*a32 - a13*a22*a31) 8（1）通过符号计算求t t y sin )(=的导数 dt dy 。（2）然后根据此结果，求- =0t dt dy 和2 π = t dt dy 。 〖目的〗 diff, limit 指令的应用。 如何理解运行结果。 〖解答〗 syms t ones表示1矩阵 zeros表示0矩阵 ones(4)表示4x4的1矩阵 zeros(4)表示4x4的0矩阵 zeros(4,5)表示4x5的矩阵 eye(10,10)表示10x10的单位矩阵 rand(4,5)表示4x5的伴随矩阵 det(a)表示计算a的行列式 inv(a)表示计算a的逆矩阵 Jordan(a)表示求a矩阵的约当标准块rank(a)表示求矩阵a的秩 [v,d]=eig(a)对角矩阵 b=a’表示求a矩阵的转置矩阵 sqrt表示求平方根 exp表示自然指数函数 log自然对数函数 abs绝对值 第一章 一、5（1） b=[97 67 34 10;-78 75 65 5;32 5 -23 -59]; >> c=[97 67;-78 75;32 5;0 -12]; >> d=[65 5;-23 -59;54 7]; >> e=b*c e = 5271 11574 -11336 664 1978 3112 （2）a=50:1:100 二、1 、x=-74; y=-27; z=(sin(x.^2+y.^2))/(sqrt(tan(abs(x+y)) )+pi) z = -0.0901 2、a=-3.0:0.1:3.0; >> b=exp(-0.3*a).*sin(a+0.3) 3、x=[2 4;-0.45 5]; y=log(x+sqrt(1+x.^2))/2 y = 0.7218 1.0474 -0.2180 1.1562 4、a*b表示a矩阵和b矩阵相乘 a.*b表示a矩阵和b矩阵单个元素相乘A(m,n)表示取a矩阵第m行，第n列 A(m,:)表示取a矩阵第m行的全部元素A(:,n)表示取a矩阵的第n列全部元素 A./B表示a矩阵除以b矩阵的对应元素， B.\A等价于A./B A.^B表示两个矩阵对应元素进行乘方运算 A.^2表示a中的每个元素的平方 A^2表示A*A 例：x=[1,2,3]; y=[4,5,6]; z=x.^y z= 1 3 2 729 指数可以是标量（如y=2）.底数也可以是标量（如x=2） 5、a=1+2i; >> b=3+4i; >> c=exp((pi*i)/6) c = 0.8660 + 0.5000i d=c+a*b/(a+b) d = 1.6353 + 1.8462i 第二章 二、4、（1） 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)]) ans = 大学社区网收集整理https://www.360docs.net/doc/f314663373.html, 评分 日期湖南商学院北津学院实验报告 课程名称MATLAB科学计算编程语言 实验名称MATLAB程序设计 专业班级信科1121班 姓名xxx 学号xxx 实验日期2012年11月5日 2012—2013学年度第一学期 一、实验目的 1.掌握利用if语句、switch语句实现选择结构的方法。 2.掌握利用for语句、while语句实现循环结构的方法。 3.熟悉利用向量运算来代替循环操作的方法并理解MATLAB程序设计的特点 4.掌握定义和调用MATLAB函数的方法。 二、实验环境 系统windows7旗舰版 处理器Intel(R)Core(TM)i7-3610M CPU @ 2.30GHz 安装内存 4.00GB （3.07GB 可用）系统类型64位操作系统运行环境 MATLAB 5.3 三、实验基本原理 利用上课所学知识解决以下问题： 1.从键盘输入一个3位数的整数，将它反向输出。如输入639，输出936。 2.输入一个百分制成绩，要求输出成绩等级A、B、C、D、E。其中90~100分为A，80~89分为B，70~79分为C，60~69分为D，60分以下为E。 要求： (1)分别用if 语句和switch 语句实现。 (2)输入百分制成绩后要判断该成绩的合理性，对不合理的成绩应输出出错信息。3.输入20个数，求其中最大数和最小数。要求分别用循环结构和调用MATLAB 的max 函数、min 函数来实现。 4.23.0ln )3.0sin(23.03.0a a e e y a a +++?=?，当a 取-3.0、-2.9、-2.8、…、2.8、2.9、 3.0时，求各点的函数值。要求分别用顺序结构和循环结构实现。 5.当n 分别取100、1000、10000时，求下列各式的值： (1)) 6...(n 1...31211122222π=+++++(2)) 2...()12)(12()2)(2(...756653443122π =??? ?????+?????????××????????××????????××n n n n 要求分别用循环结构和向量运算（使用sum 函数）来实现。 6.建立5×6矩阵，要求输出矩阵第n 行元素。当n 值超过矩阵的行数时，自动转为输出矩阵最后一行元素，并给出出错信息。 7已知,o999 ) 20()30()40(f f f y += (1)当)5ln(10)(2 ++=n n n f 时，y 的值是多小。 (2)当+×+×+×=433221)(n f …+)1(+×n n 时，y 的值是多小。 8.先用函数的递归调用定义一个函数文件求 ∑=n i m i 1，然后调用该函数文件求 第1章 MATLAB概论 1.1与其他计算机语言相比较，MATLAB语言突出的特点是什么？ MATLAB具有功能强大、使用方便、输入简捷、库函数丰富、开放性强等特点。 1.2 MATLAB系统由那些部分组成？ MATLAB系统主要由开发环境、MATLAB数学函数库、MATLAB语言、图形功能和应用程序接口五个部分组成。 1.3 安装MATLAB时，在选择组件窗口中哪些部分必须勾选，没有勾选的部分以后如何补安装？ 在安装MATLAB时，安装内容由选择组件窗口中个复选框是否被勾选来决定，可以根据自己的需要选择安装内容，但基本平台（即MATLAB选项）必须安装。第一次安装没有选择的内容在补安装时只需按照安装的过程进行，只是在选择组件时只勾选要补装的组件或工具箱即可。 1.4 MATLAB操作桌面有几个窗口？如何使某个窗口脱离桌面成为独立窗口？又如何将脱离出去的窗口重新放置到桌面上？ 在MATLAB操作桌面上有五个窗口，在每个窗口的右上角有两个小按钮，一个是关闭窗口的Close按钮，一个是可以使窗口成为独立窗口的Undock按钮，点击Undock按钮就可以使该窗口脱离桌面成为独立窗口，在独立窗口的view菜单中选择Dock ……菜单项就可以将独立的窗口重新防止的桌面上。 1.5 如何启动M文件编辑/调试器？ 在操作桌面上选择“建立新文件”或“打开文件”操作时，M文件编辑/调试器将被启动。在命令窗口中键入edit命令时也可以启动M文件编辑/调试器。 1.6 存储在工作空间中的数组能编辑吗？如何操作？ 存储在工作空间的数组可以通过数组编辑器进行编辑：在工作空间浏览器中双击要编辑的数组名打开数组编辑器，再选中要修改的数据单元，输入修改内容即可。 1.7 命令历史窗口除了可以观察前面键入的命令外，还有什么用途？ 命令历史窗口除了用于查询以前键入的命令外，还可以直接执行命令历史窗口中选定的内容、将选定的内容拷贝到剪贴板中、将选定内容直接拷贝到M文件中。 1.8 如何设置当前目录和搜索路径，在当前目录上的文件和在搜索路径上的文件有什么区别？ 当前目录可以在当前目录浏览器窗口左上方的输入栏中设置，搜索路径可以通过选择操作桌面的file 菜单中的Set Path菜单项来完成。在没有特别说明的情况下，只有当前目录和搜索路径上的函数和文件能够被MATLAB运行和调用，如果在当前目录上有与搜索路径上相同文件名的文件时则优先执行当前目录上的文件，如果没有特别说明，数据文件将存储在当前目录上。 1.9 在MATLAB中有几种获得帮助的途径？ 第1章MATLAB概论 1.1与其他计算机语言相比较，MA TLAB语言突出的特点是什么？ MATLAB具有功能强大、使用方便、输入简捷、库函数丰富、开放性强等特点。 1.2 MA TLAB系统由那些部分组成？ MATLAB系统主要由开发环境、MA TLAB数学函数库、MATLAB语言、图形功能和应用程序接口五个部分组成。 1.3 安装MATLAB时，在选择组件窗口中哪些部分必须勾选，没有勾选的部分以后如何补安装？ 在安装MA TLAB时，安装内容由选择组件窗口中个复选框是否被勾选来决定，可以根据自己的需要选择安装内容，但基本平台（即MATLAB选项）必须安装。第一次安装没有选择的内容在补安装时只需按照安装的过程进行，只是在选择组件时只勾选要补装的组件或工具箱即可。 1.4 MA TLAB操作桌面有几个窗口？如何使某个窗口脱离桌面成为独立窗口？又如何将脱离出去的窗口重新放置到桌面上？ 在MA TLAB操作桌面上有五个窗口，在每个窗口的右上角有两个小按钮，一个是关闭窗口的Close按钮，一个是可以使窗口成为独立窗口的Undock按钮，点击Undock按钮就可以使该窗口脱离桌面成为独立窗口，在独立窗口的view菜单中选择Dock ……菜单项就可以将独立的窗口重新防止的桌面上。 1.5 如何启动M文件编辑/调试器？ 在操作桌面上选择“建立新文件”或“打开文件”操作时，M文件编辑/调试器将被启动。在命令窗口中键入edit命令时也可以启动M文件编辑/调试器。 1.6 存储在工作空间中的数组能编辑吗？如何操作？ 存储在工作空间的数组可以通过数组编辑器进行编辑：在工作空间浏览器中双击要编辑的数组名打开数组编辑器，再选中要修改的数据单元，输入修改内容即可。 1.7 命令历史窗口除了可以观察前面键入的命令外，还有什么用途？ 命令历史窗口除了用于查询以前键入的命令外，还可以直接执行命令历史窗口中选定的内容、将选定的内容拷贝到剪贴板中、将选定内容直接拷贝到M文件中。 Matlab 基础练习题 常量、变量、表达式 1、 MATLAB 中，下面哪些变量名是合法的？（ ） （A ）_num （B ）num_ （C ）num- （D ）-num 2、 在MA TLAB 中，要给出一个复数z 的模，应该使用（ ）函数。 （A ）mod(z) （B ）abs(z) （C ）double(z) （D ）angle(z) 3、 下面属于MATLAB 的预定义特殊变量的是？（ ） （A ）eps （B ）none （C ）zero （D ）exp 4、 判断：在MA TLAB 的内存工作区中，存放一个英文字符 'a' 需要占用1个字节，存放 一个中文字符‘啊’需要占用2个字节。（ 错，都是2个字节 ） 5、 判断：MA TLAB 中，i 和j （ 对 ） 6、 判断：MA TLAB 中，pi 代表圆周率，它等于3.14。（ 错，后面还有很多位小数 ） 7、 在MA TLAB 中，若想计算的5 1)3.0sin(21+= πy 值，那么应该在MA TLAB 的指令窗中 输入的MA TLAB 指令是__y1=2*sin(0.3*pi)/(1+sqrt(5))_。 8、 在MA TLAB 中，a = 1，b = i ，则a 占_8__个字节，b 占_16_个字节，c 占________字 节。 9、 在MA TLAB 中，inf 的含义是__无穷大__，nan 的含义是__非数（结果不定）___。 数组 1、 在MA TLAB 中，X 是一个一维数值数组，现在要把数组X 中的所有元素按原来次序 的逆序排列输出，应该使用下面的（ ）指令。 （A ）X[end:1] （B ）X[end:-1:1] （C ）X (end:-1:1) （D ）X(end:1) 2、 在MA TLAB 中，A 是一个字二维数组，要获取A 的行数和列数，应该使用的MATLAB 的命令是（ ）。 （A ）class(A) （B ）sizeof(A) （C ）size(A) （D ）isa(A) 3、 在MATLAB 中，用指令x=1:9生成数组x 。现在要把x 数组的第二和第七个元素都 赋值为0，应该在指令窗中输入（ ） （A ）x([2 7])=(0 0) （B ）x([2,7])=[0,0] （C ）x[(2,7)]=[0 0] （D ）x[(2 7)]=(0 0) 4、 在MA TLAB 中，依次执行以下指令：clear;A=ones(3,4); A(:)=[-6:5];这时， 若在指令窗中输入指令b=A(:,2)'，那么，MATLAB 输出的结果应该是（ ） （A ）b = -3 -2 -1 （B ）b = -2 -1 0 1 （C ）b = -5 -1 3 （D ）b = -5 -2 1 4 5、 在MA TLAB 中，A = 1:9，现在执行如下指令L1 = ~(A>5)，则MATLAB 的执行结果应 该是L1 =___ 1 1 1 1 1 0 0 0 0___。 第二章 1 求下列表达式的值。 （1） w=sqrt(2)*(1+0.34245*10^(-6)) (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 (3) a=3.32; b=-7.9; y=2*pi*a^(2)*[(1-pi/4)*b-(0.8333-pi/4)*a] (4) t=[2,1-3*i;5,-0.65]; z=1/2*exp(2*t)*log(t+sqrt(1+t^(2))) 2 求下列表达式 A=[-1,5,-4;0,7,8;3,61,7]; B=[8,3,-1;2,5,3;-3,2,0]; (1) A+6*B A^2-B+eye (2) A*B A.*B B.*A (3) A/B B\A (4) [A,B] [A([1,3],:);B^2] 3 根据已知，完成下列操作 （1） A=[23,10,-0.778,0;41,-45,65,5;32,5,0,32;6,-9.54,54,3.14]; K=find(A>10&A<25); A(K) （2） A=[23,10,-0.778,0;41,-45,65,5;32,5,0,32;6,-9.54,54,3.14]; B=A(1:3,:) C=A(:,1:2) D=A(2:4,3:4) E=B*C （3） E f=input('输入一个数：','s'); f(end :-1:1) 2 用if语句 score=input('请输入成绩:'); if score>=90&&score<=100 disp('A'); elseif score>=80&&score<=89 disp('B'); elseif score>=70&&score<=79 disp('C'); elseif score>=60&&score<=69; disp('D'); elseif score<60&&score>=0; disp('E'); else disp('出错'); end 用switch语句 score=input('请输入成绩:'); switch fix(score/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'); otherwise disp('出错'); end 第四章1题 1) X=0:10; Y=x-x.^3/6; P lot(x,y) 习题二 1.如何理解“矩阵是MATLAB最基本的数据对象”？ 答：因为向量可以看成是仅有一行或一列的矩阵，单个数据（标量）可以看成是仅含一个元素的矩阵，故向量和单个数据都可以作为矩阵的特例来处理。 因此，矩阵是MATLAB最基本、最重要的数据对象。 2.设A和B是两个同维同大小的矩阵，问： (1)A*B和A.*B的值是否相等？ 答：不相等。 (2)A./B和B.\A的值是否相等？ 答：相等。 (3)A/B和B\A的值是否相等？ 答：不相等。 (4)A/B和B\A所代表的数学含义是什么？ 答：A/B等效于B的逆右乘A矩阵，即A*inv(B)，而B\A等效于B矩阵的逆左乘A矩阵，即inv(B)*A。 3.写出完成下列操作的命令。 (1)将矩阵A第2~5行中第1, 3, 5列元素赋给矩阵B。 答：B=A(2:5,1:2:5); 或B=A(2:5,[1 3 5]) (2)删除矩阵A的第7号元素。 答：A(7)=[] (3)将矩阵A的每个元素值加30。 答：A=A+30; (4)求矩阵A的大小和维数。 答：size(A); ndims(A); (5)将向量t的0元素用机器零来代替。 答：t(find(t==0))=eps; (6)将含有12个元素的向量x转换成34 矩阵。 答：reshape(x,3,4); (7)求一个字符串的ASCII码。 答：abs(‘123’); 或double(‘123’); (8)求一个ASCII码所对应的字符。 答：char(49); 4. 下列命令执行后，L1、L2、L3、L4的值分别是多少？ A=1:9;B=10-A;... L1=A==B; L2=A<=5; L3=A>3&A<7; L4=find(A>3&A<7); 答：L1的值为[0, 0, 0, 0, 1, 0, 0, 0, 0] L2的值为[1, 1, 1, 1, 1, 0, 0, 0, 0] L3的值为[0, 0, 0, 1, 1, 1, 0, 0, 0] L4的值为[4, 5, 6] 5. 已知 23100.7780414565532503269.5454 3.14A -????-??=????-?? 完成下列操作： (1) 取出A 的前3行构成矩阵B ，前两列构成矩阵C ，右下角32?子矩阵构 成矩阵D ，B 与C 的乘积构成矩阵E 。 答：B=A(1:3,:); C=A(:,1:2); D=A(2:4,3:4); E=B*C; (2) 分别求E 第4章数值运算 习题 4 及解答 1 根据题给的模拟实际测量数据的一组t和)(t y试用数值差分diff 或数值梯度gradient指令计算)(t y'曲线绘制 y'，然后把)(t y和)(t 在同一张图上，观察数值求导的后果。（模拟数据从prob_data401.mat获得） 〖目的〗 ●强调：要非常慎用数值导数计算。 ●练习mat数据文件中数据的获取。 ●实验数据求导的后果 ●把两条曲线绘制在同一图上的一种方法。 〖解答〗 （1）从数据文件获得数据的指令 假如prob_data401.mat文件在当前目录或搜索路径上 clear load prob_data401.mat （2）用diff求导的指令 dt=t(2)-t(1); yc=diff(y)/dt; %注意yc的长度将比y短1 plot(t,y,'b',t(2:end),yc,'r') （3）用gradent求导的指令（图形与上相似） dt=t(2)-t(1); yc=gradient(y)/dt; plot(t,y,'b',t,yc,'r') grid on 〖说明〗 ● 不到万不得已，不要进行数值求导。 ● 假若一定要计算数值导数，自变量增量dt 要取得比原有数据相对误差高1、2个量级 以上。 ● 求导会使数据中原有的噪声放大。 2 采用数值计算方法，画出dt t t x y x ? =0 sin )(在]10 ,0[区间曲线，并计算)5.4(y 。 〖提示〗 ● 指定区间内的积分函数可用cumtrapz 指令给出。 ● )5.4(y 在计算要求不太高的地方可用find 指令算得。 〖目的〗 ● 指定区间内的积分函数的数值计算法和cumtrapz 指令。 ● find 指令的应用。 〖解答〗 dt=1e-4; t=0:dt:10; t=t+(t==0)*eps; f=sin(t)./t; s=cumtrapz(f)*dt; plot(t,s,'LineWidth',3) ii=find(t==4.5); s45=s(ii) s45 = 第三章1. (1)A=eye(3) (2)A=100+100*rand(5,6) (3)A=1+sqrt(0.2)*randn(10,50) (4)B=ones(size(A)) (5)A+30*eye(size(A)) (6)B=diag(diag(A)) 2. B=rot90(A) C=rot90(A,-1) 3. B=inv(A) ;A的逆矩阵 C=det(A) ;A的行列式的值 D=A*B E=B*A D=E 因此A与A-1是互逆的。 4. A=[4 2 -1;3 -1 2;12 3 0]; b=[2;10;8]; x=inv(A)*b x = -6.0000 26.6667 27.3333 5. (1) diag(A) ；主对角线元素 ans = 1 1 5 9 triu(A) ；上三角阵 ans = 1 -1 2 3 0 1 -4 2 0 0 5 2 0 0 0 9 tril(A) ；下三角阵 ans = 1 0 0 0 5 1 0 0 3 0 5 0 11 15 0 9 rank(A) ；秩 ans = 4 norm(A) ；范数 ans = 21.3005 cond(A) ；条件数 ans = 11.1739 trace(A) ；迹 ans = 16 （2）略 6. A=[1 1 0.5;1 1 0.25;0.5 0.25 2] A = 1.0000 1.0000 0.5000 1.0000 1.0000 0.2500 0.5000 0.2500 2.0000 [V,D]=eig(A) V = 0.7212 0.4443 0.5315 -0.6863 0.5621 0.4615 -0.0937 -0.6976 0.7103 D = -0.0166 0 0 0 1.4801 0 0 0 2.5365 MATLAB程序设计 用MATLAB语言编写的程序，称为M文件。M文件可以根据调用方式的不同分为两类：命令文件(Script File)和函数文件(Function File)。 例3-1 分别建立命令文件和函数文件，将华氏温度f转换为摄氏温度c。 程序1：首先建立命令文件并以文件名f2c.m存盘。 clear; %清除工作空间中的变量 f=input('Input Fahrenheit temperature：'); c=5*(f-32)/9 然后在MATLAB的命令窗口中输入f2c，将会执行该命令文件，执行情况为： Input Fahrenheit temperature：73 c =22.7778 程序2：首先建立函数文件f2c.m。 function c=f2c(f) c=5*(f-32)/9 然后在MATLAB的命令窗口调用该函数文件。 clear; y=input('Input Fahrenheit temperature：'); x=f2c(y) 输出情况为： Input Fahrenheit temperature：70 c =21.1111 x =21.1111 3.1.2 M文件的建立与打开 M文件是一个文本文件，它可以用任何编辑程序来建立和编辑，而一般常用且最为方便的是使用MATLAB提供的文本编辑器。 1．建立新的M文件 为建立新的M文件，启动MATLAB文本编辑器有3种方法： (1) 菜单操作。从MATLAB主窗口的File菜单中选择New菜单项，再选择M-file命令，屏幕上将出现MATLAB 文本编辑器窗口。 (2) 命令操作。在MATLAB命令窗口输入命令edit，启动MATLAB文本编辑器后，输入M文件的内容并存盘。 (3) 命令按钮操作。单击MATLAB主窗口工具栏上的New M-File命令按钮，启动MATLAB文本编辑器后，输入M文件的内容并存盘。 2．打开已有的M文件 打开已有的M文件，也有3种方法： (1) 菜单操作。从MATLAB主窗口的File菜单中选择Open命令，则屏幕出现Open对话框，在Open对话框中选中所需打开的M文件。在文档窗口可以对打开的M文件进行编辑修改，编辑完成后，将M文件存盘。 (2) 命令操作。在MATLAB命令窗口输入命令：edit 文件名，则打开指定的M文件。 (3) 命令按钮操作。单击MATLAB主窗口工具栏上的Open File命令按钮，再从弹出的对话框中选择所需打开的M文件。 3.2 程序控制结构 3.2.1 顺序结构 1．数据的输入 从键盘输入数据，则可以使用input函数来进行，该函数的调用格式为： A=input(提示信息，选项)； 其中提示信息为一个字符串，用于提示用户输入什么样的数据。 如果在input函数调用时采用's'选项，则允许用户输入一个字符串。例如，想输入一个人的姓名，可采用命令： xm=input('What''s your name?','s'); 2．数据的输出 MATLAB提供的命令窗口输出函数主要有disp函数，其调用格式为 matlab练习题和答案 控制系统仿真实验 Matlab部分实验结果 目录 实验一 MATLAB基本操 作 ............................................................................................ 1 实验二 Matlab编程 .................................................................................................... 5 实验三Matlab底层图形控制 (6) 实验四控制系统古典分析.............................................................................................12 实验五控制系统现代分析 . (15) 实验六 PID控制器的设 计 ...........................................................................................19 实验七系统状态空间设计.............................................................................................23 实验九直流双闭环调速系统仿真 . (25) 实验一 MATLAB基本操作 1 用MATLAB可以识别的格式输入下面两个矩阵 1233，，1443678，i，，,,2357,,,,2335542，i,,,, A,1357B,,,2675342， i,,3239,,,,189543，，,,1894，， 再求出它们的乘积矩阵C，并将C矩阵的右下角2×3子矩阵赋给D矩阵。赋值完成后，调 用相应的命令查看MATLAB工作空间的占用情况。 A=[1,2,3,3;2,3,5,7;1,3,5,7;3,2,3,9;1,8,9,4]; 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+a 2（ 1） n=input(' 请输入成绩 '); switch n case num2cell(90:100) 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'; elseif n>=80&n<=89 p='B'; elseif n>=70&n<=79 p='C'; elseif n>=60&n<=69 p='D'; else p='E'; end price=p （3） try n; catch price='erroe' end 3 n=[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:20 if n(m)>a a=n(m); elseif n(m)=0 disp(A([n],:)); elseif n<0 disp(lasterr); else disp(A([6],:)); disp(lasterr); end 7（ 1） f=[]; 矩量法m ａtla ｂ程序设计实例： Ha ｌlen 方程求对称振子天线 一、条件与计算目标 已知: 对称振子天线长为Ｌ,半径为a ，且天线长度与波长得关系为,，设,半径a=0、０００0001，因此波数为。 目标： 用H ａll ｅn 方程算出半波振子、全波振子以及不同值得对应参数值。 求:（1）电流分布 （2）E 面方向图 (二维）,H 面方向图（二维）,半波振子空间方向性图(三维） 二、对称振子放置图 图1 半波振子得电流 分布 半波振子天线平行于z 轴放置，在ｘ轴与ｙ轴上得分量都为零，坐标选取方式有两种形式，一般选取图1得空间放置方 式。图１给出了天线得电流分布情况，由图可知，当天线很细时,电流分布近似正弦分布。 三、Ha ｌｌen 方程 得解题思路 ()()()()2 1 ' ' ' ' 12,cos sin sin 'z z i z z z z i z k z G z z dz c kz c kz E k z z dz j ωμ'++=-?? 对于中心馈电得偶极子，Hallen 方程为 ()22'1222 ('),'cos sin sin ,2L L i L L V i z G z z dz c kz c kz k z z j η + -- ++= <<+? 脉冲函数展开与点选配，得到 ()1121 ,''cos sin sin ,1,2,,2n n N z i n m m m m z n V I G z z dz c kz c kz k z m N j η +''=++= =???∑? 上式可以写成 矩阵形式为 四、结果与分析 （1）电流分布 习题2.1 画出下列常见曲线的图形 y （1）立方抛物线3x 命令：syms x y; ezplot('x.^(1/3)') (2)高斯曲线y=e^(-X^2); 命令：clear syms x y; ezplot('exp(-x*x)') (3)笛卡尔曲线 命令：>> clear >> syms x y; >> a=1; >> ezplot(x^3+y^3-3*a*x*y) (4)蔓叶线 命令：>> clear >> syms x y; >> a=1 ezplot(y^2-(x^3)/(a-x)) (5)摆线:()()t sin- = , = - b y 1 t x cos t a 命令：>> clear >> t=0:0.1:2*pi; >> x=t-sin(t); >>y=2*(1-cos(t)); >> plot(x,y) 7螺旋线 命令：>> clear >> t=0:0.1:2*pi; >> x=cos(t); >> y=sin(t); >> z=t; >>plot3(x,y,z) (8)阿基米德螺线 命令：clear >> theta=0:0.1:2*pi; >> rho1=(theta); >> subplot(1,2,1),polar(theta,rho1) (9) 对数螺线 命令：clear theta=0:0.1:2*pi; rho1=exp(theta); subplot(1,2,1),polar(theta,rho1) (12)心形线 命令：>> clear >> theta=0:0.1:2*pi; >> rho1=1+cos(theta); >> subplot(1,2,1),polar(theta,rho1) 练习2.2 1. 求出下列极限值 （1）n n n n 3 lim 3 +∞ → 命令：>>syms n >>limit((n^3+3^n)^(1/n)) ans = 3 （2）)121(lim n n n n ++-+∞ → 命令：>>syms n >>limit((n+2)^(1/2)-2*(n+1)^(1/2)+n^(1/2),n,inf) ans = 0 （3）x x x 2cot lim 0 → 命令：syms x ；MATLAB程序设计教程(第二版)课后答案(可编辑修改word版)
MATLAB程序设计教程课后答案
matlab课后习题解答第二章
matlab课后答案解析完整版
MATLAB程序设计教程(第二版)第三章实验报告下载
matlab课后习题及答案详解
(完整版)matlab习题及答案(2)
matlab基础练习题(带答案)
Matlab程序设计教程(第二版)刘卫国课后参考答案解析
matlab课后习题与答案
matlab课后习题答案
matlab程序设计与应用(第二版)第三章部分课后答案
MATLAB程序设计及经典例题解析3
matlab练习题和答案
MATLAB程序设计教程课后答案.doc
矩量法matlab程序设计实例
matlab课后习题答案 (附图)