From 476705d693c7122d34f9b049fa79b935405c9b49 Mon Sep 17 00:00:00 2001 From: prashantsinalkar Date: Tue, 14 Apr 2020 10:19:27 +0530 Subject: Initial commit --- .../12-Spectral_Estimation.ipynb | 101 +++++++++++++++++++++ 1 file changed, 101 insertions(+) create mode 100644 Digital_Signal_Processing_by_S_Salivahanan/12-Spectral_Estimation.ipynb (limited to 'Digital_Signal_Processing_by_S_Salivahanan/12-Spectral_Estimation.ipynb') 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 +} -- cgit