控制系统仿真.doc

《控制系统仿真.doc》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、__________________________________________________《控制系统计算机仿真—课程设计作业》姓名：专业：学号：1.构造矩阵（1）试将与的结果相加，并找出相加后新矩阵中绝对值大于10的元素。（2）组合成一个4´3的矩阵，第一列为按列顺序排列的A矩阵元素，第二列为按列顺序排列的B矩阵元素，第三列为按列顺序排列的C矩阵元素。（1）matlab：A=[-71;8-3];B=[42;57];C=[59;62];D=A.*B+A*Bnum=find(abs(D)>10);D(num)结果：D=-51

2、-557-26ans=-5157-26（2）matlab：NEW=[A(1,:),A(2,:);B(1,:),B(2,:);C(1,:),C(2,:)]’结果：NEW=-745129856____________________________________________________________________________________________________-3722.绘制函数曲线，要求写出程序代码(1)在区间均匀的取50个点，构成向量X。(2)在同一窗口绘制曲线y1=sin(2*t-0.3); y2=

3、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-');holdonplot(t,y1,'o',t,y2,'*');holdoff3.写出生成下图所示波形的MATLAB程序。图中三个波形均为余弦波，x范围为[pi/2~7*pi/2]。要求它的正半波被置零；且在和处被削顶。程序：x=linspace(pi/2,7*pi/2,10

4、0);y1=cos(x);figure（1）plot（x,y1）y1(find(y1>0))=0;figure(2)plot(x,y1)n=find((x>2*pi/3&x<4*pi/3)

5、(x>8*pi/3&x<10*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①用MATLAB语言分四个区域分别绘制的曲线，并且对图形标题及横纵坐标轴进行标注。②另建一个窗口，不分区，用不同颜色、线型绘出四条曲线

6、，并标注图例注解。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('ºá×ø±êx=[-2*pi,2*pi]');ylabel('×Ý×ø±êy=

7、sin(x)')title('y=sin(x)ÔÚ[-2*pi,2*pi]Ö®¼äµÄÇúÏßͼ')subplot(2,2,2)plot(x,y2)xlabel('ºá×ø±êx=[-2*pi,2*pi]');ylabel('×Ý×ø±êy=cos(x)')title('y=cos(x)ÔÚ[-2*pi,2*pi]Ö®¼äµÄÇúÏßͼ')subplot(2,2,3);plot(x,y3);xlabel('ºá×ø±êx=[-2*pi,2*pi]');ylabel('×Ý×ø±êy=sin(2x)')title('y=sin(

8、2x)ÔÚ[-2*pi,2*pi]Ö®¼äµÄÇúÏßͼ')subplot(2,2,4);plot(x,y4);xlabel('ºá×ø±êx=[-2*pi,2*pi]');ylabel('×Ý×ø±êy=cos(2x)')title('y=cos(2x)ÔÚ[-2*pi,2*pi]Ö®¼äµÄÇúÏßͼ')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.请

9、分析并修改下面的程序，使用矩阵或数组运算的顺序结构完成嵌套语句的相同功能。A=[123;456;789];[rc]=size(A);fori=1:1:rforj=1:1:cif(A(i,j)>8

10、A(i,j)<2)A(i,j)=0;endenden