ex0506.m

来自「不错的matlab学习资料。其中包含有matlab程序如何转化为c」· M 代码 · 共 11 行

function y=Ex0506(zeta)
% EX0506    使用if结构的二阶系统时域响应
x=0:0.1:20;
if (zeta>0)&(zeta<1)
    y=1-1/sqrt(1-zeta^2)*exp(-zeta*x).*sin(sqrt(1-zeta^2)*x+acos(zeta));
elseif zeta==1
    y=1-exp(-x).*(1+x);
end
plot(x,y)