求阶跃响应的典型指标.docVIP

%% 求阶跃响应的典型指标function main_GetPerformanceOfStepResponseglobal gTolerancegTolerance = 0.05; % 调整时间的偏差容许范围‘W=’) b=input(‘b=’)g = tf(W^2, [1, 2*W*b, b^2]);t = 0:0.01:15;y = step(g,t); G=zpk(g);G.p{1} figure(2) hold on [mag,phase,w]=bode(g); [gm,pm,wcp,wcg]=margin(mag,phase,w) margin(mag,phase,w) Wb=W*sqrt((1-2*b^2)+sqrt((1-2*b^2)^2+1)) if(b=0b=(sqrt(2)/2)) Wr=W*(sqrt(1-2*b^2) Mr=1/(2*b*sqrt(1-b^2)) else Mr=0 end %% 计算阶跃响应的指标stepvalue = 1;[OverShoot, RiseTime, PeakTime, AdjustTime, SteadyStateError] = GetPerformanceOfStepResponse(t, y, stepvalue); % 绘图figureplot(t,y)grid on line([PeakTime, PeakTime], [0, (1 + OverShoot/100)*stepvalue], color, r)text(PeakTime, stepvalue*0.05, sprintf(峰值时间%.2f,PeakTime))text(PeakTime, (1 + OverShoot/100 + 0.05)*stepvalue, sprintf(超调量%.2f%%,OverShoot)) line([RiseTime, RiseTime], [0, stepvalue], color, r)text(RiseTime, -stepvalue*0.05, sprintf(上升时间%.2f,RiseTime)) line([AdjustTime, AdjustTime], [0, stepvalue*(1 + gTolerance)], color, r)text(AdjustTime, stepvalue*0.05, sprintf(调整时间%.2f,AdjustTime)) line([AdjustTime t(end)], stepvalue*[(1 - gTolerance), (1 - gTolerance)], color, r, linestyle, --)text(AdjustTime, stepvalue*(1 - gTolerance-0.05), sprintf(容许范围%.2f, 1 - gTolerance))line([AdjustTime t(end)], stepvalue*[(1 + gTolerance), (1 + gTolerance)], color, r, linestyle, --)text(AdjustTime, stepvalue*(1 + gTolerance+0.05), sprintf(容许范围%.2f, 1 + gTolerance))text(t(end)*0.9, stepvalue*1.05, sprintf(稳态误差%f, SteadyStateError)) end %% 求阶跃响应的典型指标function [OverShoot, RiseTime, PeakTime, AdjustTime, SteadyStateError] = GetPerformanceOfStepResponse(t, y, stepvalue)global gTolerance % 超调量和峰值时间[OSValue, OSIndex] = max(y);OverShoot = (OSValue - stepvalue)/stepvalue*100;PeakTime = t(OSIndex); % 上升时间index = find(y = stepvalue, 1, first);RiseTime = t(index); % 调整时间和稳态误差index1 = find(y = stepvalue*(1 - gTolerance), 1, last); % 容许范围由全局变量指定index2 = find(y = stepvalue*(1 + gTolerance), 1, last); if isempty(index2) % 如果没有超调量，此值为空??