From b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b Mon Sep 17 00:00:00 2001 From: priyanka Date: Wed, 24 Jun 2015 15:03:17 +0530 Subject: initial commit / add all books --- 46/CH5/EX5.2/Example5_2.sce | 35 +++++++++++++++++++++++++++++++++++ 46/CH5/EX5.2/figure5_2.jpg | Bin 0 -> 88249 bytes 2 files changed, 35 insertions(+) create mode 100755 46/CH5/EX5.2/Example5_2.sce create mode 100755 46/CH5/EX5.2/figure5_2.jpg (limited to '46/CH5/EX5.2') diff --git a/46/CH5/EX5.2/Example5_2.sce b/46/CH5/EX5.2/Example5_2.sce new file mode 100755 index 000000000..475584914 --- /dev/null +++ b/46/CH5/EX5.2/Example5_2.sce @@ -0,0 +1,35 @@ +//Example 5.2 +clear +clc +tau=0.1;//min +xs=100;//Fahrenheit +ys=100;//Fahrenheit +A=2;//Fahrenheit +f=10/%pi;//cycles/min +w=2*%pi*f;//rad/min +//From Eq.(5.25), the amplitude of the response and the phase angle are calculated;thus +disp('Fahrenheit',A/sqrt((tau*w)^2+1),'A/sqrt((tau*w)^2+1)=') +phi=atan(-w*tau);//radians +phi=phi*180/%pi;//degrees +disp('degrees',phi,'phase lag=') +t=0:0.01:1; +//From Eq.(5.19), the input of the thermometer is therefore +disp("X(t)=2*sin(20*t)"); +//or +xt=xs+2*sin(20*t); +//The response of the thermometer is therefore +disp("Y(t)=0.8944*sin(20*t-63.4349)") +//or +yt=ys+0.8944*sin(20*t-63.4349); +Lag=phi/(360*f);//min +Lag=abs(Lag);//min +disp('min',Lag,'Lag=') +clf; +plot(t,yt) +plot(t,xt) +xlabel('time') +ylabel('x(t) , y(t)') +title('x(t) , y(t) Vs time') +xgrid + + diff --git a/46/CH5/EX5.2/figure5_2.jpg b/46/CH5/EX5.2/figure5_2.jpg new file mode 100755 index 000000000..a8b463c04 Binary files /dev/null and b/46/CH5/EX5.2/figure5_2.jpg differ -- cgit