summaryrefslogtreecommitdiff
path: root/1340/CH8/EX8.2/8_2.sce
blob: b523ca26aeb775a9ae7e65b5f8bcafaba24bce0c (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
clc;
s = %s;
f = (s+3)/(s*(s+1)*(s+2)*(s+4));
disp(f,"G(s)=");
evans(f);
xgrid(2);
printf("\n  Poles at 0,-1,-2 and -4, zero at -3  \n");
x = denom(f);
y = numer(f);

xr = roots(x);
yr = roots(y);

[kx,lx]=size(xr);
[ky,ly]=size(yr);

sumX=0, sumY=0;
for j  = 1:kx
    sumX = sumX + xr(j);
end
for j  = 1:ky
    sumY = sumY + yr(j);
end
disp((sumX-sumY)/(kx-ky),"Real axis intercept:");
printf("\n angle = (2k+1)*180/(sum of finite poles-sum of finte zeroes)");
printf("\n for k=0")
    theta = 180/(kx-ky);
    disp(theta,"angles made by the asymptote to the real axis for :");
printf("1 finite zero and 3 infinite zeroes");