diff options
| author | priyanka | 2015-06-24 15:03:17 +0530 |
|---|---|---|
| committer | priyanka | 2015-06-24 15:03:17 +0530 |
| commit | b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (patch) | |
| tree | ab291cffc65280e58ac82470ba63fbcca7805165 /758/CH10 | |
| download | Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.gz Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.bz2 Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.zip | |
initial commit / add all books
Diffstat (limited to '758/CH10')
| -rwxr-xr-x | 758/CH10/EX10.2/Ex_10_2.sce | 26 | ||||
| -rwxr-xr-x | 758/CH10/EX10.3/Ex_10_3.sce | 16 | ||||
| -rwxr-xr-x | 758/CH10/EX10.4/Ex_10_4.sce | 7 | ||||
| -rwxr-xr-x | 758/CH10/EX10.5/Ex_10_5.sce | 52 |
4 files changed, 101 insertions, 0 deletions
diff --git a/758/CH10/EX10.2/Ex_10_2.sce b/758/CH10/EX10.2/Ex_10_2.sce new file mode 100755 index 000000000..a27ae701a --- /dev/null +++ b/758/CH10/EX10.2/Ex_10_2.sce @@ -0,0 +1,26 @@ +//Example 10.2
+
+clc;clear;close;
+z=poly(0,'z');
+H=0.5*z/(z-0.5);
+B=8;
+pn=2^(-2*B)/12; //Noise power
+X=H*horner(H,1/z)/z;
+r=roots(denom(X));
+rl=length(r);
+rc=coeff(denom(X))
+q1=[];q2=[];
+for n=1:rl //Loop to separate poles inside the unit circle
+ if (abs(r(n))<1) then
+ q1=[q1 r(n)];
+ else
+ q2=[q2 r(n)];
+ end
+end
+P=numer(X)/rc(length(rc));
+Q1=poly(q1,'z');
+Q2=poly(q2,'z');
+I=residu(P,Q1,Q2); //Residue Calculation
+po=pn*I; //Output Noise power
+disp(pn,'Input Noise power');
+disp(po,'Output Noise power');
\ No newline at end of file diff --git a/758/CH10/EX10.3/Ex_10_3.sce b/758/CH10/EX10.3/Ex_10_3.sce new file mode 100755 index 000000000..0f7fe6c52 --- /dev/null +++ b/758/CH10/EX10.3/Ex_10_3.sce @@ -0,0 +1,16 @@ +//Example 10.3
+
+clc; clear;
+//y(n)=0.9y(n-1)+x(n)
+//Input x(n)=0
+n=-1;y=12; //Initial Condition y(-1)=12
+flag=1;
+while n<8
+ n=n+1;
+ y=[y 0.9*y(n+1)];
+ yr=round(y);
+end
+disp(n,'n=');
+disp(y,'y(n)-exact');
+disp(yr,'y(n)-rounded');
+disp([-yr(n+2) yr(n+2)],'Deadband interval ')
\ No newline at end of file diff --git a/758/CH10/EX10.4/Ex_10_4.sce b/758/CH10/EX10.4/Ex_10_4.sce new file mode 100755 index 000000000..0216a2f65 --- /dev/null +++ b/758/CH10/EX10.4/Ex_10_4.sce @@ -0,0 +1,7 @@ +//Example 10.4
+
+clc; clear;
+//y(n)=0.9y(n-1)+x(n)
+a=0.9;
+l=ceil(0.5/(1-abs(a)));
+disp([-l l],'Deadband interval ')
\ No newline at end of file diff --git a/758/CH10/EX10.5/Ex_10_5.sce b/758/CH10/EX10.5/Ex_10_5.sce new file mode 100755 index 000000000..5374722ce --- /dev/null +++ b/758/CH10/EX10.5/Ex_10_5.sce @@ -0,0 +1,52 @@ +//Example 10.5
+
+clc;clear;close;
+x=poly(0,'x'); //x=2^-2B
+z=poly(0,'z');
+H1=1/(1-0.9/z);
+H2=1/(1-0.8/z);
+H=H1*H2;
+pn=x/12; //Input Noise power
+
+//Calculation of output noise for H1(z)
+X1=H*horner(H,1/z)/z;
+r1=roots(denom(X1));
+rc1=coeff(denom(X1));
+q1=[];s1=[];
+for n=1:length(r1) //Loop to separate poles inside the unit circle
+ if (abs(r1(n))<1) then
+ q1=[q1 r1(n)];
+ else
+ s1=[s1 r1(n)];
+ end
+end
+P1=numer(X1)/rc1(length(rc1));
+Q1=poly(q1,'z');
+S1=poly(s1,'z');
+I1=abs(residu(P1,Q1,S1)); //Residue Calculation
+po1=pn*I1; //Output Noise power
+
+//Calculation of output noise for H2(z)
+X2=H2*horner(H2,1/z)/z;
+r2=roots(denom(X2));
+rc2=coeff(denom(X2));
+q2=[];s2=[];
+for n=1:length(r2) //Loop to separate poles inside the unit circle
+ if (abs(r2(n))<1) then
+ q2=[q2 r2(n)];
+ else
+ s2=[s2 r2(n)];
+ end
+end
+P2=numer(X2)/rc2(length(rc2));
+Q2=poly(q2,'z');
+S2=poly(s2,'z');
+I2=abs(residu(P2,Q2,S2)); //Residue Calculation
+po2=pn*I2; //Output Noise power
+
+po=po1+po2;
+disp(pn,'Input Noise power');
+disp(I1,'I1=');disp(I2,'I2=');
+disp(po1,'Output Noise power for H1(z)');
+disp(po2,'Output Noise power for H2(z)');
+disp(po,'Total Output Noise power');
\ No newline at end of file |