summaryrefslogtreecommitdiff
path: root/803
diff options
context:
space:
mode:
Diffstat (limited to '803')
-rw-r--r--803/CH1/EX1.3/ex1_3.sce9
-rw-r--r--803/CH2/EX2.1/ex2_1.sce7
-rw-r--r--803/CH2/EX2.2/ex2_2.sce14
-rw-r--r--803/CH3/EX3.2/ex3_2.sce7
-rw-r--r--803/CH3/EX3.3/ex3_3.sce12
-rw-r--r--803/CH3/EX3.4/ex3_4.sce12
-rw-r--r--803/CH5/EX5.1/ex5_1.sce10
-rw-r--r--803/CH5/EX5.2/ex5_2.sce13
-rw-r--r--803/CH5/EX5.3/ex5_3.sce22
-rw-r--r--803/CH5/EX5.4/ex5_4.sce22
-rw-r--r--803/CH5/EX5.5/ex5_5.sce10
-rw-r--r--803/CH5/EX5.6/ex5_6.sce8
-rw-r--r--803/CH6/EX6.1/ex6_1.sce32
-rw-r--r--803/CH6/EX6.2/ex6_2.sce22
-rw-r--r--803/CH6/EX6.3/ex6_3.sce16
-rw-r--r--803/CH6/EX6.4/ex6_4.sce7
-rw-r--r--803/CH7/EX7.1/ex7_1.sce21
-rw-r--r--803/CH7/EX7.4/ex7_4.sce6
-rw-r--r--803/CH8/EX8.2/ex8_2.sce10
-rw-r--r--803/CH8/EX8.4/ex8_4.pngbin0 -> 20789 bytes
-rw-r--r--803/CH8/EX8.4/ex8_4.sce18
-rw-r--r--803/CH8/EX8.5/ex8_5.sce21
22 files changed, 299 insertions, 0 deletions
diff --git a/803/CH1/EX1.3/ex1_3.sce b/803/CH1/EX1.3/ex1_3.sce
new file mode 100644
index 000000000..ff62b3596
--- /dev/null
+++ b/803/CH1/EX1.3/ex1_3.sce
@@ -0,0 +1,9 @@
+clc
+fk=2.5; //Desired resolution between frequencies
+fd=1.75*10^3;//maximum frequency
+dmin=1/fk;//minimum record length
+tmax=1/(2*fd);//maximum time between sampling
+N=1/(tmax*fk);//minimum number of sampling points
+disp("Seconds",dmin,"Minimum record length");
+disp("Seconds",tmax,"Maximum time betweeen sampling points");
+disp(N/2,"Minimum number of sampling points"); \ No newline at end of file
diff --git a/803/CH2/EX2.1/ex2_1.sce b/803/CH2/EX2.1/ex2_1.sce
new file mode 100644
index 000000000..a462c73d8
--- /dev/null
+++ b/803/CH2/EX2.1/ex2_1.sce
@@ -0,0 +1,7 @@
+inputl=-16;
+spuresp=70;..//spurious response
+A0=-22;..//input level
+N=4;..//fouth order
+IP4=inputl+(spuresp/(N-1));
+DelS=(IP4-A0)*(N-1);..//distortion product
+disp("dB",DelS,"Distortion Product"); \ No newline at end of file
diff --git a/803/CH2/EX2.2/ex2_2.sce b/803/CH2/EX2.2/ex2_2.sce
new file mode 100644
index 000000000..2895ec30f
--- /dev/null
+++ b/803/CH2/EX2.2/ex2_2.sce
@@ -0,0 +1,14 @@
+clc
+SNRdB=6;..//signal to noise ratio
+SNR=10^(SNRdB/10);
+Rg=50;..//resistance
+Ft=14.87;..//noise factor
+Bn=12000;..//bandwidth
+k=1.38*10^-23;..//boltzmann constant
+T=290;..//room temperature
+C=18.2;
+ev=sqrt(k*Ft*Bn*SNR*Rg*T);
+T1=(Ft-1)*T+(C+273);
+ev1=sqrt(k*Ft*Bn*SNR*Rg*T1);
+disp("Volts",ev,"Sensitivity at room temperature");
+disp("Volts",ev1,"Sensitivity at an average value of 18.2 deg C"); \ No newline at end of file
diff --git a/803/CH3/EX3.2/ex3_2.sce b/803/CH3/EX3.2/ex3_2.sce
new file mode 100644
index 000000000..5f53801de
--- /dev/null
+++ b/803/CH3/EX3.2/ex3_2.sce
@@ -0,0 +1,7 @@
+clc
+vr=20*1.852*(1000/3600);..//target radial speed
+f=10^9;..//base frequency
+c=3*(10^8);
+lamda=(c/f);
+fd=2*(vr/lamda);..//doppler speed
+disp("Hz",fd,"Doppler Shift is")
diff --git a/803/CH3/EX3.3/ex3_3.sce b/803/CH3/EX3.3/ex3_3.sce
new file mode 100644
index 000000000..8fd9d4093
--- /dev/null
+++ b/803/CH3/EX3.3/ex3_3.sce
@@ -0,0 +1,12 @@
+Ts=10^(-6);..//pulse width
+lamda=0.1;..//transmitter's duty cycle
+DH=3;..//horizontal dimension
+DV=0.5;..//vertical dimension
+DelRrmin=lamda/(2*Ts);..//resolution in range rate
+DelRamin=lamda/(Ts^2);..//resolution in acceleration range
+DelFIH=(lamda/DH)*(180/%pi);..//azimuth beamwidth
+DelFIV=(lamda/DV)*(180/%pi);..//elevation beamwidth
+disp("m/s",DelRrmin,"The nominal resolution in range rate is");
+disp("m/s^2",DelRamin,"The nominal resolution in acceleration range is");
+disp("Degree",DelFIH,"The azimuth beamwidth of the antenna is");
+disp("Degree",DelFIV,"The elevation beamwidth of the antenna is"); \ No newline at end of file
diff --git a/803/CH3/EX3.4/ex3_4.sce b/803/CH3/EX3.4/ex3_4.sce
new file mode 100644
index 000000000..7b675a415
--- /dev/null
+++ b/803/CH3/EX3.4/ex3_4.sce
@@ -0,0 +1,12 @@
+clc
+Va=500;..//radar carrying platform velocity
+lambda=0.3;
+dutcyc=10^(-6);
+PRF=10000;
+kp=1.25;..//radar platform dependent factor
+Vcomp=kp*Va;..//compression velocity
+DR=360*dutcyc*(Vcomp/lambda);..//derotation rate
+Vt=(lambda/2)*PRF;..//target velocity
+disp("m/s",Vcomp,"The compensation velocity is");
+disp("degree/range-cell",DR,"The derotation rate is");
+disp("m/s",Vt,"The maximum unambiguous target velocity is");
diff --git a/803/CH5/EX5.1/ex5_1.sce b/803/CH5/EX5.1/ex5_1.sce
new file mode 100644
index 000000000..bd98e888a
--- /dev/null
+++ b/803/CH5/EX5.1/ex5_1.sce
@@ -0,0 +1,10 @@
+clc
+Pr=10^-14;..//receiver sensitivity
+R=2.5*10^5;
+sigma=15;
+lambda=0.12;
+Gr=2000;..//receiver gain
+Gt=2000;..//transmitter gain
+pi=3.14;
+Ppeak=2.985*(Pr*(4*%pi)^3*R^4/Gt*Gr*sigma*lambda^2);
+disp("W",Ppeak,"Minimum peak transmitter power"); \ No newline at end of file
diff --git a/803/CH5/EX5.2/ex5_2.sce b/803/CH5/EX5.2/ex5_2.sce
new file mode 100644
index 000000000..1bb77279c
--- /dev/null
+++ b/803/CH5/EX5.2/ex5_2.sce
@@ -0,0 +1,13 @@
+clc
+thetaH=17.3;..//beamwidth in E plane
+thetaE=17.3;..//beamwidth in H plane
+eta=0.5;..//antenna's efficiency
+c=3*10^8;
+f=4.5*10^9;
+lambda=c/f;
+A=(78*lambda/thetaH);..//aperture dimension in H plane
+B=(54*lambda/thetaE);..//aperture dimension in E plane
+G=10*log10(eta*4*%pi*A*B/lambda^2);
+disp("m",A,"Aperture dimension in H plane");
+disp("m",B,"Aperture dimension in E plane");
+disp("dB",G,"Gain"); \ No newline at end of file
diff --git a/803/CH5/EX5.3/ex5_3.sce b/803/CH5/EX5.3/ex5_3.sce
new file mode 100644
index 000000000..89f3b90f9
--- /dev/null
+++ b/803/CH5/EX5.3/ex5_3.sce
@@ -0,0 +1,22 @@
+clc
+Ltot=123.03;
+Fn=3.16;..//noise factor
+Gt=10;..//transmitter gain
+Pt=1500;..//transmitter peak power
+To=296.7;
+R=10^5;
+k=1.38*10^-23;..//boltzmann constant
+sigma=1.5;
+Ae=8;
+Bn=10^3;..//bandwidth
+pi=3.14;
+F=1;
+w=1.67;
+angle=(4*%pi/Gt)*180/%pi;
+Beamwidth=sqrt(angle);..//elevation beamwidth
+SNR=log10(Pt*Gt*Ae*sigma*F^4/(4*%pi)^2*k*To*Bn*Fn*R^4*Ltot);..//signal to noise ratio
+T=2*%pi/w;..//time frame
+disp("dB",SNR,"Signal received");
+disp("s",T,"Time frame");
+disp("degree",Beamwidth,"Elevation Beamwidth");
+
diff --git a/803/CH5/EX5.4/ex5_4.sce b/803/CH5/EX5.4/ex5_4.sce
new file mode 100644
index 000000000..f9ec4d371
--- /dev/null
+++ b/803/CH5/EX5.4/ex5_4.sce
@@ -0,0 +1,22 @@
+clc
+Ts=1;
+n=3;
+Pd=0.9;..//detection probability
+Gt=10^5;..//transmitter gain
+k=1.38*10^-23;..//boltzmann constant
+tau=10^-6;
+To=305.8;
+Ae=8;
+Pfa=10^-6;..//probability of false alarm
+Pt=10^5;..//transmitter power
+F=1;..//noise factor
+sigma=3.2;
+Bn=1/tau;..//bandwidth
+pi=3.14;
+Fn=3.162;
+e=2.72;
+Ltot=67.608;
+alpha=0.75*(1+0.667*e^(-n/3));
+SNR=(alpha*log10(log(2)/Pfa)/n^(0.667)*(log10(1/Pd))^0.667);
+R=(tau*Pt*Gt*Ae*sigma*F^4/(4*%pi)^2*k*To*Bn*Fn*SNR*Ltot)^0.25;
+disp("m",R,"Maximum Range");
diff --git a/803/CH5/EX5.5/ex5_5.sce b/803/CH5/EX5.5/ex5_5.sce
new file mode 100644
index 000000000..d9134c6b0
--- /dev/null
+++ b/803/CH5/EX5.5/ex5_5.sce
@@ -0,0 +1,10 @@
+clc
+Bn=10^9;..//bandwidth
+T=300;..//room temperature
+k=1.38*10^-23;..//boltzmann constant
+h=6.6256*10^-34;
+f=10^9;..//frequency
+Nthermal=(10*log10(k*T*Bn));..//thermal noise power
+Ni=(10*log10(h*f*Bn));..//quantum noise power
+disp("dB",Nthermal,"Thermal noise power");
+disp("dB",Ni,"Quantum noise power"); \ No newline at end of file
diff --git a/803/CH5/EX5.6/ex5_6.sce b/803/CH5/EX5.6/ex5_6.sce
new file mode 100644
index 000000000..8e0f4399f
--- /dev/null
+++ b/803/CH5/EX5.6/ex5_6.sce
@@ -0,0 +1,8 @@
+clc
+Pt=100;..//transmitter power
+Gt=1000;..//transmitter gain
+f=3*10^9;..//frequency
+R=10^5;..//resistance
+pi=3.14;
+Smin=(Pt*(0.3*Gt/(4*%pi*f*R))^2);..//Received power
+disp("W",Smin,"Received power"); \ No newline at end of file
diff --git a/803/CH6/EX6.1/ex6_1.sce b/803/CH6/EX6.1/ex6_1.sce
new file mode 100644
index 000000000..3e250b662
--- /dev/null
+++ b/803/CH6/EX6.1/ex6_1.sce
@@ -0,0 +1,32 @@
+clc
+f=1.2*(10^6);...............................//Assigning values to parameters
+Yn=56;
+deltalat=37.45;
+x=23.44*sind(0.9856*(Yn-80.7));
+deltah=acosd(tand(x)/tand(deltalat));
+y=acosd((sind(deltalat)*sind(x))+(cosd(deltalat)*cosd(x)*cosd(deltah)));
+disp("degree",x,"Solar declination");
+disp("degree",deltah,"hour angle of the sun");
+disp("degree",y,"Solar zenith angle");
+m=4.8*(10^(-26));................//mean molecular mass of air
+k=1.38*(10^(-23));...............//Boltzmann constant
+g=9.8;...........................//gravitational constant
+e=1.6*(10^(-19));....//electron charge
+me=9.11*(10^(-31));
+epsilon=8.85*(10^(-12));
+hmax=[100 200 300 120 250];
+T=[341 1360 1710 341 1540];
+f=1.2*(10^6);
+Nm=[1.5*10^11 3*10^11 12.5*10^11 0.8*1010 4*1011];
+h=[122 256 335 132 276];
+for i=1:5
+ H(i)=(g * m * [h(i)-hmax(i)])/ (k * T(i));
+ No(i)=(Nm(i)*(secd(x)^0.5));
+ Ne(i)=((No(i))*exp(0.5*(1-(H(i)-secd(x)*exp(-H(i))))));
+ fc(i)=(1/(2*3.14))*(sqrt(((e^2)*(Ne(i)))/(epsilon*me)));
+ n(i)=sqrt(1-((fc(i)^2)/(f^2)));
+ disp(Ne(i),"Electron density of each layer:");
+ disp(fc(i),"Critical frequency of each layer:");
+ disp(n(i),"Refractive index of each layer:");
+end
+
diff --git a/803/CH6/EX6.2/ex6_2.sce b/803/CH6/EX6.2/ex6_2.sce
new file mode 100644
index 000000000..dc8c19a90
--- /dev/null
+++ b/803/CH6/EX6.2/ex6_2.sce
@@ -0,0 +1,22 @@
+clc
+deltalat=-5;..//geographic latitude
+ag=9;..//apparent elevation angle
+pi=3.14;
+x=6378.4*(10^3);..//equatorial radius of earth
+f=15*(10^6);
+r=[6493 6593 6693 6793];
+fc=[3.04*(10^6) 4.38*(10^6) 5.86*(10^6)];
+for i=1:4
+ j=1:3
+ n(j)=sqrt(1-((fc(j)^2)/f^2));..//refractive index
+ a(i)=acosd((x*cosd(ag))/r(i));..//apparent elevation angles
+ phi(i)=sind((r(i)*cosd(a(i)))/(r(i+1)));
+ theta(i)=(%pi/2)-a(i)-phi(i);
+ R012=sqrt((r(1)^2)+(r(4)^2)-(2*r(1)*r(4)*cosd(sum(theta(i)))));
+ R(i)=[r(i+1)*sind(theta(i))/cosd(a(i))];
+ dela=a(i)-acosd((r(4)/R012)*sind(sum(theta(i))));
+ delR=sum(R(i)/n(j))-R012;
+ disp(n(j),"Refractive index of each layer");
+ disp("degree",dela,"The refraction angle error");
+ disp("km",delR,"The range error");
+end; \ No newline at end of file
diff --git a/803/CH6/EX6.3/ex6_3.sce b/803/CH6/EX6.3/ex6_3.sce
new file mode 100644
index 000000000..80b2582f6
--- /dev/null
+++ b/803/CH6/EX6.3/ex6_3.sce
@@ -0,0 +1,16 @@
+clc
+f=12*(10^6);..........................//Assigning values to parameters
+c=3*(10^8);
+Vt=85*(10^3);
+alphag=7.6;
+fc=5.6*(10^6);
+y=6377;
+ro=6527*(10^3);
+x=15;................................//orientation of target velocity
+ho=150;
+n=sqrt(1-((fc^2)/f^2));
+b=acosd(y*cosd(alphag)/(n*(ro+ho)));
+delf=((-2*f*Vt*sind(x)*b)/c)
+disp(n,"Refractive index");
+disp(b,"a");
+disp("Hz",delf,"The doppler frequency error is");
diff --git a/803/CH6/EX6.4/ex6_4.sce b/803/CH6/EX6.4/ex6_4.sce
new file mode 100644
index 000000000..5f78705e3
--- /dev/null
+++ b/803/CH6/EX6.4/ex6_4.sce
@@ -0,0 +1,7 @@
+clc
+T=[341 1360 1710 341 1540];
+N=[57.6*10^5 81.5*10^5 166.5*10^5 423 947];
+for i=1:5
+ lambda=69*sqrt(T(i)/N(i));
+ disp("mm",lambda,"The debye length is");
+end \ No newline at end of file
diff --git a/803/CH7/EX7.1/ex7_1.sce b/803/CH7/EX7.1/ex7_1.sce
new file mode 100644
index 000000000..8ad0e68be
--- /dev/null
+++ b/803/CH7/EX7.1/ex7_1.sce
@@ -0,0 +1,21 @@
+clc
+D1=127;................................//Array aperture
+D2=67;
+D3=335;
+f1=3*(10^6);..//frequency at 3MHz
+f2=30*(10^6);..//frequency at 30MHz
+c=3*(10^8);
+lambda1=c/f1;
+lambda2=c/f2;
+R1=2*(D1^3) ./lambda1;
+R2=2*(D2^3) ./lambda1;
+R3=2*(D3^3) ./lambda1;
+R4=2*(D1^3) ./lambda2;
+R5=2*(D2^3) ./lambda2;
+R6=2*(D3^3) ./lambda2;
+disp("m",R1,"Range at 3MHz");
+disp("m",R2,"Range at 3MHz");
+disp("m",R3,"Range at 3MHz");
+disp("m",R4,"Range at 30MHz");
+disp("m",R5,"Range at 30MHz");
+disp("m",R6,"Range at 30MHz"); \ No newline at end of file
diff --git a/803/CH7/EX7.4/ex7_4.sce b/803/CH7/EX7.4/ex7_4.sce
new file mode 100644
index 000000000..7d801be41
--- /dev/null
+++ b/803/CH7/EX7.4/ex7_4.sce
@@ -0,0 +1,6 @@
+clc
+N=8;................................//order of 8
+P1=0.95;
+P2=0.05;
+P8=((1-P2)^8)*100;
+disp("%",P8,"Probability of finding 8kHz channel is");
diff --git a/803/CH8/EX8.2/ex8_2.sce b/803/CH8/EX8.2/ex8_2.sce
new file mode 100644
index 000000000..9724692fc
--- /dev/null
+++ b/803/CH8/EX8.2/ex8_2.sce
@@ -0,0 +1,10 @@
+clc
+R=[1/2 2/5 2/3];.........................//probability of drawing a red ball
+a=1/3;..//probability of drawing any container
+Den=sum(R)*a;
+b=R(1)*a;
+c=R(2)*a;
+d=R(3)*a;
+disp(b/Den,"ans1");
+disp(c/Den,"ans2");
+disp(d/Den,"ans3");
diff --git a/803/CH8/EX8.4/ex8_4.png b/803/CH8/EX8.4/ex8_4.png
new file mode 100644
index 000000000..888d60b31
--- /dev/null
+++ b/803/CH8/EX8.4/ex8_4.png
Binary files differ
diff --git a/803/CH8/EX8.4/ex8_4.sce b/803/CH8/EX8.4/ex8_4.sce
new file mode 100644
index 000000000..274ea86a4
--- /dev/null
+++ b/803/CH8/EX8.4/ex8_4.sce
@@ -0,0 +1,18 @@
+clc
+b=[0 1 2 3 4 5 6];..........................//breakdown/day
+f=[340 121 53 30 12 4 0];...................//frequency
+Den=sum(f);
+ans=0;
+for i=1:7
+ ans=ans+[f(i) * b(i)];
+end
+Mean=ans/Den;...............................//calculation of expectance
+disp(Mean,"mean");
+
+k=0:1:6;
+y=(exp(-Mean)*(Mean^k)). /factorial(k);
+disp(y);
+plot2d(k,y);......................//Poisson distribution
+title("Poisson Distribution",k,y);
+xlabel("Breakdown/day");
+ylabel("Probability");
diff --git a/803/CH8/EX8.5/ex8_5.sce b/803/CH8/EX8.5/ex8_5.sce
new file mode 100644
index 000000000..c5bd3544c
--- /dev/null
+++ b/803/CH8/EX8.5/ex8_5.sce
@@ -0,0 +1,21 @@
+clc
+n=10;...................................//total pulses selected
+p=0.8;..................................//probability of pulses hitting the dish
+q=0.2;..................................//probability of miss
+add=0;
+for k=2;
+ x(k)=((factorial(n)*(p^k)*((1-p)^(n-k)))/(factorial(k)*factorial(n-k)));
+ disp(x(k),"Exactly two pulses missing the target");
+end;
+for k=0:1
+ x(k)=((factorial(n)*(p^k)*((1-p)^(n-k)))/(factorial(k)*factorial(n-k)));
+ add=add+x(k);
+end;
+ y(k)=1-add;
+ disp(y(k),"Two or more pulses missing the target");
+for k=6:10
+ x(k)=((factorial(n)*(p^k)*((1-p)^(n-k)))/(factorial(k)*factorial(n-(k)));
+ y(k)=sum(x(k));
+ disp(y(k),"More than 5 pulses missing the target");
+
+end;
generated by cgit (git 2.34.1) at 2024-12-19 18:37:30 +0000