Курсовая работа: Визначення динамічних похибок вимірювань

dsolve({1.1*diff(diff(y(t),t),t)+2.5*diff(y(t),t)+

+3.1*sin(2.1*t)=0.5*Heaviside(t),y(0)=0, D(y)(0)=0});

y2(t) := 34100/115861*sin(21/10*t)+775000/2433081*cos(21/10*t)+

+1/5*Heaviside(t)*t+11/125*exp(-25/11*t)*Heaviside(t)-

-11/125*Heaviside(t)+157542/579305*exp(-25/11*t)-62/105;

plot({y2(t)},t=0..20);

dsolve({1.1*diff(diff(y(t),t),t)+2.5*diff(y(t),t)+

+3.1*sin(2.1*t)=0.5*Dirac(t),y(0)=0, D(y)(0)=0});

y3(t) := -1/5*exp(-25/11*t)*Heaviside(t)+1/5*Heaviside(t)+

+34100/115861*sin(21/10*t)+775000/2433081*cos(21/10*t)+

+273403/579305*exp(-25/11*t)-83/105;

plot({y3(t)},t=0..20);

S(f):=(4.41-(2*Pi*f)^2)/(2.2*(2*Pi*f)^4-9.3*(2*Pi*f)^2+6.51-

-5*sqrt(-1)*(2*Pi*f)^3+22.05*2*Pi*f*sqrt(-1));

S(f):=(4.41-(2*Pi*f)^2)*((2.2*(2*Pi*f)^4-9.3*(2*Pi*f)^2+6.51) +

sqrt(-1)*(5*(2*Pi*f)^3+22.05*2*Pi*f))/ ((2.2*(2*Pi*f)^4-

-9.3*(2*Pi*f)^2+6.51)^2 -( sqrt(-1)*(5*(2*Pi*f)^3+

+22.05*2*Pi*f))^2);

A(f):=(4.41-(2*Pi*f)^2)*(2.2*(2*Pi*f)^4-9.3*(2*Pi*f)^2+

+6.51) / ((2.2*(2*Pi*f)^4-9.3*(2*Pi*f)^2+6.51)^2 -( sqrt(-