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 /67 | |
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 '67')
123 files changed, 1556 insertions, 0 deletions
diff --git a/67/CH1/EX1.1.a/example11a.sce b/67/CH1/EX1.1.a/example11a.sce new file mode 100755 index 000000000..592ad32ba --- /dev/null +++ b/67/CH1/EX1.1.a/example11a.sce @@ -0,0 +1,7 @@ +//Example 1.1a
+//Determine whether the given signal is periodic or not
+clc;
+t=0:1/100:1
+x=sin(15*%pi*t);
+plot(x);
+disp('ploting the signal and showing that it is periodic with period=2pi/15pi');
\ No newline at end of file diff --git a/67/CH1/EX1.1.b/example11b.sce b/67/CH1/EX1.1.b/example11b.sce new file mode 100755 index 000000000..aa10ba752 --- /dev/null +++ b/67/CH1/EX1.1.b/example11b.sce @@ -0,0 +1,7 @@ +//Example 1.1b
+//Determine whether the given signal is periodic or not
+clc;
+t=0:1/100:5
+x=sin(sqrt(2)*%pi*t);
+plot(x);
+disp('ploting the signal and showing that itis periodic with period=2pi/sqrt(2)pi');
\ No newline at end of file diff --git a/67/CH1/EX1.14.a/example114a.sce b/67/CH1/EX1.14.a/example114a.sce new file mode 100755 index 000000000..a03938de5 --- /dev/null +++ b/67/CH1/EX1.14.a/example114a.sce @@ -0,0 +1,22 @@ +//Example 1.14a
+clc;
+x=[1,2,3,4,0,4,3,2,1]
+t=-length(x)/2:length(x)/2
+count=0
+mid=ceil(length(x)/2)
+y=zeros(1,length(x))
+y(mid+1:$)=x($:-1:mid+1)
+for t=-1:-1:-mid
+ y(t+1+mid)=x(-t)
+end
+for i=1:length(x)
+ if(y(i)==x(i))
+ count=count+1
+ end
+end
+if(count==length(x))
+ disp('THE GIVEN SYSTEM IS CAUSAL')
+else
+ disp('Since it depends on future values')
+ disp('THE GIVEN SYSTEM IS NON CAUSAL')
+end
\ No newline at end of file diff --git a/67/CH1/EX1.18.a/example118a.sce b/67/CH1/EX1.18.a/example118a.sce new file mode 100755 index 000000000..3c487489a --- /dev/null +++ b/67/CH1/EX1.18.a/example118a.sce @@ -0,0 +1,15 @@ +//Example 1.18a
+clc;
+t0=1;
+T=10;
+for t=1:T
+ x(t)=2*%pi*t/T;
+ y(t)=sin(x(t));
+end
+inputshift=sin(x(T-t0));
+outputshift=y(T-t0);
+if(inputshift==outputshift)
+ disp('THE GIVEN SYSTEM IS TIME INVARIANT')
+else
+ disp('THE GIVEN SYSTEM IS TIME VARIANT');
+end
\ No newline at end of file diff --git a/67/CH1/EX1.18.b/example118b.sce b/67/CH1/EX1.18.b/example118b.sce new file mode 100755 index 000000000..fad9bee4c --- /dev/null +++ b/67/CH1/EX1.18.b/example118b.sce @@ -0,0 +1,15 @@ +//Example 1.18b
+clc;
+t0=2;
+T=10;
+for t=1:T
+ x(t)=t;
+ y(t)=t*x(t);
+end
+inputshift=x(T-t0);
+outputshift=y(T-t0);
+if(inputshift==outputshift)
+ disp('THE GIVEN SYSTEM IS TIME INVARIANT')
+else
+ disp('THE GIVEN SYSTEM IS TIME VARIANT');
+end
\ No newline at end of file diff --git a/67/CH1/EX1.18.c/example118c.sce b/67/CH1/EX1.18.c/example118c.sce new file mode 100755 index 000000000..9f8bbc2f3 --- /dev/null +++ b/67/CH1/EX1.18.c/example118c.sce @@ -0,0 +1,15 @@ +//Example 1.18c
+clc;
+t0=2;
+T=10;
+for t=1:T
+ x(t)=t;
+ y(t)=x(t)*cos(200*%pi*t);
+end
+inputshift=x(T-t0);
+outputshift=y(T-t0);
+if(inputshift==outputshift)
+ disp('THE GIVEN SYSTEM IS TIME INVARIANT')
+else
+ disp('THE GIVEN SYSTEM IS TIME VARIANT');
+end
\ No newline at end of file diff --git a/67/CH1/EX1.19.a/example119a.sce b/67/CH1/EX1.19.a/example119a.sce new file mode 100755 index 000000000..f26983dbc --- /dev/null +++ b/67/CH1/EX1.19.a/example119a.sce @@ -0,0 +1,30 @@ +//Example 1.19a
+clc;
+x1=[1,1,1,1]
+x2=[2,2,2,2]
+a=1
+b=1
+for t=1:length(x1)
+ x3(t)=a*x1(t)+b*x2(t)
+end
+for t=1:length(x1)
+ y1(t)=t*x1(t)
+ y2(t)=t*x2(t)
+ y3(t)=t*x3(t)
+end
+for t=1:length(y1)
+ z(t)=a*y1(t)+b*y2(t)
+end
+count=0
+for n=1:length(y1)
+ if(y3(t)==z(t))
+ count=count+1;
+ end
+end
+if(count==length(y3))
+ disp('It satisfy the superposition principle');
+ disp('THE GIVEN SYSTEM IS LINEAR ');
+else
+ disp('It does not satisfy superposition principle ');
+ disp('THE GIVEN SYSTEM IS NON LINEAR');
+end
\ No newline at end of file diff --git a/67/CH1/EX1.19.b/example119b.sce b/67/CH1/EX1.19.b/example119b.sce new file mode 100755 index 000000000..d953773c2 --- /dev/null +++ b/67/CH1/EX1.19.b/example119b.sce @@ -0,0 +1,30 @@ +//Example 1.19b
+clc;
+x1=[1,1,1,1]
+x2=[2,2,2,2]
+a=1
+b=1
+for t=1:length(x1)
+ x3(t)=a*x1(t)+b*x2(t)
+end
+for t=1:length(x1)
+ y1(t)=x1(t)^2
+ y2(t)=x2(t)^2
+ y3(t)=x3(t)^2
+end
+for t=1:length(y1)
+ z(t)=a*y1(t)+b*y2(t)
+end
+count=0
+for n=1:length(y1)
+ if(y3(t)==z(t))
+ count=count+1;
+ end
+end
+if(count==length(y3))
+disp('It satisfy the superposition principle');
+disp('THE GIVEN SYSTEM IS LINEAR ');
+else
+ disp('It does not satisfy superposition principle ');
+ disp('THE GIVEN SYSTEM IS NON LINEAR');
+end
\ No newline at end of file diff --git a/67/CH1/EX1.20.b/example120b.sce b/67/CH1/EX1.20.b/example120b.sce new file mode 100755 index 000000000..f16107c5c --- /dev/null +++ b/67/CH1/EX1.20.b/example120b.sce @@ -0,0 +1,30 @@ +//Example 1.20b
+clc;
+x1=[1,1,1,1]
+x2=[2,2,2,2]
+a=1
+b=1
+for t=1:length(x1)
+ x3(t)=a*x1(t)+b*x2(t)
+end
+for t=1:length(x1)
+ y1(t)=x1(t)^2
+ y2(t)=x2(t)^2
+ y3(t)=x3(t)^2
+end
+for t=1:length(y1)
+ z(t)=a*y1(t)+b*y2(t)
+end
+count=0
+for n=1:length(y1)
+ if(y3(t)==z(t))
+ count=count+1;
+ end
+end
+if(count==length(y3))
+disp('It satisfy the superposition principle');
+disp('THE GIVEN SYSTEM IS LINEAR ');
+else
+ disp('It does not satisfy superposition principle ');
+ disp('THE GIVEN SYSTEM IS NON LINEAR');
+end
\ No newline at end of file diff --git a/67/CH1/EX1.21/example121.sce b/67/CH1/EX1.21/example121.sce new file mode 100755 index 000000000..8ce87a563 --- /dev/null +++ b/67/CH1/EX1.21/example121.sce @@ -0,0 +1,32 @@ +//Example 1.21
+clc;
+x1=[1,1,1,1]
+x2=[2,2,2,2]
+a1=1;
+b1=1;
+a=7;
+b=5;
+for t=1:length(x1)
+ x3(t)=a1*x1(t)+b1*x2(t)
+end
+for t=1:length(x1)
+ y1(t)=a*x1(t)+b
+ y2(t)=a*x2(t)+b
+ y3(t)=a*x3(t)+b
+end
+for t=1:length(y1)
+ z(t)=a1*y1(t)+b1*y2(t)
+end
+count=0
+for n=1:length(y1)
+ if(y3(t)==z(t))
+ count=count+1;
+ end
+end
+if(count==length(y3))
+disp('It satisfy the superposition principle');
+disp('THE GIVEN SYSTEM IS LINEAR ');
+else
+ disp('It does not satisfy superposition principle ');
+ disp('THE GIVEN SYSTEM IS NON LINEAR');
+end
\ No newline at end of file diff --git a/67/CH1/EX1.22/example122.sce b/67/CH1/EX1.22/example122.sce new file mode 100755 index 000000000..492abd3f6 --- /dev/null +++ b/67/CH1/EX1.22/example122.sce @@ -0,0 +1,31 @@ +//Example 1.22
+clc;
+x1=[1,1,1,1]
+x2=[2,2,2,2]
+a1=1
+b1=1
+Wc=%pi
+for t=1:length(x1)
+ x3(t)=a1*x1(t)+b1*x2(t)
+end
+for t=1:length(x1)
+ y1(t)=x1(t)*cos(Wc*t)
+ y2(t)=x2(t)*cos(Wc*t)
+ y3(t)=x3(t)*cos(Wc*t)
+end
+for t=1:length(y1)
+ z(t)=a1*y1(t)+b1*y2(t)
+end
+count=0
+for n=1:length(y1)
+ if(y3(t)==z(t))
+ count=count+1;
+ end
+end
+if(count==length(y3))
+disp('It satisfy the superposition principle');
+disp('THE GIVEN SYSTEM IS LINEAR ');
+else
+ disp('It does not satisfy superposition principle ');
+ disp('THE GIVEN SYSTEM IS NON LINEAR');
+end
\ No newline at end of file diff --git a/67/CH1/EX1.25/example125.sce b/67/CH1/EX1.25/example125.sce new file mode 100755 index 000000000..c627ac2c9 --- /dev/null +++ b/67/CH1/EX1.25/example125.sce @@ -0,0 +1,32 @@ +//Example 1.25
+clc;
+x1=[1,1,1,1]
+x2=[2,2,2,2]
+a1=1;
+b1=1;
+a=7;
+b=3;
+for t=1:length(x1)
+ x3(t)=a1*x1(t)+b1*x2(t)
+end
+for t=1:length(x1)
+ y1(t)=a*x1(t)+b
+ y2(t)=a*x2(t)+b
+ y3(t)=a*x3(t)+b
+end
+for t=1:length(y1)
+ z(t)=a1*y1(t)+b1*y2(t)
+end
+count=0
+for n=1:length(y1)
+ if(y3(t)==z(t))
+ count=count+1;
+ end
+end
+if(count==length(y3))
+disp('It satisfy the superposition principle');
+disp('THE GIVEN SYSTEM IS LINEAR ');
+else
+ disp('It does not satisfy superposition principle ');
+ disp('THE GIVEN SYSTEM IS NON LINEAR');
+end
\ No newline at end of file diff --git a/67/CH1/EX1.27/example127.sce b/67/CH1/EX1.27/example127.sce new file mode 100755 index 000000000..7af90767f --- /dev/null +++ b/67/CH1/EX1.27/example127.sce @@ -0,0 +1,6 @@ +//Example 1.27
+//Energy of the signal x(t)=Aexp(-a*t).u(t)
+clc;
+A=2;
+a=0.5;
+E=integrate('(A*exp(-a*t))^2','t',0,100);//Energy of the given signal
\ No newline at end of file diff --git a/67/CH1/EX1.28/example128.sce b/67/CH1/EX1.28/example128.sce new file mode 100755 index 000000000..e6ed32d2b --- /dev/null +++ b/67/CH1/EX1.28/example128.sce @@ -0,0 +1,5 @@ +//Example 1.28
+//Power of the signal x(t)=A
+clc;
+A=2;
+P=(integrate('A^2','t',0,100))/(2*100)
\ No newline at end of file diff --git a/67/CH1/EX1.30/example130.sce b/67/CH1/EX1.30/example130.sce new file mode 100755 index 000000000..278883820 --- /dev/null +++ b/67/CH1/EX1.30/example130.sce @@ -0,0 +1,4 @@ +//Example 1.30
+//Detremine the energy of the signal x(n)=0.5^n for n >0
+clc;
+E=integrate('(0.5^n)','n',0,1000);
\ No newline at end of file diff --git a/67/CH1/EX1.31.a/example131a.sce b/67/CH1/EX1.31.a/example131a.sce new file mode 100755 index 000000000..02666132b --- /dev/null +++ b/67/CH1/EX1.31.a/example131a.sce @@ -0,0 +1,6 @@ +//Example 1.31a
+//Determine whether the given signal is periodic or not
+clc;
+n=0:1/100:10
+x=sin(6*%pi*n/7);
+plot(x)//plotting the signal and showing it is periodic with period 2pi/(6pi/7);
\ No newline at end of file diff --git a/67/CH1/EX1.31.b/example131b.sce b/67/CH1/EX1.31.b/example131b.sce new file mode 100755 index 000000000..d0c043dac --- /dev/null +++ b/67/CH1/EX1.31.b/example131b.sce @@ -0,0 +1,6 @@ +//Example 1.31b
+//Determine whether the given signal is periodic or not
+clc;
+n=0:1/1000:100
+x=sin(n/8);
+plot(x);//plotting the signal and showing that it is periodic with period 16pi
\ No newline at end of file diff --git a/67/CH1/EX1.33/example133.sce b/67/CH1/EX1.33/example133.sce new file mode 100755 index 000000000..fe99d37d7 --- /dev/null +++ b/67/CH1/EX1.33/example133.sce @@ -0,0 +1,14 @@ +//Example 1.33
+//Find the power of the signal x(t)=Acos(Wot+theeta)
+clc;
+A=10;
+T=4;
+Wo=(2*%pi)/T;
+for i=1:T
+ x(i)=A*cos(Wo*i);
+end
+p=0;
+for i=1:T
+ p=p+(abs(x(i)^2))/T;
+end
+disp(p,'The power of the given signal is =');
\ No newline at end of file diff --git a/67/CH1/EX1.34/example134.sce b/67/CH1/EX1.34/example134.sce new file mode 100755 index 000000000..b705a4866 --- /dev/null +++ b/67/CH1/EX1.34/example134.sce @@ -0,0 +1,10 @@ +//Example 1.34
+//Find energy of x(t)=8exp(2+i4pi)t
+clc;
+E=0;
+for t=1:100
+ x(t)=8*exp((2+(%i*4*%pi))*t);
+end
+for t=1:100
+ E=E+x(t)^2;
+end
\ No newline at end of file diff --git a/67/CH1/EX1.39.a/example139a.sce b/67/CH1/EX1.39.a/example139a.sce new file mode 100755 index 000000000..f561a72bf --- /dev/null +++ b/67/CH1/EX1.39.a/example139a.sce @@ -0,0 +1,6 @@ +//Example 1.39a
+//Sketch the signal x(t)=u(t)
+clc;
+t=0:1/100:10
+x=1;
+plot(t,x);
\ No newline at end of file diff --git a/67/CH1/EX1.39.b/example139b.sce b/67/CH1/EX1.39.b/example139b.sce new file mode 100755 index 000000000..c12228730 --- /dev/null +++ b/67/CH1/EX1.39.b/example139b.sce @@ -0,0 +1,6 @@ +//Example 1.39b
+//Sketch the signal x(t)=tu(t)
+clc;
+t=0:1/100:10
+x=t
+plot(t,x)
\ No newline at end of file diff --git a/67/CH1/EX1.4/example14.sce b/67/CH1/EX1.4/example14.sce new file mode 100755 index 000000000..6f795d28b --- /dev/null +++ b/67/CH1/EX1.4/example14.sce @@ -0,0 +1,9 @@ +//Example 1.4
+//Sketch the signal x(t)=Asin(t)
+clc;
+A=0.5;
+t=0:1/100:10
+x=A*sin(t);
+plot(x);
+//since it is a periodic signal so it is power signal
+P=(integrate('((0.5)^2)*(sin(t)^2)','t',0,2*%pi))/(2*%pi);
\ No newline at end of file diff --git a/67/CH1/EX1.43/example143.sce b/67/CH1/EX1.43/example143.sce new file mode 100755 index 000000000..595fa3c2d --- /dev/null +++ b/67/CH1/EX1.43/example143.sce @@ -0,0 +1,30 @@ +//Example 1.43
+clc;
+x1=[1,1,1,1]
+x2=[2,2,2,2]
+a=1
+b=1
+for t=1:length(x1)
+ x3(t)=a*x1(t)+b*x2(t)
+end
+for t=1:length(x1)
+ y1(t)=x1(t)^2
+ y2(t)=x2(t)^2
+ y3(t)=x3(t)^2
+end
+for t=1:length(y1)
+ z(t)=a*y1(t)+b*y2(t)
+end
+count=0
+for n=1:length(y1)
+ if(y3(t)==z(t))
+ count=count+1;
+ end
+end
+if(count==length(y3))
+disp('It satisfy the superposition principle');
+disp('THE GIVEN SYSTEM IS LINEAR ');
+else
+ disp('It does not satisfy superposition principle ');
+ disp('THE GIVEN SYSTEM IS NON LINEAR');
+end
\ No newline at end of file diff --git a/67/CH1/EX1.47/example147.sce b/67/CH1/EX1.47/example147.sce new file mode 100755 index 000000000..6bc4172e9 --- /dev/null +++ b/67/CH1/EX1.47/example147.sce @@ -0,0 +1,17 @@ +//Example 1.47
+clc;
+k0=2;
+n0=2;
+N=10;
+x=[1,2,3,4,5,6,7,8,9,10];
+y=zeros(1,length(x));
+for n=1:length(x)/k0
+ y(n)=x(k0*n);
+end
+inputshift=x(N-n0);
+outputshift=y(N-n0);
+if(inputshift==outputshift)
+ disp('THE GIVEN SYSTEM IS TIME INVARIANT')
+else
+ disp('THE GIVEN SYSTEM IS TIME VARIANT');
+end
\ No newline at end of file diff --git a/67/CH1/EX1.49.a/example149a.sce b/67/CH1/EX1.49.a/example149a.sce new file mode 100755 index 000000000..15e9af5d5 --- /dev/null +++ b/67/CH1/EX1.49.a/example149a.sce @@ -0,0 +1,7 @@ +//Example 1.49a
+//Determine whether the signal x(n)=sin(7/9*pi*(n^2)+1)
+clc;
+n=0:1/100:5
+x=sin((7/9)*%pi*(n^2)+1)
+plot(x);
+disp('this shows that signal is NOT periodic');
\ No newline at end of file diff --git a/67/CH1/EX1.49.b/example149b.sce b/67/CH1/EX1.49.b/example149b.sce new file mode 100755 index 000000000..a5ac9be74 --- /dev/null +++ b/67/CH1/EX1.49.b/example149b.sce @@ -0,0 +1,9 @@ +//Example 1.49b
+//Determine whether the signal x(n)=cos(pi*n/2)cos(pi*n/4)
+clc;
+n=0:1/100:100
+x0=cos((%pi*n/2)+(%pi*n/4))
+x1=cos((%pi*n/2)-(%pi*n/4))
+x=(x0+x1)/2;
+plot(x);
+disp('plot shows that this is a periodic signal');
\ No newline at end of file diff --git a/67/CH1/EX1.5/example15.sce b/67/CH1/EX1.5/example15.sce new file mode 100755 index 000000000..1385a0e37 --- /dev/null +++ b/67/CH1/EX1.5/example15.sce @@ -0,0 +1,10 @@ +//Example 1.5
+//Sketch the signal x(t)=A[u(t+a)-u(t-a)]
+clc;
+A=1;
+a=2;
+t=-a:a
+x=1;
+plot(t,x)
+//this signal is a finite duration signal so it is energy signal
+E=integrate('1','t',-a,a);
\ No newline at end of file diff --git a/67/CH1/EX1.6/example16.sce b/67/CH1/EX1.6/example16.sce new file mode 100755 index 000000000..25fef5a15 --- /dev/null +++ b/67/CH1/EX1.6/example16.sce @@ -0,0 +1,8 @@ +//Example 1.6
+//Sketch the signal x(t)=exp(-a*t)
+clc;
+t=0:1/100:10;
+x=exp(-0.5*t);
+plot(x)
+E=integrate('(exp(-0.5*t)^2)','t',0,10)
+//Energy of the signal
\ No newline at end of file diff --git a/67/CH1/EX1.8/example18.sce b/67/CH1/EX1.8/example18.sce new file mode 100755 index 000000000..765ffa26a --- /dev/null +++ b/67/CH1/EX1.8/example18.sce @@ -0,0 +1,13 @@ +//Example 1.8
+//Find the power of the signal x(t)=Acos(Wot+theeta)
+clc;
+A=20;
+Wo=(2*%pi)/4;
+for i=1:50
+ x(i)=A*cos(Wo*i);
+end
+p=0;
+for i=1:4
+ p=p+(abs(x(i)^2))/4;
+end
+disp(p,'The power of the given signal is =');
\ No newline at end of file diff --git a/67/CH2/EX2.1/example21.sce b/67/CH2/EX2.1/example21.sce new file mode 100755 index 000000000..5d6412394 --- /dev/null +++ b/67/CH2/EX2.1/example21.sce @@ -0,0 +1,18 @@ +//Example 2.1
+clc;
+t=-8:1/100:8;
+for i=1:length(t)
+ x(i)=exp(-t(i)^2);
+ h(i)=3*t(i)^2;
+end
+y=convol(x,h);
+figure
+plot2d(t,h);
+title('Impulse responce');
+figure
+plot2d(t,x);
+title('Input signal');
+figure
+t2=-16:1/100:16
+plot2d(t2,y);
+title('Output signal');
\ No newline at end of file diff --git a/67/CH2/EX2.17.a/example217a.sce b/67/CH2/EX2.17.a/example217a.sce new file mode 100755 index 000000000..3407e0d66 --- /dev/null +++ b/67/CH2/EX2.17.a/example217a.sce @@ -0,0 +1,5 @@ +//Example 2.17a
+clc;
+disp(' y[n]=3x[n-2]+3x[n+2] ');
+disp('THE GIVEN SYSTEM IS NON-CAUSAL');
+disp('Since the value of output depends on future input');
\ No newline at end of file diff --git a/67/CH2/EX2.17.b/example217b.sce b/67/CH2/EX2.17.b/example217b.sce new file mode 100755 index 000000000..71a43610f --- /dev/null +++ b/67/CH2/EX2.17.b/example217b.sce @@ -0,0 +1,5 @@ +//Example 2.17b
+clc;
+disp(' y[n]=x[n-1]+a*x[n-2] ');
+disp('THE GIVEN SYSTEM IS CAUSAL');
+disp('Since the value of output doesnot depends on future input');
\ No newline at end of file diff --git a/67/CH2/EX2.17.c/example217c.sce b/67/CH2/EX2.17.c/example217c.sce new file mode 100755 index 000000000..ef8b994c6 --- /dev/null +++ b/67/CH2/EX2.17.c/example217c.sce @@ -0,0 +1,5 @@ +//Example 2.17c
+clc;
+disp(' y[n]=x[-n] ');
+disp('THE GIVEN SYSTEM IS NON-CAUSAL');
+disp('Since the value of output depends on future input');
\ No newline at end of file diff --git a/67/CH2/EX2.19.a/example219a.sce b/67/CH2/EX2.19.a/example219a.sce new file mode 100755 index 000000000..5471fe7f3 --- /dev/null +++ b/67/CH2/EX2.19.a/example219a.sce @@ -0,0 +1,32 @@ +//Example 2.19a
+clc;
+x1=[1,1,1,1]
+x2=[2,2,2,2]
+a=1
+b=1
+a1=0.5
+b1=0.5
+for n=1:length(x1)
+ x3(n)=a*x1(n)+b*x2(n)
+end
+for n=1:length(x1)
+ y1(n)=a1*n*x1(n)+b1
+ y2(n)=a1*n*x2(n)+b1
+ y3(n)=a1*n*x3(n)+b1
+end
+for n=1:length(y1)
+ z(n)=a*y1(n)+b*y2(n)
+end
+count=0
+for n=1:length(y1)
+ if(y3(n)==z(n))
+ count=count+1;
+ end
+end
+if(count==length(y3))
+disp('It satisfy the superposition principle');
+disp('THE GIVEN SYSTEM IS LINEAR ');
+else
+ disp('It does not satisfy superposition principle ');
+ disp('THE GIVEN SYSTEM IS NON LINEAR');
+end
\ No newline at end of file diff --git a/67/CH2/EX2.19.b/example219b.sce b/67/CH2/EX2.19.b/example219b.sce new file mode 100755 index 000000000..f13e4d4e1 --- /dev/null +++ b/67/CH2/EX2.19.b/example219b.sce @@ -0,0 +1,30 @@ +//Example 2.19b
+clc;
+x1=[1,1,1,1]
+x2=[2,2,2,2]
+a=1
+b=1
+for n=1:length(x1)
+ x3(n)=a*x1(n)+b*x2(n)
+end
+for n=1:length(x1)
+ y1(n)=exp(x1(n))
+ y2(n)=exp(x2(n))
+ y3(n)=exp(x3(n))
+end
+for n=1:length(y1)
+ z(n)=a*y1(n)+b*y2(n)
+end
+count=0
+for n=1:length(y1)
+ if(y3(n)==z(n))
+ count=count+1;
+ end
+end
+if(count==length(y3))
+disp('It satisfy the superposition principle');
+disp('THE GIVEN SYSTEM IS LINEAR ');
+else
+ disp('It does not satisfy superposition principle ');
+ disp('THE GIVEN SYSTEM IS NON LINEAR');
+end
\ No newline at end of file diff --git a/67/CH2/EX2.2/example22.sce b/67/CH2/EX2.2/example22.sce new file mode 100755 index 000000000..7c27cc5b5 --- /dev/null +++ b/67/CH2/EX2.2/example22.sce @@ -0,0 +1,24 @@ +//Example 2.2
+clc;
+t=-8:1/100:8;
+for i=1:length(t)
+ if t(i)<0 then
+ x(i)=0;
+ h(i)=0;
+ else
+ x(i)=exp(-3.*t(i));
+ h(i)=1;
+ end
+end
+t1=t+1;
+y=convol(x,h);
+figure
+plot2d(t1,h);
+title('Impulse responce');
+figure
+plot2d(t,x);
+title('Input signal');
+figure
+t2=-16:1/100:16
+plot2d(t2,y);
+title('Output signal');
\ No newline at end of file diff --git a/67/CH2/EX2.21.a/example221a.sce b/67/CH2/EX2.21.a/example221a.sce new file mode 100755 index 000000000..48d1b3db4 --- /dev/null +++ b/67/CH2/EX2.21.a/example221a.sce @@ -0,0 +1,30 @@ +//Example 2.21
+clc;
+x1=[1,1,1,1]
+x2=[2,2,2,2]
+a=1
+b=1
+for t=1:length(x1)
+ x3(t)=a*x1(t)+b*x2(t)
+end
+for t=1:length(x1)
+ y1(t)=5*sin(x1(t))
+ y2(t)=5*sin(x2(t))
+ y3(t)=5*sin(x3(t))
+end
+for t=1:length(y1)
+ z(t)=a*y1(t)+b*y2(t)
+end
+count=0
+for n=1:length(y1)
+ if(y3(t)==z(t))
+ count=count+1;
+ end
+end
+if(count==length(y3))
+disp('It satisfy the superposition principle');
+disp('THE GIVEN SYSTEM IS LINEAR ');
+else
+ disp('It does not satisfy superposition principle ');
+ disp('THE GIVEN SYSTEM IS NON LINEAR');
+end
\ No newline at end of file diff --git a/67/CH2/EX2.21.b/example221b.sce b/67/CH2/EX2.21.b/example221b.sce new file mode 100755 index 000000000..e0660a3b5 --- /dev/null +++ b/67/CH2/EX2.21.b/example221b.sce @@ -0,0 +1,30 @@ +//Example 2.21b
+clc;
+x1=[1,1,1,1]
+x2=[2,2,2,2]
+a=1
+b=1
+for t=1:length(x1)
+ x3(t)=a*x1(t)+b*x2(t)
+end
+for t=1:length(x1)
+ y1(t)=7*x1(t)+5
+ y2(t)=7*x2(t)+5
+ y3(t)=7*x3(t)+5
+end
+for t=1:length(y1)
+ z(t)=a*y1(t)+b*y2(t)
+end
+count=0
+for n=1:length(y1)
+ if(y3(t)==z(t))
+ count=count+1;
+ end
+end
+if(count==length(y3))
+disp('It satisfy the superposition principle');
+disp('THE GIVEN SYSTEM IS LINEAR ');
+else
+ disp('It does not satisfy superposition principle ');
+ disp('THE GIVEN SYSTEM IS NON LINEAR');
+end
\ No newline at end of file diff --git a/67/CH2/EX2.25/example225.sce b/67/CH2/EX2.25/example225.sce new file mode 100755 index 000000000..e4215c392 --- /dev/null +++ b/67/CH2/EX2.25/example225.sce @@ -0,0 +1,30 @@ +//Example 2.25
+clc;
+x1=[1,1,1,1]
+x2=[2,2,2,2]
+a=1
+b=1
+for n=1:length(x1)
+ x3(n)=a*x1(n)+b*x2(n)
+end
+for n=1:length(x1)
+ y1(n)=x1(n)^2
+ y2(n)=x2(n)^2
+ y3(n)=x3(n)^2
+end
+for n=1:length(y1)
+ z(n)=a*y1(n)+b*y2(n)
+end
+count=0
+for n=1:length(y1)
+ if(y3(n)==z(n))
+ count=count+1;
+ end
+end
+if(count==length(y3))
+disp('It satisfy the superposition principle');
+disp('THE GIVEN SYSTEM IS LINEAR ');
+else
+ disp('It does not satisfy superposition principle ');
+ disp('THE GIVEN SYSTEM IS NON LINEAR');
+end
\ No newline at end of file diff --git a/67/CH2/EX2.3/example23.sce b/67/CH2/EX2.3/example23.sce new file mode 100755 index 000000000..b6bc33c30 --- /dev/null +++ b/67/CH2/EX2.3/example23.sce @@ -0,0 +1,25 @@ +//Example 2.3
+clc;
+R=100;
+L=100;
+t=-8:1/100:8;
+for i=1:length(t)
+ if t(i)<0 then
+ x(i)=0;
+ h(i)=0;
+ else
+ h(i)=(R/L)*exp(-(R/L).*t(i));
+ x(i)=1;
+ end
+end
+y=convol(x,h);
+figure
+plot2d(t,h);
+title('Impulse responce');
+figure
+plot2d(t,x);
+title('Input signal');
+figure
+t2=-16:1/100:16
+plot2d(t2,y);
+title('Output signal');
\ No newline at end of file diff --git a/67/CH2/EX2.4/example24.sce b/67/CH2/EX2.4/example24.sce new file mode 100755 index 000000000..375b871e3 --- /dev/null +++ b/67/CH2/EX2.4/example24.sce @@ -0,0 +1,18 @@ +//Example 2.4
+clc;
+t=-8:1/100:8;
+for i=1:length(t)
+ x(i)=3*cos(2.*t(i));
+ h(i)=exp(-abs(t(i)));
+end
+y=convol(x,h);
+figure
+plot2d(t,h);
+title('Impulse responce');
+figure
+plot2d(t,x);
+title('Input signal');
+figure
+t2=-16:1/100:16
+plot2d(t2,y);
+title('Output signal');
\ No newline at end of file diff --git a/67/CH2/EX2.5/example25.sce b/67/CH2/EX2.5/example25.sce new file mode 100755 index 000000000..05b0eadc4 --- /dev/null +++ b/67/CH2/EX2.5/example25.sce @@ -0,0 +1,27 @@ +//Example 2.5 +clc; +Max_Limit=10; +h=ones(1,Max_Limit); +N2=0:length(h)-1; +a=0.5;//constant a>0 +for t=1:Max_Limit +x(t)=exp(-a*(t-1)); +end +N1=0:length(x)-1; +y=convol(x,h)-1; +N=0:length(x)+length(h)-2; +figure +a=gca(); +plot2d(N2,h) +xtitle('Impulse Response','t','h(t)'); +a.thickness=2; +figure +a=gca(); +plot2d(N1,x) +xtitle('Input Response','t','x(t)'); +a.thickness=2; +figure +a=gca (); +plot2d(N(1:Max_Limit),y(1:Max_Limit)) +xtitle('Output Response','t','y(t)'); +a.thickness=2;
\ No newline at end of file diff --git a/67/CH2/EX2.59/example259.sce b/67/CH2/EX2.59/example259.sce new file mode 100755 index 000000000..2661f9b75 --- /dev/null +++ b/67/CH2/EX2.59/example259.sce @@ -0,0 +1,32 @@ +//Example 2.59
+clc;
+n=-8:1:8;
+for i=1:length(n)
+ if n(i)<0 then
+ x(i)=0;
+ h(i)=0;
+ else
+ x(i)=1;
+ h(i)=2^n(i);
+ end
+end
+y=convol(x,h);
+figure
+a=gca();
+plot2d3(n,h);
+a.x_location='origin';
+a.y_location='origin';
+title('Impulse responce');
+figure
+a=gca();
+plot2d3(n,x);
+a.x_location='origin';
+a.y_location='origin';
+title('Input signal');
+figure
+a=gca();
+n1=-16:1:16
+plot2d3(n1,y);
+a.x_location='origin';
+a.y_location='origin';
+title('Output signal');
\ No newline at end of file diff --git a/67/CH2/EX2.6/example26.sce b/67/CH2/EX2.6/example26.sce new file mode 100755 index 000000000..c698b7ffa --- /dev/null +++ b/67/CH2/EX2.6/example26.sce @@ -0,0 +1,24 @@ +//Example 2.6
+clc;
+t=-8:1/100:8;
+for i=1:length(t)
+ if t(i)<0 then
+ x(i)=exp(2.*t(i));
+ h(i)=0;
+ else
+ x(i)=0;
+ h(i)=1;
+ end
+end
+t1=t+3;
+y=convol(x,h);
+figure
+plot2d(t1,h);
+title('Impulse responce');
+figure
+plot2d(t,x);
+title('Input signal');
+figure
+t2=-16:1/100:16
+plot2d(t2,y);
+title('Output signal');
\ No newline at end of file diff --git a/67/CH2/EX2.7/example27.sce b/67/CH2/EX2.7/example27.sce new file mode 100755 index 000000000..9014344d7 --- /dev/null +++ b/67/CH2/EX2.7/example27.sce @@ -0,0 +1,18 @@ +//Example 2.7
+clc;
+n=-8:1:8;
+for i=1:length(n)
+ x(i)=exp(-n(i)^2);
+ h(i)=3.*n(i)^2;
+end
+y=convol(x,h);
+figure
+plot2d3(n,h);
+title('Impulse responce');
+figure
+plot2d3(n,x);
+title('Input signal');
+figure
+n1=-16:1:16
+plot2d3(n1,y);
+title('Output signal');
\ No newline at end of file diff --git a/67/CH2/EX2.8/example28.sce b/67/CH2/EX2.8/example28.sce new file mode 100755 index 000000000..8e348364a --- /dev/null +++ b/67/CH2/EX2.8/example28.sce @@ -0,0 +1,23 @@ +//Example 2.8
+clc;
+n=-8:1:8;
+for i=1:length(n)
+ if n(i)<0 then
+ x(i)=2^n(i);
+ h(i)=0;
+ else
+ x(i)=0;
+ h(i)=1;
+ end
+end
+y=convol(x,h);
+figure
+plot2d3(n,h);
+title('Impulse responce');
+figure
+plot2d3(n,x);
+title('Input signal');
+figure
+n1=-16:1:16
+plot2d3(n1,y);
+title('Output signal');
\ No newline at end of file diff --git a/67/CH3/EX3.15/example315.sce b/67/CH3/EX3.15/example315.sce new file mode 100755 index 000000000..094c51ecf --- /dev/null +++ b/67/CH3/EX3.15/example315.sce @@ -0,0 +1,20 @@ +//Example 3.15
+clc;
+clear;
+T1=2;
+T=4*T1;
+Wo=2*%pi/T;
+W=[-Wo,0,Wo];
+ak=(2*%pi*Wo*T1/%pi)/sqrt(-1);
+XW=[-ak,0,ak];
+ak1=(2*%pi*Wo*T1/%pi);
+XW1=[ak1,0,ak1];
+figure
+a=gca();
+a.y_location="origin";
+a.x_location="origin";
+plot2d3('gnn',W,XW1,2);
+poly1=a.children(1).children(1);
+poly1.thickness=3;
+xlabel('W');
+title('CTFT of cos(Wot)');
\ No newline at end of file diff --git a/67/CH3/EX3.31/example331.sce b/67/CH3/EX3.31/example331.sce new file mode 100755 index 000000000..56b3b546a --- /dev/null +++ b/67/CH3/EX3.31/example331.sce @@ -0,0 +1,42 @@ +//Example 3.31
+clc ;
+clear;
+R=10^3;
+C=10^-3;
+A=1/(R*C);
+Dt=0.005;
+t=0:Dt:10;
+xt=A*exp(-A*t);
+Wmax=2*%pi*1;
+K=4;
+k=0:(K/1000):K;
+W=k*Wmax/K;
+XW=xt*exp(-sqrt(-1)*t'*W)*Dt;
+XW_Mag=abs(XW);
+W=[-mtlb_fliplr(W),W(2:1001)];
+XW_Mag=[mtlb_fliplr(XW_Mag),XW_Mag(2:1001)];
+[XW_Phase,db]=phasemag(XW);
+XW_Phase=[-mtlb_fliplr(XW_Phase),XW_Phase(2:1001)];
+figure
+a=gca();
+a.y_location="origin";
+plot(t,xt);
+xlabel('t in sec.');
+ylabel('x(t)');
+title('Continuous Time Signal');
+figure
+subplot(2,1,1);
+a=gca();
+a.y_location="origin";
+plot(W,XW_Mag);
+xlabel('Frequency in Radians/Seconds>W');
+ylabel('abs(X(jW))');
+title('Magnitude Response (CTFT)');
+subplot(2,1,2);
+a=gca();
+a.y_location="origin";
+a.x_location="origin";
+plot(W,XW_Phase*%pi/180) ;
+xlabel(' Frequency in Radians/Seconds> W' );
+ylabel('<X(jW)')
+title('Phase Response (CTFT) in Radians');
\ No newline at end of file diff --git a/67/CH3/EX3.8/example38.sce b/67/CH3/EX3.8/example38.sce new file mode 100755 index 000000000..78121960d --- /dev/null +++ b/67/CH3/EX3.8/example38.sce @@ -0,0 +1,43 @@ +clc ;
+close ;
+// Analog S i g n a l
+A =1; // Ampl i tude
+Dt = 0.005;
+t = 0: Dt :10;
+xt = exp(-A*t);
+Wmax = 2* %pi *1; // Analog Fr equency = 1Hz
+K = 4;
+k = 0:( K /1000) :K;
+W = k* Wmax /K;
+XW = xt* exp (- sqrt ( -1)*t'*W) * Dt;
+XW_Mag = abs(XW);
+W = [-mtlb_fliplr(W),W(2:1001)]; // Omega f rom Wmax to Wmax
+XW_Mag = [mtlb_fliplr(XW_Mag),XW_Mag(2:1001)];
+[ XW_Phase ,db] = phasemag (XW);
+XW_Phase =[-mtlb_fliplr(XW_Phase),XW_Phase(2:1001)];
+//Plotting Continuous Time Signal
+figure
+a = gca ();
+a.y_location ="origin";
+plot (t,xt);
+xlabel ( ' t in sec. ' );
+ylabel ( ' x ( t ) ' )
+title ( ' Continuous Time Signal ' )
+figure
+// Pl o t t i n g Magni tude Re spons e o f CTS
+subplot (2 ,1 ,1);
+a = gca ();
+a.y_location ="origin";
+plot (W, XW_Mag );
+xlabel ( ' Fr equency i n Radians / Seconds> W' );
+ylabel ( ' abs (X(jW) ) ' )
+title ( 'Magni tude Re spons e (CTFT) ' )
+// Pl o t t i n g Phase Reponse o f CTS
+subplot (2 ,1 ,2);
+a = gca ();
+a.y_location = "origin";
+a.x_location = "origin";
+plot (W, XW_Phase *%pi /180) ;
+xlabel ( ' Fr equency in Radians / Seconds> W' );
+ylabel ( '<X(jW) ' )
+title ( ' Phase Re spons e (CTFT) i n Radians ' )
\ No newline at end of file diff --git a/67/CH3/EX3.9/example39.sce b/67/CH3/EX3.9/example39.sce new file mode 100755 index 000000000..5aed6e28d --- /dev/null +++ b/67/CH3/EX3.9/example39.sce @@ -0,0 +1,32 @@ +//Example 3.9
+clc;
+clear;
+A=1;
+Dt=0.005;
+T1=4;
+t=-T1/2:Dt:T1/2;
+for i=1:length(t)
+xt(i)=A;
+end
+Wmax=2*%pi*1;
+K=4;
+k=0:(K/1000):K;
+W=k*Wmax/K;
+xt=xt';
+XW=xt*exp(-sqrt(-1)*t'*W)*Dt;
+XW_Mag=real(XW);
+W=[-mtlb_fliplr(W),W(2:1001)];
+XW_Mag=[mtlb_fliplr(XW_Mag),XW_Mag(2:1001)];
+subplot(2,1,1);
+a=gca();
+a.data_bounds=[-4,0;4,2];
+a.y_location="origin";
+plot(t,xt);
+xlabel('t in sec.');
+title('Continous Time Signal x(t)');
+subplot(2,1,2);
+a=gca();
+a.y_location="origin";
+plot(W,XW_Mag);
+xlabel('Frequency in Radians/Seconds');
+title('Continuous time Fourier Transform X(jW)');
\ No newline at end of file diff --git a/67/CH4/EX4.1/example41.sce b/67/CH4/EX4.1/example41.sce new file mode 100755 index 000000000..114edcfa9 --- /dev/null +++ b/67/CH4/EX4.1/example41.sce @@ -0,0 +1,6 @@ +//Example 4.1
+//Find the DTFT of (a^n)u[n],for |a|<1
+clc;
+syms w a n;
+x=a^n;
+X=symsum(x*exp(-%i*w*n),n,0,%inf);
\ No newline at end of file diff --git a/67/CH4/EX4.10/example410.sce b/67/CH4/EX4.10/example410.sce new file mode 100755 index 000000000..71c39fc42 --- /dev/null +++ b/67/CH4/EX4.10/example410.sce @@ -0,0 +1,9 @@ +//Example 4.10
+//Find DTFT of x[n]=cos(Won) with Wo=(2*pi/5)
+clc;
+syms w n;
+x1=exp(%i*(2*%pi*n/5));
+x2=exp(-%i*(2*%pi*n/5));
+X1=symsum(x1*exp(-%i*w*n),n,0,%inf);
+X2=symsum(x2*exp(-%i*w*n),n,0,%inf);
+X=(X1+X2)/2;
\ No newline at end of file diff --git a/67/CH4/EX4.12/example412.sce b/67/CH4/EX4.12/example412.sce new file mode 100755 index 000000000..ceddd6d2f --- /dev/null +++ b/67/CH4/EX4.12/example412.sce @@ -0,0 +1,6 @@ +//Example 4.12
+//Find DTFT of x[n]=u[n]
+clc;
+syms w n;
+x=1;
+X=symsum(x*exp(-%i*w*n),n,0,%inf);
\ No newline at end of file diff --git a/67/CH4/EX4.16/example416.sce b/67/CH4/EX4.16/example416.sce new file mode 100755 index 000000000..ede67c797 --- /dev/null +++ b/67/CH4/EX4.16/example416.sce @@ -0,0 +1,6 @@ +//Example 4.16
+//Find DTFT of x[n]=(a^n)u[n],for 0<a<1
+clc;
+syms w a n;
+x=a^n;
+X=symsum(x*exp(-%i*w*n),n,0,%inf);
\ No newline at end of file diff --git a/67/CH4/EX4.2/example42.sce b/67/CH4/EX4.2/example42.sce new file mode 100755 index 000000000..a71bf5e84 --- /dev/null +++ b/67/CH4/EX4.2/example42.sce @@ -0,0 +1,6 @@ +//Example 4.2
+//Find DTFT of x[n]=(a^n)u[-(n+1)]
+clc;
+syms w a n;
+x=a^n;
+X=symsum(x*exp(-%i*w*n),n,-%inf,-1);
\ No newline at end of file diff --git a/67/CH4/EX4.22/example422.sce b/67/CH4/EX4.22/example422.sce new file mode 100755 index 000000000..32f85afd2 --- /dev/null +++ b/67/CH4/EX4.22/example422.sce @@ -0,0 +1,6 @@ +//Example 4.22
+//Find DTFT of x[n]=((1/2)^(n-1))u[n-1]
+clc;
+syms w n;
+x=(1/2)^(n-1);
+X=symsum(x*exp(-%i*w*n),n,1,%inf);
\ No newline at end of file diff --git a/67/CH4/EX4.3/example43.sce b/67/CH4/EX4.3/example43.sce new file mode 100755 index 000000000..06dbbdf17 --- /dev/null +++ b/67/CH4/EX4.3/example43.sce @@ -0,0 +1,6 @@ +//Example 4.3
+//Find DTFT of unit impluse
+clc;
+syms w n;
+x=1;
+X=symsum(x*exp(-%i*w*n),n,0,0);
\ No newline at end of file diff --git a/67/CH4/EX4.5/example45.sce b/67/CH4/EX4.5/example45.sce new file mode 100755 index 000000000..cb74a3fbc --- /dev/null +++ b/67/CH4/EX4.5/example45.sce @@ -0,0 +1,10 @@ +//Example 4.5
+//Find DTFT of x[n]=a^|n| for -1<a<1
+clc;
+syms w a n;
+x1=a^n;
+x2=a^(-n);
+X1=symsum(x1*exp(-%i*w*n),n,0,%inf);
+X2=symsum(x2*exp(-%i*w*n),n,-%inf,-1);
+X=X1+X2;
+disp(X,'X(e^jw)=');
\ No newline at end of file diff --git a/67/CH4/EX4.9/example49.sce b/67/CH4/EX4.9/example49.sce new file mode 100755 index 000000000..fd4317be9 --- /dev/null +++ b/67/CH4/EX4.9/example49.sce @@ -0,0 +1,12 @@ +//Example 4.9
+clc;
+syms w a n;
+x=a^n;
+pi=22/7;
+X=symsum(x*exp(-%i*w*n),n,0,%inf);
+n1=0:10;
+a=0.5;
+x1=a^n1;
+plot2d3(n1,x1);
+xtitle('Discrete Time Signal','n','x[n]');
+a.thickness=2;
\ No newline at end of file diff --git a/67/CH5/EX5.1/example51.sce b/67/CH5/EX5.1/example51.sce new file mode 100755 index 000000000..96efb1bfe --- /dev/null +++ b/67/CH5/EX5.1/example51.sce @@ -0,0 +1,7 @@ +//Example 5.1
+//Obtain the Bode plot
+clc;
+s=%s;
+H=syslin('c',2*10^4/(s^2+100*s+10^4));
+bode(H,0.1,10000);
+funcprot(0);
\ No newline at end of file diff --git a/67/CH5/EX5.2/example52.sce b/67/CH5/EX5.2/example52.sce new file mode 100755 index 000000000..b0ad39a86 --- /dev/null +++ b/67/CH5/EX5.2/example52.sce @@ -0,0 +1,6 @@ +//Example 5.2
+//Obtain the Bode plot
+clc;
+s=%s;
+H=syslin('c',100*(1+s)/((10+s)*(100+s)));
+bode(H,0.01,2000);
\ No newline at end of file diff --git a/67/CH6/EX6.1/example61.sce b/67/CH6/EX6.1/example61.sce new file mode 100755 index 000000000..5db3959ff --- /dev/null +++ b/67/CH6/EX6.1/example61.sce @@ -0,0 +1,22 @@ +//Example 6.1
+clc;
+disp('x(t)=3cos(50%pi*t)+10sin(300%pi*t)-cos(100%pi*t)');
+w1=50*%pi;
+w2=300*%pi;
+w3=100*%pi;
+f1=w1/(2*%pi);
+f2=w2/(2*%pi);
+f3=w3/(2*%pi);
+if f1>f2 then
+ if f1>f3 then
+ disp(2*f1,'Nyquist rate=');
+ else
+ disp(2*f3,'Nyquist rate=');
+ end
+else
+ if f2>f3 then
+ disp(2*f2,'Nyquist rate=');
+ else
+ disp(2*f3,'Nyquist rate=');
+ end
+end
\ No newline at end of file diff --git a/67/CH6/EX6.2/example62.sce b/67/CH6/EX6.2/example62.sce new file mode 100755 index 000000000..1d4dc07ee --- /dev/null +++ b/67/CH6/EX6.2/example62.sce @@ -0,0 +1,15 @@ +//Example 6.2
+clc;
+disp('x(t)=(1/2%pi)cos(4000%pi*t)cos(1000%pi*t)');
+w1=5000*%pi;
+w2=3000*%pi;
+f1=w1/(2*%pi);
+f2=w2/(2*%pi);
+if f1>f2 then
+ nyquist_rate=2*f1;
+else
+ nyquist_rate=2*f2;
+end
+nyquist_interval=1/nyquist_rate;
+disp(nyquist_rate,'Nyquist rate=');
+disp(nyquist_interval,'Nyquist interval in seconds');
\ No newline at end of file diff --git a/67/CH6/EX6.26/example626.sce b/67/CH6/EX6.26/example626.sce new file mode 100755 index 000000000..779034aaf --- /dev/null +++ b/67/CH6/EX6.26/example626.sce @@ -0,0 +1,6 @@ +//Example 6.26
+//Find laplace transform x(t)=2e^(-3t)u(t)-e^(-2t)u(t)
+clc;
+syms t;
+x=2*%e^(-3*t)-%e^(-2*t);
+X=laplace(x);
\ No newline at end of file diff --git a/67/CH6/EX6.27.a/example627a.sce b/67/CH6/EX6.27.a/example627a.sce new file mode 100755 index 000000000..4ab352365 --- /dev/null +++ b/67/CH6/EX6.27.a/example627a.sce @@ -0,0 +1,6 @@ +//Example 6.27a
+//Laplace transform of x(t)=t^3+3*t^2-6*t+4
+clc;
+syms t;
+x=t^3+3*t^2-6*t+4;
+X=laplace(x);
\ No newline at end of file diff --git a/67/CH6/EX6.27.b/example627b.sce b/67/CH6/EX6.27.b/example627b.sce new file mode 100755 index 000000000..580c86850 --- /dev/null +++ b/67/CH6/EX6.27.b/example627b.sce @@ -0,0 +1,6 @@ +//Example 6.27b
+//x(t)=(cos(3t))^3
+clc;
+syms t;
+x=(cos(3*t))^3;
+X=laplace(x);
\ No newline at end of file diff --git a/67/CH6/EX6.27.c/example627c.sce b/67/CH6/EX6.27.c/example627c.sce new file mode 100755 index 000000000..a1a32beab --- /dev/null +++ b/67/CH6/EX6.27.c/example627c.sce @@ -0,0 +1,5 @@ +//Example 6.27c
+clc;
+syms a b t;
+x=sin(a*t)*cos(b*t);
+X=laplace(x);
\ No newline at end of file diff --git a/67/CH6/EX6.27.d/example627d.sce b/67/CH6/EX6.27.d/example627d.sce new file mode 100755 index 000000000..b41d846fb --- /dev/null +++ b/67/CH6/EX6.27.d/example627d.sce @@ -0,0 +1,5 @@ +//Example 6.27d
+clc;
+syms t a;
+x=t*sin(a*t);
+X=laplace(x);
\ No newline at end of file diff --git a/67/CH6/EX6.27.e/example627e.sce b/67/CH6/EX6.27.e/example627e.sce new file mode 100755 index 000000000..88fdaa2a3 --- /dev/null +++ b/67/CH6/EX6.27.e/example627e.sce @@ -0,0 +1,6 @@ +//Example 6.27e
+clc;
+syms t s;
+x1=1-%e^t;
+X1=laplace(x1);
+X=integ(X1,s,s,%inf);
\ No newline at end of file diff --git a/67/CH6/EX6.4/example64.sce b/67/CH6/EX6.4/example64.sce new file mode 100755 index 000000000..5b5aa88ee --- /dev/null +++ b/67/CH6/EX6.4/example64.sce @@ -0,0 +1,22 @@ +//Example 6.4
+clc;
+disp('x(t)=6cos(50%pi*t)+20sin(300%pi*t)-10cos(100%pi*t)');
+w1=50*%pi;
+w2=300*%pi;
+w3=100*%pi;
+f1=w1/(2*%pi);
+f2=w2/(2*%pi);
+f3=w3/(2*%pi);
+if f1>f2 then
+ if f1>f3 then
+ disp(2*f1,'Nyquist rate=');
+ else
+ disp(2*f3,'Nyquist rate=');
+ end
+else
+ if f2>f3 then
+ disp(2*f2,'Nyquist rate=');
+ else
+ disp(2*f3,'Nyquist rate=');
+ end
+end
\ No newline at end of file diff --git a/67/CH6/EX6.48/example648.sce b/67/CH6/EX6.48/example648.sce new file mode 100755 index 000000000..e96643028 --- /dev/null +++ b/67/CH6/EX6.48/example648.sce @@ -0,0 +1,8 @@ +//Example 6.48
+clc;
+syms t s;
+x=1+%e^(-3*t)-%e^(-t);
+X=laplace(x);
+H=1/((s+1)*(s^2+s+1));
+Y=X*H;
+y=ilaplace(Y);
\ No newline at end of file diff --git a/67/CH7/EX7.1/example71.sce b/67/CH7/EX7.1/example71.sce new file mode 100755 index 000000000..1a471ad25 --- /dev/null +++ b/67/CH7/EX7.1/example71.sce @@ -0,0 +1,6 @@ +//Example 7.1
+clc;
+syms a z n;
+x=a^n;
+X=symsum(x*(z^-n),n,0,%inf);
+disp(X,'X(z)=');
\ No newline at end of file diff --git a/67/CH7/EX7.10/example710.sce b/67/CH7/EX7.10/example710.sce new file mode 100755 index 000000000..53b4ba2f3 --- /dev/null +++ b/67/CH7/EX7.10/example710.sce @@ -0,0 +1,7 @@ +//Example 7.10
+//Z-transform of (2^n)u[n-2]
+clc;
+syms n z;
+x=2^n;
+X=symsum(x*(z^-n),n,2,%inf);
+disp(X,'X(z)=');
\ No newline at end of file diff --git a/67/CH7/EX7.11.a/example711a.sce b/67/CH7/EX7.11.a/example711a.sce new file mode 100755 index 000000000..45927fe7c --- /dev/null +++ b/67/CH7/EX7.11.a/example711a.sce @@ -0,0 +1,10 @@ +//Example 7.11a
+//Z transform of (a^n)cos(Wo*n)
+clc;
+syms Wo n z a;
+x1=(a^n)*exp(sqrt(-1)*Wo*n);
+X1=symsum(x1*(z^-n),n,0,%inf);
+x2=(a^n)*exp(-sqrt(-1)*Wo*n);
+X2=symsum(x2*(z^-n),n,0,%inf);
+X=(X1+X2)/2;
+disp(X,'X(z)=');
\ No newline at end of file diff --git a/67/CH7/EX7.11.b/example711b.sce b/67/CH7/EX7.11.b/example711b.sce new file mode 100755 index 000000000..2df600250 --- /dev/null +++ b/67/CH7/EX7.11.b/example711b.sce @@ -0,0 +1,10 @@ +//Example 7.11b
+//Z transform of (a^n)sin(Wo*n)
+clc;
+syms Wo n z a;
+x1=(a^n)*exp(sqrt(-1)*Wo*n);
+X1=symsum(x1*(z^-n),n,0,%inf);
+x2=(a^n)*exp(-sqrt(-1)*Wo*n);
+X2=symsum(x2*(z^-n),n,0,%inf);
+X=(X1-X2)/(2*%i);
+disp(X,'X(z)=');
\ No newline at end of file diff --git a/67/CH7/EX7.12/example712.sce b/67/CH7/EX7.12/example712.sce new file mode 100755 index 000000000..d14d7b7bf --- /dev/null +++ b/67/CH7/EX7.12/example712.sce @@ -0,0 +1,9 @@ +//Example 7.12
+//Ztransform of x[n]=(n^2)*u[n] done by Diffrentiation property
+clc;
+syms z n;
+x=1;
+X1=symsum(x*(z^-n),n,0,%inf);
+X2=(-z)*(diff(X1,z));
+X=(-z)*(diff(X2,z));
+disp(X,'X(z)=');
\ No newline at end of file diff --git a/67/CH7/EX7.13/example713.sce b/67/CH7/EX7.13/example713.sce new file mode 100755 index 000000000..a5d221976 --- /dev/null +++ b/67/CH7/EX7.13/example713.sce @@ -0,0 +1,11 @@ +//Example 7.13
+//Diffrentiation property is used here
+clc;
+syms n z;
+x1=((-1/2)^n);
+x2=(1/4)^-n;
+X1=symsum(x1*(z^-n),n,0,%inf);
+X3=(-z)*diff(X1,z);
+X2=symsum(x2*(z^-n),n,-%inf,0);
+X=X3*X2;
+disp(X,'X(z)=');
\ No newline at end of file diff --git a/67/CH7/EX7.18/example718.sce b/67/CH7/EX7.18/example718.sce new file mode 100755 index 000000000..9edff0f02 --- /dev/null +++ b/67/CH7/EX7.18/example718.sce @@ -0,0 +1,20 @@ +//Example 7.18
+//Determine the input x[n] if h[n]=[1,2,3]and y[n]=[1,1,2,-1,3]
+clc;
+clear;
+function[za]=ztransfer(sequence,n)
+ z=poly(0,'z','r')
+ za=sequence*(1/z)^n'
+endfunction
+z=poly(0,'z');
+h=[1,2,3];
+n1=0:length(h)-1;
+H=ztransfer(h,n1);
+y=[1,1,2,-1,3];
+n2=0:length(y)-1;
+Y=ztransfer(y,n2);
+X=Y/H;
+funcprot(0);
+funcprot(0);
+x=ldiv(z^2-z+1,z^2,3);
+disp(x,'x[n]=');
\ No newline at end of file diff --git a/67/CH7/EX7.19/example719.sce b/67/CH7/EX7.19/example719.sce new file mode 100755 index 000000000..aefe77b58 --- /dev/null +++ b/67/CH7/EX7.19/example719.sce @@ -0,0 +1,7 @@ +//Example 7.19
+//Find the inverse Z-transform using long division method
+clc;
+clear;
+z=poly(0,'z');
+x=ldiv(z^2,z^2-(3/2)*z+(1/2),4);
+disp(x,'x[n]=');
\ No newline at end of file diff --git a/67/CH7/EX7.2/example72.sce b/67/CH7/EX7.2/example72.sce new file mode 100755 index 000000000..e3fb3f75a --- /dev/null +++ b/67/CH7/EX7.2/example72.sce @@ -0,0 +1,6 @@ +//Example 7.2
+clc;
+syms n z;
+x=1;
+X=symsum(x*(z^-n),n,0,0);
+disp(X,'X(z)=');
\ No newline at end of file diff --git a/67/CH7/EX7.24.a/example724a.sce b/67/CH7/EX7.24.a/example724a.sce new file mode 100755 index 000000000..96207561a --- /dev/null +++ b/67/CH7/EX7.24.a/example724a.sce @@ -0,0 +1,7 @@ +//Example 7.24a
+//Inverse Z-transform using long division method
+clc;
+clear;
+z=poly(0,'z');
+x=ldiv(2*z^3+3*z^2,(z+1)*(z+0.5)*(z-0.25),4);
+disp(x,'x[n]=');
\ No newline at end of file diff --git a/67/CH7/EX7.3/example73.sce b/67/CH7/EX7.3/example73.sce new file mode 100755 index 000000000..54a9d4699 --- /dev/null +++ b/67/CH7/EX7.3/example73.sce @@ -0,0 +1,5 @@ +clc;
+syms n;
+x=ones(1);
+X=symsum(x*(z^-n),n,0,%inf);
+disp(X,'X(z)=');
\ No newline at end of file diff --git a/67/CH7/EX7.36/example736.sce b/67/CH7/EX7.36/example736.sce new file mode 100755 index 000000000..6feabd184 --- /dev/null +++ b/67/CH7/EX7.36/example736.sce @@ -0,0 +1,8 @@ +//Example 7.36
+clc;
+syms z n;
+x1=2^n;
+X1=symsum(x1*(z^-n),n,0,%inf);
+x2=3^n;
+X2=symsum(x2*(z^-n),n,0,%inf);
+X=3*X1-4*X2;
\ No newline at end of file diff --git a/67/CH7/EX7.37/example737.sce b/67/CH7/EX7.37/example737.sce new file mode 100755 index 000000000..dffbf3118 --- /dev/null +++ b/67/CH7/EX7.37/example737.sce @@ -0,0 +1,5 @@ +//Example 7.37
+clc;
+syms z n;
+x=(1/2)^n;
+X=symsum(x*(z^-n),n,0,%inf);
\ No newline at end of file diff --git a/67/CH7/EX7.38/example738.sce b/67/CH7/EX7.38/example738.sce new file mode 100755 index 000000000..3a26fe714 --- /dev/null +++ b/67/CH7/EX7.38/example738.sce @@ -0,0 +1,5 @@ +//Example 7.38
+clc;
+syms a z n;
+x=-(a^n);
+X=symsum(x*(z^-n),n,-%inf,-1);
\ No newline at end of file diff --git a/67/CH7/EX7.39/example739.sce b/67/CH7/EX7.39/example739.sce new file mode 100755 index 000000000..4eca0d008 --- /dev/null +++ b/67/CH7/EX7.39/example739.sce @@ -0,0 +1,6 @@ +//Example 7.39
+clc;
+syms z n a;
+x1=(a^n);
+X1=symsum(x1*(z^-n),n,0,%inf);
+X=(-z)*(diff(X1,z));
\ No newline at end of file diff --git a/67/CH7/EX7.42.a/example742a.sce b/67/CH7/EX7.42.a/example742a.sce new file mode 100755 index 000000000..d9a394a6f --- /dev/null +++ b/67/CH7/EX7.42.a/example742a.sce @@ -0,0 +1,10 @@ +//Example 7.42a
+clc;
+function[za]=ztransfer(sequence,n)
+ z=poly(0,'z','r')
+ za=sequence*(1/z)^n'
+endfunction
+x=[1,2,3,4,5,0,7];
+n1=0:length(x)-1;
+X=ztransfer(x,n1);
+funcprot(0);
\ No newline at end of file diff --git a/67/CH7/EX7.42.b/example742b.sce b/67/CH7/EX7.42.b/example742b.sce new file mode 100755 index 000000000..bed0b83a9 --- /dev/null +++ b/67/CH7/EX7.42.b/example742b.sce @@ -0,0 +1,10 @@ +//Example 7.42b
+clc;
+function[za]=ztransfer(sequence,n)
+ z=poly(0,'z','r')
+ za=sequence*(1/z)^n'
+endfunction
+x=[1,2,3,4,5,0,7];
+n1=-3:length(x)-4;
+X=ztransfer(x,n1);
+funcprot(0);
\ No newline at end of file diff --git a/67/CH7/EX7.43/example743.sce b/67/CH7/EX7.43/example743.sce new file mode 100755 index 000000000..1fabb7985 --- /dev/null +++ b/67/CH7/EX7.43/example743.sce @@ -0,0 +1,8 @@ +//Example 7.43
+clc;
+syms z n;
+x1=(-1/3)^n;
+x2=(1/2)^n;
+X1=symsum(x1*(z^-n),n,0,%inf);
+X2=symsum(x2*(z^-n),n,-%inf,-1);
+X=X1-X2;
\ No newline at end of file diff --git a/67/CH7/EX7.48.a/example748a.sce b/67/CH7/EX7.48.a/example748a.sce new file mode 100755 index 000000000..d9fab14af --- /dev/null +++ b/67/CH7/EX7.48.a/example748a.sce @@ -0,0 +1,8 @@ +//Example 7.48a
+clc;
+syms z n;
+x1=2^n;
+X1=symsum(x1*(z^-n),n,0,%inf);
+x2=(1/2)^n;
+X2=symsum(x2*(z^-n),n,0,%inf);
+X=X1+(3*X2);
\ No newline at end of file diff --git a/67/CH7/EX7.48.b/example748b.sce b/67/CH7/EX7.48.b/example748b.sce new file mode 100755 index 000000000..2035aa56e --- /dev/null +++ b/67/CH7/EX7.48.b/example748b.sce @@ -0,0 +1,8 @@ +//Example 7.48b
+clc;
+syms z n;
+x1=(-1/2)^n;
+X1=symsum(x1*(z^-n),n,0,%inf);
+x2=(3)^n;
+X2=symsum(x2*(z^-n),n,-%inf,-1);
+X=(3*X1)-(2*X2);
\ No newline at end of file diff --git a/67/CH7/EX7.5/example75.sce b/67/CH7/EX7.5/example75.sce new file mode 100755 index 000000000..46d3707d6 --- /dev/null +++ b/67/CH7/EX7.5/example75.sce @@ -0,0 +1,9 @@ +//Example 7.5
+clc;
+syms Wo n z;
+x1=exp(sqrt(-1)*Wo*n);
+X1=symsum(x1*(z^-n),n,0,%inf);
+x2=exp(-sqrt(-1)*Wo*n);
+X2=symsum(x2*(z^-n),n,0,%inf);
+X=(X1+X2)/2;
+disp(X,'X(z)=');
\ No newline at end of file diff --git a/67/CH7/EX7.50/example750.sce b/67/CH7/EX7.50/example750.sce new file mode 100755 index 000000000..7fa77ad27 --- /dev/null +++ b/67/CH7/EX7.50/example750.sce @@ -0,0 +1,16 @@ +//Example 7.50 +clc; +syms n z; +X1=0; +X2=0; +for i=0:2:4 + x1=(1/2)^i; + X1=X1+x1*z^-i; +end +for i=1:2:5 + x2=(1/3)^i; + X2=X2+x2*z^-i; +end +x3=(2)^n; +X3=symsum(x3*(z^-n),n,-%inf,1); +X=X1+X2+X3;
\ No newline at end of file diff --git a/67/CH7/EX7.52/example752.sce b/67/CH7/EX7.52/example752.sce new file mode 100755 index 000000000..df139447f --- /dev/null +++ b/67/CH7/EX7.52/example752.sce @@ -0,0 +1,6 @@ +//Example 7.52
+//Z transform of x[n]=(2^n)u[n-2]
+clc;
+syms z n;
+x=2^n;
+X=symsum(x*(z^-n),n,2,%inf);
\ No newline at end of file diff --git a/67/CH7/EX7.54/example754.sce b/67/CH7/EX7.54/example754.sce new file mode 100755 index 000000000..553a9ead9 --- /dev/null +++ b/67/CH7/EX7.54/example754.sce @@ -0,0 +1,9 @@ +//Example 7.54
+clc;
+clear;
+z=poly(0,'z');
+X=[2;3*z^-1;4*z^-2];
+n=0:2;
+ZI=z^n';
+x=numer(X.*ZI);
+disp(x,'x[n]=');
\ No newline at end of file diff --git a/67/CH7/EX7.56/example756.sce b/67/CH7/EX7.56/example756.sce new file mode 100755 index 000000000..37044533b --- /dev/null +++ b/67/CH7/EX7.56/example756.sce @@ -0,0 +1,20 @@ +//Example 7.56
+//Determine the input x[n] if h[n]=[1,2,3]and y[n]=[1,1,2,-1,3]
+clc;
+clear;
+function[za]=ztransfer(sequence,n)
+ z=poly(0,'z','r')
+ za=sequence*(1/z)^n'
+endfunction
+z=poly(0,'z');
+h=[1,2,3];
+n1=0:length(h)-1;
+H=ztransfer(h,n1);
+y=[1,1,2,-1,3];
+n2=0:length(y)-1;
+Y=ztransfer(y,n2);
+X=Y/H;
+funcprot(0);
+funcprot(0);
+x=ldiv(1-z+z^2,z^2,3);
+disp(x,'x[n]=');
\ No newline at end of file diff --git a/67/CH7/EX7.59.a/example759a.sce b/67/CH7/EX7.59.a/example759a.sce new file mode 100755 index 000000000..ee872dd32 --- /dev/null +++ b/67/CH7/EX7.59.a/example759a.sce @@ -0,0 +1,6 @@ +//Example 7.59a
+//Z transform of x[n]=-(a^n)u[-n-1]
+clc;
+syms a n z;
+x=-(a^n);
+X=symsum(x*(z^-n),n,-%inf,-1);
\ No newline at end of file diff --git a/67/CH7/EX7.59.b/example759b.sce b/67/CH7/EX7.59.b/example759b.sce new file mode 100755 index 000000000..5258b88ec --- /dev/null +++ b/67/CH7/EX7.59.b/example759b.sce @@ -0,0 +1,6 @@ +//Example 7.59b
+//Z transform of x[n]=(a^-n)u[-n-1]
+clc;
+syms a n z;
+x=(a^-n);
+X=symsum(x*(z^-n),n,-%inf,-1);
\ No newline at end of file diff --git a/67/CH7/EX7.6/example76.sce b/67/CH7/EX7.6/example76.sce new file mode 100755 index 000000000..a99b45cd6 --- /dev/null +++ b/67/CH7/EX7.6/example76.sce @@ -0,0 +1,6 @@ +//Example 7.6
+clc;
+syms n z;
+x=1;
+X=symsum(x*(z^-n),n,-%inf,0);
+disp(X,'X(z)=');
diff --git a/67/CH7/EX7.61.a/example761a.sce b/67/CH7/EX7.61.a/example761a.sce new file mode 100755 index 000000000..9eda77966 --- /dev/null +++ b/67/CH7/EX7.61.a/example761a.sce @@ -0,0 +1,8 @@ +//Example 7.61a
+clc;
+syms z n;
+x1=(1/2)^n;
+X1=symsum(x1*(z^-n),n,0,%inf);
+x2=(1/3)^n;
+X2=symsum(x2*(z^-n),n,0,%inf);
+X=X1+X2;
\ No newline at end of file diff --git a/67/CH7/EX7.61.b/example761b.sce b/67/CH7/EX7.61.b/example761b.sce new file mode 100755 index 000000000..bc6872be6 --- /dev/null +++ b/67/CH7/EX7.61.b/example761b.sce @@ -0,0 +1,8 @@ +//Example 7.61b
+clc;
+syms z n;
+x1=(1/3)^n;
+X1=symsum(x1*(z^-n),n,0,%inf);
+x2=(1/2)^n;
+X2=symsum(x2*(z^-n),n,-%inf,-1);
+X=X1+X2;
\ No newline at end of file diff --git a/67/CH7/EX7.61.c/example761c.sce b/67/CH7/EX7.61.c/example761c.sce new file mode 100755 index 000000000..af90ecf4a --- /dev/null +++ b/67/CH7/EX7.61.c/example761c.sce @@ -0,0 +1,8 @@ +//Example 7.61c
+clc;
+syms z n;
+x1=(1/2)^n;
+X1=symsum(x1*(z^-n),n,0,%inf);
+x2=(1/3)^n;
+X2=symsum(x2*(z^-n),n,-%inf,-1);
+X=X1+X2;
\ No newline at end of file diff --git a/67/CH7/EX7.65/example765.sce b/67/CH7/EX7.65/example765.sce new file mode 100755 index 000000000..41ee2f1ed --- /dev/null +++ b/67/CH7/EX7.65/example765.sce @@ -0,0 +1,8 @@ +//Example 7.65
+clc;
+syms z n;
+h1=(1/2)^n;
+H1=symsum(h1*(z^-n),n,0,%inf);
+h2=(-1/4)^n;
+H2=symsum(h2*(z^-n),n,0,%inf);
+H=(H1+H2)/2;
\ No newline at end of file diff --git a/67/CH7/EX7.68.a/example768a.sce b/67/CH7/EX7.68.a/example768a.sce new file mode 100755 index 000000000..08b6765c6 --- /dev/null +++ b/67/CH7/EX7.68.a/example768a.sce @@ -0,0 +1,6 @@ +//Example 7.68a
+//Z transform of x[n]=u[n]
+clc;
+syms n z;
+x=1;
+X=symsum(x*(z^-n),n,0,%inf);
\ No newline at end of file diff --git a/67/CH7/EX7.68.b/example768b.sce b/67/CH7/EX7.68.b/example768b.sce new file mode 100755 index 000000000..b8837d04f --- /dev/null +++ b/67/CH7/EX7.68.b/example768b.sce @@ -0,0 +1,6 @@ +//Example 7.68b
+//Z transform of x[n]=-u[-n-1]
+clc;
+syms n z;
+x=-1;
+X=symsum(x*(z^-n),n,-%inf,-1);
\ No newline at end of file diff --git a/67/CH7/EX7.7/example77.sce b/67/CH7/EX7.7/example77.sce new file mode 100755 index 000000000..9f4aaa626 --- /dev/null +++ b/67/CH7/EX7.7/example77.sce @@ -0,0 +1,17 @@ +//Example 7.7 +clc; +syms n z; +X1=0; +X2=0; +for i=0:2:4 + x1=(1/2)^i; + X1=X1+x1*z^-i; +end +for i=1:2:5 + x2=(1/3)^i; + X2=X2+x2*z^-i; +end +x3=2^n; +X3=symsum(x3*(z^-n),n,-%inf,1); +X=X1+X2+X3; +disp(X,'X(z)=');
\ No newline at end of file diff --git a/67/CH8/EX8.1/example81.sce b/67/CH8/EX8.1/example81.sce new file mode 100755 index 000000000..356a23826 --- /dev/null +++ b/67/CH8/EX8.1/example81.sce @@ -0,0 +1,18 @@ +//Example 8.1
+//Determine the convolution of the two finite duration sequence
+clc;
+x=[1,1,1];
+n1=-1:1;
+h=[1,1,1];
+n2=-1:1;
+y=convol(x,h);
+n=-2:1:2;
+disp(y,'y[n]=');
+a = gca ();
+a.y_location ="origin";
+a.x_location ="origin";
+plot2d3(n,round(y),5);
+poly1=a.children(1).children (1);
+poly1.thickness=2;
+xtitle('Plot of sequence y[n]','n','y[n]');
+funcprot(0);
\ No newline at end of file diff --git a/67/CH8/EX8.10/example810.sce b/67/CH8/EX8.10/example810.sce new file mode 100755 index 000000000..06da8b82d --- /dev/null +++ b/67/CH8/EX8.10/example810.sce @@ -0,0 +1,7 @@ +//Example 8.10
+//Find the IDFT of the following sequence
+clc;
+i=sqrt(-1);
+X=[3,2+i,1,2-i];
+x=fft(X,1);
+disp(x,'x[n]=');
\ No newline at end of file diff --git a/67/CH8/EX8.11/example811.sce b/67/CH8/EX8.11/example811.sce new file mode 100755 index 000000000..43ea8d057 --- /dev/null +++ b/67/CH8/EX8.11/example811.sce @@ -0,0 +1,6 @@ +//Example 8.11
+//For the given x[n] determine X[k] using FFT algorithm
+clc;
+x=[1,2,3,4,4,3,2,1];
+X=fft(x,-1);
+disp(X,'X[k]=');
\ No newline at end of file diff --git a/67/CH8/EX8.12/example812.sce b/67/CH8/EX8.12/example812.sce new file mode 100755 index 000000000..c1c7f0d67 --- /dev/null +++ b/67/CH8/EX8.12/example812.sce @@ -0,0 +1,6 @@ +//Example 8.12
+//For the given x[n] determine X[k] using FFT algorithm
+clc;
+x=[0,1,2,3,4,5,6,7];
+X=fft(x,-1);
+disp(X,'X[k]=');
\ No newline at end of file diff --git a/67/CH8/EX8.13/example813.sce b/67/CH8/EX8.13/example813.sce new file mode 100755 index 000000000..4b2059d0f --- /dev/null +++ b/67/CH8/EX8.13/example813.sce @@ -0,0 +1,6 @@ +//Example 8.13
+//Find the DFT of the following sequence
+clc;
+h=[1/3,1/3,1/3];
+H=fft(h,-1);
+disp(H,'H[k]=');
\ No newline at end of file diff --git a/67/CH8/EX8.14/example814.sce b/67/CH8/EX8.14/example814.sce new file mode 100755 index 000000000..247762739 --- /dev/null +++ b/67/CH8/EX8.14/example814.sce @@ -0,0 +1,6 @@ +//Example 8.14
+//Find the DFT of the following sequence
+clc;
+h=[1/3,1/3,1/3];
+H=fft(h,-1);
+disp(H,'H[k]=');
\ No newline at end of file diff --git a/67/CH8/EX8.15/example815.sce b/67/CH8/EX8.15/example815.sce new file mode 100755 index 000000000..647a3af32 --- /dev/null +++ b/67/CH8/EX8.15/example815.sce @@ -0,0 +1,9 @@ +//Example 8.15
+//Obtain the DFT of x[n]=(a^n).u[n]
+clc;
+a=0.5;
+for i=0:1:7
+ x(i+1)=a.^i;
+end
+X=fft(x,-1);
+disp(X,'X[k]=');
\ No newline at end of file diff --git a/67/CH8/EX8.16/example816.sce b/67/CH8/EX8.16/example816.sce new file mode 100755 index 000000000..2869a1824 --- /dev/null +++ b/67/CH8/EX8.16/example816.sce @@ -0,0 +1,6 @@ +//Example 8.16
+//Program to find DFT of sequence x=[1,1,1,1,1,1,1,0]
+clc ;
+x=[1,1,1,1,1,1,1,0];
+X=fft(x,-1);
+disp(X,'X[k]=');
\ No newline at end of file diff --git a/67/CH8/EX8.17/example817.sce b/67/CH8/EX8.17/example817.sce new file mode 100755 index 000000000..81e72ff1c --- /dev/null +++ b/67/CH8/EX8.17/example817.sce @@ -0,0 +1,6 @@ +//Example 8.17
+//Determine the DFT of the following sequence
+clc;
+x=[0,0,1,1,1,1,1,0,0,0];
+X=fft(x,-1);
+disp(X,'X[k]=');
\ No newline at end of file diff --git a/67/CH8/EX8.2/example82.sce b/67/CH8/EX8.2/example82.sce new file mode 100755 index 000000000..740e9cd29 --- /dev/null +++ b/67/CH8/EX8.2/example82.sce @@ -0,0 +1,7 @@ +//Example 8.2
+//Find the response of an FIR filter with impulse response h[n]=[1,2,4] //to the input sequence x[n]=[1,2]
+clc;
+x=[1,2];
+h=[1,2,4];
+Y=convol(x,h);
+disp(Y,'y[n]=');
\ No newline at end of file diff --git a/67/CH8/EX8.3/example8_3.sce b/67/CH8/EX8.3/example8_3.sce new file mode 100755 index 000000000..9494b3590 --- /dev/null +++ b/67/CH8/EX8.3/example8_3.sce @@ -0,0 +1,7 @@ +//Example 8.3
+//Compute DFT of x(n)={1,1,0,0} and IDFT of y(n)={1,0,1,0}
+clc;
+x=[1,1,0,0];
+Y=[1,0,1,0];
+X=fft(x,-1);
+y=fft(Y,1);
\ No newline at end of file diff --git a/67/CH8/EX8.4/example84.sce b/67/CH8/EX8.4/example84.sce new file mode 100755 index 000000000..8962598a3 --- /dev/null +++ b/67/CH8/EX8.4/example84.sce @@ -0,0 +1,6 @@ +//Example 8.4
+//Compute DFT of the following sequence
+clc;
+x=[0.25,0.25,0.25];
+X=fft(x,-1);
+disp(X,'X[k]=');
\ No newline at end of file diff --git a/67/CH8/EX8.5/example85.sce b/67/CH8/EX8.5/example85.sce new file mode 100755 index 000000000..9d058fe45 --- /dev/null +++ b/67/CH8/EX8.5/example85.sce @@ -0,0 +1,7 @@ +//Example 8.5
+//Find the DFT of the following sequence
+clc;
+x=[0.2,0.2,0.2];
+n=-1:1;
+X=fft(x,-1);
+disp(X,'X[k]=');
\ No newline at end of file diff --git a/67/CH8/EX8.6/example86.sce b/67/CH8/EX8.6/example86.sce new file mode 100755 index 000000000..fa434746b --- /dev/null +++ b/67/CH8/EX8.6/example86.sce @@ -0,0 +1,6 @@ +//Example 8.6
+//Determine the DFT of the following sequence
+clc;
+x=[1,1,2,2,3,3];
+X=fft(x,-1);
+disp(X,'X[k]=');
\ No newline at end of file diff --git a/67/CH8/EX8.7/example87.sce b/67/CH8/EX8.7/example87.sce new file mode 100755 index 000000000..3c9befd72 --- /dev/null +++ b/67/CH8/EX8.7/example87.sce @@ -0,0 +1,8 @@ +//Example 8.7
+//DFT of x[n]=a.^n
+clc;
+a=0.5;//Say for a=0.5
+n=0:4;
+x=a.^n;
+X=fft(x,-1);
+disp(X,'X[k]=');
\ No newline at end of file diff --git a/67/CH8/EX8.8/example88.sce b/67/CH8/EX8.8/example88.sce new file mode 100755 index 000000000..1891ea27c --- /dev/null +++ b/67/CH8/EX8.8/example88.sce @@ -0,0 +1,8 @@ +//Example 8.8
+//Compute 4-point DFT of the sequence x[n]=cos(n*pi/4)
+clc;
+n=0:3;
+pi=22/7;
+x=cos(n*pi/4);
+X=fft(x,-1);
+disp(X,'X[k]=');
\ No newline at end of file diff --git a/67/CH8/EX8.9/example89.sce b/67/CH8/EX8.9/example89.sce new file mode 100755 index 000000000..01f10f784 --- /dev/null +++ b/67/CH8/EX8.9/example89.sce @@ -0,0 +1,6 @@ +//Example 8.9
+//Computing IDFT of the following sequence
+clc;
+X=[1,2,3,4];
+x=fft(X,1);
+disp(x,'x[n]=');
\ No newline at end of file |