Merge pull request #1 from prashantsinalkar/masterHEADmaster

Initial commit

1 files changed, 101 insertions, 0 deletions

diff --git a/Digital_Signal_Processing_by_S_Salivahanan/12-Spectral_Estimation.ipynb b/Digital_Signal_Processing_by_S_Salivahanan/12-Spectral_Estimation.ipynb
new file mode 100644
index 0000000..c527bc1
--- /dev/null
+++ b/Digital_Signal_Processing_by_S_Salivahanan/12-Spectral_Estimation.ipynb
@@ -0,0 +1,101 @@
+{
+"cells": [
+ {
+		   "cell_type": "markdown",
+	   "metadata": {},
+	   "source": [
+       "# Chapter 12: Spectral Estimation"
+	   ]
+	},
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 12.2: Power_Spectrum.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//Example 12.2\n",
+"clc;clear;close;\n",
+"N=8;n=0:N-1;\n",
+"f1=0.6;f2=0.62;\n",
+"x=cos(2*%pi*f1*n)+cos(2*%pi*f2*n);\n",
+"L1=8;\n",
+"for k=0:L1-1\n",
+"    P1(k+1)=1/N*abs(x*(cos(%pi*n*k/L1)-%i*sin(%pi*n*k/L1))')^2\n",
+"end\n",
+"L2=16;\n",
+"for k=0:L2-1\n",
+"    P2(k+1)=1/N*abs(x*(cos(%pi*n*k/L2)-%i*sin(%pi*n*k/L2))')^2;\n",
+"end\n",
+"L3=32;\n",
+"for k=0:L3-1\n",
+"    P3(k+1)=1/N*abs(x*(cos(%pi*n*k/L3)-%i*sin(%pi*n*k/L3))')^2;\n",
+"end\n",
+"subplot(311);\n",
+"plot2d3(0:L1-1,P1);title('L=8');\n",
+"subplot(312);\n",
+"plot2d3(0:L2-1,P2);title('L=16');\n",
+"subplot(313);\n",
+"plot2d3(0:L3-1,P3);title('L=32');"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 12.4: Frequency_resolution.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//Example 12.4\n",
+"clc;clear;close;\n",
+"N=1000;\n",
+"Q=10;\n",
+"disp(N,'Length of sample sequence N=',Q,'Quality factor Q=');\n",
+"f_bart=Q/(1.11*N);\n",
+"f_w=Q/(1.39*N);\n",
+"f_bt=Q/(2.34*N);\n",
+"disp(f_bart,'Bartlett Frequency resolution =');\n",
+"disp(f_w,'Welch Frequency resolution =');\n",
+"disp(f_bt,'Blackman Turkey Frequency resolution =');"
+   ]
+   }
+],
+"metadata": {
+		  "kernelspec": {
+		   "display_name": "Scilab",
+		   "language": "scilab",
+		   "name": "scilab"
+		  },
+		  "language_info": {
+		   "file_extension": ".sce",
+		   "help_links": [
+			{
+			 "text": "MetaKernel Magics",
+			 "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md"
+			}
+		   ],
+		   "mimetype": "text/x-octave",
+		   "name": "scilab",
+		   "version": "0.7.1"
+		  }
+		 },
+		 "nbformat": 4,
+		 "nbformat_minor": 0
+}