%�@SŽšƒvƒ�ƒtƒ@ƒCƒ‹ 2006/10/10 T.Ogawa

clear
format;

Fref0 = 0.25;   %
Fref1 = 25;     %

Tstop1  = 1;
Tstop2  = 2;

Tsin = 0.5;        %‘�‰ÁŽžŠÔ[sec]

dTref = 1e-3;      %ƒTƒ“ƒvƒŠƒ“ƒOŽžŠÔ[sec]

Fref = [];
Tref = [];

for Tsim = 0:dTref:Tstop1
    Fref = [Fref
            Fref0  ];            
    Tref = [Tref
            Tsim];
end

for Tsim = (Tstop1+dTref):dTref:(Tstop1+Tsin)
    
    Fref = [Fref
             Fref0 + (Fref1-Fref0)/Tsin*((Tsim-Tstop1)-Tsin/(2*pi)*sin(2*pi/Tsin*(Tsim-Tstop1)))  ];    
    Tref = [Tref
            Tsim];      
end

for Tsim = (Tstop1+Tsin+dTref):dTref:(Tstop2)
    Fref = [Fref
            Fref1 ];
    Tref = [Tref
            Tsim];
end

 
%ƒOƒ‰ƒt 
plot(Tref,Fref,'b-');
grid on


