test4_2_22.m

两点间的头镇试验源码

% 实验4第二题的第2小题，绘制的kesai范围改动一下，kesai=0.2:0.3:2 
% 这是第二方法
clear,clc
t=linspace(0,18,200);
xi=[0.2:0.3:2];

% color_string={'r-', 'b:' 'k-.' 'g--'}; 
color_string={'b', 'g' 'r' 'c' 'm' 'y' 'k'};
% 这里用到了元胞的概念，还没讲，定义元胞用大括号，元素可以为MATLAB许可的任意数据类型

for n=1:length(xi)
    temp=1-xi(n)^2;
    beta=sqrt(temp);
    theta=atan(beta/xi(n));
    y=1-1/beta*exp(-xi(n)*t).*sin(beta*t+theta);
    plot(t,y,color_string{n}),hold on
    
    switch xi(n)
        case 0.2
            t1=2.5;
            y1=1-1/beta*exp(-xi(n)*t1).*sin(beta*t1+theta);
            text(t1+1.6,y1,'\fontsize{12}{\xi}=0.2')
        case 2.0
            t1=2.5;
            y1=1-1/beta*exp(-xi(n)*t1).*sin(beta*t1+theta);
            text(t1+0.1 ,y1,'\fontsize{12}{\xi}=0.8')
    end
%     pause
end

% 不作要求掌握，可自行模仿
title('\fontsize{12}二阶系统y=1-1/{\beta}*exp(-{\xi}t)sin({\beta}t+{\theta})的归一化响应')
xlabel('\fontsize{14} \itt \rm\rightarrow')   % \it 其后用斜体， \rm 其后用正体
ylabel('\fontsize{14} \ity \rm\rightarrow')