资源描述:
《MATLAB仿真课后习题.doc》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、.第一章习题3.请指出以下的变量名(函数名、M文件名)中,哪些是合法的?Abc2004xlil-1wu_2004a&bqst.u_xyz解:合法的变量名有:Abcwu_20044.指令窗操作(1)求[12+2×(7-4)]÷32的运算结果解:>>[12+2*(7-4)]/3^2ans=2(2)输入矩阵A=[1,2,3;4,5,6;7,8,9],观察输出。解:>>A=[1,2,3;4,5,6;7,8,9]A=123456789(3)输入以下指令,观察运算结果;clear;x=-8:0.5:8;y=x';X=ones(size(y
2、))*x;Y=y*ones(size(x));R=sqrt(X.^2+Y.^2)+eps;Z=sin(R)./R;mesh(X,Y,Z);colormap(hot)xlabel('x'),ylabel('y'),zlabel('z')解:Word文档.7.指令行编辑(1)依次键入以下字符并运行:y1=2*sin(0.3*pi)/(1+sqrt(5))解:>>y1=2*sin(0.3*pi)/(1+sqrt(5))y1=0.5000(2)通过反复按键盘的箭头键,实现指令回调和编辑,进行新的计算;y2=2*cos(0.3*pi)/
3、(1+sqrt(5))解:>>y2=2*cos(0.3*pi)/(1+sqrt(5))y2=0.363311.编写题4中(3)的M脚本文件,并运行之。解:Word文档.第二章习题1.在指令窗中键入x=1:0.2:2和y=2:0.2:1,观察所生成的数组。解:>>x=1:0.2:2x=1.00001.20001.40001.60001.80002.0000>>y=2:0.2:1y=Emptymatrix:1-by-02.要求在[0,2π]上产生50个等距采样数据的一维数组,试用两种不同的指令实现。解:y1=0:2*pi/49:2
4、*piy2=linspace(0,2*pi,50)3.计算e-2tsint,其中t为[0,2π]上生成的10个等距采样的数组。解:>>t=linspace(0,2*pi,10);x=exp(-2*t).*sin(t)x=00.15910.06030.01310.0013-0.0003-0.0002-0.0001-0.0000-0.00004.已知A=,B=,计算矩阵A、B乘积和点乘.解:>>A=[1,2;3,4];B=[5,6;7,8];x=A*Bx=19224350>>x=A.*Bx=51221325.已知A=,B=,计算A
5、&B,A
6、B,~A,A==B,A>B.解:>>A=[0,2,3,4;1,3,5,0];B=[1,0,5,3;1,5,0,5];a1=A&BWord文档.a2=A
7、Ba3=~Aa4=(A==B)a5=(A>B)a1=00111100a2=11111111a3=10000001a4=00001000a5=010100107.将题5中的A阵用串转换函数转换为串B,再size指令查看A、B的结构,有何不同?解:>>A=[0,2,3,4;1,3,5,0]B=num2str(A)size(A)size(B)A=02341350B=0234
8、1350ans=24ans=210Word文档.第三章习题1.已知系统的响应函数为,其中,要求用不同线型或颜色,在同一张图上绘制ε取值分别为0.2、0.4、0.6、0.8时,系统在t∈[0,18]区间内的响应曲线,并要求用ε=0.2和ε=0.8对他们相应的两条曲线进行文字标志。解:clccloseallclearallt=0:0.02:18;xi=[0.2,0.4,0.6,0.8]';sxi=sqrt(1-xi.^2);sita=atan(sxi./xi);y=1-exp(-xi*t).*sin(sxi*t+sita*ones
9、(1,901))./(sxi*ones(1,901))plot(t,y(1),'r-',t,y(2),'b*',t,y(3),'g+',t,y(4),'k.')text(4.2,1.4,'xi=0.2')text(3.8,0.9,'xi=0.8')Word文档.2.用plot3、mesh、surf指令绘制三维图(x,y范围自定)。解:clc;closeall;clearall;x=-5:0.1:5;y=-5:0.1:5;[X,Y]=meshgrid(x,y);a=sqrt((1-X).^2+Y.^2);b=sqrt((1+
10、X).^2+Y.^2);Z=1./(a+b);a1=sqrt((1-x).^2+y.^2);b1=sqrt((1+x).^2+y.^2);z=1./(a1+b1);subplot(1,3,1),plot3(x,y,z),xlabel('x'),ylabel('y'),zlab