Assignment 4
Assignment 4
Ans. CODE:
clc
clear
t = 0:.01:3;
Dd = Dd0/sqrt(1-z^2)*exp(-z*wn*t).*sin(wd*t + th);
Dw = -wn*Dd0/sqrt(1-z^2)*exp(-z*wn*t).*sin(wd*t);
f = f0 + Dw/(2*pi); % Frequency in Hz
D = [0; 0];
f = f0 + Dw/(2*pi); % Frequency in Hz
OUTPUT:
E = 1.3500
Pmax = 2.0769
d0 = 0.2931
Ps = 1.9884
wn = 6.1405
z = 0.2131
wd = 5.9995
fd = 0.9549
tau = 0.764
Q2.
Ans. CODE:
Pm = 0.8; E= 1.17; V = 1.0; X1 = 0.65; X2 = inf; X3 = 0.65;
if exist('Pm')~=1
if exist('E')~=1
if exist('V')~=1
if exist('X1')~=1
if exist('X2')~=1
if exist('X3')~=1
delta = 0:.01:pi;
cosdc = (Pm*(dmax-d0)+Pe3max*cos(dmax)-Pe2max*cos(d0))/(Pe3max-Pe2max);
if abs(cosdc) > 1
return
else, end
dc=acos(cosdc);
if dc > dmax
return
else, end
x=(d0:.1:dc);
y=Pe2max*sin(x*pi/180);
y1=Pe2max*sin(d0*pi/180);
y2=Pe2max*sin(dc*pi/180);
x=[d0 x dc];
y=[Pm y Pm];
xx=dc:.1:dmax;
h=Pe3max*sin(xx*pi/180);
xx=[dc xx dmax];
hh=[Pm h Pm];
delta=delta*180/pi;
if X2 == inf
if H ~= 0
d0r=d0*pi/180; dcr=dc*pi/180;
tc = sqrt(2*H*(dcr-d0r)/(pi*60*Pm));
else, end
else, end
%clc
else, end
h = figure; figure(h);
fill(x,y,'m')
hold;
fill(xx,hh,'c')
plot(delta, Pe1,'-', delta, Pe2,'r-', delta, Pe3,'g-', Pmx, Pmy,'b-', x0,y0, xc,yc, xm,ym), grid
hold off;
OUTPUT:
For this case tc can be found from analytical formula.