仿真的三次作业和三次上机实验报告,百度云链接: http://pan.baidu.com/s/1eRsWZ2i 密码: 9s8v注意不要直接抄袭喔!如果老师发现所有学生的作业都一样的话你们懂的
当时考试的时候总结的,考题基本【或者肯定】就从这里出……
采用四阶龙格库塔法求如下二阶系统的阶跃响应。
解:
T=10;%仿真时长b=0.5;w=5; %%%%%%%%%%%采用矩阵方法的四阶龙格库塔法%%%%%dt=0.05;%积分步长N=T/dt;%积分总步数Y1=[0;0];YY=Y1;A=[0 1;-w^2 -2*b*w];%%%%状态观测矩阵%%%%B=[0;1];C=[1 0];U=w^2;for i=1:N; K1=A*YY+B*U; K2=A*(YY+dt/2*K1)+B*U; K3=A*(YY+dt/2*K2)+B*U; K4=A*(YY+dt*K3)+B*U; YY=YY+dt/6*(K1+2*K2+2*K3+K4); Y1=[Y1 YY];end;figureplot([0:dt:T],C*Y1);grid on
求脉冲响应:
%把U定义在for中for i=1:N;if i < 1/dt U = w^2; else U = 0; end K1=A*YY+B*U; K2=A*(YY+dt/2*K1)+B*U; K3=A*(YY+dt/2*K2)+B*U; K4=A*(YY+dt*K3)+B*U; YY=YY+dt/6*(K1+2*K2+2*K3+K4); Y1=[Y1 YY];end;
ode45解:
function r=second_order(t,x);b=0.5;w=5;r=[x(2);-2*b*w*x(2)-w^2*x(1)+w^2];endt=[0 T];x0=[0 0];[t y]=ode45('second_order',t,x0);figure;hold onplot(t,y(:,1));
采用四阶龙格库塔法求高阶系统单位阶跃响应曲线的数值解,并从数值解中得到系统的上升时间、峰值时间、调节时间等时域指标。
解:
t=40;%仿真时长b=0.5;w=10;T = 5; %%%%%%%%%%%采用矩阵方法的四阶龙格库塔法%%%%%dt=0.05;%积分步长N=