ex0513.m

matlab课本程序

function Ex0513()
% EX0513    参数传递绘制二阶系统时域响应
z1=0.3;
[x1,y1]=ex0502(z1);
plot(x1,y1)
hold on
z1=0.5;
[x2,y2]=ex0502(z1);
plot(x2,y2)
z1=0.707;
[x3,y3]=ex0502(z1);
plot(x3,y3);

function [x,y]=ex0502(zeta)
%子函数，画二阶系统时域曲线
x=0:0.1:20;
y=1-1/sqrt(1-zeta^2)*exp(-zeta*x).*sin...
    (sqrt(1-zeta^2)*x+acos(zeta));