MATLAB教程及实训第2版曹弋8
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教程及实训第2版曹弋8
2. 脉冲响应 系统的脉冲响应使用 impulse 函数命令绘制，命 令格式与 step 函数相同。
3. 斜坡响应和加速度响应 斜坡响应 =阶跃响应 *1/s 加速度响应 =阶跃响应 *1/s^2
例8-10 使用step 和impulse 函数绘制系统
第8章 线性控制系统的分析
8.1 控制系统的数学模型 8.2 时域分析的MATLAB 实现 8.3 频域分析的MATLAB 实现 8.4 根轨迹分析的MATLAB 实现 8.5 稳定性分析的MATLAB 实现 8.6 稳态误差分析的MATLAB 实现 8.7状态分析的MATLAB 实现 8.8线性定常系统分析与设计的图形工具 8.9综合实例介绍
2. 时域分析的性能指标
在自动控制原理中，时域分析常用的系统性 能指标有超调量σp 、上升时间tr 、峰值时间tp 和过渡时间ts ，通过性能指标来分析系统暂态 性能的稳定性。
例：
>> detap=exp(-pi*zeta/sqrt(1-zeta^2))*100 %计算超调量
>> tr=(pi-acos(zeta))/(wn*sqrt(1-zeta^2)) %计算上升时间
G=zpk(z,p,k,Ts) % 由零点、极点和增益创建模型
G=zpk(z,p,k,'Property1',v1,'Porperty2',v2,…) %创建模型并设置属性
例：创建 G(s) ?
2(s ? 0.5)
(s ? 0.1 ? j)(s ? 0.1 ? j)
>> z=-0.5;
>> p=[-0.1+j -0.1-j];
(完整word版)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。
MATLAB实用教程第二版教学设计
MATLAB实用教程第二版教学设计介绍MATLAB是一种非常流行的数学软件，广泛用于科学、工程和技术中的计算和数据分析。
许多大学的计算机、工程和科学专业都要求学生掌握MATLAB的基础知识。
然而，对于初学者来说，MATLAB的学习可能会感到有些困难，需要一定的指导。
本教学设计旨在提供一些实用的教程来帮助初学者更好地掌握MATLAB的使用。
教学目标通过本教学设计，学生将能够：•理解MATLAB的基础概念，包括变量、矩阵和函数；•学习如何在MATLAB中进行基本的计算、数据操作和可视化；•学习如何编写MATLAB脚本；•掌握MATLAB编程语言的基本语法和结构；•了解MATLAB在科学和工程领域的应用。
教学内容本教学设计的内容包括以下几个方面：第一部分：MATLAB基础1.介绍MATLAB的基本概念和特点，包括什么是MATLAB，为什么要使用MATLAB，MATLAB的应用领域等。
2.MATLAB编程环境的介绍，包括工作空间、命令窗口、编辑器和变量窗口。
3.MATLAB的基本语法规则，包括变量命名、运算符、数据类型等。
4.数组和矩阵的操作，包括如何创建矩阵、矩阵的运算和元素操作等。
第二部分：MATLAB图形化1.MATLAB图形化环境的介绍，包括图形窗口、图形对象和属性编辑器。
2.基本画图函数的介绍，包括plot、scatter、bar和histogram等。
3.如何修改图形对象的属性，包括颜色、样式、标签和注释等。
4.图像处理的常用函数，包括imread、imshow、imwrite等。
第三部分：MATLAB脚本1.脚本的概念和用途，什么是脚本，如何创建脚本文件等。
2.脚本文件的操作，包括如何编写和运行脚本、调试脚本等。
3.如何使用函数和脚本文件，包括函数调用、函数参数传递等。
4.脚本文件的文件输入输出，包括文件的打开、读取、写入和关闭等。
教学方法本课程将采用以下教学方法：1.课堂教学：讲解MATLAB的基础知识和编程语言的基本语法。
matlab实用教程第二版课程设计
Matlab实用教程第二版课程设计一、简介本文档旨在介绍《Matlab实用教程》第二版课程设计的具体内容和步骤。
本课程设计以Matlab编程语言为基础，通过实践操作和代码编写，提高学生对Matlab 的掌握和应用能力，培养学生解决实际问题的能力。
二、课程设计目标1.掌握Matlab基本语法和编程技巧；2.熟悉Matlab的常用工具箱，并能灵活运用；3.能够使用Matlab解决实际问题，并能进行结果可视化；4.提高学生对数学建模和算法设计的能力。
三、课程设计内容1. 预备知识学习在课程开始前，需要学生自主学习以下预备知识：•Matlab基本操作和编程语法；•常用工具箱（Signal Processing Toolbox、Optimization Toolbox 等）；•常用函数（如fft、solve等）。
2. 课程设计步骤第一步：课程时间和人员安排确定课程设计的时间和参与人员，分配学生的任务和角色。
建议将学生分成3-5人的小组，由组长统筹安排。
第二步：题目选择和讨论老师提供一些现实问题或者学术问题，让学生根据自己感兴趣和所学知识选择题目。
然后学生在小组内进行讨论和思考，确定问题的解决方法和方案。
第三步：代码实现和调试利用Matlab编写代码实现问题的求解和仿真，并适时进行调试和修改。
建议学生在代码注释中写明代码功能和实现方法，方便以后的查看和维护。
第四步：结果可视化和分析将实现的算法结果进行可视化处理，利用Matlab的图形绘制工具进行绘图。
并进行数据分析和结果解释。
建议学生在报告中给出相关的图片和图表，增强报告的可读性。
第五步：课程设计报告撰写将实验过程和结果撰写成课程设计报告。
建议包括以下内容：•题目背景和研究目的；•算法设计和实现过程；•实验结果和分析；•总结和反思。
四、总结本课程设计旨在通过实践操作和代码编写，提高学生对Matlab的掌握和应用能力，培养学生解决实际问题的能力。
希望学生能够认真对待本次课程设计，并将所学的知识应用到实际科研和工作中。
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数学实验第二版课后练习题含答案课后练习题MATLAB数学实验第二版的课后练习题如下：第一章课后练习题1.编写MATLAB程序，计算并输出下列公式的结果：y = \\frac{1}{\\sqrt{2\\pi\\sigma^2}} e^{-\\frac{(x-\\mu)^2}{2\\sigma^2}}其中，x, $\\mu$, $\\sigma$ 分别由用户输入。
要求输出结果精确至小数点后两位。
答案如下：x=input('请输入 x 的值:');mu=input('请输入 mu 的值:');sigma=input('请输入 sigma 的值:');y=1/sqrt(2*pi*sigma^2) *exp(-(x-mu)^2/ (2*sigma^2));fprintf('y = %.2f\', y);2.编写MATLAB程序，求解下列方程的解：4x + y = 11\\\\x + 2y = 7答案如下：A= [4,1;1,2];B= [11;7];X=inv(A) *B;fprintf('x = %.2f, y = %.2f\', X(1), X(2));第二章课后练习题1.编写MATLAB程序，计算下列多项式的值：P(x) = x^4 - 2x^3 + 3x^2 - x + 1其中，x 由用户输入。
要求输出结果精确至小数点后两位。
答案如下：x=input('请输入 x 的值:');y=x^4-2*x^3+3*x^2-x+1;fprintf('P(%.2f) = %.2f\', x, y);2.编写MATLAB程序，绘制下列函数的图像：f(x) = \\begin{cases} x + 1, & x < 0 \\\\ x^2, & 0 \\leq x < 1 \\\\ 2x - 1, & x \\geq 1 \\end{cases}答案如下：x=-2:0.01:2;y1=x+1;y2=x.^2.* ((x>=0) & (x<1));y3=2*x-1;plot(x,y1,x,y2,x,y3);legend('y1 = x + 1','y2 = x^2','y3 = 2x - 1');总结本文提供了《MATLAB数学实验第二版》的部分课后练习题及其答案。
Matlab实训2
最大值和最小值
▪ 求向量的最大（小）值
▪ Y=max(X)、Y=min(X)
11
11
>> A(:,[2 3]) ans =
23 67
>> A(:,[2 3])=Sa
A= 1114 5118
第八页，共77页。
2.2 矩阵（数组）操作函数
用于矩阵（数组）操作的常用函数
函数 size length ndims numel disp cat reshape repmat fliplr flipud flipdim find
▪常用于寻找数组中所有大于某值的元素的问题
例子：找出数组A= -4 –2 0 2 4 中所有绝对值大于3的元素
-3 –5 1 3 5
>> A=[-4 -2 0 2 4;-3 -5 1 3 5]; %生成数组A
>> L=abs(A)
L=
10001
01001
>> islogical(L);
s= 235
>> A(s) ans =
523
>> Sa=[10 20 30] Sa =
10 20 30 >> A(s)=Sa
A= 1 20 30 4 10 6 7 8
>> Sa=[10 20 30]' Sa =
10 20 30 >> A(s)=Sa A=
1 20 30 4
10 6 7 8
第七页，共77页。
第二页，共77页。
A(1:4,5)
A(:,5) A(:,end) A(17:20) A(17:20)'
矩阵元素的访问与赋值（续）
最新Matlab数学实验第二版课后习题答案
Matlab数学实验第二版课后习题答案胡良剑第一章%Page20,ex1(5) 等于[exp(1),exp(2);exp(3),exp(4)](7) 3=1*3, 8=2*4(8) a为各列最小值，b为最小值所在的行号(10) 1>=4,false, 2>=3,false, 3>=2, ture, 4>=1,ture(11) 答案表明：编址第2元素满足不等式(30>=20)和编址第4元素满足不等式(40>=10)(12) 答案表明：编址第2行第1列元素满足不等式(30>=20)和编址第2行第2列元素满足不等式(40>=10)%Page20, ex2(1)a, b, c的值尽管都是1，但数据类型分别为数值，字符，逻辑，注意a与c相等，但他们不等于b(2)double(fun)输出的分别是字符a,b,s,(,x,)的ASCII码%Page20,ex3>> r=2;p=0.5;n=12;>> T=log(r)/n/log(1+0.01*p)T =11.5813%Page20,ex4>> x=-2:0.05:2;f=x.^4-2.^x;>> [fmin,min_index]=min(f)fmin =-1.3907 %最小值min_index =54 %最小值点编址>> x(min_index)ans =0.6500 %最小值点>> [f1,x1_index]=min(abs(f)) %求近似根--绝对值最小的点f1 =0.0328x1_index =24>> x(x1_index)ans =-0.8500>> x(x1_index)=[];f=x.^4-2.^x; %删去绝对值最小的点以求函数绝对值次小的点>> [f2,x2_index]=min(abs(f)) %求另一近似根--函数绝对值次小的点f2 =0.0630x2_index =65>> x(x2_index)ans =1.2500%Page20,ex5>> z=magic(10)z =92 99 1 8 15 67 74 51 58 4098 80 7 14 16 73 55 57 64 414 81 88 20 22 54 56 63 70 4785 87 19 21 3 60 62 69 71 2886 93 25 2 9 61 68 75 52 3417 24 76 83 90 42 49 26 33 6523 5 82 89 91 48 30 32 39 6679 6 13 95 97 29 31 38 45 7210 12 94 96 78 35 37 44 46 5311 18 100 77 84 36 43 50 27 59>> sum(z)ans =505 505 505 505 505 505 505 505 505 505>> sum(diag(z))ans =505>> z(:,2)/sqrt(3)ans =57.157746.188046.765450.229553.693613.85642.88683.46416.928210.3923>> z(8,:)=z(8,:)+z(3,:)z =92 99 1 8 15 67 74 51 58 4098 80 7 14 16 73 55 57 64 414 81 88 20 22 54 56 63 70 4785 87 19 21 3 60 62 69 71 2886 93 25 2 9 61 68 75 52 3417 24 76 83 90 42 49 26 33 6523 5 82 89 91 48 30 32 39 6683 87 101 115 119 83 87 101 115 11910 12 94 96 78 35 37 44 46 5311 18 100 77 84 36 43 50 27 59 1小时单位负责人接到报告后，应当于1小时内向事故发生地县级以上人民政府安全生产监督管理部门和负有安全生产监督管理职责的有关部门报告。
数字信号处理(第2版)教学课件第8章 MATLAB仿真实验
系统的稳定性是指对任意有界的输入信号，系统都能 得到有界的系统响应，或者系统的单位脉冲响应满足绝对 可和的条件。系统的稳定性由其差分方程的系数决定。
实际中检查系统是否稳定，不可能检查系统对所有有界 的输入信号、输出是否都是有界输出，或者检查系统的单位 脉冲响应满足绝对可和的条件。可行的方法是在系统的输入 端加入单位阶跃序列，如果系统的输出趋近一个常数（包括 零），就可以断定系统是稳定的。系统的稳态输出是指当n→ ∞时系统的输出。如果系统稳定，信号加入系统后，系统输出 的开始一段称为暂态效应，随着n的加大，幅度趋于稳定，达 到稳态输出。
（2）频域采样理论的验证。 给定长度为26的三角波序列x(n) 编写程序，分别对频谱函数 X (e j ) FT[x(n)] 在区间 [0, 2π] 上等间隔采样32点和16点，得到 X32 (k) 和 X16 (k) ，再分别对 X32 (k)
和 ①X16分(k)别进画行出32X点(ej和 )、16X点32 (IkF)F和T，X1得6 (k到) 的x幅32 (度n)谱和。x16 (n) 。要求：
4. 参考程序
（1）内容1参考程序，实验结果。 （2）内容2参考程序，实验结果。 （3）内容3参考程序，实验结果。
5.实验结果
图8-1 调用filter解差分方程仿真结果
5.实验结果
图8-2 稳定性分析方面的仿真结果
5.实验结果
图8-3 稳定性分析仿真结果
实验二 时域采样与频域采样
1. 实验目的
y(n) 0.5y(n 1) 0.25y(n 2) x(n) 2x(n 1) x(n 3)
MATLAB基础与应用教程(第2版)
思考与实验
3.1脚本 3.2程序控制结构
3.3函数 3.4程序调试
3.5程序性能分 析与优化
思考与实验
1
4.1二维曲线 绘制
2
4.2其他二维 图形绘制
3
4.3三维图形 绘制
4
4.4图形修饰 处理
5
4.5图像处理
4.6交互式绘图 工具
思考与实验
1
5.1特殊矩阵 的生成
2
5.2矩阵分析
3
5.3矩阵分解
5
思考与实验
11.1 Simulink概 述
11.2系统模型的构 建
11.3系统的仿真与 分析
11.4子系统
11.5 S函数的 设计与应用
思考与实验
12.2文件操作
12.1 MATLAB与 Excel的接口
12.3 MATLAB编译 器
12.4 MATLAB 与其他语言程 序的应用接口
思考与实验
13.1 MATLAB在优化 问题中的应用
13.2 MATLAB在控制 系统中的应用
13.3 MATLAB在信号 处理中的应用
13.4 MATLAB在数字 图像处理中的应用
13.5 MATLAB 在经济和金融 领域中的应用
思考与实验
作者介绍
这是《MATLAB基础与应用教程（第2版）》的读书笔记模板，暂无该书作者的介绍。
目录分析
1.2 MATLAB的工作 环境
1.1 MATLAB的发展 与基本功能
1.3 MATLAB的基本 操作
1.4 MATLAB的 帮助系统
思考与实验
2.1数值数据及操作 2.2变量及其操作
2.3 MATLAB数组 2.4 MATLAB运算
MATLAB程序设计及应用(第二版)课后实验答案(最新整理)
阵和对角阵，试通过数值计算验证。
22E R RS A OS +⎡⎤=⎢⎥⎣⎦解: M 文件如下；由ans,所以22E R RS A O S +⎡⎤=⎢⎥⎣⎦2. 产生5阶希尔伯特矩阵H 和5阶帕斯卡矩阵P ，且求其行列式的值Hh 和Hp 以及它们的条件数Th 和Tp ，判断哪个矩阵性能更好。
为什么？解：M 文件如下：因为它们的条件数Th>>Tp,所以pascal 矩阵性能更好。
3. 建立一个5×5矩阵，求它的行列式值、迹、秩和范数。
解： M 文件如下：4. 已知2961820512885A -⎡⎤⎢⎥=⎢⎥⎢⎥-⎣⎦求A 的特征值及特征向量，并分析其数学意义。
解：M 文件如图：数学意义：V的3个列向量是A的特征向量，D的主对角线上3个是A的特征值，特别的，的3个列向量分别是D的3个特征值的特征向量。
5. 下面是一个线性方程组：111⎡⎤输出结果：由结果，X和X2的值一样，这表示b的微小变化对方程解也影响较小，而A的条件数算得较小，所以数值稳定性较好，A是较好的矩阵。
6. 建立A矩阵，试比较sqrtm(A)和sqrt(A)，分析它们的区别。
解：M文件如下：分析结果知：sqrtm(A)是类似A的数值平方根（这可由b1*b1=A的结果看出），而sqrt(A)则是对A中的每个元素开根号，两则区别就在于此。
实验三 选择结构程序设计1. 求分段函数的值。
2226035605231x x x x y x x x x x x x ⎧+-<≠-⎪=-+≤<≠≠⎨⎪--⎩且且及其他用if 语句实现，分别输出x=-5.0,-3.0,1.0,2.0,2.5,3.0,5.0时的y 值。
解：M文件如下：2. 输入一个百分制成绩，要求输出成绩等级A 、B 、C 、D 、E 。
其中90分~100分为A ，80分~89分为B ，79分~79分为C ，60分~69分为D ，60分以下为E 。
要求：(1) 分别用if 语句和switch 语句实现。
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)。
0.5z z2 ? 1.5z ? 0.5
2. 零极点增益模型
G(s)
?
K
(s ? z1 )(s ? (s ? p1 )(s ?
z2 )? (s ? p2 )? (s ?
zm ) pn )
G(z) ?
K
(z ? z1 )(z ? (z ? p1 )(z ?
z2 )? (z ? zm ) p2 )? (z ? pn )
>> [num,den]=zp2tf(z,p,k);
>> G11=tf(num,den) % 转换为传递函数模型
8.1.2 系统的模型参数
2. 获取模型的参数
（1）获取模型参数的函数
tfdata 、zpkdata 和ssdata 、dssdata 函数分别用来
获取传递函数模型、零极点增益模型和状态方程 模型的参数，这些函数名都是在创建模型的函数 名后面加“ data”。
（2）get 函数
value=get(sys,'property') 统的属性
%获取当前系
8.1.3系统模型的连接和简化
1. 串联环节 G=series(G1,G2,outputs1,inputs1) 也可以直接使用G=G1*G2 u(t)
G1(s)
y(t) G2(s)
2. 并联环节
G=parallel(G1,G2,in1,in2,out1,out2)
G=zpk(z,p,k,Ts) % 由零点、极点和增益创建模型
G=zpk(z,p,k,'Property1',v1,'Porperty2',v2,…) %创建模型并设置属性
例：创建 G(s) ?
2(s ? 0.5)
(s ? 0.1 ? j)(s ? 0.1 ? j)
>> z=-0.5;
>> p=[-0.1+j -0.1-j];
第8章 线性控制系统的分析
8.1 控制系统的数学模型 8.2 时域分析的MATLAB 实现 8.3 频域分析的MATLAB 实现 8.4 根轨迹分析的MATLAB 实现 8.5 稳定性分析的MATLAB 实现 8.6 稳态误差分析的MATLAB 实现 8.7状态分析的MATLAB 实现 8.8线性定常系统分析与设计的图形工具 8.9综合实例介绍
例8-3 创建连续系统的零极点增益模型，并转换为传 递函数和状态空间模型，零极点增益模型为
G(s) ?
2(s ? 0.5)
(s ? 0.1 ? j)(s ? 0.1 ? j)
>> z=-0.5;
>> p=[-0.1+j -0.1-j];
>> k=2;
>> G=zpk(z,p,k) %创建zpk 模型
>> k=2;
>> G=zpk(z,p,k)
%创建zpk 模型3. 状方程模型? x?? Ax ? Bu
? ?
y
?
Cx
?
Du
? x(k ? 1) ? Ax(k) ? Bu(k)
? ?
y(k
?
1)
?
Cx(k) ?
Du(k)
G=ss(a,b,c,d,Ts) %由a、b、c、d参数创 建模型
G=ss(a,b,c,d,Ts,'Property1',v1,'Porperty2
',v2,…)
%创建模型并设置属性
4. 模型的转换
ss2tf
传递函数模型 (tf)
tf2ss
zp2tf
tf2zp
状态空间模型 (ss)
zp2ss ss2zp
零极点增益模型 (zpk)
5. 连续系统与离散系统模型的转换
MATLAB 控制工具箱提供了 c2d 、d2c 和d2d 函数 实现连续系统和离散系统的相互转换。
（4）列出系统总的输入和输出端的编号，使用 inputs 列出 输入端编号，outputs 列出输出端的编号；
也可以直接使用：
G1(s)
G=G1+G2
u(t)
y(t)
G2(s)
3. 反馈环节
G=feedback(G1,G2,feedin,feedout,sign) %
计算反馈模型
u(t)
y(t) G1(s)
±
G2(s)
说明：G1和G2模型必须都是连续系统或采样 周期相同的离散系统；sign 表示反馈符号， 当sign 省略或=-1时为负反馈；feedin 和 feedout 分别是G2的输入端口和G1的输出 端口，可省略，当省略时G1与G2端口正好 对应连接。
8.1 控制系统的数学模型 8.1.1 创建系统的模型并相互转换
1. 传递函数模型
G(s) ?
b1s m ? b2 s m?1 ? ? ? bm s ? bm?1 s n ? a1s n?1 ? ? ? an?1s ? an
G(z) ?
b1z m ? b2z m?1 ? ? ? bm z ? bm?1 z n ? a1z n?1 ? ? ? an?1z ? an
sys=tf(num,den,Ts) % 由分子分母得出传递函数
sys=tf(num,den,Ts,‘Property1',v1,‘Porperty2',v2,
…)
%创建传递函数并设置属性
例8-1 创建连续二阶系统，已知传递函数模
型为
5
G(s) ?
s2 ? 2s ? 2
>> num1=5; >> den1=[1 2 2]; >> sys1=tf(num1,den1) % 创建传递函数
G(s) ?
5
e ? 2s
s2 ? 2s ? 2
>> sys1t=tf(num1,den1,'inputdelay',2) % 创建带延迟环节的传递函数
>> num2=[0.5 0];
>> den2=[1 -1.5 0.5]; >> sys2=tf(num2,den2,-1)
%创建脉冲传递函数
G(z) ?
（2）获取模型尺寸的函数
n=ndims(sys)
%获取模型的维数
d=size(sys,n)
%获取模型的参数
d=size(sys,'order') %获取模型的阶数
3. 使用get和set函数
（1）set 函数
set(sys,'property1',value1,'property2',val ue2,…) %设置系统属性
4. 复杂模型的连接
（1）对框图中的每个环节进行编号并建立它们的对象模型， 环节是指一条单独的通路；
（2）建立无连接的状态空间模型，使用append 函数实现， append 的命令格式如下：
G=append(G1,G2,G3,…)
（3）写出系统的联接矩阵Q
Q是联接矩阵，第一列是各环节的编号，其后各列是与该环 节连接的输入通路编号，如果是负连接则加负号；