《精》控制系统计算机第二章作业.docVIP

第二章作业 探控0901 范项媛 4201090119 第1题： 程序： A=[1,2,3,4;4,3,2,1;2,3,4,1;3,2,4,1] 运行结果： 程序： B=[1+4j,2+3j,3+2j,4+1j;4+1j,3+2j,2+3j,1+4j;2+3j,3+2j,4+1j,1+4j;3+2j,2+3j,4+1j,1+4j] 运行结果： 程序： A(5,6)=0 运行结果： 第18题： [x,y]=meshgrid(-1:0.01: 1,-1:0.01: 1); %生成网格数据 z=sin(x.*y); %描绘出函数z=sin（x*y） subplot(1,2,1); %制定原图性的位置 mesh(x,y,z) %绘制三维网格图 z(find(x.^2+y.^2=0.5^2))=NaN; %剪切x2+y2=0.52部分 subplot(1,2,2);mesh(x,y,z) %制定剪切后图形的位置，并画出三维网格图 第23题. 步骤如下： （1）先描绘出目标函数 function dx=Lorenz(t,x) %虽不显含时间，还应写出占位 y=[-8*x(1)/3+x(2)*x(3);-10*x(2)+10*x(3);-x(1)*x(2)+28*x(2)-x(3)]; %对应方程的三行，直接将参数代入 （2）再求微分方程的数值解 x0=[0,0,10^-3]; %微分方程的初值 [t,y]=ode45(‘Lorenz’,[0,100],x0); %求解微分方程 plot(t,y); %绘制各个状态变量的时间响应 figure; plot3(y(:,1),y(:,2),y(:,3));%绘制相空间图形 plot(y(:,1),y(:,2));%绘制投影 plot(y(:,2),y(:,3)); plot(y(:,1),y(:,3)); grid， %绘制相空间图形 第26题： 步骤如下： （1）先建一个M文件，描述出目标函数和约束函数； 程序如下： function [c,ce]=c2f1a(x) ce=[]; c=[x(1)^2+x(2)^2-9;x(1)+x(2)-1]; （2）再调用非线性最优化问题求解函数可以得出一下结果 x0=[0;0];A=[1 1];B=1; x=fmincon(x(1)^2+x(2),x0,A,B,[],[],[],[],@c2f1a) 结果如下： 用图解法验证得出的解：已知限制x1和x2的条件是在一个圆心为（0,0），半径为3的圆内，用直线x1+x2≤1去截这个圆，则剩下的部分是直线下方与圆所包围的地方。观察图形可知，因为x1站的份额比x2大，所以当x1=0，x2=3时， 可使 表达式 （-（x1）2-x2）的值最小。 图形如下：