[机电系统仿真技术试题.docVIP

机电系统仿真技术试题 姓名：刘丽欢 学号：080202061009 专业：机械电子工程 电话 一、用MATLB指令求解下列各题 求的特征值和特征向量。 A=[1,6;3,5]; [V,D]=eig(A) V = -0.9125 -0.6676 0.4092 -0.7445 D = -1.6904 0 0 7.6904 求的特征多项式。 B=[2,5;7,3] B = 2 5 7 3 C=poly(B) C = -5.0000 -29.0000 该矩阵的特征多项式为 ，求其约旦标准形。 E=[1 1 2;0 1 3;0 0 2]; C=jordan(E) C = 2 0 0 0 1 1 0 0 1 E的约旦标准形是 求递推方程的通解。 f=maple(rsolve(f(n)=-3*f(n-1)-4*f(n-2),f(n))); disp(f(n)=),disp(f) function work22(n,i,j) syms a b A=ones(n) for a=1:n for b=1:n if (a+b)==(i+j) A(a,b)=0; end end end disp(A) f(n)= (3/7*i*f(1)*7^(1/2)+f(1)+1/7*i*f(0)*7^(1/2)+3*f(0))*(-8/(3-i*7^(1/2)))^n/(3-i*7^(1/2))-1/7*i*7^(1/2)*(3*f(1)+i*f(1)*7^(1/2)+f(0)+3*i*f(0)*7^(1/2))*(-8/(3+i*7^(1/2)))^n/(3+i*7^(1/2)) 求定积分 syms x; f=exp(-x^2/2); a=int(f,-15,15); an=a/(2*pi)^0.5 an = 2251799813685248/5644425081792261*erf(15/2*2^(1/2))*2^(1/2)*pi^(1/2) 完成下列各个程序段 1.求1000以内的质数。 for n=1:1000 a=primes(n); end a %即可列出所有1000内的质数阶的全1阵，使其（i，j）位置的反斜行变为0。 function c=ccc(n,i,j) if (n=1)||(i1)||(in)||(j1)||(jn) error(Input must be n1 1=i,j=n) end a=1*ones(n,n) if i+j=n+1 while i~=1 i=i-1; j=j+1; end; while j~=1 a(i,j)=0; i=i+1; j=j-1; a(i,j)=0; end; else while j~=n i=i-1; j=j+1; end; while i~=n a(i,j)=0; i=i+1; j=j-1; a(i,j)=0; end; end c=a 已知开环系统 ，绘制系统的Bode图，并判断闭环系统的稳定性，求闭环系统的单位脉冲响应。 num=50; den=[1 3 -10]; figure(1); margin(num,den); grid; [numc,denc]=cloop(num,den); figure(2); impulse(numc,denc); grid; Bode图 单位脉冲响应 已知开环系统的传递函数，绘制闭环系统的根轨迹，并分析其稳定性。 clear num=[1 2]; den=[1 8 22 24 9]; rlocus(num,den); 根轨迹 系统要求用鼠标输入要标示增益和相应得闭环极点的根轨迹上的点，执行后得到如下结果： Select a point in the graphics window selected_point =0.0083 - 3.1522i k =55.6381 p = -5.9866 0.0023 + 3.1553i 0.0023 - 3.1553i -2.0180 这里利用rlocfind（）函数来找出根轨迹与虚轴的交点，并求得交点处k