《控制系统仿真》.docxVIP

《控制系统计算机仿真一课程设计作业》 姓名： 专业： 学号: 1.构造矩阵 一-7 1 ] A = ,B = f 21 ,C』5 9〕 !8 ? 7 - !6 2 一 试将A.* B与A* B的结果相加，并找出相加后新矩阵中绝对值大于 10的元素。 组合成一个4x3的矩阵，第一列为按列顺序排列的A矩阵元素，第二列为按列顺序排列的B矩 阵元素，第三列为按列顺序排列的 C矩阵元素。 (1) matlab： A=[-7 1;8 -3]; B=[4 2;5 7]; C=[5 9;6 2]; D=A.*B+A*B num=find(abs(D)10); D(num) 结果：D = -51 -5 57 -26 ans = -51 57 -26 (2) matlab: NEW=[A(1,:),A(2,:);B(1,:),B(2,:);C(1,:),C(2,:)] 结果：NEW = -7 4 5 1 2 9 8 5 6 -3 7 2 绘制函数曲线，要求写出程序代码 在区间[0,2兀]均匀的取50个点，构成向量 X。 在同一窗口绘制曲线 y1=sin(2*t-0.3); y2=3cos(t+0.5);要求y1曲线为红色点划线, 标记点为圆圈；y2为蓝色虚线，标记点为星号。 程序代码： t=linspace(0,2*pi,50); y1=sin(2*t-0.3); y2=3*cos(t+0.5); plot(t,y1,r.,t,y2,b-); hold on plot(t,y1,o,t,y2,*); hold off 写出生成下图所示波形的 MATLAB程序。图中三个波形均为余弦波， x范围为[pi/2 ~ y1=cos(x); figure (1) plot (x,y1 ) y1(find(y10))=0; figure(2) plot(x,y1) n=find((x2*pi/3x4*pi/3)|(x8*pi/3x10*pi/3)); y1(n)=cos(2*pi/3); figure(3) plot(x,y1) 4 对于 x= [— 2兀，2 兀],y1=sinx、y2=cosx、y3=sin2x、y4=cos2x 用MATLA卵言分四个区域分别绘制的曲线，并且对图形标题及横纵坐标轴进行标注。 另建一个窗口，不分区，用不同颜色、线型绘出四条曲线，并标注图例注解。 Matlab : x=linspace(-2*pi,2*pi,1000); y1=sin(x); y2=cos(x); y3=sin(2*x); y4=cos(2*x); figure(1) subplot(2,2,1) plot(x,y1) xlabel( oa x ? + x=[-2*pi,2*pi]); ylabel( x yx ?±y=sin(x)) title(y=sin(x) ? u[-2*pi,2*pi] ???? ??u ?? i ?) subplot(2,2,2) plot(x,y2) xlabel( oa x ?±x=[-2*pi,2*pi]); ylabel( x yx ?±y=cos(x)) title(y=cos(x) ?u [-2*pi,2*pi] ????卬??U ?? 1 ?) subplot(2,2,3); plot(x,y3); xlabel( oa x ?±x=[-2*pi,2*pi]); ylabel( x yx ?±y=sin(2x)) title(y=sin(2x) ?u [-2*pi,2*pi] ????卬??U ?? 1 ?) subplot(2,2,4); plot(x,y4); xlabel( oa x ?±x=[-2*pi,2*pi]); ylabel( x yx ?±y=cos(2x)) title(y=cos(2x) ?u [-2*pi,2*pi] ????卬??U ?? 1 ?) figure(2) plot(x,y1,r-,x,y2,b.,x,y3,ko,x,y4,g*) legend(y=sin(x),y=cos(x),y=sin(2x ￡ ?,y=cos(2x)) 5.请分析并修改下面的程序，使用矩阵或数组运算的顺序结构完成嵌套 能。 A=[1 2 3; 4 5 6; 7 8 9]; [r c]=size(A); for i=1:1:r for j=1:1:c if (A(i,j)8 | A(i,j)2) A(i,j)=0; end end end 修改后的matlab : A=[1 2 3; 4 5 6; 7 8 9]; A(find(A8|A2))=0; A for语句的相同功 6.假设 a 是这样一组数组：a(n)=a(n—1)2 —a(n—2)2+a(n—3)2 ,且 a( 1 =) a , =( 2 @